[ncl-talk] Fwd: Conactenate or joint GFS grib2 files

Rabah Hachelaf hachelaf at sca.uqam.ca
Fri Jun 30 16:10:38 MDT 2017


Hello,

Why we can not concatenate or joint GFS 0.25 grib2 file
like in this example ?!
Any suggestions to solve this issue ?

Thank you
Rabah



grib_files = systemfunc("ls "+gfs_dir+"/gfs*") + ".grib2"
a = addfiles(grib_files,"r")

ListSetType (a, "join") or ListSetType (a, "cat")
rain =a[:]->APCP_P8_L1_GLL0_acc
printVarSummary(rain)

--------------------------------------------------
Variable: rain
Type: float
Total Size: 4152960 bytes
            1038240 values
Number of Dimensions: 2
Dimensions and sizes:   [lat_0 | 721] x [lon_0 | 1440]
Coordinates:
            lat_0: [90..-90]
            lon_0: [ 0..359.75]
Number Of Attributes: 15
  center :      US National Weather Service - NCEP (WMC)
  production_status :   Operational products
  long_name :   Total precipitation
  units :       kg m-2
  _FillValue :  1e+20
  grid_type :   Latitude/longitude
  parameter_discipline_and_category :   Meteorological products, Moisture
  parameter_template_discipline_category_number :       ( 8, 0, 1, 8 )
  level_type :  Ground or water surface
  level :        0
  type_of_statistical_processing :      Accumulation
  statistical_process_duration :        initial time to forecast time
  forecast_time :       6
  forecast_time_units : hours
  initial_time :        06/23/2017 (00:00)


-- 
------------------------------
Cordialement,
Best regards,
Rabah Hachelaf
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170630/0a707f83/attachment.html 


More information about the ncl-talk mailing list