[ncl-talk] time information in WGRIB2 and ncl

Tue Nov 14 13:18:13 MST 2017

Yes: NCL 6.3.0 had a problem as evidenced by the warning messgaes.

 Copyright (C) 1995-2015 - All Rights Reserved
 University Corporation for Atmospheric Research
 NCAR Command Language Version 6.3.0

warning: Entry (62001) not found in code table file
/Users/shea/devel/ncl/lib/ncarg/grib2_codetables/ncep/4/4.3.table
warning: Entry (62006) not found in code table file
/Users/shea/devel/ncl/lib/ncarg/grib2_codetables/ncep/4/4.3.table
warning: Entry (62008) not found in code table file
/Users/shea/devel/ncl/lib/ncarg/grib2_codetables/ncep/4/4.3.table
warning: Entry (62009) not found in code table file
/Users/shea/devel/ncl/lib/ncarg/grib2_codetables/ncep/4/4.3.table
warning: Entry (62010) not found in code table file
/Users/shea/devel/ncl/lib/ncarg/grib2_codetables/ncep/4/4.3.table
warning: Entry (62000) not found in code table file
/Users/shea/devel/ncl/lib/ncarg/grib2_codetables/ncep/4/4.3.table

===
There *should* be 41 forecast hours.
==
Other  issues were addressed and NCL 6.5.0 yields all species and the 41
forecast hours

 %> ncl_filedump ngac.t00z.aod_550nm.grib2 -ftime -v forecast_time0

Copyright (C) 1995-2017 - All Rights Reserved
 University Corporation for Atmospheric Research
 NCAR Command Language Version *6.5.0* -14Nov2017_0149
 The use of this software is governed by a License Agreement.
 See http://www.ncl.ucar.edu/ for more details.

Variable: f
Type: file
filename:       ngac.t00z.aod_550nm
path:   ngac.t00z.aod_550nm.grib2
   file global attributes:
   dimensions:
      forecast_time0 = 41
      lat_0 = 181
      lon_0 = 360
   variables:
      float AOTK_P48_L10_GLL0_A62000 ( forecast_time0, lat_0, lon_0 )

[SNIP]

     integer forecast_time0 ( forecast_time0 )
         long_name :    Forecast offset from initial time
         units :        hours

Variables:

Variable No. 0

*forecast_time0:       0     3     6     9    12    15    18    21    24
27    30    33    36    39    42    45    48    51    54    57    60
63    66    69    72    75    78    81    84    87    90    93    96
99   102   105   108   111   114   117   120*

On Tue, Nov 14, 2017 at 12:23 PM, Debasish Hazra <debasish.hazra5 at gmail.com>
wrote:

> Thanks Dennis. When I used ncl Verison 6.3 ncl_filedump shows :
>
> Variable: f
> Type: file
> filename:       ngac.t00z.aod_550nm
> path:   ngac.t00z.aod_550nm.grib2
>    file global attributes:
>    dimensions:
>       lat_0 = 181
>       lon_0 = 360
>    variables:
>       float AOTK_P48_L10_GLL0 ( lat_0, lon_0 )
>          generating_process_type :      62000
>          center :       US National Weather Service - NCEP (WMC)
>          production_status :    Operational products
>          long_name :    Aerosol optical thickness
>          units :        Numeric
>          _FillValue :   1e+20
>          grid_type :    Latitude/longitude
>          parameter_discipline_and_category :    Meteorological products,
> Atmospheric chemical constituents
>          parameter_template_discipline_category_number :        ( 48, 0,
> 20, 102 )
>          level_type :   Entire Atmosphere
>          level :         0
>          forecast_time :        7
>          forecast_time_units :  minutes
>          initial_time : 11/13/2017 (00:00)
>
> but when I used version 6.4 it shows forecast time as dimension :
>
> Variable: f
> Type: file
> filename:       ngac.t00z.aod_550nm
> path:   ngac.t00z.aod_550nm.grib2
>    file global attributes:
>    dimensions:
>       forecast_time0 = 41
>       lat_0 = 181
>       lon_0 = 360
>    variables:
>       float AOTK_P48_L10_GLL0 ( forecast_time0, lat_0, lon_0 )
>          generating_process_type :      62000
>          center :       US National Weather Service - NCEP (WMC)
>          production_status :    Operational products
>          long_name :    Aerosol optical thickness
>          units :        Numeric
>          _FillValue :   1e+20
>          grid_type :    Latitude/longitude
>          parameter_discipline_and_category :    Meteorological products,
> Atmospheric chemical constituents
>          parameter_template_discipline_category_number :        ( 48, 0,
> 20, 102 )
>          level_type :   Entire Atmosphere
>          level :         0
>          initial_time : 11/13/2017 (00:00)
>          aerosol_type : 62000
>
>
> Debasish.
>
>
>
> On Tue, Nov 14, 2017 at 2:05 PM, Dennis Shea <shea at ucar.edu> wrote:
>
>> NCL: 6.5.0 -14Nov2017_0149
>>
>> (available on yellowstone)
>>
>> Produced the attached dump.
>>
>> It has 5 variables
>>
>> float AOTK_P48_L10_GLL0_A62000
>> float AOTK_P48_L10_GLL0_A62001
>> float AOTK_P48_L10_GLL0_A62006
>> float AOTK_P48_L10_GLL0_A62008
>> float AOTK_P48_L10_GLL0_A62009
>> float AOTK_P48_L10_GLL0_A62010
>>
>> On Tue, Nov 14, 2017 at 9:30 AM, Debasish Hazra <
>> debasish.hazra5 at gmail.com> wrote:
>>
>>> Hi,
>>>
>>> I was trying to read forecast aerosol field using ncl and compare that
>>> with WGRIB2 and could not figure out how forecast time can be used to plot
>>> fields in different time steps. .I have uploaded the file in ftp
>>> (ngac.t00z.aod.550nm.grib2)
>>>
>>> WGRIB2 of the file shows AOD of 5 species at every 3 hourly interval
>>> upto 5days.
>>> But ncl_filedump shows no time information of 5 species :
>>>
>>> AOTK_P48_L10_GLL0 ( lat_0, lon_0 ).
>>>
>>> Any help on this.
>>>
>>> Thanks.
>>> Debasish
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171114/81906ac0/attachment.html>