[gradsusr] Problem reading sub-hourly WRF grib files (NCEP code table for 15/30 min records)

[Wesley Ebisuzaki - NOAA Federal]

Ivan,

"UPP follows NCEP code table for timing grib records - for 15 and 30
minute records use codes 13 and 14 respectively instead of simple
minutes. This seems to confuse GrADS (or g2ctl/gribmap, I'm not sure)."

"Yes it is fixed. I'm currently need to output 10 min intervals. So, for
full hour, 10m, 20m, 40m and 50m I get OK data. But for 30m I do not."

These quotes suggests that it is a problem with gribmap expecting
30 minute_time_units  and not being able to handle 1 30_minute_time_units.

Suggestion #1

    use alt_g2ctl/alt_gmp

    Both are based on wgrib2 which uses english rather than code table
numbers.

Suggestion #2

    convert forecast time from 1 (30 minutes) to 30 (minutes)

     wgrib2 IN.grb -if ":30 min fcst:" -set_ftime "30 min fcst" -fi -grib
OUT.grb

Wesley






On Tue, Nov 7, 2017 at 1:21 PM, Jeff Duda <jeffduda319 at gmail.com> wrote:

> What result do you get when you do
>
> wgrib2 -T (grib file)
>
> ?
>
> Look after the "D=", as the represents the time in the GRIB2 file.
>
> Jeff
>
> On Tue, Nov 7, 2017 at 12:10 PM, Ivan Toman <ivtoman at inet.hr> wrote:
>
>> Hello,
>>
>> Yes it is fixed. I'm currently need to output 10 min intervals. So, for
>> full hour, 10m, 20m, 40m and 50m I get OK data. But for 30m I do not.
>>
>> Thank you!
>> Ivan
>>
>>
>>
>>
>> On 11/07/2017 05:50 PM, Jeff Duda wrote:
>>
>> What is your time interval? If it's not fixed, that's one problem. If it
>> is fixed, then I don't see why regular templating wouldn't work unless the
>> time being written to the grib2 file is not right. But there are ways to
>> fix that.
>>
>> Jeff Duda
>>
>> On Tue, Nov 7, 2017 at 10:12 AM, Ivan Toman <ivtoman at inet.hr> wrote:
>>
>>> Hello,
>>>
>>> I'm having difficulties reading sub-hourly time records from WRF grib
>>> files in GrADS.
>>>
>>> UPP follows NCEP code table for timing grib records - for 15 and 30
>>> minute records use codes 13 and 14 respectively instead of simple
>>> minutes. This seems to confuse GrADS (or g2ctl/gribmap, I'm not sure).
>>>
>>> What I get as result is that I can read any sub-hourly record as long as
>>> it is not 15 or 30 minute record. For those, I get undefined grid
>>> instead of data.
>>>
>>> I use this workflow for postprocessing: wfrout >(UPP)> grib1 >(cnvgrib)>
>>> grib2 >(g2ctl,gribmap)> GrADS
>>>
>>> Does anybody know what is going on there?
>>>
>>> Thank you in advance
>>>
>>> Ivan Toman
>>>
>>> _______________________________________________
>>> gradsusr mailing list
>>> gradsusr at gradsusr.org
>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>
>>
>>
>>
>> --
>> Jeff Duda
>> Post-doctoral research fellow
>> University of Oklahoma School of Meteorology
>>
>>
>> _______________________________________________
>> gradsusr mailing listgradsusr at gradsusr.orghttp://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>
>
> --
> Jeff Duda
> Post-doctoral research fellow
> University of Oklahoma School of Meteorology
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20171107/351f9e0a/attachment-0001.html