begin name
twycross
end name

begin reactions
v1 {1.0}S = {40000.0}a0
v10 {40000.0}a4 = {20.0}c5
v11 {20.0}c5 = {40000.0}a5
v12 {40000.0}a5 = {20.0}c6
v13 {20.0}c6 = {40000.0}a6
v14 {40000.0}a6 = {20.0}c7
v15 {20.0}c7 = {40000.0}a7
v16 {40000.0}a7 = {20.0}c8
v17 {20.0}c8 = {40000.0}a8
v18 {40000.0}a8 = {20.0}c9
v19 {20.0}c9 = {40000.0}a9
v2 {40000.0}a0 = {20.0}c1
v20 {40000.0}a9 = {20.0}c10
v21 {20.0}c10 = {40000.0}a10
v22 {40000.0}a10 = {20.0}c11
v23 {20.0}c11 = {40000.0}a11
v24 {40000.0}a11 = {20.0}c12
v25 {20.0}c12 = {40000.0}a12
v26 {40000.0}a12 = {20.0}c13
v27 {20.0}c13 = {40000.0}a13
v28 {40000.0}a13 = {20.0}c14
v29 {20.0}c14 = {40000.0}a14
v3 {20.0}c1 = {40000.0}a1
v30 {40000.0}a14 = {20.0}c15
v31 {20.0}c15 = {40000.0}a15
v32 {40000.0}a15 = {20.0}c16
v33 {20.0}c16 = {40000.0}a16
v34 {40000.0}a16 = {20.0}c17
v35 {20.0}c17 = {40000.0}a17
v36 {40000.0}a17 = {20.0}c18
v37 {20.0}c18 = {40000.0}a18
v38 {40000.0}a18 = {20.0}c19
v39 {20.0}c19 = {40000.0}a19
v4 {40000.0}a1 = {20.0}c2
v40 {40000.0}a19 = {20.0}c20
v41 {20.0}c20 = {40000.0}a20
v42 {40000.0}a20 = {1.0}Fvar
v5 {20.0}c2 = {40000.0}a2
v6 {40000.0}a2 = {20.0}c3
v7 {20.0}c3 = {40000.0}a3
v8 {40000.0}a3 = {20.0}c4
v9 {20.0}c4 = {40000.0}a4
end reactions

begin rate equations
v1 = (2*Dvar*kVol*(-a0[t] + S[t]))/(lambda + Ls)
v10 = kVol*Pdiff*(Aprota*a4[t] - Aprotc*c5[t])
v11 = kVol*(Pdiff*(-(Aprota*a5[t]) + Aprotc*c5[t]) + PPin*(-(Ba*a5[t]) + Bc*c5[t]))
v12 = kVol*Pdiff*(Aprota*a5[t] - Aprotc*c6[t])
v13 = kVol*(Pdiff*(-(Aprota*a6[t]) + Aprotc*c6[t]) + PPin*(-(Ba*a6[t]) + Bc*c6[t]))
v14 = kVol*Pdiff*(Aprota*a6[t] - Aprotc*c7[t])
v15 = kVol*(Pdiff*(-(Aprota*a7[t]) + Aprotc*c7[t]) + PPin*(-(Ba*a7[t]) + Bc*c7[t]))
v16 = kVol*Pdiff*(Aprota*a7[t] - Aprotc*c8[t])
v17 = kVol*(Pdiff*(-(Aprota*a8[t]) + Aprotc*c8[t]) + PPin*(-(Ba*a8[t]) + Bc*c8[t]))
v18 = kVol*Pdiff*(Aprota*a8[t] - Aprotc*c9[t])
v19 = kVol*(Pdiff*(-(Aprota*a9[t]) + Aprotc*c10[t]) + PPin*(-(Ba*a9[t]) + Bc*c9[t]))
v2 = kVol*Pdiff*(Aprota*a0[t] - Aprotc*c1[t])
v20 = kVol*Pdiff*(Aprota*a9[t] - Aprotc*c10[t])
v21 = kVol*(Pdiff*(-(Aprota*a10[t]) + Aprotc*c10[t]) + PPin*(-(Ba*a10[t]) + Bc*c10[t]))
v22 = kVol*Pdiff*(Aprota*a10[t] - Aprotc*c11[t])
v23 = kVol*(Pdiff*(-(Aprota*a11[t]) + Aprotc*c11[t]) + PPin*(-(Ba*a11[t]) + Bc*c11[t]))
v24 = kVol*Pdiff*(Aprota*a11[t] - Aprotc*c12[t])
v25 = kVol*(Pdiff*(-(Aprota*a12[t]) + Aprotc*c12[t]) + PPin*(-(Ba*a12[t]) + Bc*c12[t]))
v26 = kVol*Pdiff*(Aprota*a12[t] - Aprotc*c13[t])
v27 = kVol*(Pdiff*(-(Aprota*a13[t]) + Aprotc*c13[t]) + PPin*(-(Ba*a13[t]) + Bc*c13[t]))
v28 = kVol*Pdiff*(Aprota*a13[t] - Aprotc*c14[t])
v29 = kVol*(Pdiff*(-(Aprota*a14[t]) + Aprotc*c14[t]) + PPin*(-(Ba*a14[t]) + Bc*c14[t]))
v3 = kVol*(Pdiff*(-(Aprota*a1[t]) + Aprotc*c1[t]) + PPin*(-(Ba*a1[t]) + Bc*c1[t]))
v30 = kVol*Pdiff*(Aprota*a14[t] - Aprotc*c15[t])
v31 = kVol*(Pdiff*(-(Aprota*a15[t]) + Aprotc*c15[t]) + PPin*(-(Ba*a15[t]) + Bc*c15[t]))
v32 = kVol*Pdiff*(Aprota*a15[t] - Aprotc*c16[t])
v33 = kVol*(Pdiff*(-(Aprota*a16[t]) + Aprotc*c16[t]) + PPin*(-(Ba*a16[t]) + Bc*c16[t]))
v34 = kVol*Pdiff*(Aprota*a16[t] - Aprotc*c17[t])
v35 = kVol*(Pdiff*(-(Aprota*a17[t]) + Aprotc*c17[t]) + PPin*(-(Ba*a17[t]) + Bc*c17[t]))
v36 = kVol*Pdiff*(Aprota*a17[t] - Aprotc*c18[t])
v37 = kVol*(Pdiff*(-(Aprota*a18[t]) + Aprotc*c18[t]) + PPin*(-(Ba*a18[t]) + Bc*c18[t]))
v38 = kVol*Pdiff*(Aprota*a18[t] - Aprotc*c19[t])
v39 = kVol*(Pdiff*(-(Aprota*a19[t]) + Aprotc*c19[t]) + PPin*(-(Ba*a19[t]) + Bc*c19[t]))
v4 = kVol*Pdiff*(Aprota*a1[t] - Aprotc*c2[t])
v40 = kVol*Pdiff*(Aprota*a19[t] - Aprotc*c20[t])
v41 = kVol*(Pdiff*(-(Aprota*a20[t]) + Aprotc*c20[t]) + PPin*(-(Ba*a20[t]) + Bc*c20[t]))
v42 = (2*Dvar*kVol*(a20[t] - Fvar[t]))/(lambda + Ls)
v5 = kVol*(Pdiff*(-(Aprota*a2[t]) + Aprotc*c2[t]) + PPin*(-(Ba*a2[t]) + Bc*c2[t]))
v6 = kVol*Pdiff*(Aprota*a2[t] - Aprotc*c3[t])
v7 = kVol*(Pdiff*(-(Aprota*a3[t]) + Aprotc*c3[t]) + PPin*(-(Ba*a3[t]) + Bc*c3[t]))
v8 = kVol*Pdiff*(Aprota*a3[t] - Aprotc*c4[t])
v9 = kVol*(Pdiff*(-(Aprota*a4[t]) + Aprotc*c4[t]) + PPin*(-(Ba*a4[t]) + Bc*c4[t]))
end rate equations

