[ncl-talk] time information in WGRIB2 and ncl

Debasish Hazra debasish.hazra5 at gmail.com
Tue Nov 14 13:28:49 MST 2017


Thanks Dennis.

Debasish

On Tue, Nov 14, 2017 at 3:18 PM, Dennis Shea <shea at ucar.edu> wrote:

> 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/35cb2bae/attachment.html>


More information about the ncl-talk mailing list