#zeilinger2



#reactions
R1:
 $pool = {1.0}cLm
  nucleus * ld * q1 * cPn
R10:
 {1.0}cTm = $pool
  nucleus * m4 * cTm / (k4 + cTm)
R11:
 $pool = {1.0}cTc
  cytoplasm * p2 * cTm
R12:
 {1.0}cTc = {1.0}cTn
  cytoplasm * r3 * cTc
R13:
 {1.0}cTn = {1.0}cTc
  nucleus * r4 * cTn
R14:
 {1.0}cTc = $pool
  cytoplasm * ((lmax - ld) * m5 + m6) * (cTc / (k5 + cTc))
R15:
 {1.0}cTn = $pool
  nucleus * ((lmax - ld) * m7 + m8) * (cTn / (k6 + cTn))
R16:
 $pool = {1.0}cXm
  nucleus * (n3 * pow(cTn, d) / (pow(g4, d) + pow(cTn, d)))
R17:
 {1.0}cXm = $pool
  nucleus * (m9 * cXm / (k7 + cXm))
R18:
 $pool = {1.0}cXc
  cytoplasm * p3 * cXm
R19:
 {1.0}cXc = {1.0}cXn
  cytoplasm * r5 * cXc
R2:
 $pool = {1.0}cLm
  nucleus * (n1 * pow(cXn, a) / (pow(g1, a) + pow(cXn, a))) * (pow(g7, h) / (pow(g7, h) + pow(cP7n, h))) * (pow(g8, i) / (pow(g8, i) + pow(cP9n, i)))
R20:
 {1.0}cXn = {1.0}cXc
  nucleus * r6 * cXn
R21:
 {1.0}cXc = $pool
  cytoplasm * m10 * cXc / (k8 + cXc)
R22:
 {1.0}cXn = $pool
  nucleus * m11 * cXn / (k9 + cXn)
R23:
 $pool = {1.0}cYm
  nucleus * (ld * q2 * cPn + (ld * n4 + n5) * pow(g5, e) / (pow(g5, e) + pow(cTn, e))) * (pow(g6, f) / (pow(g6, f) + pow(cLn, f)))
R24:
 {1.0}cYm = $pool
  nucleus * m12 * cYm / (k10 + cYm)
R25:
 $pool = {1.0}cYc
  cytoplasm * p4 * cYm
R26:
 {1.0}cYc = {1.0}cYn
  cytoplasm * r7 * cYc
R27:
 {1.0}cYn = {1.0}cYc
  nucleus * r8 * cYn
R28:
 {1.0}cYc = $pool
  cytoplasm * m13 * cYc / (k11 + cYc)
R29:
 {1.0}cYn = $pool
  nucleus * m14 * cYn / (k12 + cYn)
R3:
 {1.0}cLm = $pool
  nucleus * m1 * cLm / (k1 + cLm)
R30:
 $pool = {1.0}cPn
  nucleus * (lmax - ld) * p5
R31a:
 {1.0}cPn = $pool
  nucleus * m15 * cPn / (k13 + cPn)
R31b:
 {1.0}cPn = $pool
  nucleus * q3 * ld * cPn
R32:
 $pool = {1.0}cP7m
  nucleus * n6 * pow(cLn, j) / (pow(g9, j) + pow(cLn, j))
R33:
 {1.0}cP7m = $pool
  nucleus * m16 * cP7m / (k14 + cP7m)
R34:
 $pool = {1.0}cP7c
  cytoplasm * p6 * cP7m
R35:
 {1.0}cP7c = {1.0}cP7n
  cytoplasm * r9 * cP7c
R36:
 {1.0}cP7n = {1.0}cP7c
  nucleus * r10 * cP7n
R37:
 {1.0}cP7c = $pool
  cytoplasm * m17 * cP7c / (k15 + cP7c)
R38:
 {1.0}cP7n = $pool
  nucleus * m18 * cP7n / (k16 + cP7n)
R39:
 $pool = {1.0}cP9m
  (ld * q4 * cPn + n7 * ld + n8) * pow(cLn, k) / (pow(g10, k) + pow(cLn, k))
R4:
 $pool = {1.0}cLc
  cytoplasm * p1 * cLm
