[ncl-talk] Question about wgt_vertical_n

[Dennis Shea]

*wgt_vertical_n*
<http://www.ncl.ucar.edu/Document/Functions/Contributed/wgt_vertical_n.shtml>
  takes the input 'x' and 'dp' and multiplies them

Xdp  = X*dp

Then, depending upon the the iopt argument, it sum and/or averages  the
values at each grid point.

  vsum = *dim_sum_n_Wrap*( Xdp, lev_dim )      ; sum vertically [ie
integrate]

so, my guess is the 'dp' must be the issue.

I assume the 'dp' were calculated via *dpres_hybrid_ccm*
<http://www.ncl.ucar.edu/Document/Functions/Built-in/dpres_hybrid_ccm.shtml>
or *dpres_plevel*
<http://www.ncl.ucar.edu/Document/Functions/Built-in/dpres_plevel.shtml>

I speculate the issue may be there.

If you can provide an example, like those provide in the examples.



On Wed, Sep 15, 2021 at 10:37 AM Joseph Clark via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> 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?
>
> Thanks!
>
> Joseph
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210916/d3668300/attachment-0001.html>