b'\n \n \n
 
\n  
sfavvvdwe2d (4.0)		index
user/fperrone/Mavvvdwe2d.c
\n 2D acoustic variable-velocity variable-density time-domain FD modeling\n

\n \n \n \n \n \n
 
\n Synopsis
       sfavvvdwe2d < Fwav.rsf vel=Fvel.rsf den=Fden.rsf sou=Fsou.rsf rec=Frec.rsf > Fdat.rsf wfl=Fwfl.rsf verb=n free=n dabc=n adj=n snap=y nb=NOP jsnap=1 dpt=n

\nThe code uses a standard second-order stencil in time.
\nThe coefficients of the spatial stencil are computed
\nby matching the transfer function of the 6-point discretized
\nfirst-derivative operator to the ideal response.
\n
\nThe code implements the linearized operator obtained from the
\nsystem of first-order PDEs parametrized in incompressibility and density
\n
\ndv/dt = - 1./rho * grad(p)
\ndp/dt = - K * div(v)
\n
\nwhere
\nrho : density
\nK : incompressibility
\ndiv : divergence operator
\ngrad : gradient operator
\np,v : pressure and particle velocity wavefields
\n
\nThe models supplied by the user are wave speed and density, the code performs
\nthe conversion internally to buoyancy (inverse density) and incompressibility.
\n
\nAuthor: Francesco Perrone
\nDate: February 2020
\n\n

\n \n \n \n \n \n
 
\n Parameters
       \n \n \n
\n  
bool adj=n [y/n]		\tAdjointness
\n
\n \n\n \n \n
\n  
bool dabc=n [y/n]		\tAbsorbing BC
\n
\n \n\n \n \n
\n  
file den= 		\tauxiliary input file name
\n
\n \n\n \n \n
\n  
bool dpt=n [y/n]		\trun dot product test
\n
\n \n\n \n \n
\n  
bool free=n [y/n]		\tFree surface
\n
\n \n\n \n \n
\n  
int jsnap=1 		\tundersampling factor for the wavefields
\n
\n \n\n \n \n
\n  
int nb=NOP 		\tthickness of the absorbing boundary: NOP is the width of the FD stencil
\n
\n \n\n \n \n
\n  
file rec= 		\tauxiliary input file name
\n
\n \n\n \n \n
\n  
bool snap=y [y/n]		\twavefield snapshots
\n
\n \n\n \n \n
\n  
file sou= 		\tauxiliary input file name
\n
\n \n\n \n \n
\n  
file vel= 		\tauxiliary input file name
\n
\n \n\n \n \n
\n  
bool verb=n [y/n]		\tVerbosity
\n
\n \n\n \n \n
\n  
file wfl= 		\tauxiliary output file name
\n
\n \n
'