#gould1



#reactions
cG_cZTL_assoc:
 {1.0}cG + {1.0}cZTL = {1.0}cZG
  def * function_4_cG_cZTL_assoc_1(D, L, cG, cZG, cZTL, def, p12, p13)
cG_degr:
 {1.0}cG = $pool
  def * function_4_cG_degr_1(cG, def, m19)
cG_m_degr:
 {1.0}cG_m = $pool
  def * function_4_cG_m_degr_1(cG_m, def, m18)
cG_m_trscr:
 $pool = {1.0}cG_m
  def * function_4_cG_m_trscr_1(L, cL, cP, cT, def, g14, g15, n, n12, o, q4)
cG_trsl:
 $pool = {1.0}cG
  def * function_4_cG_trsl_1(cG_m, def, p11)
cL_degr:
 {1.0}cL = $pool
  def * function_4_cL_degr_1(c, cL, def, g3, m3, p3)
cL_m_degr:
 {1.0}cL_m = $pool
  def * function_4_cL_m_degr_1(D, L, cL_m, def, m1, m2)
cL_m_trscr:
 $pool = {1.0}cL_m
  def * function_4_cL_m_trscr_1(L, a, b, cNI, cP, cP7, cP9, cTm, def, g1, g2, n0, n1, q1)
cL_modif:
 $pool = {1.0}cLm
  def * function_4_cL_modif_1(c, cL, def, g3, p3)
cL_trsl:
 $pool = {1.0}cL
  def * function_4_cL_trsl_1(D, L, cL_m, def, p1, p2)
cLm_degr:
 {1.0}cLm = $pool
  def * function_4_cLm_degr_1(cLm, def, m4)
cNI_degr:
 {1.0}cNI = $pool
  def * function_4_cNI_degr_1(D, L, cNI, def, m17, m24)
cNI_m_degr:
 {1.0}cNI_m = $pool
  def * function_4_cNI_m_degr_1(cNI_m, def, m16)
cNI_m_trscr:
 $pool = {1.0}cNI_m
  def * function_4_cNI_m_trscr_1(cLm, cP7, def, g12, g13, l, m, n10, n11)
cNI_trsl:
 $pool = {1.0}cNI
  def * function_4_cNI_trsl_1(cNI_m, def, p10)
cP7_degr:
 {1.0}cP7 = $pool
  def * function_4_cP7_degr_1(D, L, cP7, def, m15, m23)
cP7_m_degr:
 {1.0}cP7_m = $pool
  def * function_4_cP7_m_degr_1(cP7_m, def, m14)
cP7_m_trscr:
 $pool = {1.0}cP7_m
  def * function_4_cP7_m_trscr_1(cL, cLm, cP9, def, g10, g11, j, k, n8, n9)
cP7_trsl:
 $pool = {1.0}cP7
  def * function_4_cP7_trsl_1(cP7_m, def, p9)
cP9_degr:
 {1.0}cP9 = $pool
  def * function_4_cP9_degr_1(D, L, cP9, def, m13, m22)
cP9_m_degr:
 {1.0}cP9_m = $pool
  def * function_4_cP9_m_degr_1(cP9_m, def, m12)
cP9_m_trscr:
 $pool = {1.0}cP9_m
  def * function_4_cP9_m_trscr_1(L, cL, cP, cT, def, g8, g9, h, i, n4, n7, q3)
cP9_trsl:
 $pool = {1.0}cP9
  def * function_4_cP9_trsl_1(cP9_m, def, p8)
cP_degr:
 {1.0}cP = $pool
  def * function_4_cP_degr_1(L, cP, def, m11)
cP_trsl:
 $pool = {1.0}cP
  def * function_4_cP_trsl_1(D, cP, def, p7)
cT_degr:
 {1.0}cT = $pool
  def * function_4_cT_degr_1(D, L, cT, cZG, cZTL, def, m6, m7, m8, p5)
cT_m_degr:
 {1.0}cT_m = $pool
  def * function_4_cT_m_degr_1(cT_m, def, m5)
cT_m_trscr:
 $pool = {1.0}cT_m
  def * function_4_cT_m_trscr_1(cL, cY, d, def, e, g4, g5, n2, n3)
cT_modif:
 $pool = {1.0}cTm
  def * function_4_cT_modif_1(cT, def, f, g6, p15)
cT_trsl:
 $pool = {1.0}cT
  def * function_4_cT_trsl_1(cT_m, def, p4)
cTm_degr:
 {1.0}cTm = $pool
  def * function_4_cTm_degr_1(D, L, cTm, def, m25, m26)
cY_degr:
 {1.0}cY = $pool
  def * function_4_cY_degr_1(cY, def, m10)
