(* Generated by JWS Online *)

(* This is an experimental feature of JWS Online. Please report any mistakes.*)

(* Note that the following notable SBML entities or features are not supported in notebook outputyet: *)
    (* Events *)
    (* Constraints *)
    (* Units and UnitDefinitions *)
    (* AlgebraicRules *)
    (* conversionFactors *)

variables = {
cG[t], cG\[LetterSpace]m[t], cL[t], cL\[LetterSpace]m[t], cLm[t], cNI[t], cNI\[LetterSpace]m[t], cP[t], cP7[t], cP7\[LetterSpace]m[t], cP9[t], cP9\[LetterSpace]m[t], cT[t], cT\[LetterSpace]m[t], cTm[t], cY[t], cY\[LetterSpace]m[t], cZG[t], cZTL[t]
};

initialValues = {
cG[0] == 0.0238, cG\[LetterSpace]m[0] == 0.119, cL[0] == 0.416, cL\[LetterSpace]m[0] == 1.0, cLm[0] == 0.054, cNI[0] == 0.044, cNI\[LetterSpace]m[0] == 0.0065, cP[0] == 0.825, cP7[0] == 0.019, cP7\[LetterSpace]m[0] == 0.075, cP9[0] == 0.056, cP9\[LetterSpace]m[0] == 0.35, cT[0] == 0.393, cT\[LetterSpace]m[0] == 0.25, cTm[0] == 0.24, cY[0] == 0.1, cY\[LetterSpace]m[0] == 0.093, cZG[0] == 0.0774, cZTL[0] == 0.323
};

rates = {
cG\[LetterSpace]cZTL\[LetterSpace]assoc, cG\[LetterSpace]degr, cG\[LetterSpace]m\[LetterSpace]degr, cG\[LetterSpace]m\[LetterSpace]trscr, cG\[LetterSpace]trsl, cL\[LetterSpace]degr, cL\[LetterSpace]m\[LetterSpace]degr, cL\[LetterSpace]m\[LetterSpace]trscr, cL\[LetterSpace]modif, cL\[LetterSpace]trsl, cLm\[LetterSpace]degr, cNI\[LetterSpace]degr, cNI\[LetterSpace]m\[LetterSpace]degr, cNI\[LetterSpace]m\[LetterSpace]trscr, cNI\[LetterSpace]trsl, cP7\[LetterSpace]degr, cP7\[LetterSpace]m\[LetterSpace]degr, cP7\[LetterSpace]m\[LetterSpace]trscr, cP7\[LetterSpace]trsl, cP9\[LetterSpace]degr, cP9\[LetterSpace]m\[LetterSpace]degr, cP9\[LetterSpace]m\[LetterSpace]trscr, cP9\[LetterSpace]trsl, cP\[LetterSpace]degr, cP\[LetterSpace]trsl, cT\[LetterSpace]degr, cT\[LetterSpace]m\[LetterSpace]degr, cT\[LetterSpace]m\[LetterSpace]trscr, cT\[LetterSpace]modif, cT\[LetterSpace]trsl, cTm\[LetterSpace]degr, cY\[LetterSpace]degr, cY\[LetterSpace]m\[LetterSpace]degr, cY\[LetterSpace]m\[LetterSpace]trscr, cY\[LetterSpace]trsl, cZG\[LetterSpace]degr, cZTL\[LetterSpace]degr, cZTL\[LetterSpace]trsl
};

