begin name
arnold1
end name

begin reactions
NADPHprod {1.0}$EXTERNAL = {1.0}NADPH
PGAcons {1.0}PGA = {1.0}$RuBP
PGAprodVc {1.0}$RuBP + {1.0}$CO2 + {2.0}NADPH = {2.0}PGA
PGAprodVo {2.0}NADPH + {1.0}$RuBP + {1.0}$CO2 = {1.5}PGA
end reactions

begin rate equations
NADPHprod = chloroplast*NADPHproduction[j, NADP, Nt]
PGAcons = chloroplast*PGAconsumption[PGA[t], Rp, Nt, kc, Epar, NADPH[t]]
PGAprodVc = chloroplast*carboxylation[kc, CO2, Kc, O2, Ko, Epar, RuBP]
PGAprodVo = chloroplast*oxygenation[phi, kc, CO2, Kc, O2, Ko, Epar, RuBP]
end rate equations

begin 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
end parameters

begin initial conditions
NADPH[0] = NADPHi
PGA[0] = PGAi
end initial conditions

begin initial values
NADPHi = 0.21
PGAi = 2.4
end initial values

begin assignment rules
NADP := Nt - NADPH[t]
jp := ((1 - f)*Ipar)/2
phi := (0.21*Kc*O2)/(CO2*Ko)
j := (-Abs[-467 + jp] + Abs[467 + jp])/130
end assignment rules

begin function definitions
oxygenation[phi_,kc_,S2_,K_,R_,KI_,Epar_,S1_] := (phi*((Epar*kc*S2)/(K*(1 + R/KI) + S2) + (kc*S1*S2)/(K*(1 + R/KI) + S2) - Abs[(Epar*kc*S2)/(K*(1 + R/KI) + S2) - (kc*S1*S2)/(K*(1 + R/KI) + S2)]))/2
PGAconsumption[S1_,Rp_,Nt_,kc_,Epar_,R_] := (Epar*kc*R*S1)/(Nt*Rp)
carboxylation[kc_,S2_,K_,R_,KI_,Epar_,S1_] := ((Epar*kc*S2)/(K*(1 + R/KI) + S2) + (kc*S1*S2)/(K*(1 + R/KI) + S2) - Abs[(Epar*kc*S2)/(K*(1 + R/KI) + S2) - (kc*S1*S2)/(K*(1 + R/KI) + S2)])/2
NADPHproduction[j_,S1_,Nt_] := (j*S1)/(2*Nt)
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units