<div dir="ltr">Hi:<div><br></div><div>I am trying to regrid WRF.nc output to rectilinear coordinates.</div><div><br></div><div>The code stops partway through.  It is strange because the code is already curvilinear, but the regrid program seems confused because my array sizes are equal (the code is attached below)</div><div><br></div><div><div>Working on time: 2016-01-23_12:00:00</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>get_src_grid_info: source lat dims = (624,624)</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>get_src_grid_info: source lon dims = (624,624)</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>get_src_grid_info: <b>source grid type is &#39;curvilinear&#39;</b></div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>curvilinear_to_SCRIP: calculating grid corners...</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>curvilinear_to_SCRIP: no lat values are at the poles, so</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>       calculating grid corners using</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>       calc_SCRIP_corners_noboundaries...</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>calc_SCRIP_corners_noboundaries</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>     min/max original lat: 26.0474/49.7817</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>     min/max original lon: -94.1364/-58.8636</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>calc_SCRIP_corners_noboundaries</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>     min/max Extlat2d: 26.0071/49.8189</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>     min/max Extlon2d: -94.2034/-58.7966</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>calc_SCRIP_corners_noboundaries</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>     min/max ExtGridCenter_lat: 26.0272/49.8003</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>     min/max ExtGridCenter_lon: -94.1699/-58.8301</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>get_dst_grid_info: destination lat dims = (624)</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>get_dst_grid_info: destination lon dims = (624)</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>get_dst_grid_info: Error: the lat/lon arrays are both 1D and the same length.</div><div>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>   This could be a rectilinear or unstructured grid. P<b>lease set &#39;grid_type&#39;</b></div><div><b>(0)<span class="gmail-Apple-tab-span" style="white-space:pre">        </span>   to &#39;rectilinear&#39; or &#39;unstructured&#39;.</b></div><div><br></div><div>Lastly, must we really have two exact versions of the same file &quot;a&quot; and &quot;src_file?&quot; The file &quot;a&quot; is the one usually processed by ncl.  In a program I received, they were named the same, but I am not sure of the implications of this.</div><div><br></div><div><span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px">The name to use for the NetCDF file that will contain the needed grid description of the destination grid. </span><b style="margin:0px;padding:0px;color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px">Do NOT set this to the original file that contains your data! This is the name of the *new* destination description file that will get written.</b><br></div><div><b style="margin:0px;padding:0px;color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px"><br></b></div><div><b style="margin:0px;padding:0px;color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px">Thank you!</b></div>-- <br><div class="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>