begin name
gould1
end name

begin reactions
cGcZTLassoc {1.0}cG + {1.0}cZTL = {1.0}cZG
cGdegr {1.0}cG = {1.0}$EXTERNAL
cGmdegr {1.0}cGm = {1.0}$EXTERNAL
cGmtrscr {1.0}$EXTERNAL = {1.0}cGm
cGtrsl {1.0}$EXTERNAL = {1.0}cG
cLdegr {1.0}cL = {1.0}$EXTERNAL
cLmdegr {1.0}cLm = {1.0}$EXTERNAL
cLmtrscr {1.0}$EXTERNAL = {1.0}cLm
cLmodif {1.0}$EXTERNAL = {1.0}cLm
cLtrsl {1.0}$EXTERNAL = {1.0}cL
cLmdegr {1.0}cLm = {1.0}$EXTERNAL
cNIdegr {1.0}cNI = {1.0}$EXTERNAL
cNImdegr {1.0}cNIm = {1.0}$EXTERNAL
cNImtrscr {1.0}$EXTERNAL = {1.0}cNIm
cNItrsl {1.0}$EXTERNAL = {1.0}cNI
cP7degr {1.0}cP7 = {1.0}$EXTERNAL
cP7mdegr {1.0}cP7m = {1.0}$EXTERNAL
cP7mtrscr {1.0}$EXTERNAL = {1.0}cP7m
cP7trsl {1.0}$EXTERNAL = {1.0}cP7
cP9degr {1.0}cP9 = {1.0}$EXTERNAL
cP9mdegr {1.0}cP9m = {1.0}$EXTERNAL
cP9mtrscr {1.0}$EXTERNAL = {1.0}cP9m
cP9trsl {1.0}$EXTERNAL = {1.0}cP9
cPdegr {1.0}cP = {1.0}$EXTERNAL
cPtrsl {1.0}$EXTERNAL = {1.0}cP
cTdegr {1.0}cT = {1.0}$EXTERNAL
cTmdegr {1.0}cTm = {1.0}$EXTERNAL
cTmtrscr {1.0}$EXTERNAL = {1.0}cTm
cTmodif {1.0}$EXTERNAL = {1.0}cTm
cTtrsl {1.0}$EXTERNAL = {1.0}cT
cTmdegr {1.0}cTm = {1.0}$EXTERNAL
cYdegr {1.0}cY = {1.0}$EXTERNAL
cYmdegr {1.0}cYm = {1.0}$EXTERNAL
cYmtrscr {1.0}$EXTERNAL = {1.0}cYm
cYtrsl {1.0}$EXTERNAL = {1.0}cY
cZGdegr {1.0}cZG = {1.0}$EXTERNAL
cZTLdegr {1.0}cZTL = {1.0}$EXTERNAL
cZTLtrsl {1.0}$EXTERNAL = {1.0}cZTL
end reactions

