20010918: interpolating pressure to height

Kevin,

See gdvint for interpolating between vertical coordinates. Type phelp gdvint.
The help says ZAGL is one of the supported coordinates.

You would have to have the surface pressure and heigh grids of course
(eg PRES @0%none and HGHT@0%none).

Steve Chiswell
unidata User Support


On Tue, 18 Sep 2001, Kevin Polston wrote:

> Hi Steve,
>
> Got a question for you.  I was speaking with a colleague at SPC and he sent me
> some of the scripts they use for severe weather forecasting. Several of the
> scripts use a height coordinate (zagl...not the hght coordinate). He e-mailed 
> me
> these scripts and he said "you have to use GEMPAK to interpolate height
> coordinates from the regular pressure coordinate grids".  How do you do this
> interpolation? I don't have a clue and am hoping you might know what to do. 
> I've
> also attached one of the scripts he sent so you can see for yourself. Thanks.
>
> Kevin Polston
> Senior Instructor, NWSTC
> Kansas City MO 64151
>
> 816-880-9314 ext 273
>
>
>
> GDFILE  eta40
> GDATTIM fall
> GLEVEL  4000!4000!1500:500!1500:500!1500:500!
> GVCORD  zagl
> PANEL   0
> SKIP    0
> SCALE   3
> GDPFUN
> sm9s(quo(sub(temp@2000,temp),2000!sm9s(quo(sub(temp@2000,temp),2000!sm9s(lav(mixr!sm9s(lav(mixr!sm9s(lav(mixr!
>
> TYPE    C/f!C!C!C!C
> CONTUR  2
> CINT    1/6/6!1/7!1/8/11!1/12/15!1/16
> LINE    18/1/1!5/1/1!4/10/3!26/1/2!6/1/2
> FINT    .5/6/10
> FLINE   0;0;0;10;10;19;19;12;12;16;16
> HILO
> HLSYM
> CLRBAR  1
> WIND
> REFVEC
> TITLE   5/-2/~^ 2-4 km AGL lapse rate (C/km)!6/-3/ 500-1500 m mean mixing 
> ratio
> (g/kg)!
> TEXT    1/22/hw
> CLEAR   yes
> GAREA   us
> PROJ    lcc
> MAP     1
> LATLON  0
> DEVICE  xw
> STNPLT
> SATFIL
> RADFIL
> LUTFIL
> STREAM
> POSN
> COLORS
> MARKER
> GRDLBL
> FILTER
>
>