function yprime = WeatherKappa(t,y,tt,Tt,Tt11,Pt,Pnormal,rho,kappa0) T = interp1(tt,Tt,t); T11 = interp1(tt,Tt11,t); P = interp1(tt,Pt,t); Sl = 0.95; Sm = 0.93; Tl = 17; Tm = 23; Pl = 5; A = 2; var_tl = var(Tt); var_tm = var(Tt); var_pl = var(Pnormal); beta_m = 0.2; beta_b = 0.85; k = 0.09; mu = 0.125; gamma = 0.001; rm = 2.325*0.344./(1+1.231*exp(-0.184*(T-20))); delta = 0.25/298.15*(T11+273.15).*exp(28094/1.987*(1/298.15-1./(T11+273.15)))./(1+exp(35362/1.987*(1/298.6-1./(T11+273.15)))); dl = 1-Sl*exp(-((T-Tl)/var_tl).^2); dm = 1-Sm*exp(-((T-Tm)/var_tm).^2); bm = 0.344./(1+1.231*exp(-0.184*(T-20))); kappa = (1+rho)*kappa0./(1+rho*exp(-((P-Pl)/var_pl).^2)); yprime = [rm.*(y(3)+y(4)+y(5))-delta.*y(1)-dl.*y(1)-kappa.*y(1).*y(1); 0.5.*delta.*y(1)-dm.*y(2); 0.5.*delta.*y(1)-bm.*beta_m.*y(3).*y(7)./(y(6)+y(7)+y(8)+A)-dm.*y(3); bm.*beta_m*y(3).*y(7)./(y(6)+y(7)+y(8)+A)-k.*y(4)-dm.*y(4); k*y(4)-dm.*y(5); -bm.*beta_b*y(5).*y(6)/(y(6)+y(7)+y(8)+A); bm.*beta_b*y(5).*y(6)/(y(6)+y(7)+y(8)+A)-mu*y(7)-gamma*y(7); gamma*y(7); mu*y(7)]; end