<div dir="ltr">Hello:<div><br></div><div>I have attached the wrf_times.file.</div><div><br></div><div>However, I added res@gsnAddCyclic = False</div><div><br></div><div>and sres@gsnAddCyclic = False, but the error remains.</div><div><br></div><div>The problem pertains not to the mapping of the terrain ("res"), but of the snowh_rate_dom variable.</div><div><br></div><div>Since, I am plotting a subdomain of terrain by setting:</div><div><br></div><div><div> res@mpFillOn = False</div><div> res@mpMaxLatF = lat_end ; specify the plot domain</div><div> res@mpMinLatF = lat_beg</div><div> res@mpMinLonF = lon_beg ;</div><div> res@mpMaxLonF = lon_end</div></div><div><br></div><div>I need to excerpt snowh_rate_dom from snowh_rate.</div><div><br></div><div>I do this as follows, and this is where the error message pertains to.</div><div><br></div><div>(Note, if you try to run the program, you have to read in the uploaded nc file, not wrf_output, and set ntimes = 0.)</div><div><br></div><div><div> i_loc = 1</div><div> j_loc = 1</div><div> GET_IJ::get_ij(xlat2d,xlon2d,lat_beg,lon_beg,i_loc,j_loc,dims2d(0),dims2d(1))</div><div> x_start = i_loc</div><div> y_start = j_loc</div><div><br></div><div> print("x_start = " + x_start)</div><div> print("y_start = " + y_start)</div><div> GET_IJ::get_ij(xlat2d,xlon2d,lat_end,lon_end,i_loc,j_loc,dims2d(0),dims2d(1))</div><div> x_end = i_loc</div><div> y_end = j_loc</div><div> print("x_end = " + x_end)</div><div> print("y_end = " + y_end)</div><div> snowh_rate_dom = snowh_rate(y_start:y_end,x_start:x_end)</div><div> printVarSummary(snowh_rate_dom)</div><div> printVarSummary(snowh_rate)</div><div><br></div><div> lat_dom = xlat2d(y_start:y_end,x_start:x_end)</div><div> lon_dom = xlon2d(y_start:y_end,x_start:x_end)</div><div> snowh_rate_dom!0 = "lat_dom"</div><div> snowh_rate_dom!1 = "lon_dom"</div><div> snowh_rate_dom@lat_dom = lat_dom</div><div> snowh_rate_dom@lon_dom = lon_dom</div></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jan 5, 2017 at 5:45 PM, Mary Haley <span dir="ltr"><<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-size:small">Barry,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I was unable to run the script because I was missing an ASCII file.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I think the problem is that NCL is trying to add a longitude cyclic point, and you don't want to do this for WRF data. What's happening is that NCL is taking your 270 x 270 data array and turning it into a 270 x 271 array, but the lat2d and lon2d are untouched (because gsnAddCyclic doesn't touch these), and hence you end up with two arrays of different sizes.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Try setting:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"> res@gsnAddCyclic = False</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"><div><div class="h5">On Thu, Jan 5, 2017 at 1:42 AM, Barry Lynn <span dir="ltr"><<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr"><div>Hello:</div><div><br></div><div>I thought that I have done what was indicated on this page, but I continue to get an error in dimensioning the second variable. I am trying to plot a subset of the first variable (no problem), and then extract the same subset from the second variable. I've done it before, but</div><div>it doesn't work here.</div><div><br></div><div>I been at this about 4 or 5 hours, tried different things, including adapting the suggestions here:</div><div><br></div><div><a href="https://www.ncl.ucar.edu/Support/talk_archives/2012/2697.html" target="_blank">https://www.ncl.ucar.edu/Suppo<wbr>rt/talk_archives/2012/2697.<wbr>html</a><br></div><div><br></div><div>I uploaded an extraction of the original WRF file, in case this is need.</div><div><br></div><div>Thank you.</div><div><br></div><div>ftp> put wrfsnow.ncl</div><div>local: wrfsnow.ncl remote: wrfsnow.ncl</div><div>229 Entering Extended Passive Mode (|||49155|).</div><div>150 Ok to send data.</div><div>100% |*****************************<wbr>******************************<wbr>*********| 7521 3.06 MiB/s 00:00 ETA</div><div>226 Transfer complete.</div><div>7521 bytes sent in 00:00 (11.93 KiB/s)</div><div>ftp> put wrfsnow.out.file</div><div>local: wrfsnow.out.file remote: wrfsnow.out.file</div><div>229 Entering Extended Passive Mode (|||49161|).</div><div>150 Ok to send data.</div><div>100% |*****************************<wbr>******************************<wbr>*********| 2566 8.21 MiB/s 00:00 ETA</div><div>226 Transfer complete.</div><div>2566 bytes sent in 00:00 (6.15 KiB/s)</div><div>ftp> put <a href="http://file_for_ncltalk.nc" target="_blank">file_for_ncltalk.nc</a></div><div>local: <a href="http://file_for_ncltalk.nc" target="_blank">file_for_ncltalk.nc</a> remote: <a href="http://file_for_ncltalk.nc" target="_blank">file_for_ncltalk.nc</a></div><div>229 Entering Extended Passive Mode (|||49179|).</div><div>150 Ok to send data.</div><div>100% |*****************************<wbr>******************************<wbr>*********| 2279 KiB 278.71 KiB/s 00:00 ETA</div><div>226 Transfer complete.</div><div>2333728 bytes sent in 00:08 (260.45 KiB/s)</div><div><br></div><div><br></div>-- <br><div class="m_-601649810148298472m_5559598412581958747gmail_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 <a href="tel:(914)%20432-3108" value="+19144323108" target="_blank">914 432 3108</a><br></div></div></div>
</div>
<br></div></div>______________________________<wbr>_________________<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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</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>