tslow.v:
	Window3d < 3000marm_vel_smooth1.H min1=1000 max1=8000 | Transp |\
	Grey newclip=1 color=j bclip=0.0001 eclip=0.000666667 title="(a)" label1="z" barlabel="slowness" \
	label2="x" screenratio=4 screenht=3.5 screenwd=5 wheretitle=b polarity= -1 \
	wantaxis1=y wantaxis2=y labelsz=3 titlesz=3 wantscalebar=1 barwidth=0.1 > /dev/null out=${R}/$@

bslow.v:
	Window3d < 3000marm_vel_smooth.ed.H min1=1000 max1=8000 | Transp |\
	Grey newclip=1 color=j bclip=0.0001 eclip=0.000666667 title="(b)" label1="z" barlabel="slowness" \
	label2="x" screenratio=4 screenht=3.5 screenwd=5 wheretitle=b polarity= -1 \
	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wantscalebar=1 barwidth=0.1 > /dev/null out=${R}/$@

dslow1.v:
	Window3d < 3000marm.Dslow.H min1=1000 max1=8000 min3=1500 | Transp |\
	Grey newclip=1 color=j bpclip=0 epclip=100 title="(a)" label1="z" barlabel="slowness difference" \
	label2="x" screenratio=4 screenht=3.5 screenwd=5 wheretitle=b \
	wantaxis1=y wantaxis2=y labelsz=3 titlesz=3 wantscalebar=1 barwidth=0.1 > /dev/null out=${R}/$@

dslow2.v:
	Window3d < 3000marm_vel_smooth1.H min1=-8400 n1=1076 min3=1500 squeeze=n > x1.H
	Window3d < 3000marm_vel_smooth.ed.H min1=-8400 n1=1076 min3=1500 squeeze=n > x2.H
	Math file1=x1.H file2=x2.H exp="file1/file2"  | Window3d min1=1000 max1=8000| Transp |\
	Grey newclip=1 color=j bclip=0.7 eclip=1.3 bias=1. allpos=1 title="(b)" label1="z" barlabel="slowness ratio" \
	label2="x" screenratio=4 screenht=3.5 screenwd=5 wheretitle=b \
	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wantscalebar=1 barwidth=0.1 > /dev/null out=${R}/$@
	
islow.v: tslow.v bslow.v
	vppen < ${R}/tslow.v xcenter=0 ycenter=-3. vpstyle=n > junk1.v
	vppen < ${R}/bslow.v xcenter=0 ycenter=-0. vpstyle=n > junk2.v
	vp_Overlay junk1.v junk2.v > ${R}/$@
	
dslow.v: dslow1.v dslow2.v
	vppen < ${R}/dslow1.v xcenter=0 ycenter=-3. vpstyle=n > junk1.v
	vppen < ${R}/dslow2.v xcenter=0 ycenter=-0. vpstyle=n > junk2.v
	vp_Overlay junk1.v junk2.v > ${R}/$@

biconv.v:
	Window3d < 3000marm.inv_bkim.conv.H min2=0 n2=1 min1=1000 max1=8000 n4=1 min3=1500 | Transp |\
	Grey pclip=95 title="(a)" label1="z" label2="x" screenratio=4 screenht=3.5 screenwd=4 \
	wantaxis1=y wantaxis2=y labelsz=3 titlesz=3 wheretitle=b gainpanel=a > /dev/null out=1.v
	Window3d < 3000marm.inv_bkim.conv.H n1=1 min1=5500 n4=1 min3=1500 min2=-192 n2=17 | Transp | Grey gainpanel=a pclip=95 title="  " \
	label1="z" label2="h" screenratio=4 screenht=3.5 screenwd=1 \
	wantaxis1=n wantaxis2=y labelsz=3 o2num=-192 d2num=192 n2tic=2 > /dev/null out=2.v
	vppen < 1.v xcenter=0 ycenter=-3. vpstyle=n > junk1.v
	vppen < 2.v xcenter=-3.5 ycenter=-3. vpstyle=n > junk2.v
	vp_Overlay junk1.v junk2.v > ${R}/$@

bwconv.v:
	Window3d < 3000marm.window3.H min2=0 n2=1 min1=1000 max1=8000 min3=1500 | Transp plane=23 |Stack|Transp |\
	Grey gainpanel=a title="(b)" label1="z" label2="x" screenratio=4 screenht=3.5 screenwd=4 pclip=95 \
	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wheretitle=b > /dev/null out=1.v
	Window3d < 3000marm.window3.H n1=1 min1=5500 min3=1500 min2=-192 n2=17 | Transp plane=23|Stack| Transp|Grey gainpanel=a pclip=95 title="  " \
	label1="z" label2="h" screenratio=4 screenht=3.5 screenwd=1 \
	wantaxis1=n wantaxis2=n labelsz=3 > /dev/null out=2.v
	vppen < 1.v xcenter=-0 ycenter=-0. vpstyle=n > junk1.v
	vppen < 2.v xcenter=-3.5 ycenter=-0. vpstyle=n > junk2.v
	vp_Overlay junk1.v junk2.v > ${R}/$@

