#zeilinger3



#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 * 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
  nucleus * ((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 = 1.2497
b = 4.2126
c = 1.4509
d = 1.3058
e = 2.4146
f = 2.1349
g1 = 9.041
g10 = 5.8418
g2 = 16.6598
g3 = 13.4112
g4 = 20.5277
g5 = 1.5987
g6 = 16.489
g7 = 0.2778
g8 = 0.9187
g9 = 20.3795
h = 1.4176
i = 2.0074
j = 1.7615
k = 3.8877
k1 = 13.0594
k10 = 11.5688
k11 = 26.9638
k12 = 21.8348
k13 = 1.2
k14 = 51.261
k15 = 32.939
k16 = 24.451
k17 = 9.8065
k18 = 25.9739
k19 = 21.6441
k2 = 30.5639
k3 = 33.514
k4 = 1.3722
k5 = 34.2078
k6 = 56.7596
k7 = 14.9114
k8 = 12.528
k9 = 15.0626
ld = 1.0
lmax = 1.0
m1 = 6.8248
m10 = 8.5523
m11 = 23.5996
m12 = 5.9504
m13 = 7.5959
m14 = 8.1796
m15 = 1.2
m16 = 9.3186
m17 = 3.6143
m18 = 6.7455
m19 = 1.9234
m2 = 9.4099
m20 = 3.7484
m21 = 0.0193
m3 = 13.7795
m4 = 12.1232
m5 = 7.2129
m6 = 9.5754
m7 = 1.1032
m8 = 2.2006
m9 = 4.2193
n1 = 3.2016
n2 = 11.6086
n3 = 2.4751
n4 = 1.7832
n5 = 7.4615
n6 = 11.0924
n7 = 0.1031
n8 = 3.5262
p1 = 0.6926
p2 = 0.5403
p3 = 6.9124
p4 = 6.0042
p5 = 0.5
p6 = 9.8416
p7 = 1.5323
q1 = 13.4334
q3 = 1.0
q4 = 6.274
r1 = 25.6818
r10 = 0.5024
r11 = 25.7542
r12 = 27.2451
r2 = 3.9781
r3 = 51.1965
r4 = 8.9147
r5 = 29.4607
r6 = 4.5034
r7 = 35.7842
r8 = 27.9229
r9 = 24.5689
cytoplasm = 1.0
nucleus = 1.0


#initial values
cLc = 0.0659
cLm = 0.1951
cLn = 0.3952
cP7c = 0.0861
cP7m = 0.0772
cP7n = 2.7182
cP9c = 0.0077
cP9m = 0.0011
cP9n = 0.0073
cPn = 0.0
cTc = 0.7852
cTm = 0.1234
cTn = 4.5209
cXc = 2.5615
cXm = 1.6104
cXn = 14.5474
cYc = 7.6795
cYm = 1.1149
cYn = 9.7016


#assignment rules