[gradsusr] ammended script
Eric Altshuler
ela at cola.iges.org
Thu Mar 14 18:22:31 EDT 2013
Hi kaikombo,
Here are a few script changes you should make:
1. remove the lines 'set x 1 144' and 'set y 1 73'
2. change 'set lon 0 360' to 'set lon -180 180'
3. move your sdfopen commands (uwnd and vwnd) outside the while loop--it's unnecessary to open these files every time through the loop
4. remove the 'reinit' command that's inside the while loop
5. remove 'disable sdfwrite', this is an invalid command
Best regards,
Eric L. Altshuler
Research Scientist
Center for Ocean-Land-Atmosphere Studies
4041 Powder Mill Road, Suite 302
Calverton, MD 20705-3106
USA
E-mail: ela at cola.iges.org
Phone: (301) 902-1257
Fax: (301) 595-9793
----- Original Message -----
From: "kombo kai" <kaikombo at yahoo.com>
To: gradsusr at gradsusr.org
Sent: Wednesday, March 13, 2013 12:39:15 AM
Subject: [gradsusr] ammended script
I have write the filed for the vorticity at each level using the NCEP u and v winds.
but when I plot these files at the lon and lat of my interest it leaves a verticle strip at lon 0
any one can help on I how to remove this strip?.
I have attached the image which I have plotted.
this is my 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
'set x 1 144'
'set y 1 73'
'vort=hcurl(uwnd,vwnd.2)'
'set sdfwrite d:/vorticity2/'lev'_vorty.nc'
'sdfwrite vort'
'reinit'
levcount=levcount+1
endwhile
'disable sdfwrite'
'reinit'
sdfopen d:/vorticity2/850_vorty.nc
set lon -70 120
set lat -45 45
set t 30
set gxout shaded
d vort
plsease find the attaced image which I get
regards
kaikombo
_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org
http://gradsusr.org/mailman/listinfo/gradsusr
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130314/6921288f/attachment-0003.html
More information about the gradsusr
mailing list