bimg1.v: biconv.v bwconv.v 
	vp_Overlay ${R}/biconv.v ${R}/bwconv.v > ${R}/$@

bimg1.pdf: bimg1.v
	pstexpen ${R}/bimg1.v ${R}/bimg1.ps color=n fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/bimg1.ps

wrconv.v:
	Window3d < 3000marm.rot.rec.sd3.H min2=0 n2=1 min3=1000 max3=8000 min1=1500 | Transp plane=23 |Stack|\
	Grey gainpanel=a title="(a)" label1="z" label2="x" screenratio=4 screenht=3.5 screenwd=4 pclip=100 \
	wantaxis1=y wantaxis2=y labelsz=3 titlesz=3 wheretitle=b > /dev/null out=1.v
	Window3d < 3000marm.rot.rec.sd3.H n3=1 min3=5500 min1=1500 min2=-192 n2=17 | Transp plane=23 |Stack|Grey gainpanel=a pclip=100 title="  " \
	label1="z" label2="h" screenratio=4 screenht=3.5 screenwd=1 \
	wantaxis1=n wantaxis2=y labelsz=3 o2num=-192 d2num=192 n2tic=2 > /dev/null out=2.v
	vppen < 1.v xcenter=-0 ycenter=-3 vpstyle=n > junk1.v
	vppen < 2.v xcenter=-3.5 ycenter=-3 vpstyle=n > junk2.v
	vp_Overlay junk1.v junk2.v > ${R}/$@

wsconv.v:
	Window3d < 3000marm.rot.sou.sd3.H min2=0 n2=1 min3=1000 max3=8000 min1=1500 |Transp plane=23 |Stack| \
	Grey gainpanel=a title="(b)" label1="z" label2="x" screenratio=4 screenht=3.5 screenwd=4 pclip=100 \
	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wheretitle=b > /dev/null out=1.v
	Window3d < 3000marm.rot.sou.sd3.H n3=1 min3=5500 min1=1500 min2=-192 n2=17 |Transp plane=23 |Stack| Grey gainpanel=a pclip=100 title="  " \
	label1="z" label2="h" screenratio=4 screenht=3.5 screenwd=1 \
	wantaxis1=n wantaxis2=n labelsz=3 > /dev/null out=2.v
	vppen < 1.v xcenter=-0 ycenter=-0. vpstyle=n > junk1.v
	vppen < 2.v xcenter=-3.5 ycenter=-0 vpstyle=n > junk2.v
	vp_Overlay junk1.v junk2.v > ${R}/$@

diconv.v:
	Window3d < 3000marm.Dimag.conv.H min2=0 n2=1 min1=1000 max1=8000 min3=1500 | Transp |\
	Grey title="(a)" label1="z" label2="x" screenratio=4 screenht=3.5 screenwd=4 pclip=100 \
	wantaxis1=y wantaxis2=y labelsz=3 titlesz=3 wheretitle=b > /dev/null out=1.v
	Window3d < 3000marm.Dimag.conv.H n1=1 min1=5500 min3=1500 min2=-192 n2=17| Transp | Grey title="  " pclip=100 \
	label1="z" label2="h" screenratio=4 screenht=3.5 screenwd=1 \
	wantaxis1=n wantaxis2=y labelsz=3 o2num=-192 d2num=192 n2tic=2 > /dev/null out=2.v
	vppen < 1.v xcenter=-0 ycenter=-6.5 vpstyle=n > junk1.v
	vppen < 2.v xcenter=-3.5 ycenter=-6.5 vpstyle=n > junk2.v
	vp_Overlay junk1.v junk2.v > ${R}/$@
	
diperm11.v:
	Window3d < 30001marm.Dimag11.perm.H min2=0 n2=1 min1=1000 max1=8000 min3=1500 | Transp |\
	Grey title="(b)" label1="z" label2="x" screenratio=4 screenht=3.5 screenwd=4 pclip=100 \
	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wheretitle=b > /dev/null out=1.v
	Window3d < 30001marm.Dimag11.perm.H n1=1 min1=5500 min3=1500 min2=-192 n2=17| Transp | Grey title="  " pclip=100 \
	label1="z" label2="h" screenratio=4 screenht=3.5 screenwd=1 \
	wantaxis1=n wantaxis2=n labelsz=3 > /dev/null out=2.v
	vppen < 1.v xcenter=-0. ycenter=-3.5 vpstyle=n > junk1.v
	vppen < 2.v xcenter=-3.5 ycenter=-3.5 vpstyle=n > junk2.v
	vp_Overlay junk1.v junk2.v > ${R}/$@
	
