Synopsis
		sfawefd3d < file_wav.rsf vel=file_vel.rsf sou=file_src.rsf rec=file_rec.rsf > file_dat.rsf wfl=file_wfl.rsf den=file_den.rsf verb=n snap=n expl=n dabc=n cden=n adj=n fsrf=n optfd=n fdorder=4 hybridbc=n sinc=n jsnap=nt jdata=1 nqz=sf_n(az) nqx=sf_n(ax) nqy=sf_n(ay) oqz=sf_o(az) oqx=sf_o(ax) oqy=sf_o(ay) dqz=sf_d(az) dqx=sf_d(ax) dqy=sf_d(ay)

2Nth order in space, 2nd order in time
with optimized FD scheme option and hybrid one-way ABC option
adj flag indicates backwards source injection, not exact adjoint propagator

Parameters
		bool adj=n [y/n] adjoint flag   bool cden=n [y/n] Constant density   bool dabc=n [y/n] Absorbing BC   string den= auxiliary input file name   float dqx=sf_d(ax) Saved wfld window dx   float dqy=sf_d(ay) Saved wfld window dy   float dqz=sf_d(az) Saved wfld window dz   bool expl=n [y/n] Multiple sources, one wvlt   int fdorder=4 spatial FD order   bool fsrf=n [y/n] Free surface flag   bool hybridbc=n [y/n] hybrid Absorbing BC   int jdata=1 # of t steps at which to save receiver data   int jsnap=nt # of t steps at which to save wavefield   int nqx=sf_n(ax) Saved wfld window nx   int nqy=sf_n(ay) Saved wfld window ny   int nqz=sf_n(az) Saved wfld window nz   bool optfd=n [y/n] optimized FD coefficients flag   float oqx=sf_o(ax) Saved wfld window ox   float oqy=sf_o(ay) Saved wfld window oy   float oqz=sf_o(az) Saved wfld window oz   file rec= auxiliary input file name   bool sinc=n [y/n] sinc source injection   bool snap=n [y/n] Wavefield snapshots flag   file sou= auxiliary input file name   file vel= auxiliary input file name   bool verb=n [y/n] Verbosity flag   file wfl= auxiliary output file name

Used In
		CWP          geo2011WideAzimuthAngleDecomposition/flatEICangle