begin name
arnold2
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 {1.0}$RuBP + {1.0}$O2 + {2.0}NADPH = {1.5}PGA
end reactions

begin rate equations
NADPHprod = chloroplast*NADPHproduction[J, NADP, Nt]
PGAcons = chloroplast*PGAconsumption[PGA[t], Rp, NADPH[t], Nt, Vcmax]
PGAprodVc = chloroplast*carboxylation[Vc, Vj]
PGAprodVo = chloroplast*oxygenation[phi, Vc, Vj]
end rate equations

begin 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
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
phi := (0.21*Kc*O2)/(CO2*Ko)
Vomax := 0.21*Vcmax
NADP := Nt - NADPH[t]
Vc := (CO2*Vcmax)/(CO2 + Kc*(1 + O2/Ko))
Vj := (CO2*J)/(4*(CO2 + 2*Gamma))
J := (0.03076923076923077*alpha*Q)/(10 - 9*teta)
end assignment rules

begin function definitions
oxygenation[phi_,Vc_,Vj_] := (phi*(Vc + Vj - Abs[Vc - Vj]))/2
carboxylation[Vc_,Vj_] := (Vc + Vj - Abs[Vc - Vj])/2
NADPHproduction[j_,S1_,Nt_] := (j*S1)/(2*Nt)
PGAconsumption[S1_,Rp_,R_,Nt_,Vc_] := (R*S1*Vc)/(Nt*Rp)
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units