begin name
banks1
end name

begin reactions
v1 {1.0}T1 + {1.0}Ttot = {1.0}$EXTERNAL
v10 {1.0}$EXTERNAL = {1.0}Tstar1 + {1.0}Ttot
v11 {1.0}$EXTERNAL = {1.0}Tstar1 + {1.0}Ttot
v12 {1.0}$EXTERNAL = {1.0}T2 + {1.0}Ttot
v13 {1.0}$EXTERNAL = {1.0}T2 + {1.0}Ttot
v14 {1.0}T2 + {1.0}Ttot = {1.0}$EXTERNAL
v15 {1.0}T2 = {1.0}Tstar2
v16 {1.0}Tstar2 = {1.0}T2
v17 {1.0}T2 + {1.0}Ttot = {1.0}$EXTERNAL
v18 {1.0}T2 + {1.0}Ttot = {1.0}$EXTERNAL
v19 {1.0}Tstar2 + {1.0}Ttot = {1.0}$EXTERNAL
v2 {1.0}T1 = {1.0}Tstar1
v20 {1.0}Tstar2 + {1.0}Ttot = {1.0}$EXTERNAL
v21 {1.0}Tstar2 + {1.0}Ttot = {1.0}$EXTERNAL
v22 {1.0}$EXTERNAL = {1.0}VI
v23 {1.0}VI = {1.0}VNI
v24 {1.0}VI = {1.0}$EXTERNAL
v25 {1.0}VI = {1.0}$EXTERNAL
v26 {1.0}$EXTERNAL = {1.0}VI
v27 {1.0}VI = {1.0}$EXTERNAL
v28 {1.0}$EXTERNAL = {1.0}VI
v29 {1.0}VNI = {1.0}$EXTERNAL
v3 {1.0}Tstar1 = {1.0}T1
v30 {1.0}$EXTERNAL = {1.0}E1
v31 {1.0}$EXTERNAL = {1.0}E1
v32 {1.0}E1 = {1.0}$EXTERNAL
v33 {1.0}E1 = {1.0}$EXTERNAL
v34 {1.0}E1 = {1.0}E2
v35 {1.0}$EXTERNAL = {1.0}E1
v36 {1.0}$EXTERNAL = {1.0}E2
v37 {1.0}E2 = {1.0}$EXTERNAL
v38 {1.0}E2 = {1.0}$EXTERNAL
v4 {1.0}T1 + {1.0}Ttot = {1.0}$EXTERNAL
v5 {1.0}$EXTERNAL = {1.0}T1 + {1.0}Ttot
v6 {1.0}$EXTERNAL = {1.0}T1 + {1.0}Ttot
v7 {1.0}Tstar1 + {1.0}Ttot = {1.0}$EXTERNAL
v8 {1.0}Tstar1 + {1.0}Ttot = {1.0}$EXTERNAL
v9 {1.0}Tstar1 = {1.0}Tstar2
end reactions

begin rate equations
v1 = d1*T1[t]
v10 = (aTS*pT*Tstar2[t]*VI[t])/(KV + VI[t])
v11 = aAS*pT*Tstar2[t]
v12 = (KS*lambdaT)/(KS + VI[t])
v13 = gammaT*T1[t]
v14 = d2*T2[t]
v15 = k2*T2[t]*VI[t]
v16 = f*k2*zeta1*T2[t]*VI[t]
v17 = (aT*T2[t]*VI[t])/(KV + VI[t])
v18 = aA*T2[t]
v19 = d2*Tstar2[t]
v2 = k1*T1[t]*VI[t]
v20 = (aTS*Tstar2[t]*VI[t])/(KV + VI[t])
v21 = aAS*Tstar2[t]
v22 = 1000*delta*NT*Tstar1[t]
v23 = 1000*delta*NT*zeta2*Tstar1[t]
v24 = c*VI[t]
v25 = 1000*k1*rho1*T1[t]*VI[t]
v26 = 1000*k1*rho1*zeta1*T1[t]*VI[t]
v27 = 1000*k2*rho2*T2[t]*VI[t]
v28 = 1000*f*k2*rho2*zeta1*T2[t]*VI[t]
v29 = c*VNI[t]
v3 = k1*zeta1*T1[t]*VI[t]
v30 = lambdaE
v31 = (bE1*E1[t]*Tstar1[t])/(Kb1 + Tstar1[t])
v32 = (dE*E1[t]*Tstar1[t])/(Kd + Tstar1[t])
v33 = deltaE1*E1[t]
v34 = (gammaE*E1[t]*(T1[t] + Tstar1[t]))/(Kgamma + T1[t] + Tstar1[t])
v35 = (aE*pE*E2[t]*VI[t])/(KV + VI[t])
v36 = (bE2*Kb2*E2[t])/(Kb2 + E2[t])
v37 = deltaE2*E2[t]
v38 = (aE*E2[t]*VI[t])/(KV + VI[t])
v4 = gammaT*T1[t]
v5 = (aT*pT*T2[t]*VI[t])/(KV + VI[t])
v6 = aA*pT*T2[t]
v7 = delta*Tstar1[t]
v8 = m*E1[t]*Tstar1[t]
v9 = gammaTS*Tstar1[t]
end rate equations

begin parameters
KS = 27900.0
KV = 1060.0
Kb1 = 0.0249
Kb2 = 87.0
Kd = 0.12
Kgamma = 1.36
NT = 11.7
aA = 0.00407
aAS = "9.58e-05"
aE = 0.0168
aT = 0.00985
aTS = 0.000449
bE1 = 0.0693
bE2 = 0.0061
c = 12.7
d1 = 0.0912
d2 = 0.0031
dE = 0.0472
delta = 0.171
deltaE1 = 0.0597
deltaE2 = 0.00145
f = 0.507
gammaE = 0.000689
gammaT = 0.000646
gammaTS = "7.98e-06"
k1 = "7.97e-06"
k2 = "1.05e-08"
lambdaE = 0.000488
lambdaT = 8.66
m = 1.1
pE = 0.395
pT = 9.85
rho1 = 1.0
rho2 = 1.0
zeta1 = 0.524
zeta2 = 0.16
default = 1.0
end parameters

begin initial conditions
E1[0] = E1i
E2[0] = E2i
T1[0] = T1i
T2[0] = T2i
Tstar1[0] = Tstar1i
Tstar2[0] = Tstar2i
Ttot[0] = Ttoti
VI[0] = VIi
VNI[0] = VNIi
end initial conditions

begin initial values
E1i = 0.0682
E2i = 0.691
T1i = 79.5
T2i = 289.0
Tstar1i = 58.2
Tstar2i = 0.0654
Ttoti = 426.7654
VIi = 150000.0
VNIi = 3570.0
end initial values

begin assignment rules
viruscopiesperml := VI[t] + VNI[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