Re: Ertel PV question

In 'shea_util.ncl', there are two unsupported functions
which you may find useful. On is for hybrid coordinates,
the other for pressure coordinates.

---
function PotVortHybrid\
             (p:numeric, u:numeric, v:numeric, t:numeric
             ,lat[*]:numeric, gridType[1]:integer, opt:integer)
;
; Compute Isentropic Potential Vorticity on hybrid levels
; CCM Processor User's Guide: May 1994: page B-19
; Original source P Rasch and B Boville
;
---
undef("PotVortIsobaric")
function PotVortIsobaric\
             (p[*]:numeric, u:numeric, v:numeric, t:numeric
             ,lat[*]:numeric, gridType[1]:integer, opt:integer)
;
; Compute Potential Vorticity on constant pressure levels
; Bluestein: Synoptic-Dynamic Meteorology in Midlatitudes
;            pg 264  Eq 4.5.93
;
---
  You can see the function documentation by
%> less $NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl
---
http://www.meted.ucar.edu/nwp/pcu3/cases/etapbomb/pvort.htm
Good luck
On 5/10/10 8:41 AM, JEFFREY JOSEPH WATERS wrote:
> To whom it may concern,
> I am currently trying to calculate Ertel's Potential Vorticity given
> output information from an NetCDf file. I have been reading some of the
> threads pertaining to PV, but I had two questions:
>
> 1) are there any codes for calculating Ertel's PV
> 2) if not, then are there codes for calculating relative vorticity on a
> constant theta surface (relative vorticity is one component of Ertel's
> PV equation)
>
> Any feedback would be much appreciated. Thank you.
>
> Jeff Waters
> Graduate Student
> Department of Meteorology
> Penn State University '11
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk