In [1]:

    

import os
from subprocess import call
from IPython.core.display import Image


    

In [4]:

    

cd /usgs/data2/etwomey/GOM/Final_Data_Sets/final_xyz/


    

    




/usgs/data2/etwomey/GOM/Final_Data_Sets/final_xyz

In [5]:

    

#!locate navd88.xyz


    

In [6]:

    

dataset={}
dataset['geodas']='GEODAS*navd88.xyz'
dataset['globalss']='global_ss*navd88.xyz'


    

In [7]:

    

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)


    

In [8]:

    

bbox


    

    
Out[8]:





'-R-71.500000/-63.000000/-39.500000/46.000000'

In [9]:

    

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

In [10]:

    

Image(filename='geodas.png')


    

    
Out[10]:

In [12]:

    

for dset, dfiles in dataset.iteritems():
    print dfiles


    

    




GEODAS*navd88.xyz
global_ss*navd88.xyz

In [ ]:

    

cat /usgs/data2/etwomey/GOM/Final_Data_Sets/final_xyz/geodas/