NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
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 > >