vibeta
Function for doing vertical integration using beta factors.
Synopsis
function vibeta(
p[*] : float,
x : float,
linlog[1] : integer,
psfc : float,
pbot[1] : float,
ptop[1] : float
)
Arguments
- p
- pressure levels (input, 1d array of at least length 3)
- x
- quantity to be integrated (input, array of any dimensionality,
last dimension must be same length as p). Missing values
should be indicated by x@_FillValue. If x@_FillValue is
not set, then the NCL default (-999) will be assumed.
- linlog
- scalar, linear (1) or log (2) interpolation
- psfc
- surface pressure (must be same dimensions as x, minus the last
dimension)
- pbot, ptop
- scalar, lower and upper limits of integration
Description
Function for doing vertical integration using Beta factors
(references: Boer (1982) Mon. Wea. Rev pp. 1801-20 and Trenberth (1992)
NCAR Tech Note 373). The return value is the same size as
psfc.
Example
begin
nlvl = 17
linlog = 1
p = (/ 1000.,925.,850.,700.,600.,500., \
400.,300.,250.,200.,150.,100., \
70.,50.,30.,20.,10. /)
psfc = 1013.
pbot = 1100.
ptop = p(nlvl-1)
plvcrt = p(nlvl-1)
x = new ( nlvl ,float)
do nl=0,nlvl-1
x(nl) = nlvl-nl
end do
vint = vibeta (p,x,linlog,psfc,pbot,ptop)
end
Reference Manual Control Panel
NG4.1 Home, Index, Examples, Glossary, Feedback, Ref Contents, Ref WhereAmI?
$Revision: 1.10 $ $Date: 1999/03/22 20:49:11 $