Hi Vanesa, Thank you for attaching the grid file. I was able to "plot" with GDPLOT2 and GDPLOT3 an empty map, meaning I did not encounter the erorrs GEMPLT -16] NIPBND - Invalid projection bounds. [DG -37] No reference grid navigation is set in GPLT. while running the program. However, the grid you are trying to plot is not found in the grid file you sent me. with GDINFO: GDFILE = gfsanl_4_20121028_0000_000.grd LSTALL = YES OUTPUT = T GDATTIM = 121028/00 GLEVEL = 715 GVCORD = hght GFUNC = TMPC [GDU 2] Did not find any matching grids. Your file DOES contain grids at 715m HGHT, but not for TMPC: GDFILE = gfsanl_4_20121028_0000_000.grd LSTALL = YES OUTPUT = T GDATTIM = all GLEVEL = 715 GVCORD = hght GFUNC = all r NUM TIME1 TIME2 LEVL1 LEVL2 VCORD PARM 41 121028/0000F000 715 HGHT PRES 42 121028/0000F000 715 HGHT UWND 43 121028/0000F000 715 HGHT VWND 44 121028/0000F000 715 HGHT THTA 45 121028/0000F000 715 HGHT OMEG 46 121028/0000F000 715 HGHT RELH I do not see any grids with GVCORD =pvbl, so I'm not sure if the grids created with GDVINT and GDPVSF are what you actually wanted. For example, when I plot THTA values on a 1.5 PV surface, the GDPLOT2 input used is GDFILE =vcord_thta_270-380.grd GDATTIM =f00 GLEVEL =15 GVCORD =pvbl PANEL =1 SKIP =0 SCALE =999 GDPFUN =thta TYPE =f CONTUR =1 CINT =10 LINE =3 FINT = 10/280/380 FLINE = 24;25;26;27;23;22;21;20;19;16;15;14 I've attached a scripts that I used in the past to generate THTA on PV surfaces and plot them. It may help you to take a look at the process and compare with your own to see if and where there are differences. Best, Michael James Unidata Support > > > Hi Michael, > > I was able, with the help of a professor, to get GDPVSF to work by using the > file generated by nagrib2 instead of using gdvint. I am now having an issue > getting this file to create a plot in gdplot3. Below is the output from the > program and attached is the file generated in GDPVSF that I used. I am not > sure if the error is due to the file I am using, or the way that I have > configued the gdplot3 parameters. I am still trying to create a plot of > theta on an isosurface of 1.5 PVU. The glevel set was the gdoutl that I used > in GDPVSF, just to clarify. Any help would be appreciated! > > > Thanks! > > > > GEMPAK-GDPLOT3>d > GDFILE = newgfsanl_4_20121028_0000_000.grd > GDATTIM = 121028/00 > GLEVEL = 715 > GVCORD = hght > PANEL = 0 > SKIP = 0 > SCALE = 999 > GDPFUN = thta > TYPE = C > CONTUR = 0 > CINT = 0 > LINE = 3 > FINT = 0 > FLINE = 10-20 > HILO = > HLSYM = > CLRBAR = > WIND = BM1 > REFVEC = > TITLE = 1 > TEXT = 1 > CLEAR = YES > GAREA = 35;-80;55;-65 > IJSKIP = > PROJ = MER > MAP = 1 > MSCALE = 0 > LATLON = > DEVICE = XW > STNPLT = > SATFIL = > RADFIL = > IMCBAR = > LUTFIL = > STREAM = > POSN = 0 > COLORS = 1 > MARKER = 0 > GRDLBL = 0 > FILTER = YES > BOXLIN = > REGION = VIEW > TXTCOL = > TXTYPE = > TXTFIL = > TXTLOC = > COLUMN = 1 > SHAPE = > INFO = > LOCI = > ANOTLN = > ANOTYP = > GEMPAK-GDPLOT3>proj=CED > GEMPAK-GDPLOT3>r > Number of times = 1 > 1 121028/0000 > > PARAMETERS FOR GDPLOT : FRAME NUMBER 1 INITIAL PLOT SETTINGS - > GDFILE=newgfsanl_4_20121028_0000_000.grd, GDATTIM=121028/0000 > GAREA=35;-80;55;-65, PROJ=CED, MAP=1, LATLON= > PANEL=0, TITLE=1 > OVERLAY TYPE LEVEL VCORD GRID FUNCTION > 1 C 715 HGHT THTA > [GEMPLT -16] NIPBND - Invalid projection bounds. > [DG -37] No reference grid navigation is set in GPLT. > Parameters requested: GDFILE,GDATTIM,GLEVEL,GVCORD,PANEL,SKIP,SCALE,GDPFUN, > TYPE,CONTUR,CINT,LINE,FINT,FLINE,HILO,HLSYM,CLRBAR,WIND,REFVEC,TITLE,TEXT, > CLEAR,GAREA,IJSKIP,PROJ,MAP,MSCALE,LATLON,DEVICE,STNPLT,SATFIL,RADFIL, > IMCBAR,LUTFIL,STREAM,POSN,COLORS,MARKER,GRDLBL,FILTER,BOXLIN,REGION,TXTCOL, > TXTYPE,TXTFIL,TXTLOC,COLUMN,SHAPE,INFO,LOCI,ANOTLN,ANOTYP. > > > > ~Vanesa Urango > > > ----- Original Message ----- > > From: "Unidata GEMPAK Support" <address@hidden> > To: address@hidden > Cc: address@hidden > Sent: Thursday, December 5, 2013 11:56:52 AM > Subject: [GEMPAK #XFB-497782]: GDPVSF GEMPAK: Error Messages > > Hi Vanesa, > > I made updates to GDPVSF earlier this year while working with another user. > If you could send me full the input for GDPVSF and GDVINT if you're using it, > I can go from there to help you with this. > > -Michael > Unidata > > > > Hello There, > > > > I am currently trying to create a Dynamic Tropopause (DT) map with plotted > > theta values. I am using the GDPVSF program in GEMPAK and I am getting an > > error message-- [DG_NFIL -51]. I have not been able to find any > > documentation online that has error message decoding for this particular > > program. Any help would be much appreciated. > > > > I have seen documentation online about this program being out of date/not > > maintained/not supported. In the case that this is true, do you know of > > any other program in GEMPAK that would be able to plot maps of theta on a > > iso-PV surface (aka, DT map). > > > > ~Vanesa Urango > > > > ------------------------------------------- > > Program Support Assistant, PASS Office > > Graduate Assistant, Meteorology Dept. > > BS Meteorology 2013 > > MS Applied Meteorology Candidate > > > > Dept. of Atmospheric Science & Chemistry > > Plymouth State University > > MSC #48, 17 High St. > > Plymouth, NH 03264 > > > > address@hidden > > (207)479-3876 > > > > "Laughter translates into any language" --Unknown > > > > > > > > > Ticket Details > =================== > Ticket ID: XFB-497782 > Department: Support GEMPAK > Priority: Normal > Status: Open > > > > Ticket Details =================== Ticket ID: XFB-497782 Department: Support GEMPAK Priority: Urgent Status: Open
Attachment:
pvbl.sh
Description: application/shellscript