<div dir="ltr"><div><div><div><div>FTR only: the beta version of NCL seems to have the same issue.<br><br></div>%> ncl_filedump air.2m.1901.grb | less<br><br>--<br><br> NCAR Command Language Version 6.4.0 -22Feb2017_0149<br><br><br>Variable: f<br>Type: file<br>filename: air.2m.1901<br>path: air.2m.1901.grb<br> file global attributes:<br> dimensions:<br> initial_time0_hours = 12<br> g0_lat_1 = 121<br> g0_lon_2 = 240<br> variables:<br> float 2T_GDS0_SFC_S123 ( initial_time0_hours, g0_lat_1, g0_lon_2 )<br> center : European Center for Medium-Range Weather Forecasts (RSMC)<br> long_name : 2 metre temperature<br> units : K<br> _FillValue : 1e+20<br> level_indicator : 1<br> gds_grid_type : 0<br> parameter_table_version : 128<br> parameter_number : 167<br> forecast_time : 0<br> forecast_time_units : hours<br> statistical_process_descriptor : average of N uninitialized analyses<br> statistical_process_duration : instantaneous (beginning at reference time at intervals of 3 hours)<br> N : <ARRAY of 12 elements><br><br> double initial_time0_hours ( initial_time0_hours )<br> long_name : initial time<br> units : hours since 1800-01-01 00:00<br><br> double initial_time0_encoded ( initial_time0_hours )<br> long_name : initial time encoded as double<br> units : yyyymmddhh.hh_frac<br><br> float g0_lat_1 ( g0_lat_1 )<br> long_name : latitude<br> GridType : Cylindrical Equidistant Projection Grid<br> units : degrees_north<br> Dj : 1.5<br> Di : 1.5<br> Lo2 : 358.5<br> La2 : -90<br> Lo1 : 0<br> La1 : 90<br><br> float g0_lon_2 ( g0_lon_2 )<br> long_name : longitude<br> GridType : Cylindrical Equidistant Projection Grid<br> units : degrees_east<br> Dj : 1.5<br> Di : 1.5<br> Lo2 : 358.5<br> La2 : -90<br> Lo1 : 0<br> La1 : 90<br><br> string initial_time0 ( initial_time0_hours )<br> long_name : Initial time of first record<br> units : mm/dd/yyyy (hh:mm)<br><br>=========================================<br></div>Note the attributes<br><br> statistical_process_descriptor : average of N uninitialized analyses<br> statistical_process_duration : instantaneous (beginning at reference time at intervals of 3 hours)<br> N : <ARRAY of 12 elements><br><br> f = addfile("air.2m.1901.grb","r")<br>T2m = f-> 2T_GDS0_SFC_S123<br></div> printVarSummary(T2m)<br><br></div> print(T2m@N)<br><br>(0) 248<br>(1) 224<br>(2) 248<br>(3) 240<br>(4) 248<br>(5) 240<br>(6) 248<br>(7) 248<br>(8) 240<br>(9) 248<br>(10) 240<br>(11) 248<br><br><div><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Feb 22, 2017 at 8:18 AM, Cathy Smith (NOAA Affiliate) <span dir="ltr"><<a href="mailto:cathy.smith@noaa.gov" target="_blank">cathy.smith@noaa.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
6.3.0 linux.<br>
<br>
file is<br>
<br>
air.2m.1901.grb<span class="HOEnZb"><font color="#888888"><br>
<br>
Cathy</font></span><div><div class="h5"><br>
<br>
<div class="m_6736075127538425323moz-cite-prefix">On 2/22/17 8:13 AM, Dennis Shea wrote:<br>
</div>
<blockquote type="cite">
<div dir="ltr">
<div>
<div>
<div>
<div>
<div>
<div>What version of NCL are you using?<br>
<br>
===<br>
<br>
Please do put sample file at:<br>
<br>
</div>
ftp <a href="http://ftp.cgd.ucar.edu" target="_blank">ftp.cgd.ucar.edu</a><br>
</div>
anonymous<br>
</div>
your_email<br>
</div>
cd incoming<br>
</div>
put ...<br>
</div>
quit<br>
<br>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Wed, Feb 22, 2017 at 8:09 AM, Cathy
Smith (NOAA Affiliate) <span dir="ltr"><<a href="mailto:cathy.smith@noaa.gov" target="_blank">cathy.smith@noaa.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi<br>
We have some grib files that have 10 ensemble members for
each month.<br>
NCL reads the data as one 1 value per month using
ncl_convert2nc.<br>
Looking at the grib1 files, the 1st 10 headers are the same,
the 2nd,<br>
and so forth (except for the starting value). Is it possible
to read<br>
these with ncl or should we try to create 120 grib files and
then<br>
process them? Maybe one can use -sed?<br>
<br>
Here is the header for the 1st 12 grib records of the 120.<br>
<br>
1:0:d=01010100:2T:kpds5=167:kp<wbr>ds6=1:kpds7=0:TR=123:P1=0:P2=<wbr>3:TimeU=1:sfc:anl:ave@3hr:type<wbr>=analysis:NAve=248<br>
2:58192:d=01010100:2T:kpds5=16<wbr>7:kpds6=1:kpds7=0:TR=123:P1=0:<wbr>P2=3:TimeU=1:sfc:anl:ave@3hr:<wbr>type=analysis:NAve=248<br>
3:116384:d=01010100:2T:kpds5=1<wbr>67:kpds6=1:kpds7=0:TR=123:P1=0<wbr>:P2=3:TimeU=1:sfc:anl:ave@3hr:<wbr>type=analysis:NAve=248<br>
4:174576:d=01010100:2T:kpds5=1<wbr>67:kpds6=1:kpds7=0:TR=123:P1=0<wbr>:P2=3:TimeU=1:sfc:anl:ave@3hr:<wbr>type=analysis:NAve=248<br>
5:232768:d=01010100:2T:kpds5=1<wbr>67:kpds6=1:kpds7=0:TR=123:P1=0<wbr>:P2=3:TimeU=1:sfc:anl:ave@3hr:<wbr>type=analysis:NAve=248<br>
6:290960:d=01010100:2T:kpds5=1<wbr>67:kpds6=1:kpds7=0:TR=123:P1=0<wbr>:P2=3:TimeU=1:sfc:anl:ave@3hr:<wbr>type=analysis:NAve=248<br>
7:349152:d=01010100:2T:kpds5=1<wbr>67:kpds6=1:kpds7=0:TR=123:P1=0<wbr>:P2=3:TimeU=1:sfc:anl:ave@3hr:<wbr>type=analysis:NAve=248<br>
8:407344:d=01010100:2T:kpds5=1<wbr>67:kpds6=1:kpds7=0:TR=123:P1=0<wbr>:P2=3:TimeU=1:sfc:anl:ave@3hr:<wbr>type=analysis:NAve=248<br>
9:465536:d=01010100:2T:kpds5=1<wbr>67:kpds6=1:kpds7=0:TR=123:P1=0<wbr>:P2=3:TimeU=1:sfc:anl:ave@3hr:<wbr>type=analysis:NAve=248<br>
10:523728:d=01010100:2T:kpds5=<wbr>167:kpds6=1:kpds7=0:TR=123:P1=<wbr>0:P2=3:TimeU=1:sfc:anl:ave@3hr<wbr>:type=analysis:NAve=248<br>
11:581920:d=01020100:2T:kpds5=<wbr>167:kpds6=1:kpds7=0:TR=123:P1=<wbr>0:P2=3:TimeU=1:sfc:anl:ave@3hr<wbr>:type=analysis:NAve=224<br>
12:640112:d=01020100:2T:kpds5=<wbr>167:kpds6=1:kpds7=0:TR=123:P1=<wbr>0:P2=3:TimeU=1:sfc:anl:ave@3hr<wbr>:type=analysis:NAve=224<br>
<br>
ncl_convert2nc produces<br>
<br>
netcdf air.2m.1901 {<br>
dimensions:<br>
initial_time0_hours = 12 ;<br>
g0_lat_1 = 121 ;<br>
g0_lon_2 = 240 ;<br>
ncl_strlen_0 = 18 ;<br>
variables:<br>
float \2T_GDS0_SFC_S123(initial_time<wbr>0_hours,
g0_lat_1, g0_lon_2) ;<br>
\2T_GDS0_SFC_S123:N = 248, 224, 248, 240,
248, 240,<br>
248, 248, 240, 248, 240, 248 ;<br>
<br>
etc.<br>
<br>
<br>
I can put a file in ftp if that helps.<br>
<br>
Cathy<br>
<br>
--<br>
------------------------------<wbr>----------------<br>
NOAA/ESRL PSD and CU CIRES<br>
<a href="tel:(303)%20497-6263" value="+13034976263" target="_blank">303-497-6263</a><br>
<a href="http://www.esrl.noaa.gov/psd/people/cathy.smith/" rel="noreferrer" target="_blank">http://www.esrl.noaa.gov/psd/p<wbr>eople/cathy.smith/</a><br>
<br>
Emails about data/webpages may get quicker responses from
emailing<br>
<a href="mailto:esrl.psd.data@noaa.gov" target="_blank">esrl.psd.data@noaa.gov</a><br>
<br>
______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
</blockquote>
</div>
<br>
</div>
</blockquote>
<br>
<pre class="m_6736075127538425323moz-signature" cols="72">--
------------------------------<wbr>----------------
NOAA/ESRL PSD and CU CIRES
<a href="tel:(303)%20497-6263" value="+13034976263" target="_blank">303-497-6263</a>
<a class="m_6736075127538425323moz-txt-link-freetext" href="http://www.esrl.noaa.gov/psd/people/cathy.smith/" target="_blank">http://www.esrl.noaa.gov/psd/<wbr>people/cathy.smith/</a>
Emails about data/webpages may get quicker responses from emailing
<a class="m_6736075127538425323moz-txt-link-abbreviated" href="mailto:esrl.psd.data@noaa.gov" target="_blank">esrl.psd.data@noaa.gov</a></pre>
</div></div></div>
</blockquote></div><br></div>