david Sandwell

July 22, 2011

These are grd-files of E-N-U velocity from a maxwell model produced in mid 2009.
The *tot.grd files of east and north have a splined velocity model added to the maxwell model.  
The spline model is done using the following script.

#
#   sample the grids of north and east velocity
#
# GPS file there too:  lon lat ve se vn sn vu su
#
#  correct the mean differences
#
rm junk*
rm GPS_resid_east.dat GPS_resid_north.dat
awk '{print ($1-360.,$2,$3+15.7,$5,$4-17.3,$6) }' < window_sort_edit_ven.dat > junk0.dat
#
#  extract the model from the grid
#
grdtrack junk0.dat -Gve.grd > junk1.dat
grdtrack junk1.dat -Gvn.grd > junk2.dat
awk '{print ($1,$2,$3-$7,1./$4)}' < junk2.dat | blockmedian -R-124/-114.5/32/39 -I.4 -F -Wi | awk '{print ($1, $2, $3, "1")}'  >  east.dat
awk '{print ($1,$2,$5-$8,1./$6)}' < junk2.dat | blockmedian -R-124/-114.5/32/39 -I.4 -F -Wi | awk '{print ($1, $2, $3, "1")}'  >  north.dat
echo `wc -l east.dat | awk '{print $1, "0"}' ` > head.dat
cat head.dat east.dat > GPS_resid_east.dat
cat head.dat north.dat > GPS_resid_north.dat
#
#  now grid the residuals and add back to the velocity
#
./surf2d GPS_resid_east.dat ve_resid_surf2d.grd .9 -125. -113.333 1167 32 40.336 834
./surf2d GPS_resid_north.dat vn_resid_surf2d.grd .9 -125. -113.333 1167 32 40.336 834
grdmath ve.grd ve_resid_surf2d.grd ADD = ve_tot.grd
grdmath vn.grd vn_resid_surf2d.grd ADD = vn_tot.grd
#