<div dir="ltr"><div><div>eg:<br></div><div><br>;***************************************************<br>; interpolate to pressure levels<br>;***************************************************<br> plevs = (/1000,990,980,965,950,930,900,870,840,800,750,700,600,500,400,200/)<br></div> fake = sigma   ; hyb<br></div> fake = 0          ; hya<br><div><div><br>data_plevs = vinth2p(data ,fake ,sigma ,plevs ,ps*100 ,1,1000,1,True)  <br></div><div>printVarSummary(data_plevs)<br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Apr 27, 2017 at 10:07 PM, Xi Chang <span dir="ltr">&lt;<a href="mailto:xi.chang01@gmail.com" target="_blank">xi.chang01@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi NCL team,<div><br></div><div>I am currently running some simulations with GFDL model + topography.</div><div>So, I aware that pressure level_interpolation is needed when topography is non-zero due to terrain-following sigma coordinates. I know that I can use this function from NCL &quot;vinth2p&quot; but I only have &quot;bk&quot;, &quot;pk&quot;, and &quot;ps&quot;. from the output file, but not the &quot;ak&quot;. Can you suggest me how to solve this issue?</div><div><br></div><div>Thanks</div><div>Dr. Chang</div></div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>