[ncl-talk] ncl_filedump question

Dennis Shea shea at ucar.edu
Thu Jul 19 21:13:42 MDT 2018


Can you please send ncl-talk the output from: 

%> uname  -a

%> gcc -V

--- yes. I think upgrading to 6.5.0 would. E a good idea. 

D

Sent from my iPhone

> On Jul 19, 2018, at 8:18 PM, Jim Means <jim at weatherextreme.com> wrote:
> 
> Thank you Dennis.  Unfortunately, I still just get the error message, so apparently it's something peculiar to my machine.  I wonder if there could possibly be something wrong in the way I've configured NCL? 
> I could update to the latest version of NCL.
> 
> Jim
> 
>> On 07/19/2018 06:53 PM, Dennis Shea wrote:
>> Jim
>> 
>> NCL 6.3.0 handles the WRF grib2 file just fine [well,  99%]. 
>> 
>> %> ncl_filedump 1807190600_wrfout_arw_d02.grb2  | less
>> 
>> There is only one initial time  and one forecast time on the file. If you want single dimension included in the dimensioning:
>> 
>> %> ncl_filedump 1807190600_wrfout_arw_d02.grb2 -itime -ftime  | less
>> 
>> %> ncl_filedump 1807190600_wrfout_arw_d02.grb2 -itime -ftime   >&!   out.630_itime_ftime.1807190600_wrfout_arw_d02
>> 
>> Within an NCL script:   https://www.ncl.ucar.edu/Document/Functions/Built-in/setfileoption.shtml
>> 
>>    setfileoption("grb","SingleElementDimensions","Initial_time")
>> or
>>    setfileoption("grb","SingleElementDimensions", (/"Initial_time","Forecast_time"/))
>> 
>> 
>> ===
>> I have also tried the recently released 6.5.0
>> 
>> You get the exact same result.
>> 
>> ===
>> NOTE:
>> 
>> 
>> [a]  Some of the descriptive attributes:
>> 
>>       float REFC_P8_L10_GLC0_max30min ( initial_time0_hours, ygrid_0, xgrid_0 )
>>          center :       US National Weather Service - NCEP (WMC)
>>          production_status :    Operational products
>>          long_name :    Composite reflectivity
>>          units :        dB
>>          _FillValue :   1e+20
>>          coordinates :  gridlat_0 gridlon_0
>>          grid_type :    Lambert Conformal can be secant or tangent, conical or bipolar
>>          parameter_discipline_and_category :              Meteorological products, Forecast radar imagery
>>          parameter_template_discipline_category_number :        ( 8, 0, 16, 5 )
>>          level_type :   Entire Atmosphere
>>          level :         0
>>          type_of_statistical_processing :       Maximum
>>          statistical_process_duration : 30 minutes (ending at forecast time)
>>          forecast_time :        780
>>          forecast_time_units :  minutes
>> 
>> [b] Warning messages
>> 
>> The variables are unpacked correctly!!
>> 
>> NCL uses look-up tables to match-up parameter id numbers with information such as lon_name, units, etc.
>> Unfortunately, GRIB is a moving target and NCL developers are not privy to information or it is just not publicly available.
>> 
>> So, for example:
>> 
>> warning: Entry (243) not found in code table file /Users/shea/devel/ncl/lib/ncarg/grib2_codetables/ncep/4/4.2.0.1.table
>> NCL unpacks the variables and creates the appropriate dimension names. It just does not have the needed lookup info.
>> 
>>       float VAR_0_1_243_P0_L1_GLC0 ( initial_time0_hours, ygrid_0, xgrid_0 )
>>          center :       US National Weather Service - NCEP (WMC)
>>          production_status :    Operational products
>>          long_name :    unknown variable name
>>          units :        unknown
>>          _FillValue :   1e+20
>>          coordinates :  gridlat_0 gridlon_0
>>          grid_type :    Lambert Conformal can be secant or tangent, conical or bipolar
>>          parameter_discipline_and_category :    Meteorological products, Moisture
>>          parameter_template_discipline_category_number :        ( 0, 0, 1, 243 )
>>          level_type :   Ground or water surface
>>          level :         0
>>          forecast_time :        780
>>          forecast_time_units :  minutes
>> 
>> HTH
>> D
>> 
>> 
>>> On Thu, Jul 19, 2018 at 2:26 PM, Jim Means <jim at weatherextreme.com> wrote:
>>> Hello,
>>> 
>>> I am trying to find out what's inside a grib2 file created by WRF, and I'm attempting to use ncl_filedump to find out the details. I'm running ncl_filedump with no options:
>>> 
>>> ncl_filedump 1807190600_wrfout_arw_d02.grb2
>>> and I get
>>> 
>>> fatal:NclGRIB2: Invalid Product Definition Template.
>>> 
>>> Variable: f
>>> Type: file
>>> (0) File Missing Value : -1
>>> 
>>> What am I doing wrong?
>>> 
>>> 
>>> Thanks,
>>> 
>>> Jim
>>> 
>>> -- 
>>> James D. Means, Ph.D.
>>> Senior Atmospheric & Climate Scientist
>>> California Office
>>> Tele: 619-495-1638 | Fax: 775-636-8430
>>> 930 Tahoe Blvd., Suite 802-560
>>> Incline Village, Nevada 89451
>>> jim at weatherextreme.com | vcard
>>> www.weatherextreme.com
>>> 
>>>  
>>> 
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>> 
>> 
> 
> 
> _______________________________________________
> 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/20180719/0b3fb13e/attachment.html>


More information about the ncl-talk mailing list