# crop the grid 

gmt grdcut Vtot.grd -GVcut.grd -R500/1500/200/1900 
gmt grdcut Utot.grd -GUcut.grd -R500/1500/200/1900 

gmt grd2xyz Vcut.grd > v.dat 
gmt grd2xyz Ucut.grd > u.dat 

paste u.dat v.dat > uv.dat 
awk '{print $2-1050,$1-400-500,$6,$3,0,0,0}' uv.dat > inp.dat 

trans_pole -2 2 285.6 50.1 276 54 < inp.dat > out.dat 

awk '{print $8-360,$9,$10}' < out.dat > ve.dat 
awk '{print $8-360,$9,$11}' < out.dat > vn.dat 

gmt surface ve.dat -Gve.grd -R-125/-113.333/32/40.335 -I.01/0.01
gmt surface vn.dat -Gvn.grd -R-125/-113.333/32/40.335 -I.01/0.01 

gmt makecpt -Crainbow -T-15/15/1 -D -Z > vn.cpt
gmt makecpt -Crainbow -T-15/15/1 -D -Z > ve.cpt

./plot_Lat_Long.com -125 -114 31.8 40.5 0.65 vn 5
./plot_Lat_Long.com -125 -114 31.8 40.5 0.65 ve 5

# clean up
rm v.dat u.dat uv.dat inp.dat out.dat ve.dat vn.dat 
#rm *.cpt 
rm Ucut.grd Vcut.grd 

gmt ps2raster ve_ll.ps -A -Tj
gmt ps2raster vn_ll.ps -A -Tj