#jamshidi

FIX: GSHex HXex KLEak LACEXT NALEak PYREXT adeex adoex inoex

#reactions
v_1:
 {1.0}ADO = {1.0}INO
  scaleada * adaVM * ADO / (adaKA + ADO)
v_10:
 {1.0}ATP = {1.0}ADP
  scaleatpase*katpase*(ATP)
v_11:
 {1.0}DPG23 = {1.0}PG3
  scaledpgase*dpgaseVM*DPG23/(dpgasek+DPG23)
v_12:
 {1.0}DPG13 = {1.0}DPG23
  scaledpgm*KDPGM*DPG13/(1.0+DPG23/dpgmk)
v_13:
 {1.0}PG2 = {1.0}PEP
  scaleen*(KEN1*PG2-KEN2*PEP)*exp((19.876-6152.2/TEMP) )
v_14:
 {1.0}G6P = {1.0}GL6P + {1.0}NADPH
  scaleg6pd*NUMG6PDH*3600/DENG6PDH/(1*10**-3)
v_15:
 {1.0}GAP = {1.0}DPG13 + {1.0}NADH
  scalegapdh*(KGA1*PI*GAP*(NADHT-NADH)-KGA2*DPG13*NADH*HI)*exp((16.203-5021.3/TEMP))
v_16:
 {1.0}GO6P = {1.0}NADPH + {1.0}RU5P
  scalegl6pd*NUMG16PD/DENG16PD*(160.84-39.173*HI+ 4.1322*HI*HI)
v_17:
 {1.0}GSH = {1.0}GSHex
  scalegshr*gshconst*GSH*(-2880+902.33*PHI-65.272*PHI*PHI)
v_18:
 {1.0}NADPH = {2.0}GSH
  scalegssgr*gssgNUM/gssgDEN*(-2880+902.33*PHI-65.272*PHI*PHI)/233.075
v_19:
 {1.0}HX + {1.0}PRPP = {1.0}IMP
  scalehgprt*hgprtVM*(PRPP/(PRPP+hgprtKPRPP))*(HX/(HX+hgprtKHX))*exp((14.543-4501.5/TEMP))
v_2:
 {1.0}ADE = {1.0}adeex
  scaleadeex*vmade*(ADE/(kmade+ADE)-adeex/(kmade+adeex))
v_20:
 {1.0}ATP = {1.0}ADP + {1.0}G6P
  scalehk*VM1HK*GLC/KMAGLC*MGATP/KMGATP*(1+2.1187*MG/KMGAMG)/N1HK
v_21:
 {1.0}HX = {1.0}HXex
  scalehxex*(hxexPM*HX + hxexVM*HX/(HX+hxexKM))*exp((14.543-4501.5/TEMP))
v_22:
 {1.0}IMP = {1.0}INO
  scaleimpase*KIASE*IMP/(1.0+INO/KINO)*exp(14.543-4501.5/TEMP)
v_23:
 {1.0}INO = {1.0}inoex
  vmino*(INO/(kmino+INO)-inoex/(kmino+inoex))
v_24:
 {1.0}KLEak = {1.0}KI
  scalekleak*(kV1 + kV2)*(-1.16+0.3*PHI)* exp((22.7283-7045.7976/TEMP))
v_25:
 {1.0}NADH + {1.0}PYR = {1.0}LAC
  scaleldh*(KLD1*PYR*NADH-KLD2*LAC*(NADHT-NADH))*exp((20.237-6267.2/TEMP))
v_26:
 {1.0}LAC = {1.0}LACEXT
  scalelex*VLAC*exp((14.543-4501.5/TEMP))
v_27:
 {1.0}NALEak = {1.0}NAI
  scalenaleak*(naV1 + naV2)*(-1.16+0.3*PHI)*exp((29.2222-9058.88/TEMP))
v_28:
 {1.0}ADP + {1.0}MG = {1.0}MGADP
  kAss*((ADP*MG)/(mgadpk+ MG) - MGADP)
v_29:
 {1.0}AMP + {1.0}MG = {1.0}MGAMP
  kAss*((AMP*MG)/(mgampk+ MG) - MGAMP)
v_3:
 {2.0}ADP = {1.0}AMP + {1.0}ATP
  scaleadk*(adkKAP1*ADP**2-adkKAP2*ATP*AMP)*exp((19.876-6152.2/TEMP))
v_30:
 {1.0}ATP + {1.0}MG = {1.0}MGATP
  kAss*((ATP*MG)/(mgatpk+ MG) - MGATP)
v_31:
 {1.0}DPG23 + {1.0}MG = {1.0}MGDPG23
  kAss*((DPG23*MG)/(mgdpgk+ MG) - MGDPG23)
v_32:
 {1.0}PYR = {1.0}PYREXT
  -scalepex*KPYR*(PYREXT-PYR)*exp((14.52096-4501.5/TEMP))
