[gradsusr] gribmap with GRIB2

[Gebhardt Christoph]

Hi Jennifer,

yes, you are right,  PDT 4.10 8and 4.6 as well) have the percentile value as meta datum which might need some treatment as the upper/lower limits in the probability PDTs.
The lines I added to gagmap.c works but it doesn't take into account the percentile values which will be needed if more than one percentile value is valid within the gribfile.

Unfortunately, I am not aware of any PDT 4.10 example on the web.
Is it possible to send you one of my GRIB2 files e.g. via e-mail (but not on GrADS User forum)? A file with two records to check different percentile values would require around 800kB.

Christoph.


Von: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] Im Auftrag von Jennifer Adams
Gesendet: Donnerstag, 9. Oktober 2014 19:35
An: GrADS Users Forum
Betreff: Re: [gradsusr] gribmap with GRIB2

Support for PDT 4.10 will probably require some additional grib codes, something similar to the handling of the probability forecasts (templates 4.5 and 4.9). Can you point me to a data sample?
-Jennifer

On Oct 9, 2014, at 3:42 AM, Gebhardt Christoph <Christoph.Gebhardt at dwd.de<mailto:Christoph.Gebhardt at dwd.de>> wrote:


Dear all,

just an information:

I tried to use gribmap with GRIB2 input and Product definition template 4.10 (Percentile forecasts at a horizontal level or in a horizontal layer on a continuous or non-continuous time interval).
The template 4.10 was not accepted by gribmap.

It seems to me that the current implementation in GrADS versions 2.0.2 and 2.1.a2 (I did not check 2.1.a1) is only partially prepared for PDT 4.10.
I could make it running with a simple change of source code in gagmap.c

When setting the reference time in the section
/* For fields that are statistically processed over a time interval
       e.g. averages, accumulations, extremes, et al. */
I changed the 'else if'-line by adding gfld->ipdtnum==10 to

else if (gfld->ipdtnum == 8 || gfld->ipdtnum==9 || gfld->ipdtnum==10 || gfld->ipdtnum==11 || gfld->ipdtnum==12)  {

and added the line

else if (gfld->ipdtnum==10)  endyr_idx=16;

Later on in the source code in gaint g2sp (Gets the statistical process used to derive a variable.), the PDT 4.10 has been correctly take into account in the standard source code of gagmap.c

Best regards,
Christoph.




=============================================================
Dr. Christoph Gebhardt
Forschung & Entwicklung // Research & Development
Referat Vorhersagbarkeit und Verifikation // Section Predictability & Verification
christoph.gebhardt at dwd.de<mailto:christoph.gebhardt at dwd.de>
Tel: ++49-69-8062 2743
Fax: ++49-69-8062 3721
................................................................................................................................................
Deutscher Wetterdienst
www.dwd.de<x-msg://16/www.dwd.de>
Frankfurter Straße 135
63067 Offenbach am Main
Germany

_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr

--
Jennifer M. Adams
Center for Ocean-Land-Atmosphere Studies (COLA)
111 Research Hall, Mail Stop 2B3
George Mason University
4400 University Drive
Fairfax, VA 22030




-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20141013/991432c0/attachment-0001.html