#
#   prepare the mask
#
#grdmath UCERF.grd 6390 MUL = UCERF_strain.grd
grdmath UCERF.grd 2400 MUL = UCERF_strain.grd
grdmath SEISM.grd 19800 MUL = SEISM_strain.grd
#
# mask topo < -2000 m
#
grdmath topo.grd -2000 GT 0. NAN = tmp1.grd
#
#  mask the parkfield area
#
grdmask parkfield_mask.xy -F -R-125/-114/31.0/42.0 -I.00833333 -N1/1/NaN -Gtmp2.grd 
grdmath tmp1.grd tmp2.grd OR = tmp.grd
#
#   prepare the strain data for a correlation analysis
#
foreach name (`cat stemlist_all`)
echo $name\_strain.grd
grdmath tmp.grd $name\_strain.grd OR = mask.grd
cp mask.grd tmp.grd
end
#
#   prepare the strain data for a correlation analysis
#
foreach name (`cat stemlist_all`)
echo $name\_strain.grd
grdmath $name\_strain.grd 0 AND mask.grd OR =  $name\_cut.grd
grd2xyz $name\_cut.grd -S -Z > $name.dat
end
#rm tmp*.grd