<div dir="ltr"><div class="gmail_default" style="font-size:small">Ce Zhu,<br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I have to be honest and say that I don&#39;t know anything these station model plots, but after looking in the wmstnm.f code, I do see that UNT only affects the temperature values, and not the precipitation values.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">You are looking at the correct Fortran file, but if it&#39;s precipitation that you want to change, then you need to look at this code:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><div class="gmail_default">C                                                                                                                </div><div class="gmail_default">C  Precipitation.                                                                                                </div><div class="gmail_default">C                                                                                                                </div><div class="gmail_default">      IF (IPRSNT(14).EQ.1 .AND. RRR.GT.0) THEN</div><div class="gmail_default">        IF (RRR .EQ. 990) THEN</div><div class="gmail_default">          CALL PLCHHQ(SYMPOS(1,14),SYMPOS(2,14),&#39;T&#39;,SIZ,0.,0.)</div><div class="gmail_default">          GO TO 75</div><div class="gmail_default">        ENDIF</div><div class="gmail_default">        IF (RRR .GT. 990) THEN</div><div class="gmail_default">          INCHES = NINT(10.*REAL(RRR-990)/25.4)</div><div class="gmail_default">        ELSE</div><div class="gmail_default">          INCHES = NINT(100.*REAL(RRR)/25.4)</div><div class="gmail_default">        ENDIF</div><div class="gmail_default">        CHR5 = &#39; &#39;</div><div class="gmail_default">        WRITE(CHR5,&#39;(I5)&#39;) INCHES</div><div class="gmail_default">        LL = WMGTLN(CHR5,LEN(CHR5),1)</div><div class="gmail_default">        CALL PLCHHQ(SYMPOS(1,14),SYMPOS(2,14),CHR5(LL:LEN(CHR5)),</div><div class="gmail_default">     +              SIZ,0.,0.)</div><div class="gmail_default">      ENDIF</div><div class="gmail_default">   75 CONTINUE</div><div class="gmail_default"><br></div><div class="gmail_default">And yes, if you change this as desired, and rebuild NCL from source code, then you should get the desired change.</div><div class="gmail_default"><br></div></div><div class="gmail_default" style="font-size:small">I recommend reading the description of NCL&#39;s station model data plotting at:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><a href="http://www.ncarg.ucar.edu/supplements/wmap/#HEADING1-151">http://www.ncarg.ucar.edu/supplements/wmap/#HEADING1-151</a></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">It states the following:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><p> The internal parameter <a href="http://www.ncarg.ucar.edu/supplements/wmap/#UNT">UNT</a> specifies whether 
imperial units or
metric units are used.  
The default is imperial units with the station model data
being coded according to the chart &quot;Explanation of the Daily Weather
Map&quot; provided by NOAA and published by the U.S. Government Printing
Office.  The coding of the metric units is in accordance with the
SYNOP Data Format (FM-12) as described at:
</p><p>
<a href="http://weather.unisys.com/wxp/Appendices/Formats/SYNOP.html">
http://weather.unisys.com/wxp/Appendices/Formats/SYNOP.html</a></p><p> One thing I don&#39;t understand, however, is that the <a href="http://weather.unisys.com">weather.unisys.com</a> page mentions precipitation in mm, but the wmstnm.f code seems to be using inches.</p><p>I don&#39;t know if this is a bug. Perhaps somebody else who is reading this email and has experience with this can shed some light?</p><p>--Mary</p><p><br></p></div><div class="gmail_default" style="font-size:small">​</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Nov 4, 2014 at 9:10 PM, Zhu CE <span dir="ltr">&lt;<a href="mailto:lovingzhuce@163.com" target="_blank">lovingzhuce@163.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="line-height:1.7;color:#000000;font-size:14px;font-family:Arial"><div style="line-height:1.7;color:#000000;font-size:14px;font-family:Arial"><br><div>Dear Mary,</div><div><br></div><div>Script I used is same as you, which I download from weather map example page.</div><div><br></div><div>Have you found that no matter &#39;UNT&#39; be set as &#39;0&#39; or &#39;1&#39;, the precipitation amount is still in <span style="line-height:1.7">imperial unit? I have showed the &#39;UNT=0&#39; and &#39;UNT=1&#39; figures in attachment. (where the precipitation is 40mm actually)</span></div><div><span style="line-height:1.7"><br></span></div><div><span style="line-height:1.7">I also read the wmstnm.f source code, and I find that the code transform the precipitation amount from WMO code into imperial unit. The source code is showed as below. My questions here are if I want to use metric unit in precipitation, can I just change the source code and rebuilt the NCL? Or are there any other ways I can try?</span></div><div><span style="line-height:1.7"><br></span></div><div><span style="line-height:1.7">Thanks a lot!</span></div><div><span style="line-height:1.7"><br></span></div><div><span style="line-height:1.7">Ce Zhu</span></div><div><span style="line-height:1.7"><br></span></div><div><span style="line-height:1.7"><br></span></div><div><div>      IF (IPRSNT(14).EQ.1 .AND. RRR.GT.0) THEN</div><div>        IF (RRR .EQ. 990) THEN</div><div>          CALL PLCHHQ(SYMPOS(1,14),SYMPOS(2,14),&#39;T&#39;,SIZ,0.,0.)       </div><div>          GO TO 75</div><div>        ENDIF</div><div>        IF (RRR .GT. 990) THEN</div><div>          INCHES = NINT(10.*REAL(RRR-990)/25.4)</div><div>        ELSE</div><div>          INCHES = NINT(100.*REAL(RRR)/25.4)</div><div>        ENDIF</div><div>        CHR5 = &#39; &#39;</div><div>        WRITE(CHR5,&#39;(I5)&#39;) INCHES</div><div>        LL = WMGTLN(CHR5,LEN(CHR5),1)</div><div>        CALL PLCHHQ(SYMPOS(1,14),SYMPOS(2,14),CHR5(LL:LEN(CHR5)),</div><div>     +              SIZ,0.,0.)       </div><div>      ENDIF</div></div><br><img src="cid:4fd597d9$1$1497e257a76$Coremail$lovingzhuce$163.com" style="width:843px;min-height:651px;border:none"><div><div class="h5"><br><br><div></div><div></div><br>At 2014-11-05 00:36:23, &quot;Mary Haley&quot; &lt;<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>&gt; wrote:<br> <blockquote style="PADDING-LEFT:1ex;MARGIN:0px 0px 0px 0.8ex;BORDER-LEFT:#ccc 1px solid"><div dir="ltr"><div class="gmail_default" style="font-size:small">Dear Ce Zhu,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">It is helpful if you include a script that illustrates the problem.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I just took a weather map example from our examples page, modified it, and tried it with &quot;UNT&quot; set to 0, and then 1. </div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I don&#39;t know anything about this procedure to be honest, but I do see that I get some different numbers when I use UNT = 0 or 1.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Please see attached script.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Nov 4, 2014 at 8:07 AM, Zhu CE <span dir="ltr">&lt;<a href="mailto:lovingzhuce@163.com" target="_blank">lovingzhuce@163.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><br><div><div style="line-height:1.7;color:#000000;font-size:14px;font-family:Arial"><div>Dear ncl-talk,</div><div><br></div><div>I tried to use <i><b>wmstnm</b></i> draw weather symbol. Function introduction says that it uses use imperial units as default, but when I set <b><i>wmsetp</i></b><i>(&quot;unt</i><i>&quot;,2)</i>, the precipitation unit is still in imperial unit. Is there something I missed? (e.g. when RRR=991, it means precipitation amount is 0.1mm, but the <i><b>wmstnm</b></i> gives zero precipitation in weather symbol.)</div><div><br></div><div>Thanks a lot.</div><div><br></div><div>Ce Zhu</div></div></div><br><br><span title="neteasefooter"><span></span></span></div></div><br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div></div></div></div><br><br><span title="neteasefooter"><span></span></span></blockquote></div><br></div>