b'\n \n \n
 
\n  
sfmpipfwi (4.0)
index
user/jsun/Mmpipfwi.c
\n Visco-acoustic Forward Modeling, FWI, and RTM based on SLS model \n

\n \n \n \n \n \n
 
\n Synopsis
       sfmpipfwi Fvel=Fv.rsf Fq=Fq.rsf Fwavelet=Fw.rsf output=Fdat.rsf output=Finv.rsf Fgrad=Fgrad.rsf output=Fimg.rsf Fsrc=Fsrc.rsf Fmwt=Fmwt.rsf media=1 function=2 verb=n nb=100 coef=0.003 acqui_type=1 ns= ds= s0= sz=5 nr=acpar->nx dr=acpar->dx r0=acpar->x0 rz=1 f0= interval=1 fhi=0.5/acpar->dt flo=0. onlygrad=n wt1=acpar->t0 wt2=acpar->t0+(acpar->nt-1)*acpar->dt woff1=acpar->r0 woff2=acpar->r0+(acpar->nr-1)*acpar->dr oreo=n waterz=51 grectx=3 grectz=3 niter= conv_error= c1=1e-4 c2=0.9 v1=0. v2=10. onlygrad=n wt1=acpar->t0 wt2=acpar->t0+(acpar->nt-1)*acpar->dt woff1=acpar->r0 woff2=acpar->r0+(acpar->nr-1)*acpar->dr oreo=n waterz=0 waterzb=0 grectx=3 grectz=3 niter= repeat=1 conv_error= c1=1e-4 c2=0.9 v1=0. v2=10.
\n

\n \n \n \n \n \n
 
\n Parameters
       \n \n \n
\n  
file Fgrad=
\tauxiliary output file name
\n
\n \n\n \n \n
\n  
file Fmwt=
\tauxiliary output file name
\n
\n \n\n \n \n
\n  
file Fq=
\tauxiliary input file name
\n
\n \n\n \n \n
\n  
file Fsrc=
\tauxiliary output file name
\n
\n \n\n \n \n
\n  
file Fvel=
\tauxiliary input file name
\n
\n \n\n \n \n
\n  
file Fwavelet=
\tauxiliary input file name
\n
\n \n\n \n \n
\n  
int acqui_type=1
\tif 1, fixed acquisition; if 2, marine acquisition; if 3, symmetric acquisition
\n
\n \n\n \n \n
\n  
float c1=1e-4
\t
\n
\n \n\n \n \n
\n  
float c2=0.9
\t
\n
\n \n\n \n \n
\n  
float coef=0.003
\tabsorbing boundary coefficient
\n
\n \n\n \n \n
\n  
float conv_error=
\tfinal convergence error
\n
\n \n\n \n \n
\n  
float dr=acpar->dx
\treceiver interval
\n
\n \n\n \n \n
\n  
float ds=
\tshot interval
\n
\n \n\n \n \n
\n  
float f0=
\treference frequency
\n
\n \n\n \n \n
\n  
float fhi=0.5/acpar->dt
\t
\n
\n \n\n \n \n
\n  
float flo=0.
\t
\n
\n \n\n \n \n
\n  
int function=2
\tif 1, forward modeling; if 2, FWI; if 3, RTM
\n
\n \n\n \n \n
\n  
int grectx=3
\tgradient smoothing radius in x
\n
\n \n\n \n \n
\n  
int grectz=3
\tgradient smoothing radius in z
\n
\n \n\n \n \n
\n  
int interval=1
\twavefield storing interval
\n
\n \n\n \n \n
\n  
int media=1
\tif 1, acoustic media; if 2, visco-acoustic media
\n
\n \n\n \n \n
\n  
int nb=100
\tboundary width
\n
\n \n\n \n \n
\n  
int niter=
\titeration number
\n
\n \n\n \n \n
\n  
int nr=acpar->nx
\tnumber of receiver
\n
\n \n\n \n \n
\n  
int ns=
\tshot number
\n
\n \n\n \n \n
\n  
bool onlygrad=n [y/n]
\tonly want gradident
\n
\n \n\n \n \n
\n  
bool oreo=n [y/n]
\tkeep oreo or keep cream
\n
\n \n\n \n \n
\n  
file output=
\tauxiliary output file name
\n
\n \n\n \n \n
\n  
float r0=acpar->x0
\treceiver origin
\n
\n \n\n \n \n
\n  
int repeat=1
\trepeat resetting alpha
\n
\n \n\n \n \n
\n  
int rz=1
\treceiver depth
\n
\n \n\n \n \n
\n  
float s0=
\tshot origin
\n
\n \n\n \n \n
\n  
int sz=5
\tsource depth
\n
\n \n\n \n \n
\n  
float v1=0.
\t
\n
\n \n\n \n \n
\n  
float v2=10.
\t
\n
\n \n\n \n \n
\n  
bool verb=n [y/n]
\tverbosity flag
\n
\n \n\n \n \n
\n  
int waterz=0
\twater layer depth
\n
\n \n\n \n \n
\n  
int waterzb=0
\twater layer depth from bottom up
\n
\n \n\n \n \n
\n  
float woff1=acpar->r0
\t
\n
\n \n\n \n \n
\n  
float woff2=acpar->r0+(acpar->nr-1)*acpar->dr
\t
\n
\n \n\n \n \n
\n  
float wt1=acpar->t0
\t
\n
\n \n\n \n \n
\n  
float wt2=acpar->t0+(acpar->nt-1)*acpar->dt
\t
\n
\n \n
'