Calculate the Schmidt stability
Schmidt stability, or the resistance to mechanical mixing due to the potential energy inherent in the stratification of the water column.
Schmidt stability was first defined by Schmidt (1928) and later modified by Hutchinson (1957). This stability index was formalized by Idso (1973) to reduce the effects of lake volume on the calculation (resulting in a mixing energy requirement per unit area).
schmidt.stability(wtr, depths, bthA, bthD, sal = 0)
wtr |
a numeric vector of water temperature in degrees C |
depths |
a numeric vector corresponding to the depths (in m) of the wtr measurements |
bthA |
a numeric vector of cross sectional areas (m^2) corresponding to bthD depths |
bthD |
a numeric vector of depths (m) which correspond to areal measures in bthA |
sal |
a numeric vector of salinity in Practical Salinity Scale units |
a numeric vector of Schmidt stability (J/m^2)
Schmidt, W., 1928. Ueber Temperatur and Stabilitaetsverhaltnisse von Seen. Geo- graphiska Annaler 10, 145-177.
Hutchinson, G.E., 1957. A Treatise on Limnology, vol. 1. John Wiley & Sons, Inc., New York.
Idso, S.B., 1973. On the concept of lake stability. Limnology and Oceanography 18, 681-683.
bthA <- c(1000,900,864,820,200,10)
bthD <- c(0,2.3,2.5,4.2,5.8,7)
wtr <- c(28,27,26.4,26,25.4,24,23.3)
depths <- c(0,1,2,3,4,5,6)
cat('Schmidt stability for input is: ')
cat(schmidt.stability(wtr, depths, bthA, bthD))Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.