[BUG]

[pplotn]

    # sfgpufwi is not working with shots in hh format. it outputs nan values in vsnaps

cmd((f"sfgenshots < {'vel.rsf'} csdgather=n fm={3} amp=1 dt={dt} ns={ns} ng={ng} nt={nt} "
f"sxbeg={sxbeg} szbeg=2 jsx={jsx} jsz=0 gxbeg={gxbeg} gzbeg={szbeg} jgx={jgx} jgz=0 --out=stdout> shots.rsf"))
cmd((f"sfgpufwi < {'smvel.rsf'} shots={'shots.rsf'} > "
f"{'vsnaps_.rsf'} grads={'grads_.rsf'} objs={'objs_.rsf'} illums={'illums_.rsf'} "
f" niter={str(Niter)} precon=n verb=y rbell=10 "))