<html><head><style type='text/css'>p { margin: 0; }</style></head><body><div style='font-family: arial,helvetica,sans-serif; font-size: 12pt; color: #000000'><div>Hi kaikombo,</div><div><br></div>Here are a few script changes you should make:<div><br></div><div>1. remove the lines 'set x 1 144' and 'set y 1 73'</div><div>2. change 'set lon 0 360' to 'set lon -180 180'</div><div>3. move your sdfopen commands (uwnd and vwnd) outside the while loop--it's unnecessary to open these files every time through the loop</div><div>4. remove the 'reinit' command that's inside the while loop</div><div>5. remove 'disable sdfwrite', this is an invalid command</div><div><br></div><div>Best regards,</div><div><br><div><span name="x"></span>Eric L. Altshuler<br>Research Scientist<br>Center for Ocean-Land-Atmosphere Studies<br>4041 Powder Mill Road, Suite 302<br>Calverton, MD 20705-3106<br>USA<br><br>E-mail: ela@cola.iges.org<br>Phone: (301) 902-1257<br>Fax: (301) 595-9793<span name="x"></span><br></div><br><hr id="zwchr"><b>From: </b>"kombo kai" <kaikombo@yahoo.com><br><b>To: </b>gradsusr@gradsusr.org<br><b>Sent: </b>Wednesday, March 13, 2013 12:39:15 AM<br><b>Subject: </b>[gradsusr] ammended script<br><br><table cellspacing="0" cellpadding="0" border="0"><tbody><tr><td valign="top" style="font: inherit;">I have write the filed for the vorticity at each level using the NCEP u and v winds.<br>but when I plot these files at the lon and lat of my interest it leaves a verticle strip at lon 0<br>any one can help on I how to remove this strip?.<br>I have attached the image which I have plotted.<br>this is my script<br>'reinit'<br>'c'<br>lev.1=925<br>lev.2=850<br>lev.3=700<br>lev.4=600<br>lev.5=500<br>lev.6=400<br>lev.7=300<br>lev.8=250<br>lev.9=200<br>lev.10=150<br>lev.11=100<br>levcount=1<br>while(levcount<=11)<br>lev=lev.levcount<br>say lev<br>'sdfopen d:/vorticity/uwnd.mon.mean.nc'<br>'sdfopen d:/vorticity/vwnd.mon.mean.nc'<br>'set lon 0 360'<br>'set lat -90 90'<br>'set t 1 778'<br>'set lev 'lev<br>'set x 1 144'<br>'set y 1 73'<br>'vort=hcurl(uwnd,vwnd.2)'<br>'set sdfwrite d:/vorticity2/'lev'_vorty.nc'<br>'sdfwrite vort'<br>'reinit'
<br>levcount=levcount+1<br>endwhile<br>'disable sdfwrite'<br>'reinit'<br><br>sdfopen d:/vorticity2/850_vorty.nc<br>set lon -70 120<br>set lat -45 45<br>set t 30<br>set gxout shaded<br>d vort<br> plsease find the attaced image which I get<br>regards<br>kaikombo<br><br></td></tr></tbody></table><br>_______________________________________________<br>gradsusr mailing list<br>gradsusr@gradsusr.org<br>http://gradsusr.org/mailman/listinfo/gradsusr<br></div></div></body></html>