Program sfmpiafdfwi2d | 2D Visco-acoustic Forward Modeling, FWI, and RTM based on SLS model 
[SYNOPSIS]
sfmpiafdfwi2d Fwfl1=Fwfl1.rsf Fwfl2=Fwfl2.rsf Fres=Fres.rsf 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 ompchunk=1 ompnth=0 inv=n onlysrc=n onlyvel=n sw=n ctr=n precsrc=n hidesrc=n nitersrc=1 ngrp=1 size=0 rectzsrc=1 rectxsrc=1 recttsrc=50 repeatsrc=1 perc=1.0f hard=0.1f 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.
[PARAMETERS]
file    | Fgrad | = |   | 	auxiliary output file name
file    | Fmwt | = |   | 	auxiliary output file name
file    | Fq | = |   | 	auxiliary input file name
file    | Fres | = |   | 	auxiliary output file name
file    | Fsrc | = |   | 	auxiliary output file name
file    | Fvel | = |   | 	auxiliary input file name
file    | Fwavelet | = |   | 	auxiliary input file name
file    | Fwfl1 | = |   | 	auxiliary output file name
file    | Fwfl2 | = |   | 	auxiliary output file name
int     | acqui_type | =1 |   | 	if 1, fixed acquisition; if 2, marine acquisition; if 3, symmetric acquisition 
float   | c1 | =1e-4 |   | 	
float   | c2 | =0.9 |   | 	
float   | coef | =0.003 |   | 	absorbing boundary coefficient 
float   | conv_error | = |   | 	final convergence error 
bool    | ctr | =n |  [y/n] | 	cross-correlation time-reversal imaging 
float   | dr | =acpar->dx |   | 	receiver interval 
float   | ds | = |   | 	shot interval 
float   | f0 | = |   | 	reference frequency 
float   | fhi | =0.5/acpar->dt |   | 	
float   | flo | =0. |   | 	
int     | function | =2 |   | 	if 1, forward modeling; if 2, FWI; if 3, RTM 
int     | grectx | =3 |   | 	gradient smoothing radius in x 
int     | grectz | =3 |   | 	gradient smoothing radius in z 
float   | hard | =0.1f |   | 	hard thresholding for masking 
bool    | hidesrc | =n |  [y/n] | 	hide source footprint in fwi 
int     | interval | =1 |   | 	wavefield storing interval 
bool    | inv | =n |  [y/n] | 	inversion flag 
int     | media | =1 |   | 	if 1, acoustic media; if 2, visco-acoustic media 
int     | nb | =100 |   | 	boundary width 
int     | ngrp | =1 |   | 	number of sub-groups of receivers 
int     | niter | = |   | 	iteration number 
int     | nitersrc | =1 |   | 	num of iter'ns for source inversion 
int     | nr | =acpar->nx |   | 	number of receiver 
int     | ns | = |   | 	shot number 
int     | ompchunk | =1 |   | 	
int     | ompnth | =0 |   | 	
bool    | onlygrad | =n |  [y/n] | 	only want gradident 
bool    | onlysrc | =n |  [y/n] | 	only invert for source (vel known), active when inv=y 
bool    | onlyvel | =n |  [y/n] | 	only invert for vel (source known), active when inv=y 
bool    | oreo | =n |  [y/n] | 	keep oreo or keep cream 
file    | output | = |   | 	auxiliary output file name
float   | perc | =1.0f |   | 	padding percentatge for swnorm 
bool    | precsrc | =n |  [y/n] | 	source inversion preconditioning 
float   | r0 | =acpar->x0 |   | 	receiver origin 
int     | recttsrc | =50 |   | 	source smoothing in t before masking 
int     | rectxsrc | =1 |   | 	source smoothing in x before masking 
int     | rectzsrc | =1 |   | 	source smoothing in z before masking 
int     | repeat | =1 |   | 	repeat resetting alpha 
int     | repeatsrc | =1 |   | 	source smoothing repeatation times 
int     | rz | =1 |   | 	receiver depth 
float   | s0 | = |   | 	shot origin 
int     | size | =0 |   | 	sliding window radius 
bool    | sw | =n |  [y/n] | 	sliding window normalization 
int     | sz | =5 |   | 	source depth 
float   | v1 | =0. |   | 	
float   | v2 | =10. |   | 	
bool    | verb | =n |  [y/n] | 	verbosity flag 
int     | waterz | =0 |   | 	water layer depth 
int     | waterzb | =0 |   | 	water layer depth from bottom up 
float   | woff1 | =acpar->r0 |   | 	
float   | woff2 | =acpar->r0+(acpar->nr-1)*acpar->dr |   | 	
float   | wt1 | =acpar->t0 |   | 	
float   | wt2 | =acpar->t0+(acpar->nt-1)*acpar->dt |   | 	
[DIRECTORY]
user/chenyk
