Stability class.

Computes stability class given net radiation, total cloud cover and wind speed.

Usage

stabilityClass(rad, tcc, ws, option = "impact")

Arguments

rad

The net radiation in W/m^2

tcc

The total cloud cover in a range from 1 to 8

ws

wind speed in m/s

option

This is to determine which specific categories to use to determine the stability class. It can be impact to comply with ARIA Impact(tm), pasquill or custom.

Value

stabilityClass returns a numeric vector with Pasquill stability classes coded as: A = 1, B = 2, ... , F = 6.

Details

stabilityClass() computes stability class according to IAEA method based on net radiation, total cloud cover tcc and wind speed. Net radiation and wind are used by day; tcc and wind are used by night.

Three different alogorithms are implemented; see source code for details.

See also

plotStabilityClass()

Examples

# Compute stability class with custom algorithm
stMeteo$cst <- stabilityClass(rad = stMeteo$rad,
                              tcc = stMeteo$tcc,
                              ws = stMeteo$ws,
                              option = "custom")