stability.correction {bigleaf} | R Documentation |

## Integrated Stability Correction Functions for Heat and Momentum

### Description

dimensionless stability functions needed to correct deviations from the exponential wind profile under non-neutral conditions.

### Usage

```
stability.correction(zeta, formulation = c("Dyer_1970", "Businger_1971"))
```

### Arguments

`zeta` |
Stability parameter zeta (-) |

`formulation` |
Formulation for the stability function. Either |

### Details

The functions give the integrated form of the universal functions. They
depend on the value of the stability parameter `\zeta`

,
which can be calculated from the function `stability.parameter`

.
The integration of the universal functions is:

`\psi = -x * zeta`

for stable atmospheric conditions (`\zeta`

>= 0), and

`\psi = 2 * log( (1 + y) / 2) `

for unstable atmospheric conditions (`\zeta`

< 0).

The different formulations differ in their value of x and y.

### Value

a data.frame with the following columns:

`psi_h` |
the value of the stability function for heat and water vapor (-) |

`psi_m` |
the value of the stability function for momentum (-) |

### References

Dyer, A.J., 1974: A review of flux-profile relationships. Boundary-Layer Meteorology 7, 363-372.

Dyer, A. J., Hicks, B.B., 1970: Flux-Gradient relationships in the constant flux layer. Quart. J. R. Meteorol. Soc. 96, 715-721.

Businger, J.A., Wyngaard, J. C., Izumi, I., Bradley, E. F., 1971: Flux-Profile relationships in the atmospheric surface layer. J. Atmospheric Sci. 28, 181-189.

Paulson, C.A., 1970: The mathematical representation of wind speed and temperature profiles in the unstable atmospheric surface layer. Journal of Applied Meteorology 9, 857-861.

Foken, T, 2008: Micrometeorology. Springer, Berlin, Germany.

### Examples

```
zeta <- seq(-2,0.5,0.05)
stability.correction(zeta)
stability.correction(zeta,formulation="Businger_1971")
```

*bigleaf*version 0.8.2 Index]