begin rate equations
cGcZTLassoc = def*function4cGcZTLassoc1[D, L, cG[t], cZG[t], cZTL[t], def, p12, p13]
cGdegr = def*function4cGdegr1[cG[t], def, m19]
cGmdegr = def*function4cGmdegr1[cGm[t], def, m18]
cGmtrscr = def*function4cGmtrscr1[L, cL[t], cP[t], cT[t], def, g14, g15, n, n12, o, q4]
cGtrsl = def*function4cGtrsl1[cGm[t], def, p11]
cLdegr = def*function4cLdegr1[c, cL[t], def, g3, m3, p3]
cLmdegr = def*function4cLmdegr1[D, L, cLm[t], def, m1, m2]
cLmtrscr = def*function4cLmtrscr1[L, a, b, cNI[t], cP[t], cP7[t], cP9[t], cTm[t], def, g1, g2, n0, n1, q1]
cLmodif = def*function4cLmodif1[c, cL[t], def, g3, p3]
cLtrsl = def*function4cLtrsl1[D, L, cLm[t], def, p1, p2]
cLmdegr = def*function4cLmdegr1[cLm[t], def, m4]
cNIdegr = def*function4cNIdegr1[D, L, cNI[t], def, m17, m24]
cNImdegr = def*function4cNImdegr1[cNIm[t], def, m16]
cNImtrscr = def*function4cNImtrscr1[cLm[t], cP7[t], def, g12, g13, l, m, n10, n11]
cNItrsl = def*function4cNItrsl1[cNIm[t], def, p10]
cP7degr = def*function4cP7degr1[D, L, cP7[t], def, m15, m23]
cP7mdegr = def*function4cP7mdegr1[cP7m[t], def, m14]
cP7mtrscr = def*function4cP7mtrscr1[cL[t], cLm[t], cP9[t], def, g10, g11, j, k, n8, n9]
cP7trsl = def*function4cP7trsl1[cP7m[t], def, p9]
cP9degr = def*function4cP9degr1[D, L, cP9[t], def, m13, m22]
cP9mdegr = def*function4cP9mdegr1[cP9m[t], def, m12]
cP9mtrscr = def*function4cP9mtrscr1[L, cL[t], cP[t], cT[t], def, g8, g9, h, i, n4, n7, q3]
cP9trsl = def*function4cP9trsl1[cP9m[t], def, p8]
cPdegr = def*function4cPdegr1[L, cP[t], def, m11]
cPtrsl = def*function4cPtrsl1[D, cP[t], def, p7]
cTdegr = def*function4cTdegr1[D, L, cT[t], cZG[t], cZTL[t], def, m6, m7, m8, p5]
cTmdegr = def*function4cTmdegr1[cTm[t], def, m5]
cTmtrscr = def*function4cTmtrscr1[cL[t], cY[t], d, def, e, g4, g5, n2, n3]
cTmodif = def*function4cTmodif1[cT[t], def, f, g6, p15]
cTtrsl = def*function4cTtrsl1[cTm[t], def, p4]
cTmdegr = def*function4cTmdegr1[D, L, cTm[t], def, m25, m26]
cYdegr = def*function4cYdegr1[cY[t], def, m10]
cYmdegr = def*function4cYmdegr1[cYm[t], def, m9]
cYmtrscr = def*function4cYmtrscr1[D, L, cL[t], cP[t], cT[t], def, g, g16, g7, n5, n6, q2, s]
cYtrsl = def*function4cYtrsl1[cYm[t], def, p6]
cZGdegr = def*function4cZGdegr1[cZG[t], def, m21]
cZTLdegr = def*function4cZTLdegr1[cZTL[t], def, m20]
cZTLtrsl = def*function4cZTLtrsl1[def, p14]
end rate equations

begin 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
parameter1 = 0.5
q1 = 0.8
q2 = 0.5
q3 = 2.9
q4 = 0.6
s = 3.0
def = 1.0
end parameters

begin initial conditions
cG[0] = cGi
cGm[0] = cGmi
cL[0] = cLi
cLm[0] = cLmi
cLm[0] = cLmi
cNI[0] = cNIi
cNIm[0] = cNImi
cP[0] = cPi
cP7[0] = cP7i
cP7m[0] = cP7mi
cP9[0] = cP9i
cP9m[0] = cP9mi
cT[0] = cTi
cTm[0] = cTmi
cTm[0] = cTmi
cY[0] = cYi
cYm[0] = cYmi
cZG[0] = cZGi
cZTL[0] = cZTLi
end initial conditions

begin initial values
cGi = 0.0238
cGmi = 0.119
cLi = 0.416
cLmi = 1.0
cLmi = 0.054
cNIi = 0.044
cNImi = 0.0065
cPi = 0.825
cP7i = 0.019
cP7mi = 0.075
cP9i = 0.056
cP9mi = 0.35
cTi = 0.393
cTmi = 0.25
cTmi = 0.24
cYi = 0.1
cYmi = 0.093
cZGi = 0.0774
cZTLi = 0.323
end initial values

begin assignment rules
D := 1 - L
L := 0.5*(1 + Tanh[2.*(-24 + t - 24*Floor[t/24])] - Tanh[2.*(-12 + t - 24*Floor[t/24])] + Tanh[2.*(t - 24*Floor[t/24])])
p1 := Ap1/E^(Ep1/(Rgas*Temp))
m25 := Am25/E^(Em25/(Rgas*Temp))
m15 := Am15/E^(Em15/(Rgas*Temp))
m6 := Am6/E^(Em6/(Rgas*Temp))
m1 := Am1/E^(Em1/(Rgas*Temp))
n0 := An0/E^(En0/(Rgas*Temp))
n12 := An12/E^(En12/(Rgas*Temp))
end assignment rules

