simulariatools (version 3.1.0)

stabilityClass: Stability class.

Description

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

Usage

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

Value

stabilityClass returns a numeric vector with Pasquill stability classes coded as: A = 1, B = 2, ... , F = 6 ranging from "very unstable" to "very stable".

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

The method used to determine the stability class. It can be iaea (default), pasquill or custom.

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 algorithms are implemented, selected by the option argument.

iaea option implements the *radiation-wind method recommended by the International Atomic Energy Agency (IAEA) and it is based on the net radiation during the day and cloud cover by night.

pasquill option is based on the original Pasquill formulation and lacks the "very weak" solar insolation present in the modified iaea version.

Eventually, the custom options is similar to iaea, with slightly different set of parameters for net radiation, wind speed and cloud cover.

Previously used option impact is the same as iaea and it is now deprecated.

See Also

turnerStabilityClass() which computes stability class with Turner method. plotStabilityClass() to produce graphical outputs with stability class.

Examples

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