diperm35.v:
	Window3d < 30002marm.Dimag35.perm.H min2=0 n2=1 min1=1000 max1=8000 min3=1500 | Transp |\
	Grey title="(c)" label1="z" label2="x" screenratio=4 screenht=3.5 screenwd=4 pclip=100 \
	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wheretitle=b > /dev/null out=1.v
	Window3d < 30002marm.Dimag35.perm.H n1=1 min1=5500 min3=1500 min2=-192 n2=17| Transp | Grey title="  " pclip=100 \
	label1="z" label2="h" screenratio=4 screenht=3.5 screenwd=1 \
	wantaxis1=n wantaxis2=n labelsz=3 > /dev/null out=2.v
	vppen < 1.v xcenter=-0 ycenter=-.5 vpstyle=n > junk1.v
	vppen < 2.v xcenter=-3.5 ycenter=-.5 vpstyle=n > junk2.v
	vp_Overlay junk1.v junk2.v > ${R}/$@

rimg1.v: wrconv.v wsconv.v 
	vp_Overlay ${R}/wrconv.v ${R}/wsconv.v > ${R}/$@

dimg1.v: diconv.v diperm11.v diperm35.v # diperm70.v
	vp_Overlay ${R}/diconv.v ${R}/diperm11.v ${R}/diperm35.v > ${R}/$@

dsconv.v:
	Window3d < 3000marm.Dslow.conv.adjoint.H min1=1000 max1=8000 min3=1500 | Scale scale_to=1|Transp |\
	Grey newclip=1 color=j bpclip=0 epclip=100 title="(a)" label1="z" barlabel=" "  \
	label2="x" screenratio=4 screenht=3.5 screenwd=5 wheretitle=b \
	wantaxis1=y wantaxis2=y labelsz=3 titlesz=3 wantscalebar=1 barwidth=0.1 > /dev/null out=${R}/$@

dsperm11.v:
	Window3d < 30001marm.Dslow11.perm.adjoint.H min1=1000 max1=8000 min3=1500 | Scale scale_to=1|Transp |\
	Grey newclip=1 color=j bpclip=0 epclip=100 title="(b)" label1="z" barlabel=" " \
	label2="x" screenratio=4 screenht=3.5 screenwd=5 wheretitle=b \
	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wantscalebar=1 barwidth=0.1 > /dev/null out=${R}/$@

dsperm35.v:
	Window3d < 30002marm.Dslow35.perm.adjoint.H min1=1000 max1=8000 min3=1500 | Scale scale_to=1|Transp |\
	Grey newclip=1 color=j bpclip=0 epclip=100 title="(c)" label1="z" barlabel=" " \
	label2="x" screenratio=4 screenht=3.5 screenwd=5 wheretitle=b \
	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wantscalebar=1 barwidth=0.1 > /dev/null out=${R}/$@

dsperm70.v:
	Window3d < 30001marm.Dslow35.perm.adjoint.H min1=1000 max1=8000 min3=1500 | Scale scale_to=1|Transp |\
	Grey newclip=1 color=j bpclip=0 epclip=100 title="(c)" label1="z" barlabel=" " \
	label2="x" screenratio=4 screenht=3.5 screenwd=5 wheretitle=b \
	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wantscalebar=1 barwidth=0.1 > /dev/null out=${R}/$@

dsadj.v: dsconv.v dsperm11.v dsperm35.v
	vppen < ${R}/dsconv.v xcenter=0 ycenter=-6.5 vpstyle=n > junk1.v
	vppen < ${R}/dsperm11.v xcenter=0 ycenter=-3.25 vpstyle=n > junk2.v
	vppen < ${R}/dsperm35.v xcenter=0 ycenter=0. vpstyle=n > junk3.v
	vp_Overlay junk1.v junk2.v junk3.v > ${R}/$@	

#dsoconv.v:
#	Window3d < 3000marm.inv_grad.conv.H min1=1000 max1=8000 min3=1500 | Transp |\
#	Grey newclip=1 color=j bpclip=0 epclip=100 title="(a)" label1="z" barlabel="  " \
#	label2="x" screenratio=4 screenht=3.5 screenwd=5 wheretitle=b \
#	wantaxis1=y wantaxis2=y labelsz=3 titlesz=3 wantscalebar=1 barwidth=0.1 > /dev/null out=${R}/$@
#
#dsoperm11.v:
#	Window3d < 3000marm.inv_grad11.perm.H min1=1000 max1=8000 min3=1500 | Transp |\
#	Grey newclip=1 color=j bpclip=0 epclip=100 title="(b)" label1="z" barlabel="  " \
#	label2="x" screenratio=4 screenht=3.5 screenwd=5 wheretitle=b \
#	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wantscalebar=1 barwidth=0.1 > /dev/null out=${R}/$@
#
#dsoperm31.v:
#	Window3d < 3000marm.inv_grad1.perm.H min1=1000 max1=8000 min3=1500 | Transp |\
#	Grey newclip=1 color=j bpclip=0 epclip=100 title="(c)" label1="z" barlabel="  " \
#	label2="x" screenratio=4 screenht=3.5 screenwd=5 wheretitle=b \
#	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wantscalebar=1 barwidth=0.1 > /dev/null out=${R}/$@
#
#dsdso.v: dsoconv.v dsoperm11.v dsoperm31.v
#	vppen < ${R}/dsoconv.v xcenter=0 ycenter=-6.5 vpstyle=n > junk1.v
#	vppen < ${R}/dsoperm11.v xcenter=0 ycenter=-3. vpstyle=n > junk2.v
#	vppen < ${R}/dsoperm31.v xcenter=0 ycenter=0.5. vpstyle=n > junk3.v
#	vp_Overlay junk1.v junk2.v junk3.v > ${R}/$@	

