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

20031030: nagrib & gdplot

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: 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