begin parameters
Dvar = "6.7e-10"
F = 96485.3399
L = 0.002
Ls = 0.002
PPin = "3.3e-06"
Pdiff = "5.6e-07"
R = 8.314472
T = 295.15
V = -0.12
ext = 0.0
kVol = 500.0
lambda = "5e-07"
lcyt = 0.0001
pHa = 5.3
pHc = 7.2
pK = 4.8
w = "1e-05"
defaultcompartment = 1.0
end parameters

begin initial conditions
Fvar[0] = Fvari
S[0] = Si
a0[0] = a0i
a1[0] = a1i
a10[0] = a10i
a11[0] = a11i
a12[0] = a12i
a13[0] = a13i
a14[0] = a14i
a15[0] = a15i
a16[0] = a16i
a17[0] = a17i
a18[0] = a18i
a19[0] = a19i
a2[0] = a2i
a20[0] = a20i
a3[0] = a3i
a4[0] = a4i
a5[0] = a5i
a6[0] = a6i
a7[0] = a7i
a8[0] = a8i
a9[0] = a9i
c1[0] = c1i
c10[0] = c10i
c11[0] = c11i
c12[0] = c12i
c13[0] = c13i
c14[0] = c14i
c15[0] = c15i
c16[0] = c16i
c17[0] = c17i
c18[0] = c18i
c19[0] = c19i
c2[0] = c2i
c20[0] = c20i
c3[0] = c3i
c4[0] = c4i
c5[0] = c5i
c6[0] = c6i
c7[0] = c7i
c8[0] = c8i
c9[0] = c9i
end initial conditions

begin initial values
Fvari = 0.0
Si = 1.0
a0i = 0.0
a1i = 0.0
a10i = 0.0
a11i = 0.0
a12i = 0.0
a13i = 0.0
a14i = 0.0
a15i = 0.0
a16i = 0.0
a17i = 0.0
a18i = 0.0
a19i = 0.0
a2i = 0.0
a20i = 0.0
a3i = 0.0
a4i = 0.0
a5i = 0.0
a6i = 0.0
a7i = 0.0
a8i = 0.0
a9i = 0.0
c1i = 0.0
c10i = 0.0
c11i = 0.0
c12i = 0.0
c13i = 0.0
c14i = 0.0
c15i = 0.0
c16i = 0.0
c17i = 0.0
c18i = 0.0
c19i = 0.0
c2i = 0.0
c20i = 0.0
c3i = 0.0
c4i = 0.0
c5i = 0.0
c6i = 0.0
c7i = 0.0
c8i = 0.0
c9i = 0.0
end initial values

begin assignment rules
Aprotc := (1 + 10^(pHc - pK))^(-1)
Aaniona := 1 - Aprota
Bc := -((Aanionc*subval)/(-1 + 2.71828^(-subval)))
Aprota := (1 + 10^(pHa - pK))^(-1)
Ba := (Aaniona*subval)/(-1 + 2.71828^subval)
subval := -((F*V)/(R*T))
Aanionc := 1 - Aprotc
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