# Re: Static stability

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Mon Jun 10 2013 - 07:49:53 MDT

Not sure what to say. As noted for a number of functions:
;-------------------------------------------------------------------
;----------------------- Unsupported: DJS --------------------------
;-------------------------------------------------------------------
undef("StaticStability")
function StaticStability (p:numeric, t:numeric, sopt)
;
; Compute Static Stability on constant pressure levels
; Bluestein: Synoptic-Dynamic Meteorology in Midlatitudes
; pg 197
; s = -T*d[log(theta)]/dp = -(T/theta)*d(theta)/dp
;
; Nomenclature
; p - pressure levels (Pa) [1D]
; t - temperature (K) [3D or 4D]
; (lev,lat,lon) , (time,lev,lat,lon)
; sopt - options
; sopt=0 use ===> -(T/theta)*d(theta)/dp
; sopt=1 use ===> -T*d[log(theta)]/dp
; same answer should be returned: mainly used for debug
;
; Note: In a statically stable atmosphere: d(theta)/dp<0, hence s>0
++++++++++++++++++++++++++++++++++++++++++++++++++++++++

As noted, the function is for data on constant pressure levels
(isobaric). You must interpolate the ARW WRF model data to the pressure
levels. I suggest you contact wrfhelp@ucar.edu to determine how to
do that interpolation.

Good luck

On 6/10/13 7:26 AM, Bogdan Antonescu wrote:
> Hi,
>
> I am trying to calculate the static stability at 500 hPa using the
> StaticStability function from shea_util.ncl.
> So far I have not figured out how to use this function for ARW WRF model
>
> Thank you very much.
>
> Bogdan
>
>
>
> _______________________________________________
> 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
Received on Mon Jun 10 07:49:44 2013

This archive was generated by hypermail 2.1.8 : Tue Jun 11 2013 - 12:03:58 MDT