as defined in Basler et al. 2016 (Agr. For. Meteorlogy) with a sigmoidal temperature response (Kramer 1994)
M1s(par, data)
raster or vector with estimated phenophase timing (in DOY)
if (FALSE) { # \dontrun{
estimate <- M1s(data = data, par = par)
} # }