Nothing
library(data.table)
data("MOD13A1")
dt <- tidy_MOD13(MOD13A1$dt)
st <- MOD13A1$st
sitename <- dt$site[1]
d <- dt[site == sitename, ] # get the first site data
sp <- st[site == sitename, ] # station point
# global parameter
IsPlot = TRUE
print = FALSE
nptperyear = 23
ypeak_min = 0.05
south = sp$lat < 0
dnew <- add_HeadTail(d, south, nptperyear) # add one year in head and tail
INPUT <- check_input(dnew$t, dnew$y, dnew$w, dnew$QC_flag,
nptperyear, south = south,
maxgap = nptperyear/4, alpha = 0.02, wmin = 0.2)
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.