<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; color: rgb(0, 0, 0); font-size: 14px; font-family: Calibri, sans-serif; ">
<div>Hi -</div>
<div><br>
</div>
<div>I'm using NCL version 6.1.2. I'm working with the CMIP5 CCSM4 PI simulation data. The surface U and V winds are not archived so I tried using the&nbsp;the vinth2p function to interpolate to either 850mb or 925 mb. However, when I do so the winds are very strong
 and seem to be offset. I checked the lat and lon coordinates and things seem to be fine. I'm not extrapolating (extrap = False), but this shouldn't be an issue when I interpolate to 850 mb. &nbsp;I also tried both the linear and log interpolation (interp = 1 or
 2).</div>
<div><br>
</div>
<div>The file size is large (13 Gb). Could this be the issue?</div>
<div><br>
</div>
<div>Please see the attached image file, and a&nbsp;subset of my script below.&nbsp;Thanks for any help.</div>
<div><br>
</div>
<div>Lisa</div>
<div><br>
</div>
<div>
<div>und= f1-&gt;U(:,:,:,:)</div>
<div>lat = f1-&gt;lat</div>
<div>lon = f1-&gt;lon</div>
<div>time = f1-&gt;time</div>
<div>P0u = f1-&gt;P0</div>
<div>ps = f4-&gt;PS</div>
<div><br>
</div>
<div>; for comparison with another model interpolate the model data to these</div>
<div>; 25 constant pressure levels</div>
<div>plevs = (/70.,100.,150.,200.,250.,300.,350.,400.,450.,500.,550., \</div>
<div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;600.,650.,700.,750.,800.,825.,850.,875.,900.,925.,950.,975., \</div>
<div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 990.,995.,1000./)</div>
<div>plevs@units = &quot;mb&quot;</div>
<div>plevs@positive = &quot;down&quot;</div>
<div><br>
</div>
<div>interp = 2 &nbsp;;1 = linear, 2 = log, 3 = loglog</div>
<div>extrap = False</div>
<div><br>
</div>
<div>hyam = f1-&gt;hyam</div>
<div>hybm = f1-&gt;hybm</div>
<div><br>
</div>
<div>ui = vinth2p(und,hyam,hybm,plevs,ps,interp,P0u,1,extrap)</div>
<div>printVarSummary(ui)</div>
<div><br>
</div>
<div>uil = ui(:,17,:,:)</div>
<div>copy_VarCoords(ui,uil)</div>
<div>copy_VarAtts(ui,uil)</div>
<div>printVarSummary(uil)</div>
</div>
</body>
</html>