<div dir="ltr"><div><div><div><div>[1] When reading an entire array, <br> Use<br> und= f1->U<br></div> Not<br> und= f1->U(:,:,:,:)<br><br></div> The former is more efficient. If you want to explicitly note that U is 4D, I'd suggest adding a commebt<br><br> und= f1->U ; 4D or (:,:,:,:)<br><br></div>[2] There is an odd mix of units for the P0, plev and ps variables.<br></div> I speculate your 'P0u' is in Pa when it should be in mb (hPa).<br><br><div>P0u = f1->P0</div>
<div>ps = f4->PS ; Pa<br><br></div><div>P0u = P0u*0.01<br></div><div>P0u@units = "mb" <br></div><br> <br><div><div><div><div><div><br></div></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Feb 27, 2015 at 12:18 PM, Lisa Murphy <span dir="ltr"><<a href="mailto:lmurphy@rsmas.miami.edu" target="_blank">lmurphy@rsmas.miami.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word;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 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. 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 subset of my script below. Thanks for any help.</div>
<div><br>
</div>
<div>Lisa</div>
<div><br>
</div>
<div>
<div>und= f1->U(:,:,:,:)</div>
<div>lat = f1->lat</div>
<div>lon = f1->lon</div>
<div>time = f1->time</div>
<div>P0u = f1->P0</div>
<div>ps = f4->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> 600.,650.,700.,750.,800.,825.,850.,875.,900.,925.,950.,975., \</div>
<div> 990.,995.,1000./)</div>
<div>plevs@units = "mb"</div>
<div>plevs@positive = "down"</div>
<div><br>
</div>
<div>interp = 2 ;1 = linear, 2 = log, 3 = loglog</div>
<div>extrap = False</div>
<div><br>
</div>
<div>hyam = f1->hyam</div>
<div>hybm = f1->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>
</div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>