[rsf-user] =?gb2312?q?Maybe=20a=20mistake=20about=20=A6=C7=20expression=20in=20progr?= =?gb2312?q?am=20Mpnmo.c=20?=

[<xia...@ya...>]

Hi fomel:
     Thanks for your attention! Now, I find there may be a mistake in your program Mpnmo.c.Because S=1+8η. I have deduce the formula about η in your program
   following your paper: Velocity-independent time-domain seismic imaging using local event slopes, Geophysics,vol.72,no.7. S = 1+(p(t-pl)-qlt)/sqrt(qp^3*l^3).
  Meanwhile, I find the expreesion about η in your program losts a minus sign "-". That is the expression should be :
     etr[it]= -((g-dh*dh)*t/(dt*(f*h+eps*dh))+dh)/(8+sqrtf(g)+eps*dh);
    I also  demonstrate the expression about η by numerical experiment.
    First, I make some nonhyperbolic event about  η following formulat:
           t=t0(1-1/S)+sqrt(t0^2+S*l^2/v^2)/S. S=1+8η
    Then, to calculate the slope by program sfdip.
    Last, to calculate η by program sfpnmo.
    I found the true η and estimated η are just opposite.
  Plese check it, maybe I am wrong.
  student qihao
   

       
---------------------------------
 雅虎免费邮箱3.5G容量，20M附件！