Nothing
sucal <-
function(days,lat,rad_mea,extraT =NULL,tmax,tmin,cc) {
i <- dayOfYear(days)
latt <- radians(lat)
if (is.null(extraT)) extraT <- extrat(lat=latt,i)$ExtraTerrestrialSolarRadiationDaily # [MJ]
Y <- rad_mea/extraT
dtemp <- sqrt(tmax-tmin)
cl <- sqrt(1-cc/8)
m <- lm(I(Y) ~ dtemp + cl + I(1/extraT) - 1)
rval <- c(m$coefficients,summary(m)$r.squared)
names(rval) <- c("Sa","Sb","Sc","Sr2")
rval
}
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.