Use real science/formulas for mosquito prevalence/oviposition layers

[jetpack]

We have the real formula/model for mosquito-related risk indicators, so we should use those.

[jetpack]

From Manuel:

#intrinsic oviposition rate adult mosquito
t_v <- function(x){
  return(-5.4+1.8*x-0.2124*x^2+0.01015*x^3-1.515*10^(-4)*x^4)
}

R_const =  8.314

# Celsius to Kelvin
kelvin <- function(x){
  return(x +273.15)
}

#extrinsic incubation period of adult mosquito
s_v <- function(x){
  t = kelvin(x)
  return(((3.3589*10^(-3)*t/298)*exp((1500/R_const)*(1/298-1/t)))/
           (1+exp((6.203*10^21)/R_const*(1/(-2.176*10^30))-1/t)))
}

# dengue suitability
Z_x <- function(x){
  return(exp(-mu_v(x)*s_v(x))/(mu_v(x)^2))
}