import os
from subprocess import call
from IPython.core.display import Image
cd /usgs/data2/etwomey/GOM/Final_Data_Sets/final_xyz/
/usgs/data2/etwomey/GOM/Final_Data_Sets/final_xyz
#!locate navd88.xyz
dataset={}
dataset['geodas']='GEODAS*navd88.xyz'
dataset['globalss']='global_ss*navd88.xyz'
len(dataset)
lon1=-71.5
lon2=-63.0
lat1=-39.5
lat2=46.0
box='%f/%f/%f/%f' % (lon1,lon2,lat1,lat2)
bbox='-R%f/%f/%f/%f' % (lon1,lon2,lat1,lat2)
bbox
'-R-71.500000/-63.000000/-39.500000/46.000000'
for dset, dfiles in dataset.iteritems():
# cmd = 'pscoast -JM16 %s -Dh -G120 -B1g0.5 -P -K -V > points.ps' % bbox
# print cmd
# os.system(cmd)
# log = open('points.ps', "w")
# call(['pscoast','-JM16',bbox,'-Dh -G120 -B1g0.5 -P -K -V'],shell=True,stdout=log)
os.system(('pscoast -JM16 -R%s -Dh -G120 -B1g0.5 -P -K -V > points.ps' % box))
cmd='cat %s | psxy -R%s -P -JM16 -O -K -Cgom4c.cpt -Sp -V >> points.ps' % (dfiles,box)
print cmd
# os.system('grdcontour /usgs/data0/bathy/GOM_061311_3sec.grd -JM16 -R%s -Cgom4d.cpt -P -A- -V -O >> points.ps' % box)
# os.system('convert -density 120 -trim points.ps %s.png' % dset)
# Image(filename=('%s.png' % dset))
cat GEODAS*navd88.xyz | psxy -R-71.500000/-63.000000/-39.500000/46.000000 -P -JM16 -O -K -Cgom4c.cpt -Sp -V >> points.ps cat global_ss*navd88.xyz | psxy -R-71.500000/-63.000000/-39.500000/46.000000 -P -JM16 -O -K -Cgom4c.cpt -Sp -V >> points.ps
Image(filename='geodas.png')
for dset, dfiles in dataset.iteritems():
print dfiles
GEODAS*navd88.xyz global_ss*navd88.xyz
cat /usgs/data2/etwomey/GOM/Final_Data_Sets/final_xyz/geodas/