[gradsusr] Potential vorticity computation

[Jayakrishnan P.R]

Really thankful. It worked. Attached herewith the working script for
someones reference in future may be.

'reinit'
 'sdfopen E:/Onset_Paper/vort_2015.nc'
   'set t 1'
  'set z 4'
* ########## Coriolis calculation
 'define cor=2*7.29*(1/100000)*sin(lat*3.1415/180)'
* ########## Relative Vorticity calculation
 'define vor=hcurl(u,v)'
* ########## Potential Temperature (theta) calculation
'define dt=t(z-1)*pow(1000/lev(z-1),0.286)-t(z+1)*pow(1000/lev(z+1),0.286)'
* ########## differential of theta/pressure
'define dp=100*(lev(z-1)-lev(z+1))'
 'define dtdp=dt/dp'
* ########### PV calculation on pressure surface
 'define pv=-9.8*(cor+vor)*dtdp'
'set gxout shaded'
 'd pv'
'cbarn'
'draw title Potential vorticity on 1Jan2015'
'printim Jan2015.jpg white'


On Mon, Aug 1, 2016 at 12:59 PM, Ricardo Hallak <ricardo.hallak at iag.usp.br>
wrote:

> Hi,
>
> not tested, however my guess is that you should use
> lat instead y
> when defining:
> 'define coriol=2*7.29e-5*sin(y*3.1415/180)'
>
> Other commands in the script seem ok.
>
> Hope it works,
> Ricardo
>
>
> ------------------------------
>
> 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.
>
> --
> Sincerely
>
> **************************************************
> Dr. Jayakrishnan P.R
> Post Doctoral Research Fellow
> C308, Institute of Ocean and Earth Sciences (IOES)
> University of Malaya
> 50603, Kuala Lumpur
> Tel: (+603) 7967 4640 / 6995
> Mobile: 012-3526094
>
> ****************************************************
>
>
>
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>


-- 
Sincerely

**************************************************
Dr. Jayakrishnan P.R
Post Doctoral Research Fellow
C308, Institute of Ocean and Earth Sciences (IOES)
University of Malaya
50603, Kuala Lumpur
Tel: (+603) 7967 4640 / 6995
Mobile: 012-3526094

****************************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20160801/bba4ab18/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: potvor_5.gs
Type: application/octet-stream
Size: 681 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20160801/bba4ab18/attachment.obj