sfpermtrue.v:
	< 3000marm_vel_smooth1.H Window3d min1=1000 max1=8000 min3=1500 n4=1 f4=15 | Transp |\
	Grey newclip=1 color=j bclip=0.0001 eclip=0.000666667 title="(a)" label1="z" barlabel="  " \
	label2="x" screenratio=4 screenht=2.5 screenwd=5 wheretitle=b \
	wantaxis1=y wantaxis2=y labelsz=3 titlesz=3 wantscalebar=1 barwidth=0.1 > /dev/null out=${R}/$@

sfperm11.v:
	Math file1=30003marm.inv_vout11.perm.sd.H exp="1/file1" | Window3d min1=1000 max1=8000 min3=1500 n4=1 f4=26 | Transp |\
	Grey newclip=1 color=j bclip=0.0001 eclip=0.000666667 title="(b)" label1="z" barlabel="  " \
	label2="x" screenratio=4 screenht=2.5 screenwd=5 wheretitle=b \
	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wantscalebar=1 barwidth=0.1 > /dev/null out=${R}/$@

sfperm35.v: # only reflector encode 35 shots
	Math file1=30003marm.inv_vout35.perm.sd.H exp="1/file1" | Window3d min1=1000 max1=8000 min3=1500 n4=1 f4=26 | Transp |\
	Grey newclip=1 color=j bclip=0.0001 eclip=0.000666667 title="(c)" label1="z" barlabel="  " \
	label2="x" screenratio=4 screenht=2.5 screenwd=5 wheretitle=b \
	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wantscalebar=1 barwidth=0.1 > /dev/null out=${R}/$@
	
sfperm70.v: # only reflector encode 70 shots
	Math file1=3000marm.inv_vout35.perm.sd.H exp="1/file1" | Window3d min1=1000 max1=8000 min3=1500 n4=1 f4=6 | Transp |\
	Grey newclip=1 color=j bclip=0.0001 eclip=0.000666667 title="(d)" label1="z" barlabel="  " \
	label2="x" screenratio=4 screenht=2.5 screenwd=5 wheretitle=b \
	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wantscalebar=1 barwidth=0.1 > /dev/null out=${R}/$@

sfperm.v: sfpermtrue.v sfperm11.v sfperm35.v sfperm70.v
	vppen < ${R}/sfpermtrue.v xcenter=0 ycenter=-7.75 vpstyle=n > junk0.v
	vppen < ${R}/sfperm11.v xcenter=0 ycenter=-5.5 vpstyle=n > junk1.v
	vppen < ${R}/sfperm35.v xcenter=0 ycenter=-3.25 vpstyle=n > junk2.v
	vppen < ${R}/sfperm70.v xcenter=0 ycenter=-1. vpstyle=n > junk3.v
	vp_Overlay junk0.v junk1.v junk2.v junk3.v > ${R}/$@	

dsfperm11.v:
	Window3d < 3000marm_vel_smooth1.H min1=1000 n1=291 squeeze=n min3=1500 > x1.H
	Math file1=30003marm.inv_vout11.perm.sd.H exp="1/file1" | Window3d min1=1000 max1=8000 min3=1500 n4=1 f4=26 squeeze=n > x2.H
	Math file1=x1.H file2=x2.H exp="file1/file2"  | Window3d | Transp |\
	Grey newclip=1 color=j bclip=0.7 eclip=1.3 bias=1. allpos=1 minval= 0.7 maxval= 1.3 title="(a)" label1="z" barlabel="slowness ratio" \
	label2="x" screenratio=4 screenht=3.5 screenwd=5 wheretitle=b \
	wantaxis1=y wantaxis2=y labelsz=3 titlesz=3 wantscalebar=1 barwidth=0.1 > /dev/null out=${R}/$@

