********** MODEL NAME
Diurnal Dynamics Clusia minor

********** MODEL NOTES
IC values for: treatment
Timepoint 0: 06:00!!!!
********** MODEL STATES
d/dt(SP) = (1/6)*r1a - 2*r2 - r6 + r4 + r5 - r10 - r_Starch + r11 + r_deg
d/dt(Suc) = r2 - r3 - r_exp
d/dt(Glc) = r3 - r5
d/dt(Frc) = r3 - r4
d/dt(Cit) = r6 - r7 + r9 + (1/6)*a*r1b
d/dt(Fum) = r7 - r8
d/dt(Mal) = r8 - r9 + r10 - r11 + (1/6)*(1-a)*r1b
d/dt(Starch) = r_Starch - r_deg

SP(0) = 0.5
Suc(0) = 20.49
Glc(0) = 3.27
Frc(0) = 12.16
Cit(0) = 57.28
Fum(0) = 0.04
Mal(0) = 84.80
Starch(0) = 5
********** MODEL PARAMETERS
k2_1 = 10.5146   
k2_2 = 5.79085   
k2_3 = 5.17537   
k2_4 = 0.314769   
k3_1 = 0.000255126   
k3_2 = 0.615159   
k3_3 = 2.32693e-08   
k3_4 = 7.71697   
k4_1 = 0.570163   
k4_2 = 0.043104   
k4_3 = 0.000391585   
k4_4 = 0.000375699   
k5_1 = 0.802507   
k5_2 = 0.0016181   
k5_3 = 0.000263783   
k5_4 = 0.323295   
k6_1 = 0.000246025   
k6_2 = 1.8728   
k6_3 = 1.11421   
k6_4 = 6.28683e-05   
k7_1 = 0.000718694   
k7_2 = 0.212993   
k7_3 = 0.161284   
k7_4 = 1.54841e-06   
k8_1 = 63.7199   
k8_2 = 12.5008   
k8_3 = 139.061   
k8_4 = 0.0100108   
k9_1 = 0.0420043   
k9_2 = 0.00269781   
k9_3 = 0.135888   
k9_4 = 0.0182448   
kexp_1 = 0.000191298   
kexp_2 = 0.363273   
kexp_3 = 1.17959e-06   
kexp_4 = 1.59018e-06   
kS_1 = 0.00266759   
kS_2 = 18.6898   
kS_3 = 0.229823   
kS_4 = 1.37968   
k10_1 = 14.3378   
k10_2 = 7.24513e-05   
k10_3 = 33.5106   
k10_4 = 20.2192   
k11_1 = 0.141333   
k11_2 = 0.202592   
k11_3 = 0.839475   
k11_4 = 0.141424   
kdeg_1 = 3.21125e-06   
kdeg_2 = 0.000785093   
kdeg_3 = 0.000272087   
kdeg_4 = 18.2978   
Km3 = 0.00222595   
Ki_f = 2.40126   
Ki_g = 14.6618   
a = 0.924921   
********** MODEL VARIABLES
%k1 = interpcsSB([0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24],[1.61,6.40,7.21,7.76,7.43,12.59,16.01,15.37,9.19,11.02,12.23,3.78,-0.10,-0.16,0.08,0.00,0.16,0.39,0.78,0.84,0.95,1.19,1.38,1.35,1.61],time)
%k1a = interpcsSB([0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24],[1.61,6.40,7.21,7.76,7.43,12.59,16.01,15.37,9.19,11.02,12.23,3.78,0,0,0.08,0.00,0.16,0.39,0.78,0.84,0.95,1.19,1.38,1.35,1.61],time)
%k1b = interpcsSB([0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24],[0,0,0,0,0,0,0,0,0,0,0,0,-0.10,-0.16,0,0,0,0,0,0,0,0,0,0,0],time)
k1a = interpcsSB([0,1,6,10,11,12,13,14,15,20,24],[1.61,6.40,16.01,12.23,3.78,0,0,0.08,0.00,0.95,1.61],time)
k1b = interpcsSB([0,1,6,10,11,12,13,14,15,20,24],[0,0,0,0,0,-0.10,-0.16,0,0,0,0],time)

k2 = interpcsSB([0,4,10,16,24],[k2_1,k2_2,k2_3,k2_4,k2_1],time)
k3 = interpcsSB([0,4,10,16,24],[k3_1,k3_2,k3_3,k3_4,k3_1],time)
k4 = interpcsSB([0,4,10,16,24],[k4_1,k4_2,k4_3,k4_4,k4_1],time)
k5 = interpcsSB([0,4,10,16,24],[k5_1,k5_2,k5_3,k5_4,k5_1],time)
k6 = interpcsSB([0,4,10,16,24],[k6_1,k6_2,k6_3,k6_4,k6_1],time)
k7 = interpcsSB([0,4,10,16,24],[k7_1,k7_2,k7_3,k7_4,k7_1],time)
k8 = interpcsSB([0,4,10,16,24],[k8_1,k8_2,k8_3,k8_4,k8_1],time)
k9 = interpcsSB([0,4,10,16,24],[k9_1,k9_2,k9_3,k9_4,k9_1],time)
k10 = interpcsSB([0,4,10,16,24],[k10_1,k10_2,k10_3,k10_4,k10_1],time)
k11 = interpcsSB([0,4,10,16,24],[k11_1,k11_2,k11_3,k11_4,k11_1],time)
k_Starch = interpcsSB([0,4,10,16,24],[kS_1,kS_2,kS_3,kS_4,kS_1],time)
k_deg = interpcsSB([0,4,10,16,24],[kdeg_1,kdeg_2,kdeg_3,kdeg_4,kdeg_1],time)
kexp = interpcsSB([0,4,10,16,24],[kexp_1,kexp_2,kexp_3,kexp_4,kexp_1],time)
********** MODEL REACTIONS
r1a = k1a
r1b = k1b
r2 = k2*SP
r3 = (k3*Suc)/(Km3*(1+Frc/Ki_f)+Suc*(1+Glc/Ki_g))
r4 = k4*Frc
r5 = k5*Glc
r6 = k6*SP
r7 = k7*Cit
r8 = k8*Fum
r9 = k9*Mal
r_exp = kexp*Suc
r_Starch = k_Starch*SP
r10 = k10*SP
r11 = k11*Mal
r_deg = k_deg*Starch
********** MODEL FUNCTIONS


********** MODEL EVENTS


********** MODEL MATLAB FUNCTIONS