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