dsfperm35.v:
	Window3d < 3000marm_vel_smooth1.H min1=1000 n1=291  squeeze=n min3=1500 > x1.H
	Math file1=30003marm.inv_vout35.perm.sd.H exp="1/file1" | Window3d min1=1000 max1=8000 min3=1500 n4=1 f4=26 squeeze=n > x2.H
	Math file1=x1.H file2=x2.H exp="file1/file2"  | Window3d | Transp |\
	Grey newclip=1 color=j bclip=0.7 eclip=1.3 bias=1. allpos=1 minval= 0.7 maxval= 1.3 title="(b)" label1="z" barlabel="slowness ratio" \
	label2="x" screenratio=4 screenht=3.5 screenwd=5 wheretitle=b \
	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wantscalebar=1 barwidth=0.1 > /dev/null out=${R}/$@

dsfperm70.v:
	Window3d < 3000marm_vel_smooth1.H min1=1000 n1=291  squeeze=n min3=1500 > x1.H
	Math file1=3000marm.inv_vout35.perm.sd.H exp="1/file1" | Window3d min1=1000 max1=8000 min3=1500 n4=1 f4=6 squeeze=n > x2.H
	Math file1=x1.H file2=x2.H exp="file1/file2"  | Window3d | Transp |\
	Grey newclip=1 color=j bclip=0.7 eclip=1.3 bias=1. allpos=1 minval= 0.7 maxval= 1.3 title="(c)" label1="z" barlabel="slowness ratio" \
	label2="x" screenratio=4 screenht=3.5 screenwd=5 wheretitle=b \
	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wantscalebar=1 barwidth=0.1 > /dev/null out=${R}/$@

dsfperm.v: dsfperm11.v dsfperm35.v dsfperm70.v
	vppen < ${R}/dsfperm11.v xcenter=0 ycenter=-6.5 vpstyle=n > junk1.v
	vppen < ${R}/dsfperm35.v xcenter=0 ycenter=-3.25 vpstyle=n > junk2.v
	vppen < ${R}/dsfperm70.v xcenter=0 ycenter=0. vpstyle=n > junk3.v
	vp_Overlay junk1.v junk2.v junk3.v > ${R}/$@	

hsfpermbkg.v:
	Window3d < 3000marm_vel_smooth1.H min1=1000 n1=291 squeeze=n min3=1500 > x1.H
	< 3000marm_vel_smooth.ed.H  Window3d min1=1000  min3=1500 n1=291 squeeze=n > x2.H
	Math file1=x1.H file2=x2.H exp="file1/file2"  | Window3d min3=2400 | Transp |\
	Histogram ninterv=25 dinterv=0.01 min=0.8 max=1.2 | Graph title="  " legend= y labelsz=6 \
	legendtitle="Background slowness" legendloc=tr legendpos=11.5,8 legendsz=3 plotfat=5 axiscol=1 plotcol=1 \
	label2="Count" label1="Slowness ratio" wheretitle=b min1=0.8 max1=1.2 dash=0 max2=5000 \
	wantaxis1=y labelsz=6 titlesz=6 screenratio=1.5> /dev/null out=${R}/$@

hsfperm11.v:
	Window3d < 3000marm_vel_smooth1.H min1=1000 n1=291 squeeze=n min3=1500 > x1.H
	Math file1=30003marm.inv_vout11.perm.sd.H exp="1/file1"|Window3d min1=1000 max1=8000 min3=1500 n4=1 f4=26 squeeze=n>x2.H
	Math file1=x1.H file2=x2.H exp="file1/file2"  | Window3d min3=2400 | Transp |\
	Histogram ninterv=25 dinterv=0.01 min=0.8 max=1.2 | Graph title=" " legend= y axiscol=1 plotcol=1 \
	legendtitle="Optim.slowness(11)" legendpos=12,7.5 legendsz=3 plotfat=5 \
	label2=" " label1="  " wheretitle=b min1=0.8 max1=1.2 dash=1 max2=5000 \
	wantaxis1=y labelsz=6 titlesz=6 screenratio=1.5> /dev/null out=${R}/$@

hsfperm35.v:
	Window3d < 3000marm_vel_smooth1.H min1=1000 n1=291  squeeze=n min3=1500 > x1.H
	Math file1=30003marm.inv_vout35.perm.sd.H exp="1/file1"|Window3d min1=1000 max1=8000 min3=1500 n4=1 f4=26 squeeze=n>x2.H
	Math file1=x1.H file2=x2.H exp="file1/file2"  | Window3d min3=2400 | Transp |\
	Histogram ninterv=25 dinterv=0.01 min=0.8 max=1.2 | Graph title=" " legend= y axiscol=1 plotcol=1 \
	legendtitle="Optim.slowness(35)" legendpos=12,7 legendsz=3 plotfat=5 \
	label2=" " label1="  " wheretitle=b min1=0.8 max1=1.2 dash=2 max2=5000 \
	wantaxis1=y labelsz=6 titlesz=6 screenratio=1.5 > /dev/null out=${R}/$@

