evaporation: Evaporation estimation

In brasmus/preciparea: Precipitation area ananlysis.

Description

Processing routines applied to ERAINT reanalysis to estiamte global evaporation and atmospheric moisture content

Usage

evaporation(fname = "~/Downloads/ERAINT_t2m.nc", R.s = 18.02)
moistflux(fname = "~/Downloads/ERAINT_Qs.nc")
columnH2O(fname = "~/Downloads/ERAINT_Q-columntotal.nc")

Arguments

fname	File name
R.s	specific gas constant for water vapour

Value

time series objects Es, Rho, Qs, and H2O.

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (fname = "~/Downloads/ERAINT_t2m.nc", R.s = 18.02) 
{
    t2m <- retrieve(fname)
    es <- C.C.eq(mask(t2m, land = TRUE))
    rho <- es/(R.s * (t2m + 273.15))
    es.50S50N <- aggregate.area(subset(es, is = list(lon = range(lon(es)), 
        lat = c(-50, 50))), FUN = "mean")
    es.90S50S <- aggregate.area(subset(es, is = list(lon = range(lon(es)), 
        lat = c(-90, -50))), FUN = "mean")
    es.50N90N <- aggregate.area(subset(es, is = list(lon = range(lon(es)), 
        lat = c(50, 90))), FUN = "mean")
    rho.50S50N <- aggregate.area(subset(rho, is = list(lon = range(lon(es)), 
        lat = c(-50, 50))))
    rho.90S50S <- aggregate.area(subset(rho, is = list(lon = range(lon(es)), 
        lat = c(-90, -50))))
    rho.50N90N <- aggregate.area(subset(rho, is = list(lon = range(lon(es)), 
        lat = c(50, 90))))
    Es <- combine.stations(es.50S50N, es.90S50S, es.50N90N)
    Rho <- combine.stations(rho.50S50N, rho.90S50S, rho.50N90N)
    save(Es, file = "Es.rda")
    save(Rho, file = "Rho.rda")
  }

brasmus/preciparea documentation built on May 7, 2019, 11:11 a.m.