v_33:
 {1.0}ATP + {1.0}F6P = {1.0}ADP + {1.0}FDP
  (scalepfk*pfkVM*F6P/(pfkKF6P+F6P)*MGATP/(pfkKMGATP+MGATP))/(1.0+LPFK)
v_34:
 {1.0}G6P = {1.0}F6P
  scalepgi*(kpgi1*G6P-kpgi2*F6P)*exp((13.042-4040.7/TEMP))
v_35:
 {1.0}ADP + {1.0}DPG13 = {1.0}ATP + {1.0}PG3
  scalepgk*(kpgk1*DPG13*ADP-kpgk2*PG3*ATP)*exp((12.560-3895.3/ TEMP))
v_36:
 {1.0}GL6P = {1.0}GO6P
  scalepglase*pglaseVM*GL6P/(pglaseKM + GL6P)*exp((14.167-4380.4/TEMP))
v_37:
 {1.0}PG3 = {1.0}PG2
  scalepgm*(KPGM1*PG3-KPGM2*PG2)*exp((16.907-5236.7/TEMP))
v_38:
 {1.0}ADP + {1.0}PEP = {1.0}PYR + {1.0}ATP
  scalepk * pkVM(PHI, TEMP) * pkDEN1 / (1 + pkL0 * (pkDEN2 / pkDEN3))
v_39:
 {1.0}INO = {1.0}HX + {1.0}R1P
  scalepnpase*(pnpaseKA1*INO*PI - pnpaseKA2*HX*R1P)*exp(14.543-4501.5/TEMP)
v_4:
 {1.0}ADO = {1.0}adoex
  scaleadoex*vmaex*(ADO/(kmaex+ADO)-adoex/(kmaex+adoex))
v_40:
 {1.0}R1P = {1.0}R5P
  scaleprm*(prmKA1*R1P - prmKA2*R5P )*exp((14.543-4501.5/TEMP))
v_41:
 {1.0}R5P + {1.0}ATP = {1.0}AMP + {1.0}PRPP
  scaleprppsyn*VMPRPPSYN*VRPRPPSYN*(R5P*ATP-PRPP*AMP/KEQPRPPSYN)/(VRPRPPSYN*KR5PPRPPSYN*ATP+VRPRPPSYN*KATPPRPPSYN*R5P+ VMPRPPSYN*KAMPPRPPSYN*PRPP/KEQPRPPSYN+VMPRPPSYN*KPRPPPRPPSYN*AMP/KEQPRPPSYN+VRPRPPSYN*R5P*ATP+VMPRPPSYN*PRPP*AMP/KEQPRPPSYN)*exp((14.543-4501.5/TEMP))
v_42:
 {3.0}NAI + {1.0}ATP = {1.0}ADP + {2.0}KI
  scalepump* FATP*NUMPUMP*.5/DENPUMP*exp((40.7298-12626.2626/TEMP))
v_43:
 {1.0}RU5P = {1.0}R5P
  scaleru5pi*ru5iVF*(RU5P - R5P/ru5iKEQ)/(ru5iK1+RU5P+ru5iK1/ru5iK2*R5P)
v_44:
 {1.0}GAP + {1.0}S7P = {1.0}E4P + {1.0}F6P
  scaletald*(NALD1*S7P*GAP-NALD2*E4P*F6P)*ETALD/DENTALD
v_45:
 {1.0}X5P + {1.0}R5P = {1.0}GAP + {1.0}S7P
  scaletki*(tkiN1*X5P*R5P-tkiN2*GAP*S7P)*tkiET/tkiD
v_46:
 {1.0}X5P + {1.0}E4P = {1.0}F6P + {1.0}GAP
  scaletkii*(NTKII1*X5P*E4P-NTKII2*GAP*F6P)*ETTKII/DENTKII
v_47:
 {1.0}DHAP = {1.0}GAP
  scaletpi*(tpiKTP1*DHAP-tpiKTP2*GAP)*exp((18.475-5726.8/TEMP))
v_48:
 {1.0}RU5P = {1.0}X5P
  scalexu5pe*xu5eVF*(RU5P - X5P/xu5eKEQ)/(xu5eK1+RU5P+xu5eK1/xu5eK2*X5P)
v_5:
 {1.0}PRPP + {1.0}ADE = {1.0}AMP
  scaleadprt*adprtVM*(ADE/(adprtKADE+ADE))*(PRPP/(PRPP+adprtKPRPP))
v_6:
 {1.0}ADO + {1.0}ATP = {1.0}ADP + {1.0}AMP
  scaleak*akVM*ATP/(ATP+akKMGATP)*ADO/(ADO+akKA)
v_7:
 {1.0}FDP = {1.0}DHAP + {1.0}GAP
  scaleald*(aldN1*FDP-aldN2*GAP*DHAP)*aldET/(aldD1+aldD2*FDP+aldD3*GAP+aldD4*DHAP+aldD5*FDP*GAP+aldD6*GAP*DHAP)
