def rr(rr, par):
par["ixrec"] = (4.0 - par["ox"]) / par["dx"]
xrec = par["ox"] + par["ixrec"] * par["dx"]
par["izrec"] = par["shft"]
zrec = par["oz"] + par["izrec"] * par["dz"]
geom.point2d(rr, xrec, zrec, "", par)
Plot(rr, wplot.rrplot2d("plotcol=2 plotfat=20", par))
def ss(ss, par):
par["ixsou"] = (0.0 - par["ox"]) / par["dx"]
xsou = par["ox"] + par["ixsou"] * par["dx"]
par["izsou"] = par["shft"]
zsou = par["oz"] + par["izsou"] * par["dz"]
geom.point2d(ss, xsou, zsou, "", par)
Plot(ss, wplot.rrplot2d("plotcol=1 plotfat=20", par))