Nothing
RotCmoistcorrection=function(P, E, S.Thick, clay, pE, fk)
{
MTSMD = -(20 + 1.3 * clay - 0.01 * (clay^2)) * (S.Thick/23) *(1/fk)
Moistdeficit = P - E * pE
ATSMD = 0
for (i in 2:length(Moistdeficit)) {
ATSMD[1] = ifelse(Moistdeficit[1] > 0, 0, Moistdeficit[1])
if (ATSMD[i - 1] + Moistdeficit[i] < 0) {
ATSMD[i] = ATSMD[i - 1] + Moistdeficit[i]
}
else (ATSMD[i] = 0)
if (ATSMD[i] <= MTSMD) {
ATSMD[i] = MTSMD
}
}
b = ifelse(ATSMD > 0.444 * MTSMD, 1, (0.2 + (1-0.2) * ((MTSMD - ATSMD)/(MTSMD - 0.444 * MTSMD))))
return(b)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.