This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.
Hi Dan, I fear than GDPVSF has not been maintained in quite some time, and it having problems due to changes in the some grid routines made around release 6.0. I am looking at it now and hopefully can get you a working version to test, and include in the next Unidata release. Thanks for sending the grid file for me to try. Michael > The goal is to produce plots of theta and pressure on a constant > potential vorticity surface (1.5 PVU in the example below). > > I have a grid (ftp://ftp1.esrl.noaa.gov/users/dgottas/Gempak/test.grd) > that contains data on both constant-pressure (native) and > potential-temperature (produced with gdvint) surfaces. The listing of > these grids, as produced by gdinfo, is attached below in TestGridInfo.txt. > > When attempting to produce grids of pressure, potential temperature, and > wind on the 1.5 PVU surface using the following settings for gdpvsf, > > GDFILE = test.grd > GDOUTF = test.grd > GFUNC = abs(pvor(pres,obs)) > GDATTIM = last > GVCORD = thta > STARTL = 400 > STOPL = 200 > DESIRE = 0.0000015 > GDOUTL = 15 > GVOUTC = PVAB > GPACK > GLIST = pres;thta;urel;vrel > PMAX = 700 > > I receive the sequential errors/warnings "DG_INIT has not been called", > "No functional grids could be found.", and "Grid size is too large.", > and no new data/grids are written to the output file. A detailed listing > of this output is contained in the attached file GdpvsfOutput.txt. > > Also, tried to produce the desired grids using the native pressure grids > as input: > > GDFILE = test.grd > GDOUTF = test.grd > GFUNC = abs(pvor(thta,obs)) > GDATTIM = last > GVCORD = pres > STARTL = 100 > STOPL = 1000 > DESIRE = 0.0000015 > GDOUTL = 15 > GVOUTC = PVAB > GPACK > GLIST = pres;thta;urel;vrel > PMAX = 700 > > and received the same errors/warnings. > > Any insight or guidance you can provide would be greatly appreciated. > > Thanks, > > Dan > > Ticket Details =================== Ticket ID: SRJ-938159 Department: Support GEMPAK Priority: Critical Status: Open