[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[GEMPAK #SRJ-938159]: trouble using gdpvsf

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.

  • Subject: [GEMPAK #SRJ-938159]: trouble using gdpvsf
  • Date: Tue, 20 Nov 2012 08:48:03 -0700
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