hsfperm70.v:
	Window3d < 3000marm_vel_smooth1.H min1=1000 n1=291  squeeze=n min3=1500 > x1.H
	Math file1=3000marm.inv_vout35.perm.sd.H exp="1/file1"|Window3d min1=1000 max1=8000 min3=1500 n4=1 f4=6 squeeze=n>x2.H
	Math file1=x1.H file2=x2.H exp="file1/file2"  | Window3d min3=2400 | Transp |\
	Histogram ninterv=25 dinterv=0.01 min=0.8 max=1.2 | Graph title=" " legend= y axiscol=1 plotcol=1 \
	legendtitle="Optim.slowness(70)" legendpos=12,6.5 legendsz=3 plotfat=5 \
	label2=" " label1="  " wheretitle=b min1=0.8 max1=1.2 dash=6 max2=5000 \
	wantaxis1=y labelsz=6 titlesz=6 screenratio=1.5 > /dev/null out=${R}/$@

hsfperm.v: hsfpermbkg.v hsfperm11.v hsfperm35.v hsfperm70.v
	vppen < ${R}/hsfpermbkg.v xcenter=0 ycenter=0 vpstyle=n > junk1.v
	vppen < ${R}/hsfperm11.v xcenter=0 ycenter=0 pstyle=n > junk2.v
	vppen < ${R}/hsfperm35.v xcenter=0 ycenter=0 vpstyle=n > junk3.v
	vppen < ${R}/hsfperm70.v xcenter=0 ycenter=0. vpstyle=n > junk4.v
	vp_Overlay junk1.v junk2.v junk3.v junk4.v > junk5.v
	vp_annotate batch=y text=hist_anno < junk5.v > ${R}/$@


################
fipermback.v:
	Window3d < 3000marm.inv_bkim.conv.H min2=0 n2=1 min1=1000 max1=8000 min3=1500| Transp |\
	Grey title="(a)" label1="z" label2="x" screenratio=4 screenht=2 screenwd=6 pclip=95 \
	wantaxis1=y wantaxis2=y labelsz=3 titlesz=3 wheretitle=b  > /dev/null out=1.v
	vppen <1.v xcenter=-1. ycenter=-7 vpstyle=n > ${R}/$@

fipermtrue.v:
	Window3d < 3000marm.bkim.vtrue.H min2=0 n2=1 min1=1000 max1=8000 min3=1500| Transp |\
	Grey title="(b)" label1="z" label2="x" screenratio=4 screenht=2 screenwd=6 pclip=95 \
	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wheretitle=b  > /dev/null out=1.v
	vppen <1.v xcenter=-1. ycenter=-5.3 vpstyle=n > ${R}/$@

fiperm11.v:
	Window3d < 3000marm.bkim.v11a.H min2=0 n2=1 min1=1000 max1=8000 min3=1500| Transp |\
	Grey title="(c)" label1="z" label2="x" screenratio=4 screenht=2 screenwd=6 pclip=95 \
	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wheretitle=b  > /dev/null out=1.v
	vppen <1.v xcenter=-1. ycenter=-3.6 vpstyle=n > ${R}/$@

fiperm35.v:
	Window3d < 3000marm.bkim.v35a.H min2=0 n2=1 min1=1000 max1=8000 min3=1500 | Transp |\
	Grey title="(d)" label1="z" label2="x" screenratio=4 screenht=2 screenwd=6 pclip=95 \
	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wheretitle=b > /dev/null out=1.v
	vppen < 1.v xcenter=-1. ycenter=-1.9 vpstyle=n > ${R}/$@

fiperm70.v:
	Window3d < 3000marm.bkim.v70.H min2=0 n2=1 min1=1000 max1=8000 min3=1500 | Transp |\
	Grey title="(e)" label1="z" label2="x" screenratio=4 screenht=2 screenwd=6 pclip=95 \
	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wheretitle=b  > /dev/null out=1.v
	vppen <1.v xcenter=-1. ycenter=-0.2 vpstyle=n > ${R}/$@

fimg1.v: fipermback.v fipermtrue.v fiperm11.v fiperm35.v fiperm70.v
	vp_Overlay ${R}/fipermback.v ${R}/fipermtrue.v ${R}/fiperm11.v ${R}/fiperm35.v ${R}/fiperm70.v > ${R}/$@

fimg1.pdf: fimg1.v
	pstexpen ${R}/fimg1.v ${R}/fimg1.ps color=y fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/fimg1.ps


