<div dir="ltr">Hi:<div><br></div><div>I think you meant:</div><div><br></div><div><div>; snowh_rate_dom = snowh_rate(y_start:y_end,x_start:x_end)</div><div> snowh_rate_dom = (/snowh_rate(y_start:y_end,x_start:x_end<b><font size="4">)/)</font></b></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_<b>dom</b>@lat2d= lat_dom</div><div> snowh_rate_<b>dom</b>@lon2d = lon_dom</div></div><div><br></div><div><font size="4">Thank you!</font> -- I had adapted this code from a rectangular grid, and that's why it didn't work (I needed to strip out the old curvelinear coordinates, as you noted). </div><div><br></div><div>I will try to attach both for future reference, after making some additional changes to the placement/size of the maps.</div><div><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jan 6, 2017 at 2:21 AM, 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 style="font-size:small">Thanks for cleaning up the script.</div><div style="font-size:small"><br></div><div style="font-size:small">I was able to get it to run and I see what the problem is.</div><div style="font-size:small"><br></div><div style="font-size:small">You are subscripting snowh_rate with:</div><div style="font-size:small"><br></div><div><span class=""><div style="font-size:small"><font face="monospace, monospace"> snowh_rate_dom = snowh_rate(y_start:y_end,x_<wbr>start:x_end)</font></div><div style="font-size:small"><br></div></span><div style="font-size:small">and then trying to attach coordinate arrays to it with:</div><span class=""><div style="font-size:small"><font face="monospace, monospace"><br></font></div><div style="font-size:small"><div><font face="monospace, monospace"> lat_dom = xlat2d(y_start:y_end,x_start:<wbr>x_end)</font></div><div><font face="monospace, monospace"> lon_dom = xlon2d(y_start:y_end,x_start:<wbr>x_end)</font></div></div></span><div><span class=""><div style="font-size:small"><font face="monospace, monospace"> snowh_rate_dom!0 = "lat_dom"</font></div><div style="font-size:small"><font face="monospace, monospace"> snowh_rate_dom!1 = "lon_dom"</font></div><div style="font-size:small"><font face="monospace, monospace"> snowh_rate_dom@lat_dom = lat_dom</font></div><div style="font-size:small"><font face="monospace, monospace"> snowh_rate_dom@lon_dom = lon_dom</font></div><div style="font-size:small"><br></div></span><div style="font-size:small">This is not going to work, because snowh_rate_dom is still a curvilinear grid (i.e. a grid defined by 2D lat / lon coordinates), so you still need to do the "@lat2d" and "@lon2d" attachments.</div><div style="font-size:small"><br></div><div style="font-size:small">Furthermore, when when you do this:<br></div><div style="font-size:small"><br></div><div><span class=""><div style="font-size:small"><font face="monospace, monospace"> snowh_rate_dom = snowh_rate(y_start:y_end,x_<wbr>start:x_end)</font></div></span><div style="font-size:small"><br>the original lat2d and lon2d attributes from snowh_rate are copied over, and you end up with a lat2d and lon2d that are the wrong size. That's why you are getting these two errors:</div><div style="font-size:small"><br></div><div><span class=""><div><font face="monospace, monospace">0)<span class="m_-9206883669308222897gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>is_valid_latlon2d_attr: Warning: The 'lat2d' attribute must either be the same dimension sizes as the data, or one element larger in both directions. Your data will most likely not be overlaid on the map correctly.</font></div></span><span class=""><div><font face="monospace, monospace">(0)<span class="m_-9206883669308222897gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>is_valid_latlon2d_attr: Warning: The 'lon2d' attribute must either be the same dimension sizes as the data, or one element larger in both directions. Your data will most likely not be overlaid on the map correctly.</font></div></span></div><div><br></div><div>You almost have the right idea. </div><div><br></div><div>First, when you subscript snowh_rate, make sure the lat2d and lon2d don't get copied over by using the special (/ and /) syntax:</div><div><br></div><div><div><font face="monospace, monospace"> snowh_rate_dom = (/snowh_rate(y_start:y_end,x_<wbr>start:x_end/)) ; enclosing a variable with (/ and /) will strip all metadata</font></div></div><div><br></div><div>Now you need to attach new lat2d and lon2d to snowh_rate_dom that are subscripted the same way. You already have this with lat_dom and lon_dom, so just use those:</div><div><br></div><div><font face="monospace, monospace">snowh_rate@lat2d = lat_dom</font></div><div><font face="monospace, monospace">snowh_rate@lon2d = lon_dom</font></div><div><br></div><div>Make sure you remove these lines:</div><div><font face="monospace, monospace"><br></font></div><div><span class=""><div><font face="monospace, monospace"> snowh_rate_dom!0 = "lat_dom"</font></div><div><font face="monospace, monospace"> snowh_rate_dom!1 = "lon_dom"</font></div><div><font face="monospace, monospace"> snowh_rate_dom@lat_dom = lat_dom</font></div><div><font face="monospace, monospace"> snowh_rate_dom@lon_dom = lon_dom</font></div><div><br></div></span><span class="HOEnZb"><font color="#888888"><div>--Mary</div><div><br></div><div><br></div></font></span></div><div><br></div><div><br></div></div></div><div style="font-size:small"><br></div></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jan 5, 2017 at 10:26 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><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Mary:<div><br></div><div>This reads the files that I previously attached.</div><div><br></div><div>I left the wrf calls, but put them in a logical if statement.</div><div><br></div><div>I did not remove the code that sets the label bars, makes nice maps.</div><div><br></div><div>If I need to do this, please let me know.</div><div><br></div><div>The code now gives these errors (below), but I think you can ignore them.</div><div>(I am not sure how to manually add this information, for this "clean" version.</div><div><br></div><div>It still gives the one about the dimensions.</div><div><br></div><div>(Thanks.)</div><div><br></div><div><div>0)<span class="m_-9206883669308222897m_-8541097583354215696gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>is_valid_latlon2d_attr: Warning: The 'lat2d' attribute must either be the same dimension sizes as the data, or one element larger in both directions. Your data will most likely not be overlaid on the map correctly.</div><div>(0)<span class="m_-9206883669308222897m_-8541097583354215696gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>check_for_y_lat_coord: Warning: Data either does not contain a valid latitude coordinate array or doesn't contain one at all.</div><div>(0)<span class="m_-9206883669308222897m_-8541097583354215696gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>A valid latitude coordinate array should have a 'units' attribute equal to one of the following values: </div><div>(0)<span class="m_-9206883669308222897m_-8541097583354215696gmail-Apple-tab-span" style="white-space:pre-wrap"> </span> 'degrees_north' 'degrees-north' 'degree_north' 'degrees north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'</div><div>(0)<span class="m_-9206883669308222897m_-8541097583354215696gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>is_valid_latlon2d_attr: Warning: The 'lon2d' attribute must either be the same dimension sizes as the data, or one element larger in both directions. Your data will most likely not be overlaid on the map correctly.</div><div>(0)<span class="m_-9206883669308222897m_-8541097583354215696gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>check_for_lon_coord: Warning: Data either does not contain a valid longitude coordinate array or doesn't contain one at all.</div><div>(0)<span class="m_-9206883669308222897m_-8541097583354215696gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>A valid longitude coordinate array should have a 'units' attribute equal to one of the following values: </div><div>(0)<span class="m_-9206883669308222897m_-8541097583354215696gmail-Apple-tab-span" style="white-space:pre-wrap"> </span> 'degrees_east' 'degrees-east' 'degree_east' 'degrees east' 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' '</div></div></div><div class="m_-9206883669308222897HOEnZb"><div class="m_-9206883669308222897h5"><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jan 5, 2017 at 7:11 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">Barry,<br><br>Offline for a bit. I still can't run your script because it's assuming that I have a WRF output file, but the "<a href="http://file_for_ncltalk.nc" target="_blank">file_for_ncltalk.nc</a>" file you gave me is not a WRF output file. I<div style="font-size:small;display:inline">t's very time-consuming for me to try to figure out how to change your script to work with this file, so please, send me a *clean* script that works with the files you've sent me so far.</div><div><div style="font-size:small;display:inline"><br></div></div><div><div style="font-size:small;display:inline">Thanks,</div></div><div><div style="font-size:small;display:inline"><br></div></div><div><div style="font-size:small;display:inline">--Mary</div></div><div><div style="font-size:small;display:inline"><br></div></div><div><div style="font-size:small;display:inline"></div>
</div></div><div class="m_-9206883669308222897m_-8541097583354215696HOEnZb"><div class="m_-9206883669308222897m_-8541097583354215696h5"><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jan 5, 2017 at 9:58 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><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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,l<wbr>at_beg,lon_beg,i_loc,j_loc,dim<wbr>s2d(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,l<wbr>at_end,lon_end,i_loc,j_loc,dim<wbr>s2d(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_sta<wbr>rt:x_end)</div><div> printVarSummary(snowh_rate_dom<wbr>)</div><div> printVarSummary(snowh_rate)</div><div><br></div><div> lat_dom = xlat2d(y_start:y_end,x_start:x<wbr>_end)</div><div> lon_dom = xlon2d(y_start:y_end,x_start:x<wbr>_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="m_-9206883669308222897m_-8541097583354215696m_-4294232288934853161HOEnZb"><div class="m_-9206883669308222897m_-8541097583354215696m_-4294232288934853161h5"><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 style="font-size:small">Barry,</div><div style="font-size:small"><br></div><div style="font-size:small">I was unable to run the script because I was missing an ASCII file.</div><div style="font-size:small"><br></div><div 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 style="font-size:small"><br></div><div style="font-size:small">Try setting:</div><div style="font-size:small"><br></div><div style="font-size:small"> res@gsnAddCyclic = False</div><div style="font-size:small"><br></div><div style="font-size:small">--Mary</div><div style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="m_-9206883669308222897m_-8541097583354215696m_-4294232288934853161m_-8642245181990319030h5">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="m_-9206883669308222897m_-8541097583354215696m_-4294232288934853161m_-8642245181990319030h5"><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.htm<wbr>l</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_-9206883669308222897m_-8541097583354215696m_-4294232288934853161m_-8642245181990319030m_-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="m_-9206883669308222897m_-8541097583354215696m_-4294232288934853161m_-8642245181990319030gmail_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 <a href="tel:(914)%20432-3108" value="+19144323108" target="_blank">914 432 3108</a><br></div></div></div>
</div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="m_-9206883669308222897m_-8541097583354215696gmail_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 <a href="tel:(914)%20432-3108" value="+19144323108" target="_blank">914 432 3108</a><br></div></div></div>
</div>
</div></div></blockquote></div><br></div>
</div></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></div>