[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[GEMPAK #WYS-922950]: RE: NMAP Loops (#1562)
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 #WYS-922950]: RE: NMAP Loops (#1562)
- Date: Mon, 26 Feb 2007 13:34:36 -0700
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
- Prev by Date: [GEMPAK #ISM-524898]: Computing max or min over subset of a grid
- Next by Date: [Support #BEM-969521]: Gempak - NAGRIB and rotated lat/lon grid
- Previous by thread: [GEMPAK #WYS-922950]: RE: NMAP Loops (#1562)
- Next by thread: [GEMPAK #LXA-807502]: gdstat
- Index(es):