fangback.v:
	Window3d < 3000marm.ang.vback.H n3=1 min3=1500 min1=1500|\
	Grey title="1500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=y wantaxis2=y labelsz=2 > /dev/null out=0.v
	Window3d < 3000marm.ang.vback.H n3=1 min3=2500 min1=1500|\
	Grey title="3500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=1.v
	Window3d < 3000marm.ang.vback.H n3=1 min3=3500 min1=1500|\
	Grey title="5500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=2.v
	Window3d < 3000marm.ang.vback.H n3=1 min3=4500 min1=1500|\
	Grey title="1500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=3.v
	Window3d < 3000marm.ang.vback.H n3=1 min3=5500 min1=1500|\
	Grey title="3500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=4.v
	Window3d < 3000marm.ang.vback.H n3=1 min3=6500 min1=1500|\
	Grey title="5500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=5.v
	Window3d < 3000marm.ang.vback.H n3=1 min3=7500 min1=1500|\
	Grey title="7500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=6.v
	vppen < 0.v xcenter=-0. ycenter=-7 vpstyle=n > junk0.v
	vppen < 1.v xcenter=-1. ycenter=-7 vpstyle=n > junk1.v
	vppen < 2.v xcenter=-2. ycenter=-7 vpstyle=n > junk2.v
	vppen < 3.v xcenter=-3. ycenter=-7 vpstyle=n > junk3.v
	vppen < 4.v xcenter=-4. ycenter=-7 vpstyle=n > junk4.v
	vppen < 5.v xcenter=-5. ycenter=-7 vpstyle=n > junk5.v
	vppen < 6.v xcenter=-6. ycenter=-7 vpstyle=n > junk6.v
	vp_Overlay junk0.v junk1.v junk2.v junk3.v junk4.v junk5.v junk6.v > $@	
	rm -f [0-7].v junk[0-7].v
	pstexpen fangback.v fangback.ps color=y fat=1 fatmult=1.5 invras=n force=  invras=y 

fangtrue.v:
	Window3d < 3000marm.ang.vtrue.H n3=1 min3=1500 min1=1500|\
	Grey title="1500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=y wantaxis2=y labelsz=2 > /dev/null out=0.v
	Window3d < 3000marm.ang.vtrue.H n3=1 min3=2500 min1=1500|\
	Grey title="3500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=1.v
	Window3d < 3000marm.ang.vtrue.H n3=1 min3=3500 min1=1500|\
	Grey title="5500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=2.v
	Window3d < 3000marm.ang.vtrue.H n3=1 min3=4500 min1=1500|\
	Grey title="1500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=3.v
	Window3d < 3000marm.ang.vtrue.H n3=1 min3=5500 min1=1500|\
	Grey title="3500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=4.v
	Window3d < 3000marm.ang.vtrue.H n3=1 min3=6500 min1=1500|\
	Grey title="5500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=5.v
	Window3d < 3000marm.ang.vtrue.H n3=1 min3=7500 min1=1500|\
	Grey title="7500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=6.v
	vppen < 0.v xcenter=-0. ycenter=-7 vpstyle=n > junk0.v
	vppen < 1.v xcenter=-1. ycenter=-7 vpstyle=n > junk1.v
	vppen < 2.v xcenter=-2. ycenter=-7 vpstyle=n > junk2.v
	vppen < 3.v xcenter=-3. ycenter=-7 vpstyle=n > junk3.v
	vppen < 4.v xcenter=-4. ycenter=-7 vpstyle=n > junk4.v
	vppen < 5.v xcenter=-5. ycenter=-7 vpstyle=n > junk5.v
	vppen < 6.v xcenter=-6. ycenter=-7 vpstyle=n > junk6.v
	vp_Overlay junk0.v junk1.v junk2.v junk3.v junk4.v junk5.v junk6.v > $@	
	rm -f [0-7].v junk[0-7].v
	pstexpen fangtrue.v fangtrue.ps color=y fat=1 fatmult=1.5 invras=n force=  invras=y 

fang11.v:
	Window3d < 3000marm.ang.v11.H n3=1 min3=1500 min1=1500|\
	Grey title="1500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=y wantaxis2=y labelsz=2 > /dev/null out=0.v
	Window3d < 3000marm.ang.v11.H n3=1 min3=2500 min1=1500|\
	Grey title="3500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=1.v
	Window3d < 3000marm.ang.v11.H n3=1 min3=3500 min1=1500|\
	Grey title="5500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=2.v
	Window3d < 3000marm.ang.v11.H n3=1 min3=4500 min1=1500|\
	Grey title="1500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=3.v
	Window3d < 3000marm.ang.v11.H n3=1 min3=5500 min1=1500|\
	Grey title="3500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=4.v
	Window3d < 3000marm.ang.v11.H n3=1 min3=6500 min1=1500|\
	Grey title="5500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=5.v
	Window3d < 3000marm.ang.v11.H n3=1 min3=7500 min1=1500|\
	Grey title="7500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=6.v
	vppen < 0.v xcenter=-0. ycenter=-7 vpstyle=n > junk0.v
	vppen < 1.v xcenter=-1. ycenter=-7 vpstyle=n > junk1.v
	vppen < 2.v xcenter=-2. ycenter=-7 vpstyle=n > junk2.v
	vppen < 3.v xcenter=-3. ycenter=-7 vpstyle=n > junk3.v
	vppen < 4.v xcenter=-4. ycenter=-7 vpstyle=n > junk4.v
	vppen < 5.v xcenter=-5. ycenter=-7 vpstyle=n > junk5.v
	vppen < 6.v xcenter=-6. ycenter=-7 vpstyle=n > junk6.v
	vp_Overlay junk0.v junk1.v junk2.v junk3.v junk4.v junk5.v junk6.v > $@	
	rm -f [0-7].v junk[0-7].v
	pstexpen fang11.v fang11.ps color=y fat=1 fatmult=1.5 invras=n force=  invras=y 

