<div dir="ltr"><div><a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/wgt_vertical_n.shtml"><b>wgt_vertical_n</b></a>   takes the input 'x' and 'dp' and multiplies them</div><div><br></div><div>Xdp  = X*dp         <br></div><div><br></div><div>Then, depending upon the the iopt argument, it sum and/or averages  the values at each grid point.</div><div><br></div><div>  vsum = <b>dim_sum_n_Wrap</b>( Xdp, lev_dim )      ; sum vertically [ie integrate]</div><div><br></div><div>so, my guess is the 'dp' must be the issue.</div><div><br></div><div>I assume the 'dp' were calculated via <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/dpres_hybrid_ccm.shtml"><b>dpres_hybrid_ccm</b></a>  or <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/dpres_plevel.shtml"><b>dpres_plevel</b></a></div><div><br></div><div>I speculate the issue may be there.  <br></div><div><br></div><div>If you can provide an example, like those provide in the examples.</div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Sep 15, 2021 at 10:37 AM Joseph Clark via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">I don't have a coding error, but a question: how does wgt_vertical_n deal with values beneath the surface. In other words, for a location with a surface pressure of 980 hPa, would the extrapolated value on the 1000 hPa surface be included in the vertical integral? Based on my testing of the script (by setting enormous values to a variable at the 1000 hPa level), this appears to be the case. Is there a justification for including subsurface values in the vertical integral or should the user avoid this by setting subsurface values to zero?<br><br><div>Thanks!<br><br>Joseph</div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>