ps2psl: Surface pressure reduction to sea-level

In SantanderMetGroup/convertR: A climate4R Package for Unit Conversion and Variable Derivation

Surface pressure reduction to sea-level

Description

Estimate the sea-level pressure from surface pressure, using temperature and surface geopotential height

Usage

ps2psl(ps, tas, zgs)

Arguments

ps	Surface pressure
tas	Near-surface air temperature
zgs	surface geopotential height. This is a static variable (i.e., constant in time). Thus, member and time dimensions must be either singletons or just dropped.

Value

A climate4R CDM grid of estimated sea-level pressure (in Pascals)

Input units

In principle, the derivation functions will work regardless of the input units, as all units are internally converted as necessary according to the specific formulae definitions. However, the units string (see getGridUnits) must be parseable. In the same vein, the input units need to be convertible to the required ones. Unit consistency is internally achieved by the function udConvertGrid.

Author(s)

J. Bedia, S. Herrera

References

The adjustment formulae of pressure values to other levels are those routinely used by the WMO. These are described in the following technical report: https://www.wmo.int/pages/prog/www/IMOP/meetings/SI/ET-Stand-1/Doc-10_Pressure-red.pdf

See Also

Other derivation: hurs2huss(), hurs2tdps(), hurs2w(), huss2hurs(), huss2pvp(), psl2ps(), rad2cc(), tas2ws(), tdps2hurs()

Other pressure: huss2pvp(), psl2ps(), tas2ws()

Examples

data("ps.iberia")
data("tas.iberia")
data("zgs.iberia")
psl.hat <- ps2psl(ps = ps.iberia, tas = tas.iberia, zgs = zgs.iberia)
require(transformeR)
# Error map:
psl.err <- gridArithmetics(psl.hat, psl.iberia, operator = "-")
## Not run: 
require(visualizeR)
err <- climatology(psl.err)
spatialPlot(err, set.min = -10000,
            at = seq(-10000, 10000, 500),
            color.theme = "Spectral", rev.colors = TRUE,
            backdrop.theme = "coastline",
            main = "Mean estimate error (Pascals)")

## End(Not run)

SantanderMetGroup/convertR documentation built on June 30, 2023, 3:03 a.m.