<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
        {font-family:Helvetica;
        panose-1:2 11 6 4 2 2 2 2 2 4;}
@font-face
        {font-family:Helvetica;
        panose-1:2 11 6 4 2 2 2 2 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:Tahoma;
        panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri","sans-serif";
        mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
span.E-MailFormatvorlage17
        {mso-style-type:personal;
        font-family:"Calibri","sans-serif";
        color:windowtext;}
span.apple-style-span
        {mso-style-name:apple-style-span;}
span.E-MailFormatvorlage19
        {mso-style-type:personal-reply;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:70.85pt 70.85pt 2.0cm 70.85pt;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="DE" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US" style="color:#1F497D">Hi Jennifer,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:#1F497D">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.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:#1F497D">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.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:#1F497D">Unfortunately, I am not aware of any PDT 4.10 example on the web.
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:#1F497D">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.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:#1F497D">Christoph.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:#1F497D"><o:p> </o:p></span></p>
<div style="border:none;border-left:solid blue 1.5pt;padding:0cm 0cm 0cm 4.0pt">
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";mso-fareast-language:DE">Von:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";mso-fareast-language:DE"> gradsusr-bounces@gradsusr.org [mailto:gradsusr-bounces@gradsusr.org]
<b>Im Auftrag von </b>Jennifer Adams<br>
<b>Gesendet:</b> Donnerstag, 9. Oktober 2014 19:35<br>
<b>An:</b> GrADS Users Forum<br>
<b>Betreff:</b> Re: [gradsusr] gribmap with GRIB2<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">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? <o:p></o:p></p>
<div>
<p class="MsoNormal">—Jennifer<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">On Oct 9, 2014, at 3:42 AM, Gebhardt Christoph <<a href="mailto:Christoph.Gebhardt@dwd.de">Christoph.Gebhardt@dwd.de</a>> wrote:<o:p></o:p></p>
</div>
<p class="MsoNormal"><br>
<br>
<o:p></o:p></p>
<div>
<p class="MsoNormal"><span lang="EN-US">Dear all,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US">just an information:</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US">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).</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US">The template 4.10 was not accepted by gribmap.</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US">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.</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US">I could make it running with a simple change of source code in gagmap.c</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US">When setting the reference time in the section
</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US">/* For fields that are statistically processed over a time interval
</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"> e.g. averages, accumulations, extremes, et al. */</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US">I changed the 'else if'-line by adding gfld->ipdtnum==10 to</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US">else if (gfld->ipdtnum == 8 || gfld->ipdtnum==9 || gfld->ipdtnum==10 || gfld->ipdtnum==11 || gfld->ipdtnum==12) {</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US">and added the line </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US">else if (gfld->ipdtnum==10) endyr_idx=16;</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US">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</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US">Best regards,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US">Christoph.</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span style="mso-fareast-language:DE">=============================================================</span><o:p></o:p></p>
<p class="MsoNormal"><span style="mso-fareast-language:DE">Dr. Christoph Gebhardt</span><o:p></o:p></p>
<p class="MsoNormal"><span style="mso-fareast-language:DE">Forschung & Entwicklung // Research & Development</span><o:p></o:p></p>
<p class="MsoNormal"><span style="mso-fareast-language:DE">Referat Vorhersagbarkeit und Verifikation // Section Predictability & Verification</span><o:p></o:p></p>
<p class="MsoNormal"><span style="mso-fareast-language:DE"><a href="mailto:christoph.gebhardt@dwd.de">christoph.gebhardt@dwd.de</a></span><o:p></o:p></p>
<p class="MsoNormal"><span style="mso-fareast-language:DE">Tel: ++49-69-8062 2743</span><o:p></o:p></p>
<p class="MsoNormal"><span style="mso-fareast-language:DE">Fax: ++49-69-8062 3721</span><o:p></o:p></p>
<p class="MsoNormal"><span style="mso-fareast-language:DE">………………………………………………………………………….…………………………………………………..</span><o:p></o:p></p>
<p class="MsoNormal"><span style="mso-fareast-language:DE">Deutscher Wetterdienst</span><o:p></o:p></p>
<p class="MsoNormal"><span style="mso-fareast-language:DE"><a href="x-msg://16/www.dwd.de">www.dwd.de</a></span><o:p></o:p></p>
<p class="MsoNormal"><span style="mso-fareast-language:DE">Frankfurter Straße 135</span><o:p></o:p></p>
<p class="MsoNormal"><span style="mso-fareast-language:DE">63067 Offenbach am Main</span><o:p></o:p></p>
<p class="MsoNormal"><span style="mso-fareast-language:DE">Germany</span><o:p></o:p></p>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman","serif";mso-fareast-language:DE">_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a><o:p></o:p></span></p>
</div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman","serif";mso-fareast-language:DE"><o:p> </o:p></span></p>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Helvetica","sans-serif";color:black;mso-fareast-language:DE">--<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-size:9.0pt;font-family:"Helvetica","sans-serif";color:black;mso-fareast-language:DE">Jennifer M. Adams<br>
Center for Ocean-Land-Atmosphere Studies (COLA)<br>
111 Research Hall, Mail Stop 2B3<br>
George Mason University<br>
4400 University Drive<br>
Fairfax, VA 22030 <o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Helvetica","sans-serif";color:black;mso-fareast-language:DE"><o:p> </o:p></span></p>
</div>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Helvetica","sans-serif";color:black;mso-fareast-language:DE"><br>
<br>
</span><span style="font-size:12.0pt;font-family:"Times New Roman","serif";mso-fareast-language:DE"><o:p></o:p></span></p>
</div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman","serif";mso-fareast-language:DE"><o:p> </o:p></span></p>
</div>
</div>
</div>
</div>
</body>
</html>