rateEquations = {
cG\[LetterSpace]cZTL\[LetterSpace]assoc -> def*function\[LetterSpace]4\[LetterSpace]cG\[LetterSpace]cZTL\[LetterSpace]assoc\[LetterSpace]1[D, L, cG[t], cZG[t], cZTL[t], def, p12, p13], cG\[LetterSpace]degr -> def*function\[LetterSpace]4\[LetterSpace]cG\[LetterSpace]degr\[LetterSpace]1[cG[t], def, m19], cG\[LetterSpace]m\[LetterSpace]degr -> def*function\[LetterSpace]4\[LetterSpace]cG\[LetterSpace]m\[LetterSpace]degr\[LetterSpace]1[cG\[LetterSpace]m[t], def, m18], cG\[LetterSpace]m\[LetterSpace]trscr -> def*function\[LetterSpace]4\[LetterSpace]cG\[LetterSpace]m\[LetterSpace]trscr\[LetterSpace]1[L, cL[t], cP[t], cT[t], def, g14, g15, n, n12, o, q4], cG\[LetterSpace]trsl -> def*function\[LetterSpace]4\[LetterSpace]cG\[LetterSpace]trsl\[LetterSpace]1[cG\[LetterSpace]m[t], def, p11], cL\[LetterSpace]degr -> def*function\[LetterSpace]4\[LetterSpace]cL\[LetterSpace]degr\[LetterSpace]1[c, cL[t], def, g3, m3, p3], cL\[LetterSpace]m\[LetterSpace]degr -> def*function\[LetterSpace]4\[LetterSpace]cL\[LetterSpace]m\[LetterSpace]degr\[LetterSpace]1[D, L, cL\[LetterSpace]m[t], def, m1, m2], cL\[LetterSpace]m\[LetterSpace]trscr -> def*function\[LetterSpace]4\[LetterSpace]cL\[LetterSpace]m\[LetterSpace]trscr\[LetterSpace]1[L, a, b, cNI[t], cP[t], cP7[t], cP9[t], cTm[t], def, g1, g2, n0, n1, q1], cL\[LetterSpace]modif -> def*function\[LetterSpace]4\[LetterSpace]cL\[LetterSpace]modif\[LetterSpace]1[c, cL[t], def, g3, p3], cL\[LetterSpace]trsl -> def*function\[LetterSpace]4\[LetterSpace]cL\[LetterSpace]trsl\[LetterSpace]1[D, L, cL\[LetterSpace]m[t], def, p1, p2], cLm\[LetterSpace]degr -> def*function\[LetterSpace]4\[LetterSpace]cLm\[LetterSpace]degr\[LetterSpace]1[cLm[t], def, m4], cNI\[LetterSpace]degr -> def*function\[LetterSpace]4\[LetterSpace]cNI\[LetterSpace]degr\[LetterSpace]1[D, L, cNI[t], def, m17, m24], cNI\[LetterSpace]m\[LetterSpace]degr -> def*function\[LetterSpace]4\[LetterSpace]cNI\[LetterSpace]m\[LetterSpace]degr\[LetterSpace]1[cNI\[LetterSpace]m[t], def, m16], cNI\[LetterSpace]m\[LetterSpace]trscr -> def*function\[LetterSpace]4\[LetterSpace]cNI\[LetterSpace]m\[LetterSpace]trscr\[LetterSpace]1[cLm[t], cP7[t], def, g12, g13, l, m, n10, n11], cNI\[LetterSpace]trsl -> def*function\[LetterSpace]4\[LetterSpace]cNI\[LetterSpace]trsl\[LetterSpace]1[cNI\[LetterSpace]m[t], def, p10], cP7\[LetterSpace]degr -> def*function\[LetterSpace]4\[LetterSpace]cP7\[LetterSpace]degr\[LetterSpace]1[D, L, cP7[t], def, m15, m23], cP7\[LetterSpace]m\[LetterSpace]degr -> def*function\[LetterSpace]4\[LetterSpace]cP7\[LetterSpace]m\[LetterSpace]degr\[LetterSpace]1[cP7\[LetterSpace]m[t], def, m14], cP7\[LetterSpace]m\[LetterSpace]trscr -> def*function\[LetterSpace]4\[LetterSpace]cP7\[LetterSpace]m\[LetterSpace]trscr\[LetterSpace]1[cL[t], cLm[t], cP9[t], def, g10, g11, j, k, n8, n9], cP7\[LetterSpace]trsl -> def*function\[LetterSpace]4\[LetterSpace]cP7\[LetterSpace]trsl\[LetterSpace]1[cP7\[LetterSpace]m[t], def, p9], cP9\[LetterSpace]degr -> def*function\[LetterSpace]4\[LetterSpace]cP9\[LetterSpace]degr\[LetterSpace]1[D, L, cP9[t], def, m13, m22], cP9\[LetterSpace]m\[LetterSpace]degr -> def*function\[LetterSpace]4\[LetterSpace]cP9\[LetterSpace]m\[LetterSpace]degr\[LetterSpace]1[cP9\[LetterSpace]m[t], def, m12], cP9\[LetterSpace]m\[LetterSpace]trscr -> def*function\[LetterSpace]4\[LetterSpace]cP9\[LetterSpace]m\[LetterSpace]trscr\[LetterSpace]1[L, cL[t], cP[t], cT[t], def, g8, g9, h, i, n4, n7, q3], cP9\[LetterSpace]trsl -> def*function\[LetterSpace]4\[LetterSpace]cP9\[LetterSpace]trsl\[LetterSpace]1[cP9\[LetterSpace]m[t], def, p8], cP\[LetterSpace]degr -> def*function\[LetterSpace]4\[LetterSpace]cP\[LetterSpace]degr\[LetterSpace]1[L, cP[t], def, m11], cP\[LetterSpace]trsl -> def*function\[LetterSpace]4\[LetterSpace]cP\[LetterSpace]trsl\[LetterSpace]1[D, cP[t], def, p7], cT\[LetterSpace]degr -> def*function\[LetterSpace]4\[LetterSpace]cT\[LetterSpace]degr\[LetterSpace]1[D, L, cT[t], cZG[t], cZTL[t], def, m6, m7, m8, p5], cT\[LetterSpace]m\[LetterSpace]degr -> def*function\[LetterSpace]4\[LetterSpace]cT\[LetterSpace]m\[LetterSpace]degr\[LetterSpace]1[cT\[LetterSpace]m[t], def, m5], cT\[LetterSpace]m\[LetterSpace]trscr -> def*function\[LetterSpace]4\[LetterSpace]cT\[LetterSpace]m\[LetterSpace]trscr\[LetterSpace]1[cL[t], cY[t], d, def, e, g4, g5, n2, n3], cT\[LetterSpace]modif -> def*function\[LetterSpace]4\[LetterSpace]cT\[LetterSpace]modif\[LetterSpace]1[cT[t], def, f, g6, p15], cT\[LetterSpace]trsl -> def*function\[LetterSpace]4\[LetterSpace]cT\[LetterSpace]trsl\[LetterSpace]1[cT\[LetterSpace]m[t], def, p4], cTm\[LetterSpace]degr -> def*function\[LetterSpace]4\[LetterSpace]cTm\[LetterSpace]degr\[LetterSpace]1[D, L, cTm[t], def, m25, m26], cY\[LetterSpace]degr -> def*function\[LetterSpace]4\[LetterSpace]cY\[LetterSpace]degr\[LetterSpace]1[cY[t], def, m10], cY\[LetterSpace]m\[LetterSpace]degr -> def*function\[LetterSpace]4\[LetterSpace]cY\[LetterSpace]m\[LetterSpace]degr\[LetterSpace]1[cY\[LetterSpace]m[t], def, m9], cY\[LetterSpace]m\[LetterSpace]trscr -> def*function\[LetterSpace]4\[LetterSpace]cY\[LetterSpace]m\[LetterSpace]trscr\[LetterSpace]1[D, L, cL[t], cP[t], cT[t], def, g, g16, g7, n5, n6, q2, s], cY\[LetterSpace]trsl -> def*function\[LetterSpace]4\[LetterSpace]cY\[LetterSpace]trsl\[LetterSpace]1[cY\[LetterSpace]m[t], def, p6], cZG\[LetterSpace]degr -> def*function\[LetterSpace]4\[LetterSpace]cZG\[LetterSpace]degr\[LetterSpace]1[cZG[t], def, m21], cZTL\[LetterSpace]degr -> def*function\[LetterSpace]4\[LetterSpace]cZTL\[LetterSpace]degr\[LetterSpace]1[cZTL[t], def, m20], cZTL\[LetterSpace]trsl -> def*function\[LetterSpace]4\[LetterSpace]cZTL\[LetterSpace]trsl\[LetterSpace]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\[LetterSpace]1 -> 0.5, q1 -> 0.8, q2 -> 0.5, q3 -> 2.9, q4 -> 0.6, s -> 3.0, def -> 1.0
};

