[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
20031030: nagrib & gdplot
- Subject: 20031030: nagrib & gdplot
- Date: Fri, 31 Oct 2003 14:46:35 -0700
>From: "Seth O McDowell" <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200310301715.h9UHFaOb005888
>After converting a file into gempak format using nagrib, I tried running gdplo
> t to show the data. The map & lat/lon grid show, but I can't get it to output
> the information.
>http://www.wku.edu/~seth.mcdowell/gribstuff.zip is both the original file & it
> converted to gempak.
>
Do the nagrib conversion this way:
GBFILE = pgbexample
INDXFL =
GDOUTF = reanalysis.gem
PROJ =
GRDAREA =
KXKY =
MAXGRD = 2000
CPYFIL = #2
GAREA = dset
OUTPUT = t
GBTBLS = ncarncep1.tbl;ncarncep1.tbl
GBDIAG =
PDSEXT = NO
GEMPAK-NAGRIB>
That will create the gempak file with:
GRID FILE: reanalysis.gem
GRID NAVIGATION:
PROJECTION: CED
ANGLES: 0.0 0.0 0.0
GRID SIZE: 144 73
LL CORNER: -90.00 0.00
UR CORNER: 90.00 -2.50
GRID ANALYSIS BLOCK:
ANALYSIS TYPE: BARNES
DELTAN: 4.000
DELTAX: -9999.000
DELTAY: -9999.000
GRID AREA: -90.00 -180.00 90.00 180.00
EXTEND AREA: -90.00 -180.00 90.00 180.00
DATA AREA: -90.00 -180.00 90.00 180.00
Number of grids in file: 140
Maximum number of grids in file: 2000
The problem you had previously with your testme.gem appears to be a confusion
with cpyfil=#2 or garea=dset.
I was able to decode, then plot with gfunc using 500mb hght as a quick test:
GDATTIM = first
GLEVEL = 500
GVCORD = pres
GFUNC = hght
GDFILE = reanalysis.gem
CINT = 60
LINE = 3
MAP = 1
TITLE = 1
DEVICE = XW
SATFIL =
RADFIL =
PROJ =
GAREA = dset
IJSKIP =
CLEAR = YES
PANEL = 0
TEXT = 1
SCALE = 999
LATLON =
HILO =
HLSYM =
CLRBAR =
CONTUR = 0
SKIP = 0
FINT = 0
FLINE = 10-20
CTYPE = C
LUTFIL =
STNPLT =
IMBAR =
GEMPAK-GDCNTR>
Steve Chiswell