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

[Support #PRI-562474]: Gempak/ER40 reanalysis data?



> 
> In a previous message to me, you wrote:
> 
> >Pete,
> >The ensemble extension comes from the gb/gbensemble.c routine. The ENHW is a 
> >whole
> >ensemble weighted median.
> 
> That is bizzare.. I don't think these are ensemble members, but rather just
> a single value, from what I read of the data set information.

Pete,

The whole ensemble weighted mean bytes are in your PDS. The grib data has the 
ensemble extended pds
past octet 28. Your other software may just not be showing all the bytes of the 
pds.
Presumably, the reanalysis uses perterbations to get a handle on variability.

Steve Chiswell
Unidata User SUpport
> 
> A short bit of output from wgrib for the U/V values shows:
> 
> 949:20033136:d=97020600:U:kpds5=131:kpds6=100:kpds7=300:TR=0:P1=0:P2=0:TimeU=1:300
>  mb:anl:ECMWF stream?:NAve=0
> 950:20054268:d=97020600:V:kpds5=132:kpds6=100:kpds7=300:TR=0:P1=0:P2=0:TimeU=1:300
>  mb:anl:ECMWF stream?:NAve=0
> 951:20075400:d=97020600:U:kpds5=131:kpds6=100:kpds7=400:TR=0:P1=0:P2=0:TimeU=1:400
>  mb:anl:ECMWF stream?:NAve=0
> 952:20096532:d=97020600:V:kpds5=132:kpds6=100:kpds7=400:TR=0:P1=0:P2=0:TimeU=1:400
>  mb:anl:ECMWF stream?:NAve=0
> 953:20117664:d=97020600:U:kpds5=131:kpds6=100:kpds7=500:TR=0:P1=0:P2=0:TimeU=1:500
>  mb:anl:ECMWF stream?:NAve=0
> 954:20138796:d=97020600:V:kpds5=132:kpds6=100:kpds7=500:TR=0:P1=0:P2=0:TimeU=1:500
>  mb:anl:ECMWF stream?:NAve=0
> 955:20159928:d=97020600:U:kpds5=131:kpds6=100:kpds7=600:TR=0:P1=0:P2=0:TimeU=1:600
>  mb:anl:ECMWF stream?:NAve=0
> 956:20181060:d=97020600:V:kpds5=132:kpds6=100:kpds7=600:TR=0:P1=0:P2=0:TimeU=1:600
>  mb:anl:ECMWF stream?:NAve=0
> 
> Why would these be interpreted as ensemble members instead of just a normal
> single U or V value?
> 
> Pete
> 
> >
> >The nagrib decoder provides a parameter PDSEXT that gives you the control to 
> >turn off the
> >appending of extensions to the parameter names (such as the  C001, P001, 
> >etc).
> >Beware, if you turn that off, you need to decode all ensemble members
> >into different files so that they don't overwrite each other (eg, when you 
> >have
> >HGHTC001, HGHTP001, HGHTN001, HGHTP002, etc., without the extension, they
> >will all be HGHT). The ens_ grid functions are designed to work that way with
> >separate files as well, so its not a bad thing, just more complicated in
> >determining the output file name.
> >
> >For dcgrib2, I have the flag to add extensions hardcoded in 
> >$NAWIPS/unidata/ldmbridge/dcgrib2/dcwppg.c
> >to "pdsextflag = TRUE;".
> >
> >I could create a command line variable for dcgrib2 to toggle that as well.
> >
> >Steve Chiswell
> >Unidata User SUpport
> >
> >
> >>
> >> I am trying to convert the ECMWF ER40 reanalysis data from GRIB
> >> to GEMPAK format, (NCAR MSS data set DS118.1) using the following
> >> command:
> >>
> >> cat U57535.feb0105 | /home/gempak/GEMPAK5.9.3/os/linux/bin/dcgrib2 -d - -m 
> >> 10000 -e GEMTBL=/home/gempak/GEMPAK5.9.3/gempak/tables 
> >> YYMMDDHH_ecmwf_er40.gem
> >>
> >> I get a gempak file, but the variables all have ENHW appended to
> >> them, for example:
> >>
> >> GEMPAK-GDINFO>l
> >> GDFILE   = 97020512_ecmwf_er40.gem
> >> LSTALL   = YES
> >> OUTPUT   = T
> >> GDATTIM  = all
> >> GLEVEL   = all
> >> GVCORD   = PRES
> >> GFUNC    = all
> >> GEMPAK-GDINFO>r
> >>
> >> GRID FILE: 97020512_ecmwf_er40.gem
> >>
> >> GRID NAVIGATION:
> >> PROJECTION:          CED
> >> GRID SIZE:          144  73
> >> LL CORNER:             -90.00      0.00
> >> UR CORNER:              90.00    357.50
> >>
> >> GRID ANALYSIS BLOCK:
> >> UNKNOWN ANALYSIS TYPE
> >>
> >> Number of grids in file:   253
> >>
> >> Maximum number of grids in file:  10000
> >>
> >> NUM       TIME1              TIME2           LEVL1 LEVL2  VCORD PARM
> >> 1     970205/1200F000                       1000         PRES PVORENHW
> >> 2     970205/1200F000                       1000         PRES HGHTENHW
> >> 3     970205/1200F000                       1000         PRES TMPKENHW
> >> 4     970205/1200F000                       1000         PRES SPFHENHW
> >> 5     970205/1200F000                       1000         PRES OMEGENHW
> >> 6     970205/1200F000                       1000         PRES VORTENHW
> >> 7     970205/1200F000                       1000         PRES DIVGENHW
> >> 8     970205/1200F000                       1000         PRES RELHENHW
> >> 9     970205/1200F000                       1000         PRES OZMRENHW
> >> 10     970205/1200F000                       1000         PRES UWNDENHW
> >> 11     970205/1200F000                       1000         PRES VWNDENHW
> >> 12     970205/1200F000                        925         PRES PVORENHW
> >> 13     970205/1200F000                        925         PRES HGHTENHW
> >> 14     970205/1200F000                        925         PRES TMPKENHW
> >> 15     970205/1200F000                        925         PRES SPFHENHW
> >> 16     970205/1200F000                        925         PRES OMEGENHW
> >> 17     970205/1200F000                        925         PRES VORTENHW
> >> 18     970205/1200F000                        925         PRES DIVGENHW
> >> 19     970205/1200F000                        925         PRES RELHENHW
> >> 20     970205/1200F000                        925         PRES OZMRENHW
> >> 21     970205/1200F000                        925         PRES UWNDENHW
> >> 22     970205/1200F000                        925         PRES VWNDENHW
> >> 23     970205/1200F000                        850         PRES PVORENHW
> >> 24     970205/1200F000                        850         PRES HGHTENHW
> >> 25     970205/1200F000                        850         PRES TMPKENHW
> >> 26     970205/1200F000                        850         PRES SPFHENHW
> >> 27     970205/1200F000                        850         PRES OMEGENHW
> >> 28     970205/1200F000                        850         PRES VORTENHW
> >> 29     970205/1200F000                        850         PRES DIVGENHW
> >> 30     970205/1200F000                        850         PRES RELHENHW
> >> 31     970205/1200F000                        850         PRES OZMRENHW
> >> 32     970205/1200F000                        850         PRES UWNDENHW
> >> 33     970205/1200F000                        850         PRES VWNDENHW
> >> 34     970205/1200F000                        775         PRES PVORENHW
> >> 35     970205/1200F000                        775         PRES HGHTENHW
> >>
> >>
> >> Any idea where the ENHW comes from, or how to make it not happen?
> >>
> >> Thanks,
> >>
> >> Pete
> >>
> >> --
> >> +>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>+<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<+
> >> ^ Pete Pokrandt                    V 1447  AOSS Bldg  1225 W Dayton St^
> >> ^ Systems Programmer               V Madison,         WI     53706    ^
> >> ^                                  V      address@hidden       ^
> >> ^ Dept of Atmos & Oceanic Sciences V (608) 262-3086 (Phone/voicemail) ^
> >> ^ University of Wisconsin-Madison  V       262-0166 (Fax)             ^
> >> +<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<+>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>+
> >>
> >>
> >
> >
> >Ticket Details
> >===================
> >Ticket ID: PRI-562474
> >Department: Support GEMPAK
> >Priority: Normal
> >Status: Closed
> >
> 
> 
> --
> +>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>+<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<+
> ^ Pete Pokrandt                    V 1447  AOSS Bldg  1225 W Dayton St^
> ^ Systems Programmer               V Madison,         WI     53706    ^
> ^                                  V      address@hidden       ^
> ^ Dept of Atmos & Oceanic Sciences V (608) 262-3086 (Phone/voicemail) ^
> ^ University of Wisconsin-Madison  V       262-0166 (Fax)             ^
> <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<+>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>+
> 
> 


Ticket Details
===================
Ticket ID: PRI-562474
Department: Support GEMPAK
Priority: Normal
Status: Closed