#arnold1

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, Nt, kc, Epar, NADPH)
PGA_prod_Vc:
 {1.0}RuBP + {1.0}CO2 + {2.0}NADPH = {2.0}PGA
  chloroplast * carboxylation(kc, CO2, Kc, O2, Ko, Epar, RuBP)
PGA_prod_Vo:
 {2.0}NADPH + {1.0}RuBP + {1.0}CO2 = {1.5}PGA
  chloroplast * oxygenation(phi, kc, CO2, Kc, O2, Ko, Epar, RuBP)


#parameters
Epar = 1.33846153846
Ipar = 1000.0
Kc = 460.0
Ko = 330.0
Nt = 0.5
Rp = 3.2
f = 0.23
kc = 2.5
CO2 = 230.0
O2 = 210.0
RuBP = 2.0
chloroplast = 1.0


#initial values
NADPH = 0.21
PGA = 2.4


#assignment rules
!F NADP = Nt - NADPH
!F jp = 1.0 / 2.0 * (1.0 - f) * Ipar
!F phi = 0.21 * (O2 / Ko) / (CO2 / Kc)
!F j = (abs(jp + 467.0) - abs(jp - 467.0)) / 2.0 / 65.0