#arnold2

FIX: CO2 O2 RuBP

#reactions
NADPH_prod:
 $pool = {1.0}NADPH
  chloroplast * NADPH_production(J, NADP, Nt)
PGA_cons:
 {1.0}PGA = {1.0}RuBP
  chloroplast * PGA_consumption(PGA, Rp, NADPH, Nt, Vcmax)
PGA_prod_Vc:
 {1.0}RuBP + {1.0}CO2 + {2.0}NADPH = {2.0}PGA
  chloroplast * carboxylation(Vc, Vj)
PGA_prod_Vo:
 {1.0}RuBP + {1.0}O2 + {2.0}NADPH = {1.5}PGA
  chloroplast * oxygenation(phi, Vc, Vj)


#parameters
Gamma = 42.8926627111176
Kc = 406.066146391504
Ko = 276.900084789612
Nt = 0.5
Q = 1000.0
Rp = 3.2
Vcmax = 2.53232284114507
alpha = 0.3
teta = 0.9
CO2 = 245.0
O2 = 210.0
RuBP = 2.0
chloroplast = 1.0


#initial values
NADPH = 0.21
PGA = 2.4


#assignment rules
!F phi = 0.21 * (O2 / Ko) / (CO2 / Kc)
!F Vomax = 0.21 * Vcmax
!F NADP = Nt - NADPH
!F Vc = Vcmax * CO2 / ((1.0 + O2 / Ko) * Kc + CO2)
!F Vj = J / 4.0 * CO2 / (CO2 + 2.0 * Gamma)
!F J = alpha * Q / (10.0 - 9.0 * teta) / (65.0 * 0.5)