[gradsusr] gribmap with GRIB2

[Gebhardt Christoph]

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
Frankfurter Straße 135
63067 Offenbach am Main
Germany

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20141009/e0d04c9f/attachment.html