cY_m_degr:
 {1.0}cY_m = $pool
  def * function_4_cY_m_degr_1(cY_m, def, m9)
cY_m_trscr:
 $pool = {1.0}cY_m
  def * function_4_cY_m_trscr_1(D, L, cL, cP, cT, def, g, g16, g7, n5, n6, q2, s)
cY_trsl:
 $pool = {1.0}cY
  def * function_4_cY_trsl_1(cY_m, def, p6)
cZG_degr:
 {1.0}cZG = $pool
  def * function_4_cZG_degr_1(cZG, def, m21)
cZTL_degr:
 {1.0}cZTL = $pool
  def * function_4_cZTL_degr_1(cZTL, def, m20)
cZTL_trsl:
 $pool = {1.0}cZTL
  def * function_4_cZTL_trsl_1(def, p14)


#parameters
Am1 = 2.3323
Am15 = 1287.0
Am25 = 0.559
Am6 = 1499400.0
An0 = 4638600.0
An12 = 3.5369
Ap1 = 559910.0
Em1 = 3.6511
Em15 = 20.7914
Em25 = 1.7253
Em6 = 38.9484
En0 = 40.5939
En12 = 1.0203
Ep1 = 35.3172
Rgas = 0.0083145
Temp = 300.15
a = 2.0
b = 3.0
c = 3.0
d = 2.5
dawn = 0.0
dawn1 = 0.0
dawn2 = 9.0
dusk = 12.0
dusk1 = 3.0
dusk2 = 12.0
e = 2.0
f = 3.0
g = 2.0
g1 = 0.1
g10 = 0.7
g11 = 0.7
g12 = 0.5
g13 = 0.6
g14 = 0.17
g15 = 0.4
g16 = 0.2
g2 = 0.28
g3 = 0.4
g4 = 0.91
g5 = 0.3
g6 = 0.3
g7 = 0.18
g8 = 0.14
g9 = 0.3
h = 2.0
i = 3.0
j = 3.0
k = 3.0
l = 2.0
m = 2.0
m10 = 0.3
m11 = 1.0
m12 = 1.0
m13 = 0.32
m14 = 0.28
m16 = 0.5
m17 = 0.3
m18 = 1.0
m19 = 0.2
m2 = 0.24
m20 = 1.2
m21 = 0.2
m22 = 2.0
m23 = 1.0
m24 = 0.405
m26 = 0.14
m3 = 0.2
m4 = 0.2
m5 = 0.3
m7 = 0.5
m8 = 0.1
m9 = 1.0
n = 1.0
n1 = 1.8
n10 = 0.18
n11 = 0.71
n2 = 0.7
n3 = 0.06
n4 = 0.0
n5 = 3.4
n6 = 1.25
n7 = 0.2
n8 = 0.42
n9 = 0.26
o = 2.0
p10 = 0.36
p11 = 0.23
p12 = 30.0
p13 = 0.4
p14 = 0.45
p15 = 0.05
p2 = 0.27
p3 = 0.1
p4 = 0.268
p5 = 1.0
p6 = 0.44
p7 = 0.3
p8 = 0.7
p9 = 0.4
parameter_1 = 0.5
q1 = 0.8
q2 = 0.5
q3 = 2.9
q4 = 0.6
s = 3.0
def = 1.0


#initial values
cG = 0.0238
cG_m = 0.119
cL = 0.416
cL_m = 1.0
cLm = 0.054
cNI = 0.044
cNI_m = 0.0065
cP = 0.825
cP7 = 0.019
cP7_m = 0.075
cP9 = 0.056
cP9_m = 0.35
cT = 0.393
cT_m = 0.25
cTm = 0.24
cY = 0.1
cY_m = 0.093
cZG = 0.0774
cZTL = 0.323


#assignment rules
!F D = 1.0 - L
!F L = 0.5 * (1.0 + tanh((time - 24.0 * floor(time / 24.0)) / 0.5) - (1.0 + tanh((time - 24.0 * floor(time / 24.0) - 12.0) / 0.5)) + 1.0 + tanh((time - 24.0 * floor(time / 24.0) - 24.0) / 0.5))
!F p1 = Ap1 * exp(-Ep1 / (Rgas * Temp))
!F m25 = Am25 * exp(-Em25 / (Rgas * Temp))
!F m15 = Am15 * exp(-Em15 / (Rgas * Temp))
!F m6 = Am6 * exp(-Em6 / (Rgas * Temp))
!F m1 = Am1 * exp(-Em1 / (Rgas * Temp))
!F n0 = An0 * exp(-En0 / (Rgas * Temp))
!F n12 = An12 * exp(-En12 / (Rgas * Temp))