<div dir="ltr">Really thankful. It worked. Attached herewith the working script for someones reference in future may be.<div><br></div><div><div>'reinit'</div><div> 'sdfopen E:/Onset_Paper/<a href="http://vort_2015.nc" target="_blank">vort_2015.nc</a>'</div><div> 'set t 1'</div><div> 'set z 4'</div><div>* ########## Coriolis calculation</div><div> 'define cor=2*7.29*(1/100000)*sin(lat*3.1415/180)'</div><div>* ########## Relative Vorticity calculation</div><div> 'define vor=hcurl(u,v)'</div><div>* ########## Potential Temperature (theta) calculation</div><div>'define dt=t(z-1)*pow(1000/lev(z-1),0.286)-t(z+1)*pow(1000/lev(z+1),0.286)'</div><div>* ########## differential of theta/pressure</div><div>'define dp=100*(lev(z-1)-lev(z+1))'</div><div> 'define dtdp=dt/dp'</div><div>* ########### PV calculation on pressure surface</div><div> 'define pv=-9.8*(cor+vor)*dtdp'</div><div>'set gxout shaded'</div><div> 'd pv'</div><div>'cbarn'</div><div>'draw title Potential vorticity on 1Jan2015'</div><div>'printim Jan2015.jpg white'</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Aug 1, 2016 at 12:59 PM, Ricardo Hallak <span dir="ltr"><<a href="mailto:ricardo.hallak@iag.usp.br" target="_blank">ricardo.hallak@iag.usp.br</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="font-family:times new roman,new york,times,serif;font-size:12pt;color:#000000">Hi,<br><br>not tested, however my guess is that you should use <br>lat instead y<br>when defining:<br>'define coriol=2*7.29e-5*sin(y*3.1415/180)'<br><br>Other commands in the script seem ok.<br><br>Hope it works,<br>Ricardo<br><br><br><hr><blockquote style="border-left:2px solid rgb(16,16,255);margin-left:5px;padding-left:5px;color:#000;font-weight:normal;font-style:normal;text-decoration:none;font-family:Helvetica,Arial,sans-serif;font-size:12pt"><div><div><div dir="ltr">Any scripts for potential vorticity computation for the layer between surface to 600 hPa. I have zonal, meridional wind and air temperature at different levels to compute the same. I tried with the scripts that are already therein gradsusr archives. Not seems to work. One such script is attached herewith. Any help in this regard will be much appreciated. Thanks. <br clear="all"><div><br></div>-- <br><div><div dir="ltr"><div style="font-size:12.8px">Sincerely</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">**************************************************</div><div style="font-size:12.8px">Dr. Jayakrishnan P.R</div><div style="font-size:12.8px">Post Doctoral Research Fellow</div><div style="font-size:12.8px">C308, Institute of Ocean and Earth Sciences (IOES)</div><div style="font-size:12.8px">University of Malaya</div><div style="font-size:12.8px">50603, Kuala Lumpur</div><div style="font-size:12.8px">Tel: <a href="tel:%28%2B603%29%207967%204640" value="+60379674640" target="_blank">(+603) 7967 4640</a> / 6995</div><div style="font-size:12.8px">Mobile: 012-3526094</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">****************************************************</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><span><font size="1"><br></font></span></div><div style="font-size:12.8px"><span><font size="1"><br></font></span></div><div style="font-size:12.8px"><br></div></div></div>
</div>
<br></div></div>_______________________________________________<br>gradsusr mailing list<br><a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br><a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><span name="x"></span><br></blockquote></div></div><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div data-smartmail="gmail_signature"><div dir="ltr"><div style="font-size:12.8px">Sincerely</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">**************************************************</div><div style="font-size:12.8px">Dr. Jayakrishnan P.R</div><div style="font-size:12.8px">Post Doctoral Research Fellow</div><div style="font-size:12.8px">C308, Institute of Ocean and Earth Sciences (IOES)</div><div style="font-size:12.8px">University of Malaya</div><div style="font-size:12.8px">50603, Kuala Lumpur</div><div style="font-size:12.8px">Tel: <a href="tel:%28%2B603%29%207967%204640" value="+60379674640" target="_blank">(+603) 7967 4640</a> / 6995</div><div style="font-size:12.8px">Mobile: 012-3526094</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">****************************************************</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><span style="color:rgb(7,55,99)"><font size="1"><br></font></span></div><div style="font-size:12.8px"><span style="color:rgb(7,55,99)"><font size="1"><br></font></span></div><div style="font-size:12.8px"><br></div></div></div>
</div></div>