begin name
montagne1
end name

begin reactions
assaal {1.0}T1 + {1.0}alpha = {1.0}alphaT1
assabl {1.0}alpha + {1.0}T2 = {1.0}alphaT2
assbcl {1.0}beta + {1.0}T3 = {1.0}betaT3
assbcr {1.0}T3 + {1.0}Inh = {1.0}T3Inh
assbcrl {1.0}beta + {1.0}T3Inh = {1.0}betaT3Inh
disaa {1.0}alphaT1alpha = {1.0}alpha + {1.0}alphaalphaT1
disab {1.0}alphaT2beta = {1.0}beta + {1.0}alphabetaT2
disbc {1.0}betaT3Inh = {1.0}Inh + {1.0}betaInhT3
exoa {1.0}alpha = {1.0}$empty
exob {1.0}beta = {1.0}$empty
exoc {1.0}Inh = {1.0}$empty
inhac {1.0}T1 + {1.0}Inh = {1.0}InhT1
inhdisplac {1.0}T1alpha + {1.0}Inh = {1.0}alpha + {1.0}InhT1
massaalr {1.0}alphaT1alpha = {1.0}alpha + {1.0}alphaT1
massaar {1.0}T1alpha = {1.0}T1 + {1.0}alpha
massaarl {1.0}alphaT1alpha = {1.0}alpha + {1.0}T1alpha
massablr {1.0}alphaT2beta = {1.0}alphaT2 + {1.0}beta
massabr {1.0}T2beta = {1.0}T2 + {1.0}beta
massabrl {1.0}alphaT2beta = {1.0}alpha + {1.0}T2beta
massbclr {1.0}betaT3Inh = {1.0}betaT3 + {1.0}Inh
minhdisplca {1.0}alpha + {1.0}InhT1 = {1.0}alphaT1 + {1.0}Inh
nickaa {1.0}alphaalphaT1 = {1.0}alphaT1alpha
nickab {1.0}alphabetaT2 = {1.0}alphaT2beta
nickbc {1.0}betaInhT3 = {1.0}betaT3Inh
polaa {1.0}alphaT1 = {1.0}alphaalphaT1
polab {1.0}alphaT2 = {1.0}alphabetaT2
polbc {1.0}betaT3 = {1.0}betaInhT3
end reactions

begin rate equations
assaal = sample*(-(k0r*alphaT1[t]) + k0d*alpha[t]*T1[t])
assabl = sample*(-(k7r*alphaT2[t]) + k7d*alpha[t]*T2[t])
assbcl = sample*(-(k14r*betaT3[t]) + k14d*beta[t]*T3[t])
assbcr = sample*(k16d*Inh[t]*T3[t] - k16r*T3Inh[t])
assbcrl = sample*(-(k17r*betaT3Inh[t]) + k17d*beta[t]*T3Inh[t])
disaa = k5d*sample*alphaT1alpha[t]
disab = k12d*sample*alphaT2beta[t]
disbc = k19d*sample*betaT3Inh[t]
exoa = k24d*sample*alpha[t]
exob = k25d*sample*beta[t]
exoc = k26d*sample*Inh[t]
inhac = sample*(-(k21r*InhT1[t]) + k21d*Inh[t]*T1[t])
inhdisplac = sample*(-(k22r*alpha[t]*InhT1[t]) + k22d*Inh[t]*T1alpha[t])
massaalr = sample*(-(k1r*alpha[t]*alphaT1[t]) + k1d*alphaT1alpha[t])
massaar = sample*(-(k2r*alpha[t]*T1[t]) + k2d*T1alpha[t])
massaarl = sample*(k3d*alphaT1alpha[t] - k3r*alpha[t]*T1alpha[t])
massablr = sample*(k8d*alphaT2beta[t] - k8r*alphaT2[t]*beta[t])
massabr = sample*(-(k9r*beta[t]*T2[t]) + k9d*T2beta[t])
massabrl = sample*(k10d*alphaT2beta[t] - k10r*alpha[t]*T2beta[t])
massbclr = sample*(k15d*betaT3Inh[t] - k15r*betaT3[t]*Inh[t])
minhdisplca = sample*(-(k23r*alphaT1[t]*Inh[t]) + k23d*alpha[t]*InhT1[t])
nickaa = k6d*sample*alphaalphaT1[t]
nickab = k13d*sample*alphabetaT2[t]
nickbc = k20d*sample*betaInhT3[t]
polaa = k4d*sample*alphaT1[t]
polab = k11d*sample*alphaT2[t]
polbc = k18d*sample*betaT3[t]
end rate equations

