<html><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:12px"><div dir="ltr" id="yui_3_16_0_1_1427149847145_16204"><span>Hi Matt, <br></span></div><div id="yui_3_16_0_1_1427149847145_17299" dir="ltr"><br><span></span></div><div id="yui_3_16_0_1_1427149847145_17450" dir="ltr"><span id="yui_3_16_0_1_1427149847145_17449">I may have encountered this problem before. If I recall correctly, wrf_user_intrp3d requires the plane to be floats, rather than integers.</span></div><div id="yui_3_16_0_1_1427149847145_17751" dir="ltr"><br><span id="yui_3_16_0_1_1427149847145_17449"></span></div><div id="yui_3_16_0_1_1427149847145_17752" dir="ltr"><span id="yui_3_16_0_1_1427149847145_17449">You could just change your script to </span></div><div id="yui_3_16_0_1_1427149847145_16205" class="qtdSeparateBR"><br><div id="yui_3_16_0_1_1427149847145_17757" dir="ltr"> tmp_plane = (/ 40,81 , 259,81 /)<br style="" class=""></div><div id="yui_3_16_0_1_1427149847145_18063" dir="ltr">plane = tofloat(tmp_plane)<br></div><div id="yui_3_16_0_1_1427149847145_18068" dir="ltr"> tc_plane = wrf_user_intrp3d(tc,z,"v",plane,0.,True)</div><div id="yui_3_16_0_1_1427149847145_18275" dir="ltr"><br></div><div id="yui_3_16_0_1_1427149847145_18276" dir="ltr">Hope this helps,</div><div id="yui_3_16_0_1_1427149847145_18277" dir="ltr">Dave.<br style="" class=""></div><div id="yui_3_16_0_1_1427149847145_18282" style="" class="" dir="ltr"><br style="" class=""></div></div> <div style="display: block;" id="yui_3_16_0_1_1427149847145_16210" class="yahoo_quoted"> <div id="yui_3_16_0_1_1427149847145_16209" style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 12px;"> <div id="yui_3_16_0_1_1427149847145_16208" style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 16px;">------------------------------<br><div id="yui_3_16_0_1_1427149847145_16217" class="y_msg_container"><br>Message: 2<br>Date: Mon, 23 Mar 2015 06:21:51 +0000<br>From: Matthew Fearon <<a id="yui_3_16_0_1_1427149847145_18067" ymailto="mailto:Matthew.Fearon@dri.edu" href="mailto:Matthew.Fearon@dri.edu">Matthew.Fearon@dri.edu</a>><br>Subject: [ncl-talk] diagonal cross section with wrf_user_intrp3d<br>To: "<a id="yui_3_16_0_1_1427149847145_17749" ymailto="mailto:wrf-help@ucar.edu" href="mailto:wrf-help@ucar.edu">wrf-help@ucar.edu</a>" <<a id="yui_3_16_0_1_1427149847145_18065" ymailto="mailto:wrf-help@ucar.edu" href="mailto:wrf-help@ucar.edu">wrf-help@ucar.edu</a>>, "<a id="yui_3_16_0_1_1427149847145_18066" ymailto="mailto:ncl-talk@ucar.edu" href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>"<br> <<a id="yui_3_16_0_1_1427149847145_18064" ymailto="mailto:ncl-talk@ucar.edu" href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>Message-ID:<br> <<a id="yui_3_16_0_1_1427149847145_16531" ymailto="mailto:CE379185CDDDEA4CB6796C413A518A638A2F3AD3@Ex-Reno-MB1.DRI.local" href="mailto:CE379185CDDDEA4CB6796C413A518A638A2F3AD3@Ex-Reno-MB1.DRI.local">CE379185CDDDEA4CB6796C413A518A638A2F3AD3@Ex-Reno-MB1.DRI.local</a>><br>Content-Type: text/plain; charset="iso-8859-1"<br><br>Dear WRF-help and NCL users,<br><br>Has anyone successfully plotted a cross section with wrf data<br>along a diagonal, say southeast to northwest or vice versa, using<br>wrf_user_intrp3d? North/south or east/west works fine, but from the<br>documentation, a diagonal seemed possible using the example below,<br>although I can't get it working?<br><br>Also, the plane is defined (/xstart,xend , ystart,yend/) for sure<br>rather than (/xstart,ystart , xend,yend/) instead?<br><br>thanks,<br>Matt<br><br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"<br><br> a = addfile<<a id="yui_3_16_0_1_1427149847145_17291" href="https://www.ncl.ucar.edu/Document/Functions/Built-in/addfile.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/Built-in/addfile.shtml</a>>("wrfout_d01_2000-01-24_12:00:00.nc","r")<br><br> time = 1<br> tc = wrf_user_getvar<<a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml</a>>(a,"tc",time) ; T [C]<br> z = wrf_user_getvar<<a id="yui_3_16_0_1_1427149847145_17244" href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml</a>>(a,"z",time) ; z on mass points<br><br> plane = (/ 40,81 , 259,81 /) ; approx. START x;y and END x;y point<br><br> ; Extract cross section from point A to point B, as defined in "plane"<br> ; And vertically interpolate to height coordinates ("z")<br> tc_plane = wrf_user_intrp3d(tc,z,"v",plane,0.,True)<br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a id="yui_3_16_0_1_1427149847145_17750" href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150323/4b9f2038/attachment-0001.html" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150323/4b9f2038/attachment-0001.html </a><br><br><br></div> </div> </div> </div> </div></body></html>