R40:
 {1.0}cP9m = $pool
  nucleus * m19 * cP9m / (k17 + cP9m)
R41:
 $pool = {1.0}cP9c
  cytoplasm * p7 * cP9m
R42:
 {1.0}cP9c = {1.0}cP9n
  cytoplasm * r11 * cP9c
R43:
 {1.0}cP9n = {1.0}cP9c
  nucleus * r12 * cP9n
R44:
 {1.0}cP9c = $pool
  cytoplasm * m20 * cP9c / (k18 + cP9c)
R45:
 {1.0}cP9n = $pool
  nucleus * m21 * cP9n / (k19 + cP9n)
R5:
 {1.0}cLc = {1.0}cLn
  cytoplasm * r1 * cLc
R6:
 {1.0}cLn = {1.0}cLc
  nucleus * r2 * cLn
R7:
 {1.0}cLc = $pool
  cytoplasm * m2 * cLc / (k2 + cLc)
R8:
 {1.0}cLn = $pool
  nucleus * m3 * cLn / (k3 + cLn)
R9:
 $pool = {1.0}cTm
  nucleus * (n2 * pow(cYn, b) / (pow(g2, b) + pow(cYn, b))) * (pow(g3, c) / (pow(g3, c) + pow(cLn, c)))


#parameters
Day_in_hours = 24.0
a = 2.2802
b = 3.1075
c = 1.6808
d = 1.0164
e = 1.4943
f = 1.9491
g1 = 16.3389
g10 = 5.6855
g2 = 16.7487
g3 = 11.5922
g4 = 11.3625
g5 = 0.5061
g6 = 7.8469
g7 = 0.4444
g8 = 11.0459
g9 = 14.5219
h = 2.2116
i = 1.1065
j = 2.5579
k = 3.3953
k1 = 22.3951
k10 = 16.1162
k11 = 48.5862
k12 = 23.2876
k13 = 1.2
k14 = 50.9418
k15 = 49.4094
k16 = 42.4837
k17 = 18.6089
k18 = 16.2407
k19 = 26.5795
k2 = 32.7881
k3 = 29.0823
k4 = 4.0551
k5 = 16.9133
k6 = 43.7049
k7 = 8.6873
k8 = 13.4324
k9 = 14.605
ld = 1.0
lmax = 1.0
m1 = 8.0568
m10 = 9.2511
m11 = 7.9066
m12 = 8.4753
m13 = 6.8544
m14 = 3.2581
m15 = 1.2
m16 = 9.531
m17 = 5.4062
m18 = 8.671
m19 = 6.1155
m2 = 10.4609
m20 = 3.4152
m21 = 0.028
m3 = 12.7853
m4 = 8.5185
m5 = 9.3024
m6 = 10.899
m7 = 0.7527
m8 = 13.7459
m9 = 2.6345
n1 = 2.3023
n2 = 7.5433
n3 = 0.6703
n4 = 1.5293
n5 = 2.6296
n6 = 11.3117
n7 = 0.0833
n8 = 2.0738
p1 = 1.2294
p2 = 1.0494
p3 = 8.583
p4 = 14.6828
p5 = 0.5
p6 = 6.7738
p7 = 10.4532
q1 = 7.9798
q2 = 2.5505
q3 = 1.0
q4 = 7.4548
r1 = 31.5166
r10 = 0.4557
r11 = 34.6266
r12 = 22.838
r2 = 9.1138
r3 = 29.4222
r4 = 33.6178
r5 = 27.818
r6 = 4.2863
r7 = 9.1917
r8 = 25.8963
r9 = 31.0318
cytoplasm = 1.0
nucleus = 1.0


#initial values
cLc = 0.0731
cLm = 0.1114
cLn = 0.2366
cP7c = 0.0266
cP7m = 0.0204
cP7n = 1.5103
cP9c = 0.734
cP9m = 0.002
cP9n = 1.1162
cPn = 0.0
cTc = 5.2235
cTm = 3.6732
cTn = 4.5333
cXc = 2.4188
cXm = 0.652
cXn = 14.7289
cYc = 49.2611
cYm = 0.2992
cYn = 17.4355


#assignment rules