psl2ps: Sea-level pressure to surface pressure

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

Sea-level pressure to surface pressure

Description

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

Usage

psl2ps(psl, tas, zgs)

Arguments

psl	Sea-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.

Details

Surface pressure is often needed for other variable derivations, but in many models is often not available, while mean sea-level pressure is usually provided.

Value

A climate4R CDM grid of estimated surface 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(), ps2psl(), rad2cc(), tas2ws(), tdps2hurs()

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

Examples

data("psl.iberia")
data("tas.iberia")
data("zgs.iberia")
ps.hat <- psl2ps(psl = psl.iberia, tas = tas.iberia, zgs = zgs.iberia)
require(transformeR)
# Error map:
data("ps.iberia")
ps.err <- gridArithmetics(ps.hat, ps.iberia, operator = "-")
## Not run: 
require(visualizeR)
err <- climatology(ps.err)
spatialPlot(err, set.max = 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.