<div dir="ltr"><div dir="ltr"><div class="gmail_default" style="font-size:small">Sorry, I forgot to mention that if you want to get contour lines for RH2, but you also want to fill in particular areas, then you need to use gsn_contour_shade. The script and image I just sent you has an example of this.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">For information and examples on gsn_contour_shade, see coneff_13.ncl at:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/coneff.shtml#ex13">http://www.ncl.ucar.edu/Applications/coneff.shtml#ex13</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div></div></div><br><div class="gmail_quote"><div dir="ltr">On Fri, Nov 16, 2018 at 10:29 AM Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div class="gmail_default" style="font-size:small">Ehsan,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Part of the problem may be that wrf_map_resources sets a bunch of resources that sets your map outlines to gray, and makes them half as thick, making them hard to see.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I've modified the same script I sent you before, but changed some resources to make the map outlines really thick and white, so you could see them better. I commented out the shapefile outlines.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">If you continue to have problems, I will need to see your script. </div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">For more information about map outlines, go to </div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/mapoutlines.shtml" target="_blank">http://www.ncl.ucar.edu/Applications/mapoutlines.shtml</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">and look at mapoutlines_7.ncl. It shows how to get different levels of map outlines in your map.</div><div class="gmail_default"><br></div><div class="gmail_default">I suggest looking at the WRF/GSN examples page more closely, as it has examples on how to customize the map outlines. In particular, look at example wrf_gsn_6.ncl:</div><div class="gmail_default"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/wrfgsn.shtml" target="_blank">http://www.ncl.ucar.edu/Applications/wrfgsn.shtml</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div><div class="gmail_default" style="font-size:small"><br></div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr">On Fri, Nov 16, 2018 at 5:25 AM Ehsan Taghizadeh <<a href="mailto:ehsantaghizadeh@yahoo.com" target="_blank">ehsantaghizadeh@yahoo.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457ydp11c48cdbyahoo-style-wrap" style="font-family:bookman old style,new york,times,serif;font-size:16px"><div></div>
<div>Dear Mary</div><div><br></div><div>I've changed a little your script, as I just want "2m RH" and "10 u,v". I've also added my shape files. However I had my two main problems, yet.</div><div>1- I couldn't have countries boundaries, in spite of adding following lines:</div><div><br></div><div><span><div style="color:rgb(0,0,0);font-family:bookman old style,new york,times,serif;font-size:16px"><span><div style="color:rgb(0,0,0);font-family:bookman old style,new york,times,serif;font-size:16px">rh_res = wrf_map_resources(f,rh_res) ; Required for setting native WRF map projection</div><div style="color:rgb(0,0,0);font-family:bookman old style,new york,times,serif;font-size:16px">(instead of: tc_res = wrf_map_resources(f,tc_res))</div><div style="color:rgb(0,0,0);font-family:bookman old style,new york,times,serif;font-size:16px"><br></div></span>;---Boundaries<span style="white-space:pre-wrap"> (</span>3 following lines Don't work???)</div><div style="color:rgb(0,0,0);font-family:bookman old style,new york,times,serif;font-size:16px"> rh_res@mpDataSetName = "Earth..4" ; 1 through 4</div><div style="color:rgb(0,0,0);font-family:bookman old style,new york,times,serif;font-size:16px"> rh_res@mpNationalLineColor = "black" ; "mpCountyLineColor" doesn't work</div><div style="color:rgb(0,0,0);font-family:bookman old style,new york,times,serif;font-size:16px"> rh_res@mpNationalLineThicknessF = 2.5 ; "mpCountyLineThicknessF" doesn't work</div><div><br></div></span></div><div><br></div><div>2- I want to draw "RH" contours for all values, and fill it just for RH>=80%.</div><div><br></div><div>Attached file is output of edited script.</div><div><br></div><div><span><div> </div></span></div><div><br></div><div>Sincerely</div><div>Ehsan</div><div><br></div>
</div><div id="m_3547087428720066326m_5592500688940257208m_2663260063975951457yahoo_quoted_3033562869" class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yahoo_quoted">
<div style="font-family:'Helvetica Neue',Helvetica,Arial,sans-serif;font-size:13px;color:#26282a">
<div>
On Thursday, November 15, 2018, 9:23:21 PM GMT+3:30, Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>> wrote:
</div>
<div><br></div>
<div><br></div>
<div><div id="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default" style="font-size:small">Ehsan,</div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default" style="font-size:small"><br clear="none"></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default" style="font-size:small">Please include ncl-talk when you have followup questions about ncl-talk, unless we specifically ask to take the conversation offline.</div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default" style="font-size:small"><br clear="none"></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default" style="font-size:small">[1] If you are attaching shapefile outlines, in general you don't also want to be drawing NCL's map outlines. So, you shouldn't need to set mpOutlineBoundarySets or mpDataSetName.</div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default" style="font-size:small"><br clear="none"></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default" style="font-size:small">If you do need to set these for some reason, then you need to set these resources for when you call one of the gsn map functions, like gsn_csm_map or gsn_csm_contour_map.</div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default" style="font-size:small"><br clear="none"></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default" style="font-size:small">[2] To just fill in the RH2 values above 80, you can give gsn_csm_contour one contour level and then provide two colors to indicate what color you want for everything less than that level, and what color you want for everything greater than or equal to it.</div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default" style="font-size:small"><br clear="none"></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default" style="font-size:small">Something like this:</div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default" style="font-size:small"><br clear="none"></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default" style="font-size:small"><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><font face="monospace, monospace"> rh_res@cnFillOn = True</font></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><font face="monospace, monospace"> rh_res@cnLevelSelectionMode = "ExplicitLevels"</font></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><font face="monospace, monospace"> rh_res@cnLevels = 80.</font></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><font face="monospace, monospace"> rh_res@cnFillColors = (/"transparent","navyblue"/)</font></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><font face="monospace, monospace"> rh_res@cnLineLabelsOn = False</font></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><font face="monospace, monospace"> rh_res@lbLabelBarOn = False</font></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><font face="monospace, monospace"><br clear="none"></font></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><font face="arial, helvetica, sans-serif">Note that this will fill in everything where RH2 >= to 80 (not RH 2 > 80). If you really need to do this for RH2 > 80, then you can use a slightly bigger value for the contour level.</font></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><font face="arial, helvetica, sans-serif"><br clear="none"></font></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><font face="arial, helvetica, sans-serif">Another way to do this is to use the "where" function:</font></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><br clear="none"></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><font face="arial, helvetica, sans-serif">rh2@_FillValue = -999</font></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><font face="arial, helvetica, sans-serif">rh2 = where(rh2.gt.80,rh2,rh2@_FillValue)</font></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><font face="arial, helvetica, sans-serif"><br clear="none"></font></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><font face="arial, helvetica, sans-serif">When you do this, however, you may notice that the RH2 contours are now blocky looking, because NCL can't give you smooth contours at locations that go up against a bunch of missing values. Also you will now get a warning:</font></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><font face="arial, helvetica, sans-serif"><br clear="none"></font></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><font face="arial, helvetica, sans-serif"></font><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default">warning:ContourPlotSetValues: Data values out of range of levels set by EXPLICITLEVELS mode</div><div><br clear="none"></div><div>because NCL is not finding any values below 80, and it's just warning you about this.</div></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><font face="arial, helvetica, sans-serif"><br clear="none"></font></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default">I've attached a sample script that shows RH2 (filled in navy blue) and the difference between using "where" and not.</div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><br clear="none"></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default">--Mary</div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default"><br clear="none"></div></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default" style="font-size:small"><br clear="none"></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default" style="font-size:small"><br clear="none"></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default" style="font-size:small"><br clear="none"></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default" style="font-size:small"><br clear="none"></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default" style="font-size:small"><br clear="none"></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_default" style="font-size:small"><br clear="none"></div></div></div></div><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491yqt5688011736" id="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491yqt81433"><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_extra"><br clear="none"><div class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_quote">On Wed, Nov 14, 2018 at 11:11 PM, Ehsan Taghizadeh <span dir="ltr"><<a rel="nofollow" shape="rect" href="mailto:ehsantaghizadeh@yahoo.com" target="_blank">ehsantaghizadeh@yahoo.com</a>></span> wrote:<br clear="none"><blockquote class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="font-family:bookman old style,new york,times,serif;font-size:16px"><div></div>
<div>Dear Mary</div><div>Thank you very much for your helpful reply. I used "<span>wrf_gsn_5.ncl"</span> from</div><div><span><a rel="nofollow" shape="rect" href="https://www.ncl.ucar.edu/Applications/wrfgsn.shtml" style="color:rgb(25,106,212);text-decoration-line:underline;font-family:Helvetica Neue,Helvetica,Arial,sans-serif" target="_blank">https://www.ncl.ucar.edu/ Applications/wrfgsn.shtml</a></span><br clear="none"></div><div>and added shapefiles to "map".</div><div>I need two <span><span style="color:rgb(0,0,0);font-family:bookman old style,new york,times,serif;font-size:16px">more </span></span>helps now, and I'll be thankful if I could have any suggestion again:</div><div><br clear="none"></div><div>1- where should I add boundaries</div><div><span></span><div>@mpOutlineBoundarySets <span style="white-space:pre-wrap"> </span>= "National"</div><div>@mpDataSetName = "Earth..4" ; 1 through 4<br clear="none"></div><div><br clear="none"></div></div><div>I've tried to add them to "res_map", but it seems doesn't work!</div><div><br clear="none"></div><div>2- I've used "2m Relative Humidity" and contoured it. How could I fill it just for RH>80% (just one color)</div><div><br clear="none"></div><div>If it helps I could attach the script. However I've attached one sample of output.</div><div><br clear="none"></div><div>Sincerely</div><div>Ehsan</div>
</div><div id="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491m_5839482278452081718yahoo_quoted_2292740708">
<div style="font-family:'Helvetica Neue',Helvetica,Arial,sans-serif;font-size:13px;color:#26282a"><span class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491">
</span><div>
On Wednesday, November 14, 2018, 8:51:23 PM GMT+3:30, Mary Haley <<a rel="nofollow" shape="rect" href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>> wrote:
</div>
<div><br clear="none"></div>
<div><br clear="none"></div>
<div><div id="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491m_5839482278452081718yiv1388694067"><div><span class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491"></span><div dir="ltr"><div dir="ltr"><div style="font-size:small">Ehsan,</div><div style="font-size:small"><br clear="none"></div><div style="font-size:small">Please see example wrf_gsn_5.ncl and wrf_nogsn_5.ncl at</div><div style="font-size:small"><br clear="none"></div><div><a rel="nofollow" shape="rect" href="https://www.ncl.ucar.edu/Applications/wrfgsn.shtml" target="_blank">https://www.ncl.ucar.edu/ Applications/wrfgsn.shtml</a><br clear="none"></div><div><br clear="none"></div><div>It shows you how to overlay several WRF plots on top of each other, using WRF plotting functions, and gsn_csm plotting functions.</div><div><br clear="none"></div><div>Adding shapefile outlines should be done to the map plot, and I believe there's an example on the same page.</div><div><br clear="none"></div><div>--Mary</div><div><br clear="none"></div></div></div><div><br clear="none"><div><span class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491">On Wed, Nov 14, 2018 at 12:49 AM, Ehsan Taghizadeh <span dir="ltr"><<a rel="nofollow" shape="rect" href="mailto:ehsantaghizadeh@yahoo.com" target="_blank">ehsantaghizadeh@yahoo.com</a>></span> wrote:<br clear="none"></span><blockquote style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div id="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491m_5839482278452081718yiv1388694067yqt98647"><div><span class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491"></span><div>Hi,</div><div>I have some basic problems about overlaying different quantities of WRF, for a few months and unfortunately, until now I couldn't solve it and I'll be thankful if I could have any help for this.</div><div>I saw "<a rel="nofollow" shape="rect" href="https://www.ncl.ucar.edu/Training/Tutorials/WRF_Users_Workshop/WRF_NCL.shtml" target="_blank">https://www.ncl.ucar.edu/ Training/Tutorials/WRF_Users_ Workshop/WRF_NCL.shtml</a>", (<span><a rel="nofollow" shape="rect" href="https://www.ncl.ucar.edu/Training/Tutorials/WRF_Users_Workshop/Scripts/wrf_demo_plot_tc_gsn_shapefiles.ncl" style="margin:0px;padding:0px;color:rgb(133,45,133);font-family:verdana,sans-serif;font-size:13.3333px" target="_blank">wrf_demo_plot_tc_gsn_ shapefiles.ncl</a></span>). It is near to what if want. I just want to add some quantities to it and overlay on it"</div><div><span>uv10 = wrf_user_getvar(f0,"uvmet10", it); <span> u at 10 m, mass point</span> </span></div><span class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491"></span><div><span><span>rh2 = wrf_user_getvar(f,"rh2",0)<span><span style="color:rgb(0,0,0);font-family:bookman old style,new york,times,serif;font-size:16px">;</span></span><span style="white-space:pre-wrap"> </span> 2m Relative Humidity (%)</span><br clear="none"></span></div><div><span><span><br clear="none"></span></span></div><div><span><span>I mean overlay of "uv10" as vector and "rh2" as contour, on "terrain" and "shapefiles".</span></span></div><div><span><span>I've written some scripts and I could send it. I want to use "gsn_csm_countour", "gsn_csm_vector", and "overlay" for related quantities.</span></span></div><div><span><span>I hear from Mary Haley "<span><span style="font-family:Helvetica Neue,Helvetica,Arial,sans-serif;font-size:small">I recommend that people use the gsn_csm_contour_map function for plotting WRF data when they need to add a lot of customization to the plot</span></span>", for that reason I pointed "gsn".</span></span></div><div><span><span><br clear="none"></span></span></div><div><span><span>I'll be thankful for any help. As I said above I have this problem for a few months, before August, when I tried to plot a cross-section of WRF quantities with overlaying them. For now, I just want vector and contour, not a cross-section.</span></span></div><div><span><span><br clear="none"></span></span></div><div><span><span>Sincerely</span></span></div><span><font color="#888888"></font></span><div><span><span>Ehsan</span></span></div><div style="font-family:bookman old style,new york,times,serif;font-size:16px"></div></div></div><span class="m_3547087428720066326m_5592500688940257208m_2663260063975951457yiv7279426491"><br clear="none">______________________________ _________________<br clear="none">
ncl-talk mailing list<br clear="none">
<a rel="nofollow" shape="rect" href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br clear="none">
List instructions, subscriber options, unsubscribe:<br clear="none">
<a rel="nofollow" shape="rect" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/ mailman/listinfo/ncl-talk</a><br clear="none">
<br clear="none"></span></blockquote></div><br clear="none"></div></div></div></div>
</div>
</div></div></blockquote></div><br clear="none"></div></div></div></div></div>
</div>
</div></div>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div>
</blockquote></div>