R/ssd2igl.R

In rsachse/limnotools: Help Tools of TU Dresden's Limnology Workgroup

`ssd2igl` <-
function(S, doy, a=0.25, b=0.5, rho=50.29) {
  xi <- 0.0172 * doy - 1.39
  S0 <- 12.3 + sin(xi) * (4.3 + (rho - 51))/6
  Rex <- 245 * (9.9 + 7.08 * sin(xi) + 0.18 * (rho - 51) * (sin(xi) - 1))
  Rex * (a + b * S/S0)
}