assignments = {
function\[LetterSpace]4\[LetterSpace]cP9\[LetterSpace]m\[LetterSpace]degr\[LetterSpace]1[cP9\[LetterSpace]m_,def_,m12_] -> (cP9\[LetterSpace]m*m12)/def, function\[LetterSpace]4\[LetterSpace]cL\[LetterSpace]degr\[LetterSpace]1[c_,cL_,def_,g3_,m3_,p3_] -> (cL*m3 + (cL^c*p3)/(cL^c + g3^c))/def, function\[LetterSpace]4\[LetterSpace]cL\[LetterSpace]m\[LetterSpace]trscr\[LetterSpace]1[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)), function\[LetterSpace]4\[LetterSpace]cL\[LetterSpace]m\[LetterSpace]degr\[LetterSpace]1[D_,L_,cL\[LetterSpace]m_,def_,m1_,m2_] -> (cL\[LetterSpace]m*(L*m1 + D*m2))/def, function\[LetterSpace]4\[LetterSpace]cL\[LetterSpace]trsl\[LetterSpace]1[D_,L_,cL\[LetterSpace]m_,def_,p1_,p2_] -> (cL\[LetterSpace]m*(L*p1 + D*p2))/def, function\[LetterSpace]4\[LetterSpace]cL\[LetterSpace]modif\[LetterSpace]1[c_,cL_,def_,g3_,p3_] -> (cL^c*p3)/(def*(cL^c + g3^c)), function\[LetterSpace]4\[LetterSpace]cLm\[LetterSpace]degr\[LetterSpace]1[cLm_,def_,m4_] -> (cLm*m4)/def, function\[LetterSpace]4\[LetterSpace]cT\[LetterSpace]m\[LetterSpace]trscr\[LetterSpace]1[cL_,cY_,d_,def_,e_,g4_,g5_,n2_,n3_] -> (g5^e*((cY^d*n2)/(cY^d + g4^d) + n3))/(def*(cL^e + g5^e)), function\[LetterSpace]4\[LetterSpace]cT\[LetterSpace]m\[LetterSpace]degr\[LetterSpace]1[cT\[LetterSpace]m_,def_,m5_] -> (cT\[LetterSpace]m*m5)/def, function\[LetterSpace]4\[LetterSpace]cT\[LetterSpace]trsl\[LetterSpace]1[cT\[LetterSpace]m_,def_,p4_] -> (cT\[LetterSpace]m*p4)/def, function\[LetterSpace]4\[LetterSpace]cT\[LetterSpace]degr\[LetterSpace]1[D_,L_,cT_,cZG_,cZTL_,def_,m6_,m7_,m8_,p5_] -> (cT*m8 + cT*(L*m6 + D*m7)*(cZG + cZTL*p5))/def, function\[LetterSpace]4\[LetterSpace]cT\[LetterSpace]modif\[LetterSpace]1[cT_,def_,f_,g6_,p15_] -> (cT^f*p15)/(def*(cT^f + g6^f)), function\[LetterSpace]4\[LetterSpace]cY\[LetterSpace]m\[LetterSpace]trscr\[LetterSpace]1[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, function\[LetterSpace]4\[LetterSpace]cY\[LetterSpace]trsl\[LetterSpace]1[cY\[LetterSpace]m_,def_,p6_] -> (cY\[LetterSpace]m*p6)/def, function\[LetterSpace]4\[LetterSpace]cY\[LetterSpace]m\[LetterSpace]degr\[LetterSpace]1[cY\[LetterSpace]m_,def_,m9_] -> (cY\[LetterSpace]m*m9)/def, function\[LetterSpace]4\[LetterSpace]cTm\[LetterSpace]degr\[LetterSpace]1[D_,L_,cTm_,def_,m25_,m26_] -> (cTm*(L*m25 + D*m26))/def, function\[LetterSpace]4\[LetterSpace]cP\[LetterSpace]degr\[LetterSpace]1[L_,cP_,def_,m11_] -> (cP*L*m11)/def, function\[LetterSpace]4\[LetterSpace]cP\[LetterSpace]trsl\[LetterSpace]1[D_,cP_,def_,p7_] -> ((1 - cP)*D*p7)/def, function\[LetterSpace]4\[LetterSpace]cY\[LetterSpace]degr\[LetterSpace]1[cY_,def_,m10_] -> (cY*m10)/def, function\[LetterSpace]4\[LetterSpace]cP9\[LetterSpace]trsl\[LetterSpace]1[cP9\[LetterSpace]m_,def_,p8_] -> (cP9\[LetterSpace]m*p8)/def, function\[LetterSpace]4\[LetterSpace]cP9\[LetterSpace]m\[LetterSpace]trscr\[LetterSpace]1[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, function\[LetterSpace]4\[LetterSpace]cP7\[LetterSpace]trsl\[LetterSpace]1[cP7\[LetterSpace]m_,def_,p9_] -> (cP7\[LetterSpace]m*p9)/def, function\[LetterSpace]4\[LetterSpace]cP7\[LetterSpace]m\[LetterSpace]degr\[LetterSpace]1[cP7\[LetterSpace]m_,def_,m14_] -> (cP7\[LetterSpace]m*m14)/def, function\[LetterSpace]4\[LetterSpace]cP7\[LetterSpace]m\[LetterSpace]trscr\[LetterSpace]1[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, function\[LetterSpace]4\[LetterSpace]cP9\[LetterSpace]degr\[LetterSpace]1[D_,L_,cP9_,def_,m13_,m22_] -> (cP9*(L*m13 + D*m22))/def, function\[LetterSpace]4\[LetterSpace]cNI\[LetterSpace]trsl\[LetterSpace]1[cNI\[LetterSpace]m_,def_,p10_] -> (cNI\[LetterSpace]m*p10)/def, function\[LetterSpace]4\[LetterSpace]cNI\[LetterSpace]m\[LetterSpace]degr\[LetterSpace]1[cNI\[LetterSpace]m_,def_,m16_] -> (cNI\[LetterSpace]m*m16)/def, function\[LetterSpace]4\[LetterSpace]cP7\[LetterSpace]degr\[LetterSpace]1[D_,L_,cP7_,def_,m15_,m23_] -> (cP7*(L*m15 + D*m23))/def, function\[LetterSpace]4\[LetterSpace]cNI\[LetterSpace]m\[LetterSpace]trscr\[LetterSpace]1[cLm_,cP7_,def_,g12_,g13_,l_,m_,n10_,n11_] -> ((cLm^l*n10)/(cLm^l + g12^l) + (cP7^m*n11)/(cP7^m + g13^m))/def, function\[LetterSpace]4\[LetterSpace]cG\[LetterSpace]cZTL\[LetterSpace]assoc\[LetterSpace]1[D_,L_,cG_,cZG_,cZTL_,def_,p12_,p13_] -> (cG*cZTL*L*p12 - cZG*D*p13)/def, function\[LetterSpace]4\[LetterSpace]cG\[LetterSpace]degr\[LetterSpace]1[cG_,def_,m19_] -> (cG*m19)/def, function\[LetterSpace]4\[LetterSpace]cG\[LetterSpace]trsl\[LetterSpace]1[cG\[LetterSpace]m_,def_,p11_] -> (cG\[LetterSpace]m*p11)/def, function\[LetterSpace]4\[LetterSpace]cZG\[LetterSpace]degr\[LetterSpace]1[cZG_,def_,m21_] -> (cZG*m21)/def, function\[LetterSpace]4\[LetterSpace]cZTL\[LetterSpace]degr\[LetterSpace]1[cZTL_,def_,m20_] -> (cZTL*m20)/def, function\[LetterSpace]4\[LetterSpace]cZTL\[LetterSpace]trsl\[LetterSpace]1[def_,p14_] -> p14/def, function\[LetterSpace]4\[LetterSpace]cG\[LetterSpace]m\[LetterSpace]degr\[LetterSpace]1[cG\[LetterSpace]m_,def_,m18_] -> (cG\[LetterSpace]m*m18)/def, function\[LetterSpace]4\[LetterSpace]cNI\[LetterSpace]degr\[LetterSpace]1[D_,L_,cNI_,def_,m17_,m24_] -> (cNI*(L*m17 + D*m24))/def, function\[LetterSpace]4\[LetterSpace]cG\[LetterSpace]m\[LetterSpace]trscr\[LetterSpace]1[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, 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))
};

events = {

};

speciesAnnotations = {
cY\[LetterSpace]m[t]->"http://identifiers.org/chebi/CHEBI:33699", cY\[LetterSpace]m[t]->"http://identifiers.org/kegg.compound/C00046"
};

reactionAnnotations = {

};

units = {
{"time" -> "", "metabolite" -> "", "extent" -> ""}
};

(* Time evolution *)
odes = {
cG'[t] == 1.0*cG\[LetterSpace]trsl -1.0*cG\[LetterSpace]degr -1.0*cG\[LetterSpace]cZTL\[LetterSpace]assoc, cG\[LetterSpace]m'[t] == 1.0*cG\[LetterSpace]m\[LetterSpace]trscr -1.0*cG\[LetterSpace]m\[LetterSpace]degr, cL'[t] == 1.0*cL\[LetterSpace]trsl -1.0*cL\[LetterSpace]degr, cL\[LetterSpace]m'[t] == 1.0*cL\[LetterSpace]m\[LetterSpace]trscr -1.0*cL\[LetterSpace]m\[LetterSpace]degr, cLm'[t] == 1.0*cL\[LetterSpace]modif -1.0*cLm\[LetterSpace]degr, cNI'[t] == 1.0*cNI\[LetterSpace]trsl -1.0*cNI\[LetterSpace]degr, cNI\[LetterSpace]m'[t] == 1.0*cNI\[LetterSpace]m\[LetterSpace]trscr -1.0*cNI\[LetterSpace]m\[LetterSpace]degr, cP'[t] == 1.0*cP\[LetterSpace]trsl -1.0*cP\[LetterSpace]degr, cP7'[t] == 1.0*cP7\[LetterSpace]trsl -1.0*cP7\[LetterSpace]degr, cP7\[LetterSpace]m'[t] == 1.0*cP7\[LetterSpace]m\[LetterSpace]trscr -1.0*cP7\[LetterSpace]m\[LetterSpace]degr, cP9'[t] == 1.0*cP9\[LetterSpace]trsl -1.0*cP9\[LetterSpace]degr, cP9\[LetterSpace]m'[t] == 1.0*cP9\[LetterSpace]m\[LetterSpace]trscr -1.0*cP9\[LetterSpace]m\[LetterSpace]degr, cT'[t] == 1.0*cT\[LetterSpace]trsl -1.0*cT\[LetterSpace]degr, cT\[LetterSpace]m'[t] == 1.0*cT\[LetterSpace]m\[LetterSpace]trscr -1.0*cT\[LetterSpace]m\[LetterSpace]degr, cTm'[t] == 1.0*cT\[LetterSpace]modif -1.0*cTm\[LetterSpace]degr, cY'[t] == 1.0*cY\[LetterSpace]trsl -1.0*cY\[LetterSpace]degr, cY\[LetterSpace]m'[t] == 1.0*cY\[LetterSpace]m\[LetterSpace]trscr -1.0*cY\[LetterSpace]m\[LetterSpace]degr, cZG'[t] == 1.0*cG\[LetterSpace]cZTL\[LetterSpace]assoc -1.0*cZG\[LetterSpace]degr, cZTL'[t] == 1.0*cZTL\[LetterSpace]trsl -1.0*cG\[LetterSpace]cZTL\[LetterSpace]assoc -1.0*cZTL\[LetterSpace]degr
};
timeCourse = NDSolve[Join[odes, initialValues]//.rateEquations//.assignments//.parameters, variables, {t, 0, 100}];

(* Steady-state solution initialized with result of time evolution *)
findRootEquations = odes /.D[_[t],t]->0;
findRootVariables = Partition[Flatten[{#, #/.timeCourse/.t->100} &/@variables],2];
steadyStateVariables = FindRoot[findRootEquations//.rateEquations//.assignments//.parameters, findRootVariables, MaxIterations->100]
fluxes = #//.assignments//.parameters/.steadyStateVariables&/@rateEquations

(* Plot the time evolution of the variables *)
plotTable=Table[Plot[variables[[i]]/.parameters/.timeCourse,{t,0,100},PlotLegends->variables[[i]],PlotRange->Full],{i,Length[variables]}]