<div dir="ltr">Hi:<div><br></div><div>WRF zoom has its own way of doing things.</div><div><br></div><div>If you deviate, you can run into trouble.</div><div><br></div><div>
<p class="gmail-p1"><span class="gmail-s1">; bounding box for January 26 2015</span></p>
<p class="gmail-p1"><span class="gmail-s1"> lats = (/ 37.50, 46.50 /)</span></p>
<p class="gmail-p1"><span class="gmail-s1"> lons = (/ -77.0, -64.50 /)</span></p><p class="gmail-p1"><span class="gmail-s1"><br></span></p><p class="gmail-p1"><span class="gmail-s1"> </span><span class="gmail-s2">mpres</span><span class="gmail-s1"> = True ; Set map options (none set here)</span></p><p class="gmail-p1"><span class="gmail-s1"> </span><span class="gmail-s2">mpres</span><span class="gmail-s1">@mpOutlineOn = True</span></p><p class="gmail-p1"><span class="gmail-s1"> </span><span class="gmail-s2">mpres</span><span class="gmail-s1">@mpOutlineBoundarySets = "USStates"</span></p><p class="gmail-p1"><span class="gmail-s1"> </span><span class="gmail-s2">mpres</span><span class="gmail-s1">@mpUSStateLineColor = "Black"</span></p><p class="gmail-p1"><span class="gmail-s1"> </span><span class="gmail-s2">mpres</span><span class="gmail-s1">@mpNationalLineColor = "Black"</span></p><p class="gmail-p1"><span class="gmail-s1"> </span><span class="gmail-s2">mpres</span><span class="gmail-s1">@mpGeophysicalLineColor = "Black"</span></p><p class="gmail-p1"><span class="gmail-s1"> </span><span class="gmail-s2">mpres</span><span class="gmail-s1">@mpUSStateLineThicknessF = 3.0</span></p><p class="gmail-p1"><span class="gmail-s1"> </span><span class="gmail-s2">mpres</span><span class="gmail-s1">@mpNationalLineThicknessF = 3.0</span></p><p class="gmail-p1"><span class="gmail-s1"> </span><span class="gmail-s2">mpres</span><span class="gmail-s1">@mpGeophysicalLineThicknessF = 3.0</span></p><p class="gmail-p1"><span class="gmail-s1"> </span><span class="gmail-s2">mpres</span><span class="gmail-s1">@mpGridAndLimbOn = False</span></p><p class="gmail-p1"><span class="gmail-s1">; </span><span class="gmail-s2">mpres</span><span class="gmail-s1">@mpOutlineBoundarySets = "AllBoundaries"</span></p><p class="gmail-p1"><span class="gmail-s1">;<b> loc(0,;) is west-east (x) ; loc(1,:) is south-north (y)</b></span></p><p class="gmail-p1"><span class="gmail-s1"><b>; subtract one since we want to use it as an index in NCL</b></span></p><p class="gmail-p1"><span class="gmail-s1"><b> loc = wrf_user_ll_to_ij(a, lons, lats, True)</b></span></p><p class="gmail-p1"><span class="gmail-s1"><b> x_start = loc(0,0) - 1</b></span></p><p class="gmail-p1"><span class="gmail-s1"><b> x_end = loc(0,1) - 1</b></span></p><p class="gmail-p1"><span class="gmail-s1"><b> y_start = loc(1,0) - 1</b></span></p><p class="gmail-p1"><span class="gmail-s1"><b> y_end = loc(1,1) - 1</b></span></p><p class="gmail-p2"><b><span class="gmail-s1"></span><br></b></p><p class="gmail-p2"><b><span class="gmail-s1"></span><br></b></p><p class="gmail-p1"><b><span class="gmail-s1"> </span><span class="gmail-s2">mpres</span><span class="gmail-s1">@ZoomIn = True</span></b></p><p class="gmail-p1"><b><span class="gmail-s1"> </span><span class="gmail-s2">mpres</span><span class="gmail-s1">@Xstart = x_start</span></b></p><p class="gmail-p1"><b><span class="gmail-s1"> </span><span class="gmail-s2">mpres</span><span class="gmail-s1">@Ystart = y_start</span></b></p><p class="gmail-p1"><b><span class="gmail-s1"> </span><span class="gmail-s2">mpres</span><span class="gmail-s1">@Xend = x_end</span></b></p><p class="gmail-p1"><span class="gmail-s1"><b>
</b></span></p><p class="gmail-p1"><b><span class="gmail-s1"> </span><span class="gmail-s2">mpres</span><span class="gmail-s1">@Yend = y_end</span></b></p><p class="gmail-p1"><b><span class="gmail-s1"><br></span></b></p><p class="gmail-p1"><span class="gmail-s1"> var_</span><span class="gmail-s2">zoom</span><span class="gmail-s1"> = avo_ave_new(y_start:y_end,x_start:x_end)</span></p><p class="gmail-p1">
</p><p class="gmail-p1"><span class="gmail-s1"> contour_avo = wrf_contour(a,wks,var_</span><span class="gmail-s2">zoom</span><span class="gmail-s1">,opts_avo)</span></p><p class="gmail-p1"><span class="gmail-s1"><br></span></p><p class="gmail-p1"><span class="gmail-s1">
</span></p><p class="gmail-p1"><span class="gmail-s1"> plot = wrf_map_overlays(a,wks,(/contour_avo,contour_z/),opts_avo,</span><span class="gmail-s2">mpres</span></p></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 1, 2017 at 2:53 PM, Luis Cesareo Cana Cascallar <span dir="ltr"><<a href="mailto:luis.cana@ulpgc.es" target="_blank">luis.cana@ulpgc.es</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word">
Dear all,
<div>I’m trying to redraw a plot with maps, contours and vectors zooming in a specific area. To do so, I’ve set the following parameters:</div>
<div><br>
</div>
<div>
<div>; mpres@mpLimitMode = "Corners" </div>
<div>; mpres@mpLeftCornerLatF = 28.0 </div>
<div>; mpres@mpLeftCornerLonF = -13.65 </div>
<div>; mpres@mpRightCornerLatF = 28.8 </div>
<div>; mpres@mpRightCornerLonF = -14.55</div>
</div>
<div><br>
</div>
<div>The map looks nice, but the variables are not drawn accordingly to this new zoomed map. NCL overlay the original data over the zoomed map because because that’s what I’m telling it to do...</div>
<div><br>
</div>
<div>I’ve taken a look to the NCL WRF tutorial about zooming in a plot, but they took another approach, locating (x,y) coordinates from the given corners and creating new zoomed variables. I think this method may be a bit outdated. I’ve also googled
in ncl-talk archives about this topic and I found an old mail (Jan 9, 2011) with a similar problem <span style="background-color:rgb(255,255,255)">but it required a mod to WRFUserARW.ncl</span><span style="background-color:rgb(255,255,255)"> in
order to add special resources. As it’s an old mail, is still valid? </span></div>
<div><br>
</div>
<div>Is it possible to set the above limits to the contour and vector plots just by setting the left and right corner as above? If not, which is the best approach to this problem?</div>
<div><br>
</div>
<div>Best,</div>
<div><br>
</div>
<div>Luis</div>
<div><br>
</div>
<div><br>
<div><span><img height="111" width="446" id="m_3370875791777812057D006734A-B335-4714-A829-EDB7F7283221" src="cid:B00E1AE2-614B-461B-9561-42FA796CA6DF@ccbb.ulpgc.es"></span>
</div>
<br>
</div>
</div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">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/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior Lecturer,</div><div><div><span style="color:rgb(136,136,136)">The Institute of the Earth Science, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">The Hebrew University of Jerusalem, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Givat Ram, Jerusalem 91904, Israel </span><br style="color:rgb(136,136,136)"></div><span style="color:rgb(136,136,136)">Tel: 972 547 231 170</span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Fax: (972)-25662581</span></div></div><div><span style="color:rgb(136,136,136)"><br></span></div><div>C.E.O, Weather It Is, LTD<br>Weather and Climate Focus<br><a href="http://weather-it-is.com" target="_blank">http://weather-it-is.com</a><br>Jerusalem, Israel<br>Local: 02 930 9525<br>Cell: 054 7 231 170<br>Int-IS: x972 2 930 9525<br>US 914 432 3108<br></div></div></div>
</div>