<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">The input coords should be monotonically increasing, I don’t think there is any need for the output coords to be structured in any particular manner. <div class=""><br class=""><div class=""><div class="">Chances are that your file has decreasing latitudes, you can flip these easily on read in. e.g.</div><div class=""><div dir="ltr" class=""><div class=""><div class="gmail_signature"><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><span lang="en-US" class=""><div class="" style="font-size: 12.8px;"><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span><font face="Andale Mono" class="">vwnd = f->V(1,{preL:preH},{lat1:lat2:-1}<wbr class="">,{lon1:lon2})</font></div></div></span></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div><div class=""><br class=""></div>oh actually, probably just change the order of lat1 and lat2, as is, you’re specifically asking for decreasing order. </div><div class=""><div dir="ltr" class=""><div class=""><div class="gmail_signature"><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><span lang="en-US" class=""><div class="" style="font-size: 12.8px;"><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span><font face="Andale Mono" class="">vwnd = f->V(1,{preL:preH},{lat2:lat1}<wbr class="">,{lon1:lon2})</font></div></div></span></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div><br class=""></div><div class="">I would also recommend you read in at least an extra grid point as the interpolation requires the surrounding points. e.g.</div><div class="">(I don’t know how big your grid is, so for simplicity add an extra degree to all dimensions )</div><div class=""><div dir="ltr" class=""><div class=""><div class="gmail_signature"><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><span lang="en-US" class=""><div class="" style="font-size: 12.8px;"><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span><font face="Andale Mono" class="">vwnd = f->V(1,{preL:preH},{lat2-1:lat1+1}<wbr class="">,{lon1-1:lon2+1})</font></div></div></span></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">If that doesn’t work, include a printVarSummary(vwnd) in your reply to ncl-talk</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">Good luck, </div><div class=""><br class=""></div><div class="">Alan</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""><div class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; border-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-stroke-width: 0px;"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">##############################<br class="">Alan Brammer,<div class="">Post-Doc Researcher</div><div class=""><br class=""></div><div class="">Department of Atmospheric and Environmental Sciences,<br class="">University at Albany, State University of New York, Albany, NY, 12222<div class=""><a href="mailto:abrammer@albany.edu" class="">abrammer@albany.edu</a><br class="">##############################</div></div></div></span></div></div>
</div>
<br class=""><div><blockquote type="cite" class=""><div class="">On 18 May 2017, at 09:14, James Russell <<a href="mailto:jorussel@ncsu.edu" class="">jorussel@ncsu.edu</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class=""><div class="gmail_signature"><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><span lang="en-US" class=""><span style="font-size:12.8px" class="">Hi all,</span><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class="">I'm attempting to plot a cross-section with height that traverses a line of lat/lon points. The data set in is currently in pressure x lat x lon coordinates. I'm using gc_latlon to create two vectors with the lat and lon coordinates I need and then attempting to use linint2_points_Wrap to interpolate to those points. This is similar to the narr_5.ncl example on the Slices page. My issue is that linint2_points_Wrap needs both the xcoord and ycoord to be monotonically increasing while I want to plot a cross-section along a line that has latitude decreasing but longitude increasing.</div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class="">My script is below. With the current configuration, because the input latitudes are monotonically decreasing, it outputs only missing values for the result.</div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class="">Can anyone help me with this? Is there a better way to go about creating the slice? Or am I missing an obvious working around in my code?</div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class="">Thanks,</div><div style="font-size:12.8px" class="">James</div><div style="font-size:12.8px" class=""><br class=""></div><div style="font-size:12.8px" class=""><div class="">lat1 = 12.</div><div class="">lat2 = 10.</div><div class="">lon1 = -15.</div><div class="">lon2 = 0.</div><div class="">preL = 1000.</div><div class="">preH = 200.</div><div class=""><br class=""></div><div class="">pathin = "/home/james/"</div><div class="">filein = "WRFNCLPP_V.nc"</div><div class=""><br class=""></div><div class="">npts = 30</div><div class=""><br class=""></div><div class="">;*****************************<wbr class="">*****************************</div><div class="">; Create great circle vectors of lat/lon points</div><div class="">;*****************************<wbr class="">*****************************</div><div class=""><br class=""></div><div class="">pts = gc_latlon(lat1,lon1,lat2,lon2,<wbr class="">npts,-2)</div><div class=""><br class=""></div><div class="">;*****************************<wbr class="">*****************************</div><div class="">; open file and read required data</div><div class="">;*****************************<wbr class="">*****************************</div><div class=""><br class=""></div><div class="">print("Opening file and reading data")</div><div class="">f = addfile(pathin+filein,"r")</div><div class="">vwnd = f->V(1,{preL:preH},{lat1:lat2}<wbr class="">,{lon1:lon2})</div><div class=""><br class=""></div><div class="">;*****************************<wbr class="">*****************************</div><div class="">; Interpolate to points on great circle</div><div class="">;*****************************<wbr class="">*****************************</div><div class=""><br class=""></div><div class="">vwnd_xs = linint2_points_Wrap(vwnd&lon,<wbr class="">vwnd&lat,vwnd, \</div><div class=""> False,pts@gclon,pts@gclat,0)</div></div></span></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</div>
_______________________________________________<br class="">ncl-talk mailing list<br class=""><a href="mailto:ncl-talk@ucar.edu" class="">ncl-talk@ucar.edu</a><br class="">List instructions, subscriber options, unsubscribe:<br class="">http://mailman.ucar.edu/mailman/listinfo/ncl-talk<br class=""></div></blockquote></div><br class=""></div></div></div></body></html>