begin function definitions
function4cP9mdegr1[cP9m_,def_,m12_] := (cP9m*m12)/def
function4cLdegr1[c_,cL_,def_,g3_,m3_,p3_] := (cL*m3 + (cL^c*p3)/(cL^c + g3^c))/def
function4cLmtrscr1[L_,a_,b_,cNI_,cP_,cP7_,cP9_,cTm_,def_,g1_,g2_,n0_,n1_,q1_] := (g1^a*(L*n0 + (cTm^b*n1)/(cTm^b + g2^b) + cP*L*q1))/(def*((cNI + cP7 + cP9)^a + g1^a))
function4cLmdegr1[D_,L_,cLm_,def_,m1_,m2_] := (cLm*(L*m1 + D*m2))/def
function4cLtrsl1[D_,L_,cLm_,def_,p1_,p2_] := (cLm*(L*p1 + D*p2))/def
function4cLmodif1[c_,cL_,def_,g3_,p3_] := (cL^c*p3)/(def*(cL^c + g3^c))
function4cLmdegr1[cLm_,def_,m4_] := (cLm*m4)/def
function4cTmtrscr1[cL_,cY_,d_,def_,e_,g4_,g5_,n2_,n3_] := (g5^e*((cY^d*n2)/(cY^d + g4^d) + n3))/(def*(cL^e + g5^e))
function4cTmdegr1[cTm_,def_,m5_] := (cTm*m5)/def
function4cTtrsl1[cTm_,def_,p4_] := (cTm*p4)/def
function4cTdegr1[D_,L_,cT_,cZG_,cZTL_,def_,m6_,m7_,m8_,p5_] := (cT*m8 + cT*(L*m6 + D*m7)*(cZG + cZTL*p5))/def
function4cTmodif1[cT_,def_,f_,g6_,p15_] := (cT^f*p15)/(def*(cT^f + g6^f))
function4cYmtrscr1[D_,L_,cL_,cP_,cT_,def_,g_,g16_,g7_,n5_,n6_,q2_,s_] := ((g16^g*g7^s*(L*n5 + D*n6))/((cL^g + g16^g)*(cT^s + g7^s)) + cP*L*q2)/def
function4cYtrsl1[cYm_,def_,p6_] := (cYm*p6)/def
function4cYmdegr1[cYm_,def_,m9_] := (cYm*m9)/def
function4cTmdegr1[D_,L_,cTm_,def_,m25_,m26_] := (cTm*(L*m25 + D*m26))/def
function4cPdegr1[L_,cP_,def_,m11_] := (cP*L*m11)/def
function4cPtrsl1[D_,cP_,def_,p7_] := ((1 - cP)*D*p7)/def
function4cYdegr1[cY_,def_,m10_] := (cY*m10)/def
function4cP9trsl1[cP9m_,def_,p8_] := (cP9m*p8)/def
function4cP9mtrscr1[L_,cL_,cP_,cT_,def_,g8_,g9_,h_,i_,n4_,n7_,q3_] := ((g8^h*(L*n4 + (cL^i*n7)/(cL^i + g9^i)))/(cT^h + g8^h) + cP*L*q3)/def
function4cP7trsl1[cP7m_,def_,p9_] := (cP7m*p9)/def
function4cP7mdegr1[cP7m_,def_,m14_] := (cP7m*m14)/def
function4cP7mtrscr1[cL_,cLm_,cP9_,def_,g10_,g11_,j_,k_,n8_,n9_] := (((cL + cLm)^j*n8)/((cL + cLm)^j + g10^j) + (cP9^k*n9)/(cP9^k + g11^k))/def
function4cP9degr1[D_,L_,cP9_,def_,m13_,m22_] := (cP9*(L*m13 + D*m22))/def
function4cNItrsl1[cNIm_,def_,p10_] := (cNIm*p10)/def
function4cNImdegr1[cNIm_,def_,m16_] := (cNIm*m16)/def
function4cP7degr1[D_,L_,cP7_,def_,m15_,m23_] := (cP7*(L*m15 + D*m23))/def
function4cNImtrscr1[cLm_,cP7_,def_,g12_,g13_,l_,m_,n10_,n11_] := ((cLm^l*n10)/(cLm^l + g12^l) + (cP7^m*n11)/(cP7^m + g13^m))/def
function4cGcZTLassoc1[D_,L_,cG_,cZG_,cZTL_,def_,p12_,p13_] := (cG*cZTL*L*p12 - cZG*D*p13)/def
function4cGdegr1[cG_,def_,m19_] := (cG*m19)/def
function4cGtrsl1[cGm_,def_,p11_] := (cGm*p11)/def
function4cZGdegr1[cZG_,def_,m21_] := (cZG*m21)/def
function4cZTLdegr1[cZTL_,def_,m20_] := (cZTL*m20)/def
function4cZTLtrsl1[def_,p14_] := p14/def
function4cGmdegr1[cGm_,def_,m18_] := (cGm*m18)/def
function4cNIdegr1[D_,L_,cNI_,def_,m17_,m24_] := (cNI*(L*m17 + D*m24))/def
function4cGmtrscr1[L_,cL_,cP_,cT_,def_,g14_,g15_,n_,n12_,o_,q4_] := ((g14^n*g15^o*L*n12)/((cT^n + g14^n)*(cL^o + g15^o)) + cP*L*q4)/def
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
cY_m http://identifiers.org/chebi/CHEBI:33699
cY_m http://identifiers.org/kegg.compound/C00046
end species annotations

begin units
end units