Hi Kathryn,
Having had a look at the input GEMPAK file, it appears that it exceeds
the maximum # of vertical levels (MAXLVL=65) that gdvint accepts. I
would assume that recompiling the gdvint program, with
$GEMPAK/source/programs/gd/gdvint/vicmn.cmn edited so that MAXLVL is
larger (or use LLMXLV, set to 500 in $GEMPAK/include/GEMPRM.PRM), would
be a safe road to take.
--Kevin
______________________________________________________________________
Kevin Tyle, Systems Administrator **********************
Dept. of Earth & Atmospheric Sciences ktyle@xxxxxxxxxxxxxxxx
University at Albany, ES-235 518-442-4578 (voice)
1400 Washington Avenue 518-442-5825 (fax)
Albany, NY 12222 **********************
______________________________________________________________________
Kathryn Shontz wrote:
Hello Gembuds-
I have been using gempak files (converted from netcdf using wrf2gem)
to plot wind, vorticity, helicity and precipitable water with no
problem. For purposes of comparison to other work in the field, I
wanted to convert my pressure files to zagl.
Not a problem with gdvint it appeared... until I tried to do heights
exceeding 8200m. My .gem files were created in pressure coordinates,
starting at 1000mb to 200mb with a dP of 10mb. I can make plots of
350mb (~8400m) with the original files but for some reason, gdvint
wont let any interpolation happen beyond this point. I tried to make
the change in height greater than my initial 100m dZ step, but even at
1000m dZ the outcome was the same... Please see the messages I got below:
Has anyone encountered a similar problem and if so, how was it
resolved? THANK YOU SO MUCH FOR THE HELP!
Kathryn Shontz
University of California, Los Angeles
Department of Atmospheric and Oceanic Science
Graduate Student
GEMPAK-GDVINT>l
GDFILE = gordon_850wind_0913_12.gem
GDOUTF = test200m.gem
GDATTIM = 060913/1200f000
GVCORD = pres/zagl
GLEVEL = 6600-12000-200
MAXGRD = 10000
GAREA = grid
VCOORD = none
GEMPAK-GDVINT>glev=7000-
12000-500
GEMPAK-GDVINT>gdoutf=test500.gem
GEMPAK-GDVINT>r
Too many input levels, top level reset to 370.
PRES found.
TMPK found.
HGHT found.
UREL found.
VREL found.
W found.
SPFH found.
TMPK found.
TMPK not interpolated to level = 8500
TMPK not interpolated to level = 9000
TMPK not interpolated to level = 9500
TMPK not interpolated to level = 10000
TMPK not interpolated to level = 10500
TMPK not interpolated to level = 11000
TMPK not interpolated to level = 11500
TMPK not interpolated to level = 12000
UREL not interpolated to level = 8500
UREL not interpolated to level = 9000
UREL not interpolated to level = 9500
UREL not interpolated to level = 10000
UREL not interpolated to level = 10500
UREL not interpolated to level = 11000
UREL not interpolated to level = 11500
UREL not interpolated to level = 12000
VREL not interpolated to level = 8500
VREL not interpolated to level = 9000
VREL not interpolated to level = 9500
VREL not interpolated to level = 10000
VREL not interpolated to level = 10500
VREL not interpolated to level = 11000
VREL not interpolated to level = 11500
VREL not interpolated to level = 12000
W not interpolated to level = 8500
W not interpolated to level = 9000
W not interpolated to level = 9500
W not interpolated to level = 10000
W not interpolated to level = 10500
W not interpolated to level = 11000
W not interpolated to level = 11500
W not interpolated to level = 12000
SPFH not interpolated to level = 8500
SPFH not interpolated to level = 9000
SPFH not interpolated to level = 9500
SPFH not interpolated to level = 10000
SPFH not interpolated to level = 10500
SPFH not interpolated to level = 11000
SPFH not interpolated to level = 11500
SPFH not interpolated to level = 12000
Interpolation is done. No errors.
Parameters requested: GDFILE,GDOUTF,GDATTIM,GVCORD,GLEVEL,MAXGRD,GAREA,
VCOORD.
------------------------------------------------------------------------
_______________________________________________
gembud mailing list
gembud@xxxxxxxxxxxxxxxx
For list information or to unsubscribe, visit: http://www.unidata.ucar.edu/mailing_lists/
