[gradsusr] Hello everyone

avi kojokro kojokroavi at gmail.com
Sun Aug 21 13:02:34 EDT 2016


I have this script (which is attached to this letter) that I use to enter
the three servers
Is there a way to modify this script, so I'll write the date and the latitude
and longitude points once at the beginning of the script?


Thanks

Avi.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20160821/27c95e4b/attachment.html 
-------------- next part --------------




'reinit'
'c'
'set display color white'
'c'
'sdfopen http://www.esrl.noaa.gov/psd/thredds/dodsC/Datasets/ncep.reanalysis.derived/pressure/air.mon.mean.nc'

'set lat 32.0'
'set lon 35.0'

'set time 00z01aug2016'
'q time'
say ' 'result
'set dignum 2'

'set lev 500'

*'set gxout print'
'd air'
tt1=subwrd(result,4)
tt=tt' 'tt1
say ' 'T500 result 

'set lev 700'

*'set gxout print'
'd air'
tt1=subwrd(result,4)
tt=tt' 'tt1
say ' 'T700 result 

'set lev 850'

*'set gxout print'
'd air'
tt1=subwrd(result,4)
tt=tt' 'tt1
say ' 'T850 result

'set lev 1000'

*'set gxout print'
'd air'
tt1=subwrd(result,4)
tt=tt' 'tt1
say ' 'T2m result



'reinit'
'c'
'set display color white'


'sdfopen http://www.esrl.noaa.gov/psd/thredds/dodsC/Datasets/ncep.reanalysis.derived/pressure/hgt.mon.mean.nc'

'set lat 32.0'
'set lon 35.0'

'set time 00z01aug2016'

'q time'
say ' 'result

'set lev 500'

*'set gxout print'
'd hgt'
tt1=subwrd(result,4)
tt=tt' 'tt1
say ' 'H500 result

'set lev 700'

*'set gxout print'
'd hgt'
tt1=subwrd(result,4)
tt=tt' 'tt1
say ' 'H700 result

'set lev 850'


*'set gxout print'
'd hgt'
tt1=subwrd(result,4)
tt=tt' 'tt1
say ' 'H850 result


*'set gxout print'
'd hgt(lev=500)-hgt(lev=1000)'
tt1=subwrd(result,4)
tt=tt' 'tt1
say ' '1000-500 result



'reinit'
'c'
'set display color white'


'sdfopen http://www.esrl.noaa.gov/psd/thredds/dodsC/Datasets/ncep.reanalysis.derived/surface/slp.mon.mean.nc'

'set lat 32.0'
'set lon 35.0'

'set time 00z01aug2016'


'q time'
say ' 'result

'set lev 0'

*'set gxout print'
'd slp'
tt1=subwrd(result,4)
tt=tt' 'tt1
say ' 'SLP result










More information about the gradsusr mailing list