stabilityClass {simulariatools}

Stability class.

Description

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.

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.

Value

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

See Also

plotStabilityClass()

Examples

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