<div dir="ltr">It work,<div><br></div><div>Thank you.</div><div><br></div><div><br></div><div><br></div><div><br></div><div class="gmail_extra"><br><div class="gmail_quote">2017-07-01 10:58 GMT-04:00 Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div><div><div><div>Hi Rabah,<br><br></div>NCL's default mode when importing dimensions of size one is to eliminate that dimension. Size one dimensions are considered 'degenerate.' GRIB formatted 'gfs' files [ eg: gfs.t00z.pgrb2.0p25.f000 , gfs.t00z.pgrb2.0p25.f006, gfs.t00z.pgrb2.0p25.f012, ...] have two 'size one' temporal dimensions: initial_time0_hours and forecast_time0:<br><br></div>initial_time0_hours: initial time of the forecast; this should be the same for each file<br></div>forecast_time0: forecast time [hours offset from the initial_time0_hours]; 0, 6, 12, ... <br><br></div>For the gfs files you are reading, the 'initial_time0_hours' will be the same for all the files but the 'forecast_time0' will will differ.<br><br></div>Try the following:<br><br>%> ncl_filedump gfs.t00z.pgrb2.0p25.f006.grb2 <wbr> | less<br>%> ncl_filedump gfs.t00z.pgrb2.0p25.f006.grb2 -ftime | less<br>%> ncl_filedump gfs.t00z.pgrb2.0p25.f006.grb2 -itime -ftime | less<br><br></div>The -itime and -ftime force NCL to not consider the single element dimension sizes as degenerate.<br><br></div><div>%> ncl_filedump -h<br></div><div><br>------<br></div>Your situation <br><br>[1]<br>grib_files = systemfunc("ls gfs.t00z.pgrb2.0p25.f*") + ".grib2"<span class=""><br>a = addfiles(grib_files,"r")<br></span>print( a[0] ) ; no 'time' dimensions<br> <wbr> ; 'time information' is included but are variable attributes<br><br>x = a[:]->.... <wbr> ; where ... is some variable on your file<br></div>print(x) <wbr> ; x will be 2-dimension ... no time dimension(s)<br></div> <wbr> ; 'time information' is included but are attributes<br><div>==============<br></div><div>Start over ...<br><br>... aggregating the 'forecast dimension'<br><br>Read 'SingleDimensions' portion of:<br> <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/setfileoption.shtml" target="_blank">http://www.ncl.ucar.edu/<wbr>Document/Functions/Built-in/<wbr>setfileoption.shtml</a><br><br>[2]<br><div>setfileoption("grb","<wbr>SingleElementDimensions","<wbr>forecast_time") <br>grib_files = systemfunc("ls gfs.t00z.pgrb2.0p25.f*") + ".grib2"<span class=""><br>a = addfiles(grib_files,"r")<br></span>print( a[0] ) <br><br>x = a[:]->...<br>print(x)<br><br></div><div>Cheers<br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote"><span class="">On Fri, Jun 30, 2017 at 4:10 PM, Rabah Hachelaf <span dir="ltr"><<a href="mailto:hachelaf@sca.uqam.ca" target="_blank">hachelaf@sca.uqam.ca</a>></span> wrote:<br></span><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr"><div class="gmail_quote"><div dir="ltr">Hello,<div><br></div><div>Why we can not concatenate or joint GFS 0.25 grib2 file </div><div>like in this example ?!</div><div>Any suggestions to solve this issue ?</div><div><br></div><div>Thank you</div><div>Rabah<br></div><div><br></div><div><br></div><div><br></div><div>grib_files = systemfunc("ls "+gfs_dir+"/gfs*") + ".grib2"</div><div>a = addfiles(grib_files,"r") </div><div><br></div><div>ListSetType (a, "join") or ListSetType (a, "cat")<br clear="all"><div>rain =a[:]->APCP_P8_L1_GLL0_acc<br></div><div>printVarSummary(rain)<br></div><div><br></div><div>------------------------------<wbr>--------------------</div><div><div>Variable: rain</div><div>Type: float</div><div>Total Size: 4152960 bytes</div><div> 1038240 values</div><div>Number of Dimensions: 2</div><div>Dimensions and sizes: [lat_0 | 721] x [lon_0 | 1440]</div><div>Coordinates:</div><div> lat_0: [90..-90]</div><div> lon_0: [ 0..359.75]</div><div>Number Of Attributes: 15</div><div> center : US National Weather Service - NCEP (WMC)</div><div> production_status : Operational products</div><div> long_name : Total precipitation</div><div> units : kg m-2</div><div> _FillValue : 1e+20</div><div> grid_type : Latitude/longitude</div><div> parameter_discipline_and_categ<wbr>ory : Meteorological products, Moisture</div><div> parameter_template_discipline_<wbr>category_number : ( 8, 0, 1, 8 )</div><div> level_type : Ground or water surface</div><div> level : 0</div><div> type_of_statistical_processing : Accumulation</div><div> statistical_process_duration : initial time to forecast time</div><div> forecast_time : 6</div><div> forecast_time_units : hours</div><div> initial_time : 06/23/2017 (00:00)</div></div><span class="m_1298647055638245156HOEnZb"><font color="#888888"><span class="m_1298647055638245156m_1963491885106915389HOEnZb"><font color="#888888"><div><br></div><div><br></div>-- <br><div class="m_1298647055638245156m_1963491885106915389m_5898273391518277278gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>------------------------------</div>Cordialement,<br>Best regards,<br>Rabah Hachelaf <br> <br><br></div></div></div></div>
</font></span></font></span></div></div>
</div><br><br clear="all"><div><br></div><div class="m_1298647055638245156m_1963491885106915389gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"> <br></div></div>
</div>
<br></div></div><span class="">______________________________<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>
<br></span></blockquote></div><br></div>
</blockquote></div><br><br clear="all"><div><br></div>
</div></div>