[gradsusr] writing vort.nc file

kombo kai kaikombo at yahoo.com
Tue Mar 12 11:40:20 EDT 2013


Hi grad user
I have the u and v wind components from the NCEP/NCAR reanalysis 1 data sets
I am writing the vorticity files (vorty.nc) at each level, 
I have tried difeent  options but the results of the files puts the data only to the 0 120 and not -70 0 I think there is problem with my script
please any one with suggestion or modification of my scripr so as to work properly as I need
here is the script

'reinit'
'c'
lev.1=925
lev.2=850
lev.3=700
lev.4=600
lev.5=500
lev.6=400
lev.7=300
lev.8=250
lev.9=200
lev.10=150
lev.11=100
levcount=1
while(levcount<=11)
lev=lev.levcount
say lev
'sdfopen d:/vorticity/uwnd.mon.mean.nc'
'sdfopen d:/vorticity/vwnd.mon.mean.nc'
'set lon 0 360'
'set lat -90 90'
'set t 1 778'
'set lev 'lev
'vort=hcurl(uwnd,vwnd.2)'
'set sdfwrite d:/vorticity2/'lev'_vorty.nc'
'sdfwrite vort'
'reinit' 
levcount=levcount+1
endwhile
'disable sdfwrite'
'reinit'
regards
kaikombo


-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130312/6594f064/attachment-0003.html 


More information about the gradsusr mailing list