v_8:
 {1.0}AMP = {1.0}ADO
  scaleampase*(KAMPAS*AMP)*exp((19.876-6152.2/TEMP))
v_9:
 {1.0}AMP = {1.0}IMP
  scaleampda*ampdaVM*AMP/(ampdaKAMP+AMP)*exp((19.876-6152.2/TEMP))


#parameters
CO2 = 1.2
CONSTHK = 1.55
ETALD = 0.00069
ETG16PD = 2.1e-06
ETG6PD = 9.3e-08
ETTKII = 0.00033
FMAXN = 2.8166
G16P = 0.16
GLC = 5.0
GSSGT = 3.32
KALD1 = 2088000.0
KALD2 = 163080.0
KALD3 = 58680.0
KALD4 = 3636000.0
KALD5 = 1764000.0
KALD6 = 216000.0
KALD7 = 61200.0
KALD8 = 284400.0
KAMPPRPPSYN = 0.275
KATPPRPPSYN = 0.17
KE = 5.0
KEN = 1.708543057
KEN1 = 50000.0
KEQALD = 1.05
KEQPRPPSYN = 28.6
KEQTKII = 10.3
KG16P = 0.334
KG16PD1 = 1200000.0
KG16PD10 = 450000.0
KG16PD11 = 300.0
KG16PD12 = 9900000.0
KG16PD2 = 410.0
KG16PD3 = 1000000000.0
KG16PD4 = 26000.0
KG16PD5 = 48.0
KG16PD6 = 30.0
KG16PD7 = 630.0
KG16PD8 = 36000.0
KG16PD9 = 800.0
KG6P = 0.334
KG6PD1 = 110000000.0
KG6PD10 = 1400000000.0
KG6PD2 = 870.0
KG6PD3 = 26000000.0
KG6PD4 = 300.0
KG6PD5 = 750.0
KG6PD6 = 2000.0
KG6PD7 = 220000.0
KG6PD8 = 1100000000.0
KG6PD9 = 10000.0
KGA1 = 115975.0
KGAPDH = 1.1294e-06
KGLC = 0.04
KINO = 25.0
KLAC = 12.8
KLD1 = 308642.2
KLDH = 44.64
KMAGLC = 0.051
KMG = 1.03
KMGAMG = 1.14
KMGATP = 1.13
KMGPG = 3.44
KNAPUMP = 6.2672
KPG = 1800.0
KPG23 = 2.7
KPGI = 0.41
KPGM = 0.1470588235
KPGM1 = 5000.0
KPRPPPRPPSYN = 0.09
KPUMP1 = 0.0082
KPUMP2 = 0.0501
KPYR = 63.8
KR5PPRPPSYN = 0.65
KSPPUMP = 0.1328
KSPUMP = 0.0617
KTKII1 = 777600.0
KTKII2 = 136800.0
KTKII3 = 122400.0
KTKII4 = 561600.0
KTKII5 = 806400.0
KTKII6 = 630000.0
KTKII7 = 144000.0
KTKII8 = 76680.0
Kglcg6p = 0.069
MGT = 2.7
NADHT = 0.089
NADPHT = 0.066
NAE = 140.0
PHE = 7.4
PI = 1.2
PNA = 3.96e-07
PSIPUMP = 0.7114
RT = 0.6
S = 1.55e-06
TEMP = 310.0
V = 8.7e-11
adaKA = 0.052
adaKII = 1.3
adkK1 = 0.081
adkK2 = 0.81
adkK3 = 22.2
adkk0 = 0.22933781
adprtKADE = 0.0023
adprtKPRPP = 0.0195
akKA = 0.0004
akKMGATP = 0.8
aldET = 0.00037
aldK1 = 38500000.0
aldK2 = 838800.0
aldK3 = 6840000.0
aldK4 = 40320000.0
aldK5 = 252000.0
aldK6 = 23076000.0
aldkeq = 0.081
ampdaKAMP = 0.8
ampdaKIMP = 2.5
ampdaVM = 0.01
dpgaseVM = 0.52
gssgET = 1.25e-07
gssgK1 = 85000000.0
gssgK10 = 50000000.0
gssgK11 = 7000.0
gssgK12 = 100000000.0
gssgK2 = 510.0
gssgK3 = 100000000.0
gssgK4 = 7200.0
gssgK5 = 810.0
gssgK6 = 1000.0
gssgK7 = 1000000.0
gssgK8 = 50000000.0
gssgK9 = 1000000.0
hgprtKHX = 0.22
hgprtKPRPP = 0.005
hgprtVM = 0.2011
hxexKM = 0.4
hxexPM = 37.8
hxexVM = 151.6
kAss = 10000.0
kFMAXK = 3.115
kKM = 4.0
kPK = 3.564e-07
kmade = 2.6
kmaex = 0.13
kmino = 0.13
kpgi1 = 60290.0
kpgk1 = 8633310.0
mgadpk = 0.81
mgampk = 22.2
mgatpk = 0.081
mgdpgk = 1.67
naKM = 21.0
pfkKAMP = 0.033
pfkKF6P = 0.1
pfkKMG = 0.44
pfkKMGATP = 0.068
pfkLO = 0.001072
pglaseKM = 0.08
pglaseVM = 9.93
pkKADP = 0.474
pkKATP = 3.39
pkKFDP = 0.005
pkKPEP = 0.225
pkL0 = 19.0
pnpaseKA1 = 1950.0
pnpaseKeq = 0.09
prmKA1 = 427.59
prmKeq = 13.3
ru5iK1 = 0.78
ru5iK2 = 2.2
ru5iKEQ = 2.57
ru5iVF = 1702.296
scaleada = 1.0
scaleadeex = 1.0
scaleadk = 1.0
scaleadoex = 1.0
scaleadprt = 1.0
scaleak = 1.0
scaleald = 1.0
scaleampase = 1.0
scaleampda = 1.0
scaleatpase = 1.0
scaledpgase = 1.0
scaledpgm = 1.0
scaleen = 1.0
scaleg6pd = 1.0
scalegapdh = 1.0
scalegl6pd = 1.0
scalegshr = 1.0
scalegssgr = 1.0
scalehgprt = 1.0
scalehk = 1.0
scalehxex = 1.0
scaleimpase = 1.0
scalekleak = 1.0
scaleldh = 1.0
scalelex = 1.0
scalenaleak = 1.0
scalepex = 1.0
scalepfk = 1.0
scalepgi = 1.0
scalepgk = 1.0
scalepglase = 1.0
scalepgm = 1.0
scalepk = 1.0
scalepnpase = 1.0
scaleprm = 1.0
scaleprppsyn = 1.0
scalepump = 1.0
scaleru5pi = 1.0
scaletald = 1.0
scaletki = 1.0
scaletkii = 1.0
scaletpi = 1.0
scalexu5pe = 1.0
tkiET = 0.00033
tkiK1 = 777600.0
tkiK2 = 136800.0
tkiK3 = 122400.0
tkiK4 = 561600.0
tkiK5 = 1184400.0
tkiK6 = 630000.0
tkiK7 = 144000.0
tkiK8 = 161280.0
tkiKEQ = 1.2
tpiEPSI = 200.0
tpiKTP1 = 31056.0
tpiKTPI = 22.0
vmade = 90.0
vmaex = 61.2
vmino = 61.2
xu5eK1 = 0.19
xu5eK2 = 1.5
xu5eKEQ = 3.0
xu5eVF = 4633.56
GSHex = 0.0
HXex = 0.0
KLEak = 0.0
LACEXT = 1.2
NALEak = 0.0
PYREXT = 0.06
adeex = 0.0013
adoex = 8e-05
inoex = 0.0002
default_compartment = 1.0


