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

20030619: NCEP/NCAR Reanalysis on cd and dcgrib



Ruth,

The grib data you have is specifying a period of time for the average
precipitation rate, and so you end up with 2 times (time1:time2).
In gdlist, to list the first time in your file, you would need to set
GDATTIM=790101/0000:790101/0000

I downloaded your file and sucessfully listed one of the fields using:
 GDATTIM  = 790801/0000:790801/0000
 GLEVEL   = 0
 GVCORD   = none
 GFUNC    = ave_pr--
 GDFILE   = 2001.cmporain.sfc.gem
 GAREA    = dset
 PROJ     = ced
 SCALE    = 999
 OUTPUT   = t
 GEMPAK-GDLIST>



Actually, I would have expected the average to be from Jan 1 to Jan 31 etc,
instead of having the same time for both time1 and time2.
I'd suggest using dcgrib2 rather than dcgrib (depricated) as
it handles the different grib time types better. See if
that gives the same start/stop times for the average (if so, its
really that way in the GRIB). But, thats just cosmetic...it 
won't affect your displays.


Steve Chiswell

>From: Ruth Platner <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200306191735.h5JHZwLd017057

>
>--------------730365F8EB909D7E997E0C18
>Content-Type: text/plain; charset=us-ascii
>Content-Transfer-Encoding: 7bit
>
>Hi,
>
>I have the NCEP/NCAR Reanalysis data on cd. I've converted many of these
>files to gempak files using dcgrib. I haven't had any problems until I
>converted the monthly average of observations for PRATE, precipitation
>rate.
>
>I don't get any errors from dcgrib, I logged the conversion verbosely
>and everything seems fine.
>
>Here's the command I used:
>
>/ljuka/ldm/gempak/bin/sol/dcgrib -d dcgrib.log -g /ljuka/ldm/gempak/gempak/tab
> les -v 2001.cmporain.sfc.gem < /cdrom/reanal01/data/obs/cmporain.sfc
>
>gdinfo will list the file information as long as gdattim=all.
>
> GDFILE   = 2001.cmporain.sfc.gem
>
> LSTALL   = YES
>
> OUTPUT   = T
>
> GDATTIM  = all
>
> GLEVEL   = 0
>
> GVCORD   = none
>
> GFUNC    = AVE_PR--
>
> GRID FILE: 2001.cmporain.sfc.gem
>
> GRID NAVIGATION:
>
>     PROJECTION:          CED
>
>     GRID SIZE:          144  73
>
>     LL CORNER:             -90.00      0.00
>
>     UR CORNER:              90.00     -2.50
>
> GRID ANALYSIS BLOCK:
>
>      UNKNOWN ANALYSIS TYPE
>
> Number of grids in file:   278
>
> Maximum number of grids in file:   1500
>
>  NUM       TIME1              TIME2           LEVL1 LEVL2  VCORD PARM
>
>    1     790101/0000       790101/0000            0         NONE AVE_PR--
>
>    2     790201/0000       790201/0000            0         NONE AVE_PR--
>
>    3     790301/0000       790301/0000            0         NONE AVE_PR--
>
>If I set gdattim=first, it can't find the data
>
> GEMPAK-GDINFO>gdattim=first
>
> GEMPAK-GDINFO>run
>
> GRID FILE: 2001.cmporain.sfc.gem
>
> GRID NAVIGATION:
>
>     PROJECTION:          CED
>
>     GRID SIZE:          144  73
>
>     LL CORNER:             -90.00      0.00
>
>     UR CORNER:              90.00     -2.50
>
> GRID ANALYSIS BLOCK:
>
>      UNKNOWN ANALYSIS TYPE
>
> Number of grids in file:   278
>
> Maximum number of grids in file:   1500
>
> [GDU 2]  Did not find any matching grids.
>
>Gdlist gives me the following:
>
>aspre monthlyrean 69% gdlist
>
> GDATTIM   Grid date/time                    790301/0000
>
> GLEVEL    Grid level                        0
>
> GVCORD    Grid vertical coordinate          none
>
> GFUNC     Scalar grid                       AVE_PR--
>
> GDFILE    Grid file                         2001.cmporain.sfc.gem
>
> GAREA     Graphics area                     20;85;60;150
>
> PROJ      Map projection/angles/margins|dr  ced
>
> SCALE     Scalar scale / vector scale       999
>
> OUTPUT    Output device/filename            T
>
> GEMPAK-GDLIST>run
>
> [DG -7]  Input grid AVE_PR-- ^790301/0000 @0 %NONE cannot be found.
>
>This is true for all the dates, not just the first date. It seems that
>the problem is with the dates, but I'm not sure.
>
>Do you have any ideas?
>
>2001.cmporain.sfc.gem is an 11MB file so I haven't attached it, but I
>put it and the dcgrib.log file at
>ftp://oldeuros.gso.uri.edu/pub/rlp/gempak/
>
>Thanks,
>
>Ruth
>
>
>
>--------------730365F8EB909D7E997E0C18
>Content-Type: text/html; charset=us-ascii
>Content-Transfer-Encoding: 7bit
>
><!doctype html public "-//w3c//dtd html 4.0 transitional//en">
><html>
>Hi,
><p>I have the NCEP/NCAR Reanalysis data on cd. I've converted many of these
>files to gempak files using dcgrib. I haven't had any problems until I
>converted the monthly average of observations for PRATE, precipitation
>rate.
><p>I don't get any errors from dcgrib, I logged the conversion verbosely
>and everything seems fine.
><p>Here's the command I used:
><pre>/ljuka/ldm/gempak/bin/sol/dcgrib -d dcgrib.log -g /ljuka/ldm/gempak/gempa
> k/tables -v 2001.cmporain.sfc.gem &lt; /cdrom/reanal01/data/obs/cmporain.sfc<
> /pre>
>gdinfo will list the file information as long as gdattim=all.
><pre>&nbsp;GDFILE&nbsp;&nbsp; = 2001.cmporain.sfc.gem</pre>
>
><pre>&nbsp;LSTALL&nbsp;&nbsp; = YES</pre>
>
><pre>&nbsp;OUTPUT&nbsp;&nbsp; = T</pre>
>
><pre>&nbsp;GDATTIM&nbsp; = all</pre>
>
><pre>&nbsp;GLEVEL&nbsp;&nbsp; = 0</pre>
>
><pre>&nbsp;GVCORD&nbsp;&nbsp; = none</pre>
>
><pre>&nbsp;GFUNC&nbsp;&nbsp;&nbsp; = AVE_PR--</pre>
>
><pre>&nbsp;GRID FILE: 2001.cmporain.sfc.gem</pre>
>
><pre>&nbsp;GRID NAVIGATION:</pre>
>
><pre>&nbsp;&nbsp;&nbsp;&nbsp; PROJECTION:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&
> nbsp;&nbsp;&nbsp; CED</pre>
>
><pre>&nbsp;&nbsp;&nbsp;&nbsp; GRID SIZE:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n
> bsp;&nbsp;&nbsp; 144&nbsp; 73</pre>
>
><pre>&nbsp;&nbsp;&nbsp;&nbsp; LL CORNER:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n
> bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -90.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00<
> /pre>
>
><pre>&nbsp;&nbsp;&nbsp;&nbsp; UR CORNER:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n
> bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 90.00&nbsp;&nbsp;&nbsp;&nbsp; -2.50<
> /pre>
>
><pre>&nbsp;GRID ANALYSIS BLOCK:</pre>
>
><pre>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; UNKNOWN ANALYSIS TYPE</pre>
>
><pre>&nbsp;Number of grids in file:&nbsp;&nbsp; 278</pre>
>
><pre>&nbsp;Maximum number of grids in file:&nbsp;&nbsp; 1500</pre>
>
><pre>&nbsp; NUM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; TIME1&nbsp;&nbsp;&nbsp;&nb
> sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; TIME2&nbsp;&nbsp;&n
> bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LEVL1 LEVL2&nbsp; VCORD PARM</
> pre>
>
><pre>&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp; 790101/0000&nbsp;&nbsp;&nbsp
> ;&nbsp;&nbsp;&nbsp; 790101/0000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs
> p;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NONE A
> VE_PR--</pre>
>
><pre>&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp; 790201/0000&nbsp;&nbsp;&nbsp
> ;&nbsp;&nbsp;&nbsp; 790201/0000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs
> p;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NONE A
> VE_PR--</pre>
>
><pre>&nbsp;&nbsp;&nbsp; 3&nbsp;&nbsp;&nbsp;&nbsp; 790301/0000&nbsp;&nbsp;&nbsp
> ;&nbsp;&nbsp;&nbsp; 790301/0000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs
> p;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NONE A
> VE_PR--</pre>
>If I set gdattim=first, it can't find the data
><pre>&nbsp;GEMPAK-GDINFO>gdattim=first</pre>
>
><pre>&nbsp;GEMPAK-GDINFO>run</pre>
>
><pre>&nbsp;GRID FILE: 2001.cmporain.sfc.gem</pre>
>
><pre>&nbsp;GRID NAVIGATION:</pre>
>
><pre>&nbsp;&nbsp;&nbsp;&nbsp; PROJECTION:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&
> nbsp;&nbsp;&nbsp; CED</pre>
>
><pre>&nbsp;&nbsp;&nbsp;&nbsp; GRID SIZE:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n
> bsp;&nbsp;&nbsp; 144&nbsp; 73</pre>
>
><pre>&nbsp;&nbsp;&nbsp;&nbsp; LL CORNER:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n
> bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -90.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00<
> /pre>
>
><pre>&nbsp;&nbsp;&nbsp;&nbsp; UR CORNER:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n
> bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 90.00&nbsp;&nbsp;&nbsp;&nbsp; -2.50<
> /pre>
>
><pre>&nbsp;GRID ANALYSIS BLOCK:</pre>
>
><pre>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; UNKNOWN ANALYSIS TYPE</pre>
>
><pre>&nbsp;Number of grids in file:&nbsp;&nbsp; 278</pre>
>
><pre>&nbsp;Maximum number of grids in file:&nbsp;&nbsp; 1500</pre>
>
><pre>&nbsp;[GDU 2]&nbsp; Did not find any matching grids.</pre>
>Gdlist gives me the following:
><pre>aspre monthlyrean 69% gdlist</pre>
>
><pre>&nbsp;GDATTIM&nbsp;&nbsp; Grid date/time&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb
> sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n
> bsp; 790301/0000</pre>
>
><pre>&nbsp;GLEVEL&nbsp;&nbsp;&nbsp; Grid level&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n
> bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&
> nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0</pre>
>
><pre>&nbsp;GVCORD&nbsp;&nbsp;&nbsp; Grid vertical coordinate&nbsp;&nbsp;&nbsp;
> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; none</pre>
>
><pre>&nbsp;GFUNC&nbsp;&nbsp;&nbsp;&nbsp; Scalar grid&nbsp;&nbsp;&nbsp;&nbsp;&n
> bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&
> nbsp;&nbsp;&nbsp;&nbsp;&nbsp; AVE_PR--</pre>
>
><pre>&nbsp;GDFILE&nbsp;&nbsp;&nbsp; Grid file&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb
> sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n
> bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2001.cmporain.sfc.gem</pre>
>
><pre>&nbsp;GAREA&nbsp;&nbsp;&nbsp;&nbsp; Graphics area&nbsp;&nbsp;&nbsp;&nbsp;
> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp
> ;&nbsp;&nbsp;&nbsp; 20;85;60;150</pre>
>
><pre>&nbsp;PROJ&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Map projection/angles/margins|dr
> &nbsp; ced</pre>
>
><pre>&nbsp;SCALE&nbsp;&nbsp;&nbsp;&nbsp; Scalar scale / vector scale&nbsp;&nbs
> p;&nbsp;&nbsp;&nbsp;&nbsp; 999</pre>
>
><pre>&nbsp;OUTPUT&nbsp;&nbsp;&nbsp; Output device/filename&nbsp;&nbsp;&nbsp;&n
> bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; T</pre>
>
><pre>&nbsp;GEMPAK-GDLIST>run</pre>
>
><pre>&nbsp;[DG -7]&nbsp; Input grid AVE_PR-- ^790301/0000 @0 %NONE cannot be f
> ound.</pre>
>This is true for all the dates, not just the first date. It seems that
>the problem is with the dates, but I'm not sure.
><p>Do you have any ideas?
><p>2001.cmporain.sfc.gem is an 11MB file so I haven't attached it, but
>I put it and the dcgrib.log file at <A HREF="ftp://oldeuros.gso.uri.edu/pub/rl
> p/gempak/">ftp://oldeuros.gso.uri.edu/pub/rlp/gempak/</A>
><p>Thanks,
><p>Ruth
><br>&nbsp;
><br>&nbsp;</html>
>
>--------------730365F8EB909D7E997E0C18--
>