A vertical integral is just a weighted sum.
vertical_integral = SUM( LH*dz )
Please read:
http://www.ncl.ucar.edu/Document/Functions/Contributed/dim_sum_wgt_Wrap.shtml
http://www.ncl.ucar.edu/Document/Functions/Built-in/dim_sum_wgt.shtml
dz = (/0.5, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0 \
,1.0, 1.0, 1.0, 1.5, 2.0, 2.0, 1.0/)
If these are km you should multiply by 1000 for meters
LHsum = dim_sum_wgt_Wrap( LH(time|:,lat|:,lon|:,lev|:), dz, 0 )
Good luck
Liang Guo wrote:
> Dear All,
>
> I have a dataset which has 4 dimension LH(time,lev,lat,lon). The lev
> dimension use height coordinate, which contains 14 levels (/ 0.0, 1.0,
> 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0, 12.0, 14.0, 16.0/).
>
> Which function I can use to calculate vertical integral. Thanks.
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sat Sep 15 2007 - 10:59:35 MDT
This archive was generated by hypermail 2.2.0 : Mon Sep 17 2007 - 14:58:21 MDT