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

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?