#shell to run rfgen_asc and to compare its output with the observed set 
#of R and T receiver functions
#inputs - $1 max freq,  $2 baz, $3 model name and $4 phase velocity, for #modeling
# $5,$6 - names of RRF and TRF 

/bin/rm .gmtcommands
gmtset ANOT_FONT_SIZE  12
gmtset LABEL_FONT_SIZE 16

#first - run the simulation
#take max frequency, baz, model name and phase velocity as input
#
echo $1 > runfile
echo $2 >> runfile
echo $3 >> runfile
echo $4 >> runfile


rfsyn < runfile > /dev/null

/bin/mv RRF.table $3.$1.$2.$4.rrftab
/bin/mv TRF.table $3.$1.$2.$4.trftab

echo syn done
#now make plots
set BOX = -R-2/15/-1.5/3
set FRAME = -JX2/2
set dirname = `pwd`

#RADIAL
psxy $FRAME $BOX -G0/255/0 -K -N -M -Y6 -X2 -P << END > datasynmod.ps
-2 0 
30 0
>
0 -2
0 2
>
5 -2
5 2
>
10 -2
10 2
>
15 -2
15 2 
END


set scaled = ` awk '$1 > -1 && $1 < 1 {print $2} ' $5 | sort -n | tail -1`

awk '$1 > -2.1 && $1 < 20.1 {print $1,0,$2}' $5 | pswiggle $FRAME $BOX -Z$scaled -W3p/0/0/255 -N -K -O >> datasynmod.ps   

set scales = `awk '$1 > -1 && $1 < 1 {print $3}' $3.$1.$2.$4.rrftab | sort -n | tail -1`

awk '$1 < 20.1 {print $1, 0, $3}' $3.$1.$2.$4.rrftab | pswiggle $FRAME $BOX -Z$scales -W1p/255/0/0 -N -O -K >> datasynmod.ps


pstext $FRAME $BOX -O -K -N << END >> datasynmod.ps
15 0.5 18 0 5 6 R
10 6 14 0 5 6 $dirname
END

echo rad done

#TRANS
psxy $FRAME $BOX -G0/255/0 -K -O -N -M -Y-0.8 -Ba5f1:"time, s":\S << END >> datasynmod.ps
-2 0 
30 0
>
0 -2
0 2
>
5 -2
5 2
>
10 -2
10 2
>
15 -2
15 2 
END


awk '$1 > -2.1 && $1 < 20.1 {print $1,0,$2}' $6 | pswiggle $FRAME $BOX -Z$scaled -W3p/0/0/255 -N -K -O >> datasynmod.ps   
awk '$1 < 20.1 {print $1, 0, $3}' $3.$1.$2.$4.trftab | pswiggle $FRAME $BOX -Z$scales -W/255/0/0 -N -O -K >> datasynmod.ps


pstext $FRAME $BOX -O -N -K << END >> datasynmod.ps
15 4.5 18 0 5 6 $3 $1 $2 $4
15 0.5 18 0 5 6 T
END

#now plot the model, S velocity and E value vs depth

/bin/cp $3 tmpmod
jp2tab tmpmod
awk '{print $1,$5}' tmpmod.tab |psxy -JX2/1.3 -R0/100/1.5/5/ -W2p -O -K -Y-2.5 -Ba25f5:"depth, km":/a0.5f0.1:"S vel, km/s":/SE >> datasynmod.ps

#awk '{print $1,$2/$5}' tmpmod.tab |psxy -JX2/1.3 -R0/100/1/2.5/ -W1/255/0/0 -O -K #-Ba25f5:"depth, km":/a0.5f0.1:"Vp/Vs ratio":/SE >> datasynmod.ps
#NOTE you can plot other values from the model, like Vp/Vs ratio

awk '{print $1,$6*100}' tmpmod.tab |psxy -JX2/1.3 -R0/100/-20/20  -W1pto -O  -Ba25f5/a10f2:"S anisotropy, %":/WN >> datasynmod.ps