fang35.v:
	Window3d < 3000marm.ang.v35.H n3=1 min3=1500 min1=1500|\
	Grey title="1500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=y wantaxis2=y labelsz=2 > /dev/null out=0.v
	Window3d < 3000marm.ang.v35.H n3=1 min3=2500 min1=1500|\
	Grey title="3500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=1.v
	Window3d < 3000marm.ang.v35.H n3=1 min3=3500 min1=1500|\
	Grey title="5500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=2.v
	Window3d < 3000marm.ang.v35.H n3=1 min3=4500 min1=1500|\
	Grey title="1500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=3.v
	Window3d < 3000marm.ang.v35.H n3=1 min3=5500 min1=1500|\
	Grey title="3500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=4.v
	Window3d < 3000marm.ang.v35.H n3=1 min3=6500 min1=1500|\
	Grey title="5500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=5.v
	Window3d < 3000marm.ang.v35.H n3=1 min3=7500 min1=1500|\
	Grey title="7500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=6.v
	vppen < 0.v xcenter=-0. ycenter=-7 vpstyle=n > junk0.v
	vppen < 1.v xcenter=-1. ycenter=-7 vpstyle=n > junk1.v
	vppen < 2.v xcenter=-2. ycenter=-7 vpstyle=n > junk2.v
	vppen < 3.v xcenter=-3. ycenter=-7 vpstyle=n > junk3.v
	vppen < 4.v xcenter=-4. ycenter=-7 vpstyle=n > junk4.v
	vppen < 5.v xcenter=-5. ycenter=-7 vpstyle=n > junk5.v
	vppen < 6.v xcenter=-6. ycenter=-7 vpstyle=n > junk6.v
	vp_Overlay junk0.v junk1.v junk2.v junk3.v junk4.v junk5.v junk6.v > ${R}/$@	
	rm -f [0-7].v junk[0-7].v
	pstexpen ${R}/fang35.v ${R}/fang35.ps color=y fat=1 fatmult=1.5 invras=n force=  invras=y 

fang70.v:
	Window3d < 3000marm.ang.v70.H n3=1 min3=1500 min1=1500|\
	Grey title="1500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=y wantaxis2=y labelsz=2 > /dev/null out=0.v
	Window3d < 3000marm.ang.v70.H n3=1 min3=2500 min1=1500|\
	Grey title="3500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=1.v
	Window3d < 3000marm.ang.v70.H n3=1 min3=3500 min1=1500|\
	Grey title="5500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=2.v
	Window3d < 3000marm.ang.v70.H n3=1 min3=4500 min1=1500|\
	Grey title="1500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=3.v
	Window3d < 3000marm.ang.v70.H n3=1 min3=5500 min1=1500|\
	Grey title="3500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=4.v
	Window3d < 3000marm.ang.v70.H n3=1 min3=6500 min1=1500|\
	Grey title="5500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=5.v
	Window3d < 3000marm.ang.v70.H n3=1 min3=7500 min1=1500|\
	Grey title="7500" pclip=95 label1="z" titlesz=3 o2num=-60 d2num=60 n2tic=2 \
	label2="angle" screenratio=4 screenht=2 screenwd=1 wheretitle=t grid2=1 gridcol=1 axiscol=1 \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=6.v
	vppen < 0.v xcenter=-0. ycenter=-7 vpstyle=n > junk0.v
	vppen < 1.v xcenter=-1. ycenter=-7 vpstyle=n > junk1.v
	vppen < 2.v xcenter=-2. ycenter=-7 vpstyle=n > junk2.v
	vppen < 3.v xcenter=-3. ycenter=-7 vpstyle=n > junk3.v
	vppen < 4.v xcenter=-4. ycenter=-7 vpstyle=n > junk4.v
	vppen < 5.v xcenter=-5. ycenter=-7 vpstyle=n > junk5.v
	vppen < 6.v xcenter=-6. ycenter=-7 vpstyle=n > junk6.v
	vp_Overlay junk0.v junk1.v junk2.v junk3.v junk4.v junk5.v junk6.v > ${R}/$@	
	rm -f [0-7].v junk[0-7].v
	pstexpen ${R}/fang70.v ${R}/fang70.ps color=y fat=1 fatmult=1.5 invras=n force=  invras=y