[GEMPAK #WYS-922950]: RE: NMAP Loops (#1562)

Charley,

The default maximun number of grids in a file created by dcgrib2 is 5000, 
however, you
can configure this on the command line using the "-m #" option (see dcgrib2 -h 
for help!).

You can use "dcgrib2 -m 24000 etc....."

There are just under 18000 grids in the 3 hourly data through 84 hours:


 Number of grids in file: 17827

 Maximum number of grids in file:  24000

Steve Chiswell
Unidata User Support

> Hi,
> 
> Thanks for the help on this issue.  I am getting closer.  I got dcgrib2
> to work using the following script under gempak 5.9.4; however, it
> appears to be stopping somewhere within the 48 hour forecast run (or
> after 5000 records, see snippet from last few lines of gdinfo for this
> gem file, below script).  Any ideas on what is causing this?
> 
> #! /bin/csh
> # Script for decoding NAM model grib data to gempak format for use with
> NMAP and GEMPAK.
> # Uses dcgrib2 and can be run from cron.
> #
> # Script created February 20, 2007 by Charley Knoderer
> #
> # Define current date
> set yy=`date -u +%y`
> set mm=`date -u +%m`
> set dd=`date -u +%d`
> set date=${yy}${mm}${dd}
> set hh=`date -u +%H`
> #
> # Define the output data file names. Naming convention:
> YYMMDDHH_nam_212.gem
> #
> set Grb12=/usr/local/ldm/data/gempak/model/sundog-grib/${date}12
> #
> # Define the output data file names. Naming convention:
> YYMMDDHH_nam_212.gem
> #
> set
> gdf12=/usr/local/ldm/data/gempak/model/sundog-grib/${date}12_nam_212.gem
> #set gdf12=/home/gempak/exe/${date}12_nam_212.gem
> #
> cd /usr/local/ldm/data/gempak/model/sundog-grib
> dcgrib2 -d logfile ${gdf12} < ${Grb12}.nam.t12z.awip3d00.tm00
> dcgrib2 -d logfile ${gdf12} < ${Grb12}.nam.t12z.awip3d06.tm00
> dcgrib2 -d logfile ${gdf12} < ${Grb12}.nam.t12z.awip3d12.tm00
> dcgrib2 -d logfile ${gdf12} < ${Grb12}.nam.t12z.awip3d18.tm00
> dcgrib2 -d logfile ${gdf12} < ${Grb12}.nam.t12z.awip3d24.tm00
> dcgrib2 -d logfile ${gdf12} < ${Grb12}.nam.t12z.awip3d30.tm00
> dcgrib2 -d logfile ${gdf12} < ${Grb12}.nam.t12z.awip3d36.tm00
> dcgrib2 -d logfile ${gdf12} < ${Grb12}.nam.t12z.awip3d42.tm00
> dcgrib2 -d logfile ${gdf12} < ${Grb12}.nam.t12z.awip3d48.tm00
> dcgrib2 -d logfile ${gdf12} < ${Grb12}.nam.t12z.awip3d54.tm00
> dcgrib2 -d logfile ${gdf12} < ${Grb12}.nam.t12z.awip3d60.tm00
> dcgrib2 -d logfile ${gdf12} < ${Grb12}.nam.t12z.awip3d66.tm00
> dcgrib2 -d logfile ${gdf12} < ${Grb12}.nam.t12z.awip3d72.tm00
> dcgrib2 -d logfile ${gdf12} < ${Grb12}.nam.t12z.awip3d78.tm00
> dcgrib2 -d logfile ${gdf12} < ${Grb12}.nam.t12z.awip3d84.tm00
> 
> mv -f ${gdf12} /usr/local/ldm/data/gempak/model/sundog-grids/.
> 
> 
> 
> Logfile:
> 
> 4971     070226/1200F048                          0         NONE WXTZ
> 4972     070226/1200F048                          0         NONE WXTR
> 4973     070226/1200F048                          0         NONE CAPE
> 4974     070226/1200F048                          0         NONE CINS
> 4975     070226/1200F048                          0         NONE PWTR
> 4976     070226/1200F048                          0         NONE CLD
> 4977     070226/1200F048                          0         NONE SWRD06
> 4978     070226/1200F048                          0         NONE LAND
> 4979     070226/1200F048                        180  0      PDLY LFT4
> 4980     070226/1200F048                        180  0      PDLY CAPE
> 4981     070226/1200F048                        180  0      PDLY CINS
> 4982     070226/1200F048                        150  120    PDLY TMPK
> 4983     070226/1200F048                        150  120    PDLY RELH
> 4984     070226/1200F048                        150  120    PDLY UREL
> 4985     070226/1200F048                        150  120    PDLY VREL
> 4986     070226/1200F048                        120  90     PDLY TMPK
> 4987     070226/1200F048                        120  90     PDLY UREL
> 4988     070226/1200F048                        120  90     PDLY VREL
> 4989     070226/1200F048                         90  60     PDLY TMPK
> 4990     070226/1200F048                         90  60     PDLY RELH
> 4991     070226/1200F048                         90  60     PDLY UREL
> 4992     070226/1200F048                         90  60     PDLY VREL
> 4993     070226/1200F048                         60  30     PDLY TMPK
> 4994     070226/1200F048                         60  30     PDLY UREL
> 4995     070226/1200F048                         60  30     PDLY VREL
> 4996     070226/1200F048                         30  0      PDLY TMPK
> 4997     070226/1200F048                         30  0      PDLY RELH
> 4998     070226/1200F048                         30  0      PDLY UREL
> 4999     070226/1200F048                         30  0      PDLY VREL
> 5000     070226/1200F048                         30  0      PDLY LIFT
> 
> Thanks,
> Charley
> 
> -----Original Message-----
> From: Unidata GEMPAK Support
> [mailto:address@hidden]
> Sent: Tuesday, February 20, 2007 11:27 AM
> To: Charley Knoderer
> Cc: address@hidden
> Subject: [GEMPAK #WYS-922950]: RE: NMAP Loops (#1562)
> 
> 
> Charley,
> 
> Make sure you are decoding using the dcgrib2 program from 5.9.4 and not
> 5.7.3. As I mentioned before, NCEP changed one byte of the grib headers
> that changed the F000 meaning to an analysis time - so the 5.7.3 version
> would be decoding these initial times as YYMMDD/HHNN rather than
> YYMMDD/HHNNf000.
> 
> You should be able to see the grid times in the file if you correctly
> fill out the entries in gdinfo rather than leaving blank.
> 
> Try
> GDATTIM = fall
> GLEVEL = all
> GVCORD = all
> GFUNC = all
> 
> Steve Chiswell
> Unidata User Support
> 
> 
> > Hi,
> >
> > Here's the listing I get off of gdinfo:
> >
> > GEMPAK-GDINFO>l
> > GDFILE   = 07021918_nam_212.gem
> > LSTALL   = YES
> > OUTPUT   =
> > GDATTIM  = all
> > GLEVEL   =
> > GVCORD   =
> > GFUNC    =
> > GEMPAK-GDINFO>r
> >
> > GRID FILE: 07021918_nam_212.gem
> >
> > GRID NAVIGATION:
> > PROJECTION:          LCC
> > ANGLES:                25.0   -95.0    25.0
> > GRID SIZE:          185 129
> > LL CORNER:              12.19   -133.46
> > UR CORNER:              57.29    -49.39
> >
> > GRID ANALYSIS BLOCK:
> > ANALYSIS TYPE:        BARNES
> > DELTAN:               4.000
> > DELTAX:           -9999.000
> > DELTAY:           -9999.000
> > GRID AREA:            11.00 -154.00   62.00  -48.00
> > EXTEND AREA:          11.00 -154.00   62.00  -48.00
> > DATA AREA:            11.00 -154.00   62.00  -48.00
> >
> > Number of grids in file:  9210
> >
> > Maximum number of grids in file:  29999
> >
> > [GDU 2]  Did not find any matching grids.
> > Parameters requested:
> > GDFILE,LSTALL,OUTPUT,GDATTIM,GLEVEL,GVCORD,GFUNC.
> >
> > To add further confusion to the matter, we have one machine running
> > gempak 5.7.3 and it displays the first time period in nmap, but only
> > the first time period.  On a separate machine, we are running gempak
> > 5.9.4 and reading the decoded grib from 5.7.3 and nmap will not read
> > these files at all.
> >
> > Any suggestions?
> >
> > Thanks,
> > Charley
> >
> > -----Original Message-----
> > From: Steve Chiswell [mailto:address@hidden]
> > Sent: Wednesday, January 31, 2007 11:39 AM
> > To: Charley Knoderer
> > Cc: address@hidden; GEMPAK support
> > Subject: Re: NMAP Loops (#1562)
> >
> >
> > Charley,
> >
> > Check your "GDINFO" listing of your decoded file. Your statement about
> 
> > the "first" time period being missing may indicate that the source of
> > your grids changed from using the PDS identifier for F000 to an
> > analysis time . The Jan 22 time frame is coincident with changes in
> > the mainframe at NWS. If that is the case, try upgrading to 5.9.4.
> >
> > Steve Chiswell
> > Unidata User Support
> >
> >
> > On Wed, 2007-01-31 at 10:29 -0800, Charley Knoderer wrote:
> > > Hi,
> > >
> > > We are having a very strange problem with NMAP not reading the
> > > decoded
> >
> > > grib files correctly.  Gempak seems to be reading them for the maps
> > > I create, with the exception of the first time period. However the
> > > loops
> >
> > > are not loading.
> > >
> > > This problem first appeared between 12 Z January 22 and 00 Z January
> 
> > > 23.  We have not made any changes to the decoding process (we use
> > > nagrib to decode the grib files) before or after this date.  This
> > > problem is apparent on using NMAP 5.7.3 and 5.9.3.
> > >
> > > Here's the script we are using:
> > >
> > > #!/bin/csh
> > > #
> > > #  CONVERT_GEMPAK.csh
> > > #
> > > #  This is a simple script to write model GRIB data to a GEMPAK #
> > > file using the NAGRIB program found in the NAWIPS distribution. #
> > >
> > ######################################################################
> > ##
> > >
> > > #  BEGIN THE SCRIPT
> > >
> > > #  Check to see if the user has the NAWIPS environmental variables
> > > set. #  Uncomment out the following line if NAWIPS is not included
> > > in the #  users environment. #
> > > #  setenv NAWIPS /usr1/nawips
> > > if ( $?NAWIPS ) then
> > >    source $NAWIPS/Nawips.cshrc
> > > else
> > >    echo "WARNING: Please Define the location of the NAWIPS
> > > distribution at the top"
> > >    echo "         of ${GMASTER_BIN}/convert_gempak and rerun the
> > > script. - EXIT"
> > >    exit 1
> > > endif
> > >
> > > set YYMMDDHH = $1
> > > set MODEL    = $2
> > > set SUFFX    = $3
> > >
> > > if ( ! -e ${GRIB_PATH}/.local_gribs || -z ${GRIB_PATH}/.local_gribs
> > > ) then
> > >    echo " "
> > >    echo "   ERROR: No grib files found - EXIT"
> > >    echo " "
> > >    exit 1
> > > endif
> > >
> > > set GRIBS = `cat ${GRIB_PATH}/.local_gribs`
> > >
> > > set GEMFILE = ${YYMMDDHH}_${MODEL}_${SUFFX}.gem
> > >
> > > if ( ! -e ${GEM_PATH} ) mkdir -p ${GEM_PATH}
> > >
> > > echo "  CONVERT: Creating GEMPAK file ${GEM_PATH}/${GEMFILE}"
> > >
> > > if ( -e gemfile ) rm -f gemfile
> > >
> > > echo "INDXFL   =           "  > gemfile
> > > echo "GDOUTF   = ${GEM_PATH}/${GEMFILE}" >> gemfile
> > > echo "PROJ     =           " >> gemfile
> > > echo "GRDAREA  =           " >> gemfile
> > > echo "KXKY     =           " >> gemfile
> > > echo "MAXGRD   = 29999     " >> gemfile
> > > echo "CPYFIL   = GDS       " >> gemfile
> > > echo "GAREA    =           " >> gemfile
> > > echo "OUTPUT   = F/${GEM_PATH}/nagrib.fil" >> gemfile
> > > echo "GBTBLS   =           " >> gemfile
> > > echo "GBDIAG   =           " >> gemfile
> > >
> > > foreach GRIB ($GRIBS)
> > >    echo "GBFILE   = ${GRIB_PATH}/${GRIB}" >> gemfile
> > >    echo "run               " >> gemfile
> > >    echo "                  " >> gemfile
> > > end
> > >
> > > echo "exit"                  >> gemfile
> > >
> > > $GEMEXE/nagrib < gemfile >& /dev/null
> > > $GEMEXE/gpend            >& /dev/null
> > >
> > > rm -f gemfile *.nts >& /dev/null
> > >
> > > #  End of the Show
> > > exit 0
> > >
> > > Anyone have any ideas?
> > >
> > > Thanks,
> > > Charley
> > >
> > > Charley Knoderer
> > > Meteorologist II
> > > Sonoma Technology, Inc.
> > > 1360 Redwood Way, Suite C
> > > Petaluma, CA  94954
> > > http://www.sonomatech.com
> > > Phone: (707) 665-9900
> > > Fax: (707) 665-9800
> > >
> > --
> > Steve Chiswell <address@hidden>
> > Unidata
> >
> >
> 
> 
> Ticket Details
> ===================
> Ticket ID: WYS-922950
> Department: Support GEMPAK
> Priority: Normal
> Status: Closed
> 
> 
> 


Ticket Details
===================
Ticket ID: WYS-922950
Department: Support GEMPAK
Priority: Normal
Status: Closed