begin parameters
k0d = 0.0294
k0r = 3.43457943925
k10d = 3.43457943925
k10r = 0.0294
k11d = 11.8408
k12d = 9.2239832
k13d = 2.60186
k14d = 0.0171
k14r = 0.610714285714
k15d = 0.00186296832954
k15r = 0.027
k16d = 0.027
k16r = 0.00186296832954
k17d = 0.0171
k17r = 0.610714285714
k18d = 17.024
k19d = 5.566848
k1d = 3.43457943925
k1r = 0.0294
k20d = 3.2064
k21d = 0.027
k21r = 0.00608108108108
k22d = 0.021546
k22r = "4.15391351351e-05"
k23d = "4.15391351351e-05"
k23r = 0.021546
k24d = 0.411
k25d = 0.485802
k26d = 1.7262
k2d = 3.43457943925
k2r = 0.0294
k3d = 3.43457943925
k3r = 0.0294
k4d = 15.2
k5d = 11.8408
k6d = 3.34
k7d = 0.0294
k7r = 3.43457943925
k8d = 0.610714285714
k8r = 0.0171
k9d = 0.610714285714
k9r = 0.0171
empty = 0.0
sample = 1.0
end parameters

begin initial conditions
Inh[0] = Inhi
InhT1[0] = InhT1i
T1[0] = T1i
T1alpha[0] = T1alphai
T2[0] = T2i
T2beta[0] = T2betai
T3[0] = T3i
T3Inh[0] = T3Inhi
alpha[0] = alphai
alphaT1[0] = alphaT1i
alphaT1alpha[0] = alphaT1alphai
alphaT2[0] = alphaT2i
alphaT2beta[0] = alphaT2betai
alphaalphaT1[0] = alphaalphaT1i
alphabetaT2[0] = alphabetaT2i
beta[0] = betai
betaInhT3[0] = betaInhT3i
betaT3[0] = betaT3i
betaT3Inh[0] = betaT3Inhi
end initial conditions

begin initial values
Inhi = 0.0
InhT1i = 0.0
T1i = 38.5
T1alphai = 0.0
T2i = 3.89
T2betai = 0.0
T3i = 38.5
T3Inhi = 0.0
alphai = 10.0
alphaT1i = 0.0
alphaT1alphai = 0.0
alphaT2i = 0.0
alphaT2betai = 0.0
alphaalphaT1i = 0.0
alphabetaT2i = 0.0
betai = 0.0
betaInhT3i = 0.0
betaT3i = 0.0
betaT3Inhi = 0.0
end initial values

begin assignment rules
alphatotal := alpha[t] + alphaT1[t] + 2*alphaT1alpha[t] + 2*alphaT2[t] + alphaT2beta[t] + T1alpha[t]
betatotal := alphaT2beta[t] + beta[t] + betaT3[t] + betaT3Inh[t] + T2beta[t]
Inhtotal := betaT3Inh[t] + Inh[t] + InhT1[t] + T3Inh[t]
fluorescence := -0.38 + 0.00093*(22*(alphaalphaT1[t] + alphabetaT2[t] + alphaT1alpha[t] + alphaT2beta[t]) + 27*(betaInhT3[t] + betaT3Inh[t]) + 11*(alphaT1[t] + alphaT2[t] + betaT3[t] + T1alpha[t] + T2beta[t]) + 16*(InhT1[t] + T3Inh[t]))
Bpconcentration := 22*(alphaalphaT1[t] + alphabetaT2[t] + alphaT1alpha[t] + alphaT2beta[t]) + 27*(betaInhT3[t] + betaT3Inh[t]) + 11*(alphaT1[t] + alphaT2[t] + betaT3[t] + T1alpha[t] + T2beta[t]) + 16*(InhT1[t] + T3Inh[t])
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units