[ncl-talk] reading ensemble grib files in NCL

[Cathy Smith (NOAA Affiliate)]

6.3.0 linux.

file is

air.2m.1901.grb

Cathy

On 2/22/17 8:13 AM, Dennis Shea wrote:
> What version of NCL are you using?
>
> ===
>
> Please do put sample file at:
>
> ftp ftp.cgd.ucar.edu <http://ftp.cgd.ucar.edu>
> anonymous
> your_email
> cd incoming
> put ...
> quit
>
>
> On Wed, Feb 22, 2017 at 8:09 AM, Cathy Smith (NOAA Affiliate) 
> <cathy.smith at noaa.gov <mailto:cathy.smith at noaa.gov>> wrote:
>
>     Hi
>     We have some grib files that have 10 ensemble members for each month.
>     NCL reads the data as one 1 value per month using ncl_convert2nc.
>     Looking at the grib1 files, the 1st 10 headers are the same, the 2nd,
>     and so forth (except for the starting value). Is it possible to read
>     these with ncl or should we try to create 120 grib files and then
>     process them? Maybe one can use -sed?
>
>     Here is the header for the 1st 12 grib records of the 120.
>
>     1:0:d=01010100:2T:kpds5=167:kpds6=1:kpds7=0:TR=123:P1=0:P2=3:TimeU=1:sfc:anl:ave at 3hr:type=analysis:NAve=248
>     2:58192:d=01010100:2T:kpds5=167:kpds6=1:kpds7=0:TR=123:P1=0:P2=3:TimeU=1:sfc:anl:ave at 3hr:type=analysis:NAve=248
>     3:116384:d=01010100:2T:kpds5=167:kpds6=1:kpds7=0:TR=123:P1=0:P2=3:TimeU=1:sfc:anl:ave at 3hr:type=analysis:NAve=248
>     4:174576:d=01010100:2T:kpds5=167:kpds6=1:kpds7=0:TR=123:P1=0:P2=3:TimeU=1:sfc:anl:ave at 3hr:type=analysis:NAve=248
>     5:232768:d=01010100:2T:kpds5=167:kpds6=1:kpds7=0:TR=123:P1=0:P2=3:TimeU=1:sfc:anl:ave at 3hr:type=analysis:NAve=248
>     6:290960:d=01010100:2T:kpds5=167:kpds6=1:kpds7=0:TR=123:P1=0:P2=3:TimeU=1:sfc:anl:ave at 3hr:type=analysis:NAve=248
>     7:349152:d=01010100:2T:kpds5=167:kpds6=1:kpds7=0:TR=123:P1=0:P2=3:TimeU=1:sfc:anl:ave at 3hr:type=analysis:NAve=248
>     8:407344:d=01010100:2T:kpds5=167:kpds6=1:kpds7=0:TR=123:P1=0:P2=3:TimeU=1:sfc:anl:ave at 3hr:type=analysis:NAve=248
>     9:465536:d=01010100:2T:kpds5=167:kpds6=1:kpds7=0:TR=123:P1=0:P2=3:TimeU=1:sfc:anl:ave at 3hr:type=analysis:NAve=248
>     10:523728:d=01010100:2T:kpds5=167:kpds6=1:kpds7=0:TR=123:P1=0:P2=3:TimeU=1:sfc:anl:ave at 3hr:type=analysis:NAve=248
>     11:581920:d=01020100:2T:kpds5=167:kpds6=1:kpds7=0:TR=123:P1=0:P2=3:TimeU=1:sfc:anl:ave at 3hr:type=analysis:NAve=224
>     12:640112:d=01020100:2T:kpds5=167:kpds6=1:kpds7=0:TR=123:P1=0:P2=3:TimeU=1:sfc:anl:ave at 3hr:type=analysis:NAve=224
>
>     ncl_convert2nc produces
>
>     netcdf air.2m.1901 {
>     dimensions:
>              initial_time0_hours = 12 ;
>              g0_lat_1 = 121 ;
>              g0_lon_2 = 240 ;
>              ncl_strlen_0 = 18 ;
>     variables:
>              float \2T_GDS0_SFC_S123(initial_time0_hours, g0_lat_1,
>     g0_lon_2) ;
>                      \2T_GDS0_SFC_S123:N = 248, 224, 248, 240, 248, 240,
>     248, 248, 240, 248, 240, 248 ;
>
>     etc.
>
>
>     I can put a file in ftp if that helps.
>
>     Cathy
>
>     --
>     ----------------------------------------------
>     NOAA/ESRL PSD and CU CIRES
>     303-497-6263
>     http://www.esrl.noaa.gov/psd/people/cathy.smith/
>     <http://www.esrl.noaa.gov/psd/people/cathy.smith/>
>
>     Emails about data/webpages may get quicker responses from emailing
>     esrl.psd.data at noaa.gov <mailto:esrl.psd.data at noaa.gov>
>
>     _______________________________________________
>     ncl-talk mailing list
>     ncl-talk at ucar.edu <mailto:ncl-talk at ucar.edu>
>     List instructions, subscriber options, unsubscribe:
>     http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>     <http://mailman.ucar.edu/mailman/listinfo/ncl-talk>
>
>

-- 
----------------------------------------------
NOAA/ESRL PSD and CU CIRES
303-497-6263
http://www.esrl.noaa.gov/psd/people/cathy.smith/

Emails about data/webpages may get quicker responses from emailing
esrl.psd.data at noaa.gov

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170222/f204ab1e/attachment.html