def vrefl(vel, par):
geom.zero2d(vel + "zero", "", par)
Flow(
vel,
vel + "zero",
"""
spike nsp=1 mag=-%g k1=%d l1=%d |
add add=%g
"""
% (par["vel"], (par["zrefl"] - par["oz"]) / par["dz"] + 1, par["nz"], par["vel"]),
)
Plot(vel, wplot.igrey2d("mean=y", par))
def vtran(vel, par):
geom.zero2d(vel + "zero", "", par)
Flow(vel, vel + "zero", 'math output="%(vel)g"' % par)
Plot(vel, wplot.igrey2d("mean=y", par))