#initial values
ADE = 0.00103
ADO = 5.19e-05
ADP = 0.492
AMP = 0.148
ATP = 2.401
DHAP = 0.147
DPG13 = 0.000228
DPG23 = 4.65
E4P = 0.00507
F6P = 0.0272
FDP = 0.01386
G6P = 0.066
GAP = 0.00669
GL6P = 0.001759
GO6P = 0.0374
GSH = 3.312
HX = 0.000226
IMP = 0.0192
INO = 4.27e-05
KI = 139.882
LAC = 1.398
MG = 0.233836
MGADP = 0.110216
MGAMP = 0.00154265
MGATP = 1.78328
MGDPG23 = 0.571129
NADH = 0.030592
NADPH = 0.0659
NAI = 14.135
PEP = 0.02087
PG2 = 0.01226
PG3 = 0.0839
PRPP = 0.0136
PYR = 0.062
R1P = 0.0117
R5P = 0.01266
RU5P = 0.00493
S7P = 0.02398
X5P = 0.01478


#assignment rules
!F NADP = NADPHT-NADPH
!F HE = 10^(-PHE)*1000
!F NAD = NADHT-NADH
!F HI = HE/(92/100*RT)
!F pfkKATP = Exp((-5.8693*10^-2+0.32796*PHI))/1000.0
!F kpgi2 = kpgi1/KPGI
!F aldN1 = aldK1*aldK3*aldK5
!F DTKII5 = KTKII1*KTKII5*(KTKII3+KTKII7)
!F DTKII6 = KTKII4*KTKII8*(KTKII2+KTKII6)
!F aldD1 = aldK5*(aldK3+aldK2)
!F aldD5 = aldK1*aldK4
!F aldD6 = aldK4*aldK6
!F tpiKTP2 = tpiKTP1*tpiKTPI
!F aldD3 = aldK2*aldK4
!F DALD3 = KALD2*KALD4*(KALD6+KALD7)
!F DALD5 = KALD1*KALD5*(KALD3+KALD7)
!F DALD7 = KALD5*KALD8*(KALD2+KALD3)
!F DALD2 = KALD5*KALD7*(KALD2+KALD3)
!F dpgasek = 0.5208*Exp(-8.0287+0.98217*PHI)
!F DENTKII = DTKII1*X5P+DTKII2*E4P+DTKII3*GAP+DTKII4*F6P+DTKII5*X5P*E4P+DTKII6*GAP*F6P+DTKII7*E4P*F6P+DTKII8*X5P*GAP
!F NALD1 = KALD1*KALD3*KALD5*KALD7
!F KEN2 = KEN1/KEN
!F prmKA2 = prmKA1/prmKeq
!F pkDEN3 = ((1+PEP/pkKPEP)^4)*((1+FDP/pkKFDP)^4)
!F KLD2 = KLD1/KLDH
!F akVM = 8.6229/(1+1.636*10^-4/HI)*Exp((19.876-6152.2/TEMP))
!F adkKAP1 = 200*0.272
!F adkKAP2 = adkKAP1/adkKADK
!F NALD2 = NALD1/KEQALD
!F DENG6PD1 = DG6PD1+(DG6PD2*NADP+DG6PD3*G6P+DG6PD4*GL6P+DG6PD5*NADPH)*1.*10^-3
!F KPGM2 = KPGM1/KPGM
!F katpase = 0.4149938/(1+HI/(3.7342*10^-3)+8.2714*10^-7/HI)*Exp((14.543-4501.5/TEMP))
!F MG6PD2 = KG6PD2*KG6PD4*KG6PD6*KG6PD8*KG6PD10*1.1631/(1.0+HI/(3.867*10^-4))*Exp((14.167-4380.4/TEMP))*ETG6PD
!F DG6PD1 = KG6PD2*KG6PD9*(KG6PD4*KG6PD6+KG6PD4*KG6PD7+KG6PD5*KG6PD7)
!F DG6PD5 = KG6PD2*KG6PD10*(KG6PD4*KG6PD6+KG6PD4*KG6PD7+KG6PD5*KG6PD7)
!F DG6PD3 = KG6PD3*KG6PD5*KG6PD7*KG6PD9
!F DG6PD4 = KG6PD2*KG6PD4*KG6PD6*KG6PD8
!F DG6PD7 = KG6PD1*KG6PD4*KG6PD6*KG6PD8
!F DG6PD8 = KG6PD3*KG6PD5*KG6PD7*KG6PD10
!F pkDEN2 = (1+ATP/pkKATP)^4
!F DG6PD11 = KG6PD3*KG6PD8*KG6PD10*(KG6PD5+KG6PD6)
!F DENG6PD3 = DG6PD9*GL6P*NADPH*1.*10^-6+DG6PD10*NADP*G6P*GL6P*1.*10^-9
!F DENG6PDH = DENG6PD1+DENG6PD2+DENG6PD3+DENG6PD4
!F MG16PD2 = KG16PD2*KG16PD4*KG16PD6*KG16PD8*KG16PD10*KG16PD12*Exp((17.896-5546.8/TEMP))*ETG16PD
!F DG6PD10 = KG6PD1*KG6PD3*KG6PD8*(KG6PD5+KG6PD6)
!F DG16PD3 = KG16PD3*KG16PD5*KG16PD7*KG16PD9*KG16PD11
!F gssgDEN3 = gssgD9*GSH*GSH*1*10^-6+(gssgD10+gssgD11)*GSH*NADP*1*10^-6
!F gssgDEN5 = (gssgD14*NADPH*GSH*GSH+gssgD15*GSSG*GSH*NADP)*1*10^-9
!F gssgDEN7 = gssgD18*GSSG*GSH*GSH*NADP*1*10^-12
!F DG16PD4 = KG16PD2*KG16PD4*KG16PD6*KG16PD8*KG16PD11
!F AG16PD1 = KG16PD1*KG16PD3*(KG16PD5*KG16PD9*KG16PD11+KG16PD6*KG16PD9*KG16PD11)
!F DG16PD6 = AG16PD1+BG16PD1+CG16PD1
!F DG16PD18 = KG16PD3*KG16PD8*KG16PD10*KG16PD12*(KG16PD5+KG16PD6)
!F gssgD12 = gssgK1*gssgK3*gssgK8*gssgK11*(gssgK5+gssgK6)
!F gssgD13 = gssgK1*gssgK3*gssgK5*gssgK7*gssgK10
!F gshconst = 1.1429*10^-3*Exp((5.0476-1567.2/TEMP))
!F tkiN1 = tkiK1*tkiK3*tkiK5*tkiK7
!F tkiN2 = tkiN1/tkiKEQ
!F tkiD1 = tkiK1*tkiK3*(tkiK6+tkiK7)
!F tkiD3 = tkiK2*tkiK4*(tkiK6+tkiK7)
!F tkiD2 = tkiK5*tkiK7*(tkiK2+tkiK3)
!F tkiD4 = tkiK6*tkiK8*(tkiK2+tkiK3)
!F tkiD6 = tkiK4*tkiK8*(tkiK2+tkiK6)
!F tkiD8 = tkiK1*tkiK4*(tkiK6+tkiK7)
!F NTKII1 = KTKII1*KTKII3*KTKII5*KTKII7
!F DG16PD17 = KG16PD1*KG16PD3*KG16PD8*KG16PD10*(KG16PD5+KG16PD6)
!F NTKII2 = NTKII1/KEQTKII
!F DTKII1 = KTKII1*KTKII3*(KTKII6+KTKII7)
!F DTKII3 = KTKII2*KTKII4*(KTKII6+KTKII7)
!F DTKII4 = KTKII6*KTKII8*(KTKII2+KTKII3)
!F DENG16PD4 = (DG16PD12*NADP*GO6P*CO2+DG16PD13*NADP*GO6P*RU5P)*1.*10^-9
!F DENG16PD6A = DG16PD16*CO2*RU5P*NADPH*1.*10^-9
!F DENG16PD6B = DG16PD17*NADP*GO6P*CO2*RU5P*1.*10^-12
!F DENG16PD7 = DG16PD18*GO6P*CO2*RU5P*NADPH*1.*10^-12
!F gssgM2 = gssgK2*gssgK4*gssgK6*gssgK8*gssgK10*gssgK12*Exp(15.136-4690.8/TEMP)*gssgET*15455.6
!F gssgD2 = gssgK1*gssgK9*gssgK11*(gssgK4*gssgK6+gssgK4*gssgK7+gssgK5*gssgK7)
!F gssgD3 = gssgK3*gssgK5*gssgK7*gssgK9*gssgK11
!F gssgD4 = gssgK2*gssgK4*gssgK6*gssgK8*gssgK11
!F DG16PD15 = KG16PD3*KG16PD5*KG16PD7*KG16PD10*KG16PD12
!F DTKII7 = KTKII5*KTKII8*(KTKII2+KTKII3)
!F DTKII8 = KTKII1*KTKII4*(KTKII6+KTKII7)
!F naV2 = FMAXN*((NAE/(naKM+NAE))-(NAI*RT/(NAI*RT+naKM)))
!F gssgD7 = gssgK1*gssgK4*gssgK6*gssgK8*gssgK11
!F gssgD8 = gssgK3*gssgK5*gssgK7*gssgK9*gssgK12
!F VMPUMP = 2.318/(1+HI/(3.7342*10^-3)+8.2714*10^-7/HI)
!F NUMPUMP = VMPUMP*(KE^2+KSPPUMP*KE*PSIPUMP*0.5)
!F DENPUMP1 = KSPUMP*KSPPUMP+2.*KSPPUMP*KE+KE^2
!F DENPUMP2 = (KNAPUMP/NAI+1)^3
!F DENPUMP3 = (KPUMP1*KSPUMP*KSPPUMP+KPUMP2*(KE^2+PSIPUMP*KSPPUMP*KE))
!F DENPUMP = DENPUMP1+DENPUMP2*DENPUMP3
!F FATP = ATP/(ATP+0.04)
!F DALD1 = KALD1*KALD3*(KALD6+KALD7)
!F DG16PD10 = KG16PD2*KG16PD4*KG16PD6*KG16PD8*KG16PD12
!F gssgD10 = gssgK2*gssgK4*gssgK6*gssgK8*gssgK12
!F pnpaseKA2 = pnpaseKA1/pnpaseKeq
!F DG16PD12 = KG16PD1*KG16PD3*KG16PD8*KG16PD11*(KG16PD5+KG16PD6)
!F DALD4 = KALD6*KALD8*(KALD2+KALD3)
!F DALD6 = KALD4*KALD8*(KALD2+KALD6)
!F DENG16PD1 = DG16PD1+(DG16PD2*NADP+DG16PD3*GO6P+DG16PD4*CO2+DG16PD5*NADPH)*1.*10^-3
!F DENG16PD3 = (DG16PD9*CO2*RU5P+DG16PD10*CO2*NADPH+DG16PD11*RU5P*NADPH)*1.*10^-6
!F gssgD15 = gssgK3*gssgK5*gssgK7*gssgK10*gssgK12
!F gssgD17 = gssgK1*gssgK3*gssgK8*gssgK10*(gssgK5+gssgK6)
!F gssgD18 = gssgK3*gssgK8*gssgK10*gssgK12*(gssgK5+gssgK6)
!F gssgDEN2 = (gssgD6*NADPH*GSSG+gssgD7*NADPH*GSH+gssgD8*GSSG*NADP)*1*10^-6
!F DALD8 = KALD1*KALD4*(KALD6+KALD7)
!F DENTALD = DALD1*S7P+DALD2*GAP+DALD3*E4P+DALD4*F6P+DALD5*S7P*GAP+DALD6*E4P*F6P+DALD7*GAP*F6P+DALD8*S7P*E4P
!F KIASE = 0.09*1.0431/(1.0+HI/(2.5642*10^-3)+1.2159*10^-6/HI)
!F VMPRPPSYN = 0.5054*1.6548/(1.0+HI/(1.0099*10^-4)+1.8835*10^-6/HI)
!F VRPRPPSYN = VMPRPPSYN/7.5
!F kV1 = kPK*SV*Ln(RT)/(RT-1)*(KE-RT*KI)
!F kV2 = kFMAXK*((KE/(kKM+KE))-(KI*RT/(KI*RT+kKM)))
!F naV1 = PNA*SV*Ln((RT))/(RT-1)*(NAE-RT*NAI)
!F SV = S/V
!F GSSG = 0.5*(GSSGT-GSH)
!F ATPF = ATP-MGATP
!F VM1HK = 7.5095/(1+HI/(2.2718*10^-4)+4.6015*10^-6/HI)*Exp(18.210-5637.9/TEMP)
!F N1HK = 1+GLC/KGLC+(1+GLC/KMAGLC)*MGATP/KMGATP+(1+GLC/KGLC)*MG/KMG+(1+GLC/KMAGLC)*MGATP*MG/KMGATP/KMGAMG+(G6P/KG6P+GLC*G6P/KGLC/Kglcg6p+G16P/KG16P+GLC/KGLC*CONSTHK)*(1+MG/KMG)+(1+GLC/KGLC)*(DPG23/KPG23+MG*DPG23/KMG/KMGPG)
!F pfkVM = 250.0*Exp((13.474-4174.8/TEMP))
!F LPFK = pfkLO*((1.0+ATPF/pfkKATP)^4*(1.0+MG/pfkKMG)^4)/((1.0+F6P/pfkKF6P)^4*(1.0+AMP/pfkKAMP)^4)
!F aldN2 = aldN1/aldkeq
!F aldD2 = aldK1*(aldK3+aldK5)
!F KGA2 = KGA1/KGAPDH
!F kpgk2 = kpgk1/KPG
!F dpgmk = 17778/10000*Exp(-144503/10000+148/100*PHI)
!F adaVM = 36.348/(1+HI/(1.995*10^-4)+3.162*10^-5/HI)*Exp((19.876-6152.2/TEMP))
!F DG6PD9 = KG6PD8*KG6PD10*(KG6PD2*KG6PD4+KG6PD2*KG6PD5+KG6PD2*KG6PD6+KG6PD4*KG6PD6)
!F DG16PD2 = KG16PD1*KG16PD9*KG16PD11*(KG16PD4*KG16PD6+KG16PD4*KG16PD7+KG16PD5*KG16PD7)
!F CG16PD1 = KG16PD1*KG16PD3*(KG16PD5*KG16PD7*KG16PD11)
!F DG16PD8 = KG16PD3*KG16PD5*KG16PD7*KG16PD9*KG16PD12
!F DG16PD11 = KG16PD2*KG16PD10*KG16PD12*(KG16PD4*KG16PD6+KG16PD4*KG16PD7+KG16PD5*KG16PD7)
!F DG16PD13 = KG16PD1*KG16PD3*KG16PD5*KG16PD7*KG16PD10
!F DG16PD14 = KG16PD1*KG16PD4*KG16PD6*KG16PD8*KG16PD10
!F gssgD6 = gssgK1*gssgK3*(gssgK5*gssgK9*gssgK11+gssgK6*gssgK9*gssgK11+gssgK7*gssgK9*gssgK11+gssgK5*gssgK7*gssgK9+gssgK5*gssgK7*gssgK11)
!F PHI = -log10(HI)+3
!F aldD4 = aldK6*(aldK3+aldK2)
!F pkDEN1 = MGADP/(pkKADP+MGADP)*PEP/(pkKPEP+PEP)
!F VLAC = -KLAC*(LACEXT-LAC)
!F KAMPAS = 1.648/(1+HI/(2.5642*10^-3)+1.2159*10^-6/HI)
!F adkKADK = adkk0*(MG+adkK1)*(MG+adkK3)/(MG+adkK2)/(MG+adkK2)
!F adprtVM = 0.202/(1+HI/(3.981*10^-5)+3.162*10^-7/HI)*Exp((23.356-7240.2/TEMP))
!F MG6PD1 = KG6PD1*KG6PD3*KG6PD5*KG6PD7*KG6PD9*1.1631/(1.0+HI/(3.867*10^-4))*Exp((14.167-4380.4/TEMP))*ETG6PD
!F DG6PD2 = KG6PD1*KG6PD9*(KG6PD4*KG6PD6+KG6PD4*KG6PD7+KG6PD5*KG6PD7)
!F DG6PD6 = KG6PD1*KG6PD3*(KG6PD5*KG6PD7+KG6PD5*KG6PD9+KG6PD6*KG6PD9+KG6PD7*KG6PD9)
!F NUMG6PDH = (MG6PD1*NADP*G6P-MG6PD2*GL6P*NADPH)*1.*10^-6
!F DENG6PD2 = (DG6PD6*NADP*G6P+DG6PD7*NADP*GL6P+DG6PD8*G6P*NADPH)*1.*10^-6
!F DENG6PD4 = DG6PD11*G6P*GL6P*NADPH*1.*10^-9
!F MG16PD1 = KG16PD1*KG16PD3*KG16PD5*KG16PD7*KG16PD9*KG16PD11*Exp((17.896-5546.8/TEMP))*ETG16PD
!F DG16PD1 = KG16PD2*KG16PD9*KG16PD11*(KG16PD4*KG16PD6+KG16PD4*KG16PD7+KG16PD5*KG16PD7)
!F DG16PD5 = KG16PD2*KG16PD9*KG16PD12*(KG16PD4*KG16PD6+KG16PD4*KG16PD7+KG16PD5*KG16PD7)
!F BG16PD1 = KG16PD1*KG16PD3*(KG16PD7*KG16PD9*KG16PD11+KG16PD5*KG16PD7*KG16PD9)
!F DG16PD7 = KG16PD1*KG16PD4*KG16PD6*KG16PD8*KG16PD11
!F DG16PD9 = KG16PD2*KG16PD4*KG16PD6*KG16PD8*KG16PD10
!F DG16PD16 = KG16PD8*KG16PD10*KG16PD12*(KG16PD2*KG16PD4+KG16PD2*KG16PD5+KG16PD2*KG16PD6+KG16PD4*KG16PD6)
!F NUMG16PD = MG16PD1*NADP*GO6P*(36.7065*10^-3)-MG16PD2*CO2*RU5P*NADPH*(36.7065*10^-6)
!F DENG16PD2 = (DG16PD6*NADP*GO6P+DG16PD7*NADP*CO2+DG16PD8*GO6P*NADPH)*1.*10^-6
!F DENG16PD5 = (DG16PD14*NADP*CO2*RU5P+DG16PD15*GO6P*RU5P*NADPH)*1.*10^-9
!F DENG16PD6 = DENG16PD6A+DENG16PD6B
!F DENG16PD = DENG16PD1+DENG16PD2+DENG16PD3+DENG16PD4+DENG16PD5+DENG16PD6+DENG16PD7
!F gssgM1 = gssgK1*gssgK3*gssgK5*gssgK7*gssgK9*gssgK11*Exp(15.136-4690.8/TEMP)*gssgET*15455.6
!F gssgD1 = gssgK2*gssgK9*gssgK11*(gssgK4*gssgK6+gssgK4*gssgK7+gssgK5*gssgK7)
!F gssgD5 = gssgK2*gssgK9*gssgK12*(gssgK4*gssgK6+gssgK4*gssgK7+gssgK5*gssgK7)
!F gssgD9 = gssgK2*gssgK4*gssgK6*gssgK8*gssgK10
!F gssgD11 = gssgK2*gssgK10*gssgK12*(gssgK4*gssgK6+gssgK4*gssgK7+gssgK5*gssgK7)
!F gssgD14 = gssgK1*gssgK4*gssgK6*gssgK8*gssgK10
!F gssgD16 = gssgK8*gssgK10*gssgK12*(gssgK2*gssgK4+gssgK2*gssgK5+gssgK2*gssgK6+gssgK4*gssgK6)
!F gssgNUM = gssgM1*NADPH*GSSG-gssgM2*GSH*GSH*NADP
!F gssgDEN1 = gssgD1+(gssgD2*NADPH+gssgD3*GSSG+gssgD4*GSH+gssgD5*NADP)*1*10^-3
!F gssgDEN4 = (gssgD12+gssgD13)*NADPH*GSSG*GSH*1*10^-9
!F gssgDEN6 = gssgD16*GSH*GSH*NADP*1*10^-9+gssgD17*NADPH*GSSG*GSH*GSH*1*10^-12
!F gssgDEN = gssgDEN1+gssgDEN2+gssgDEN3+gssgDEN4+gssgDEN5+gssgDEN6+gssgDEN7
!F tkiD5 = tkiK1*tkiK5*(tkiK3+tkiK7)
!F tkiD7 = tkiK5*tkiK8*(tkiK2+tkiK3)
!F tkiD = tkiD1*X5P+tkiD2*R5P+tkiD3*GAP+tkiD4*S7P+tkiD5*X5P*R5P+tkiD6*GAP*S7P+tkiD7*R5P*S7P+tkiD8*X5P*GAP
!F DTKII2 = KTKII5*KTKII7*(KTKII2+KTKII3)
!F KDPGM = piecewise(0.0,PHI<6.85,2.75*10.0^5.0*exp(16.907-5236.7/TEMP))