#locke1



#reactions
four12:
 {1.0}cLc = $pool
  compartment * (m2 * cLc / (k2 + cLc))
four15:
 {1.0}cTc = $pool
  compartment * ((1 - LD(time, dayLength)) * m5 + m6) * (cTc / (k5 + cTc))
four18:
 {1.0}cXc = $pool
  compartment * (m10 * cXc / (k8 + cXc))
four21:
 {1.0}cYc = $pool
  compartment * (m13 * cYc / (k11 + cYc))
one11:
 $pool = {1.0}cLm
  compartment * LD(time, dayLength) * q1 * cPn
one12:
 $pool = {1.0}cLc
  compartment * p1 * cLm
one13:
 {1.0}cLc = {1.0}cLn
  compartment * r1 * cLc
one14:
 $pool = {1.0}cTm
  compartment * (n2 * pow(cYn, b) / (pow(g2, b) + pow(cYn, b))) * (pow(g3, c) / (pow(g3, c) + pow(cLn, c)))
one15:
 $pool = {1.0}cTc
  compartment * p2 * cTm
one16:
 {1.0}cTc = {1.0}cTn
  compartment * r3 * cTc
one17:
 $pool = {1.0}cXm
  compartment * (n3 * pow(cTn, d) / (pow(g4, d) + pow(cTn, d)))
one18:
 $pool = {1.0}cXc
  compartment * p3 * cXm
one19:
 {1.0}cXc = {1.0}cXn
  compartment * r5 * cXc
one20:
 $pool = {1.0}cYm
  compartment * (LD(time, dayLength) * q2 * cPn + (LD(time, dayLength) * n4 + n5) * g5**e / (g5**e + cTn**e)) * (g6**f / (g6**f + cLn**f))
one21:
 $pool = {1.0}cYc
  compartment * p4 * cYm
one22:
 {1.0}cYc = {1.0}cYn
  compartment * r7 * cYc
one23:
 $pool = {1.0}cPn
  compartment * (1 - LD(time, dayLength)) * p5
three11:
 {1.0}cLm = $pool
  compartment * (m1 * cLm / (k1 + cLm))
three12:
 {1.0}cLn = {1.0}cLc
  compartment * r2 * cLn
three13:
 {1.0}cLn = $pool
  compartment * (m3 * cLn / (k3 + cLn))
three15:
 {1.0}cTn = {1.0}cTc
  compartment * r4 * cTn
three16:
 {1.0}cTn = $pool
  compartment * ((1 - LD(time, dayLength)) * m7 + m8) * (cTn / (k6 + cTn))
three18:
 {1.0}cXn = {1.0}cXc
  compartment * r6 * cXn
three19:
 {1.0}cXn = $pool
  compartment * (m11 * cXn / (k9 + cXn))
three21:
 {1.0}cYn = {1.0}cYc
  compartment * r8 * cYn
three22:
 {1.0}cYn = $pool
  compartment * (m14 * cYn / (k12 + cYn))
three23:
 {1.0}cPn = $pool
  compartment * q3 * LD(time, dayLength) * cPn
two11:
 $pool = {1.0}cLm
  compartment * (n1 * pow(cXn, a) / (pow(g1, a) + pow(cXn, a)))
two14:
 {1.0}cTm = $pool
  compartment * (m4 * cTm / (k4 + cTm))
two17:
 {1.0}cXm = $pool
  compartment * (m9 * cXm / (k7 + cXm))
two20:
 {1.0}cYm = $pool
  compartment * (m12 * cYm / (k10 + cYm))
two23:
 {1.0}cPn = $pool
  compartment * (m15 * cPn / (k13 + cPn))


#parameters
CP = 24.0
Lmax = 1.0
Lmin = 0.0
a = 3.3064
b = 1.0258
c = 1.0258
d = 1.4422
dayLength = 12.0
e = 3.6064
f = 1.0237
g1 = 0.876738488
g2 = 0.036805783
g3 = 0.26593318
g4 = 0.538811228
g5 = 1.17803247
g6 = 0.064455137
k1 = 1.817
k10 = 1.7303
k11 = 1.8258
k12 = 1.8066
k13 = 1.2
k2 = 1.5644
k3 = 1.2765
k4 = 2.5734
k5 = 2.7454
k6 = 0.4033
k7 = 6.5585
k8 = 0.6632
k9 = 17.1111
m1 = 1.5283
m10 = 0.2179
m11 = 3.3442
m12 = 4.297
m13 = 0.1347
m14 = 0.6114
m15 = 1.2
m2 = 20.44
m3 = 3.6888
m4 = 3.8231
m5 = 0.0013
m6 = 3.1741
m7 = 0.0492
m8 = 4.0424
m9 = 10.1132
n1 = 5.1694
n2 = 3.0087
n3 = 0.2431
n4 = 0.0857
n5 = 0.1649
p1 = 0.8295
p2 = 4.324
p3 = 2.147
p4 = 0.2485
p5 = 0.5
q1 = 2.4514
q2 = 2.40178
q3 = 1.0
r1 = 16.8363
r2 = 0.1687
r3 = 0.3166
r4 = 2.1509
r5 = 1.0352
r6 = 3.3017
r7 = 2.2123
r8 = 0.2002
two11_Fch_0 = 8.0
four12_Fch_0 = 8.0
compartment = 1.0


#initial values
cLc = 0.015
cLm = 0.539
cLn = 0.0855
cPn = 0.835
cTc = 8.7
cTm = 0.456
cTn = 0.4
cXc = 1.34
cXm = 0.08
cXn = 0.4
cYc = 0.004
cYm = 0.037
cYn = 0.013


#assignment rules