#banks1



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


#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


#initial values
E1 = 0.0682
E2 = 0.691
T1 = 79.5
T2 = 289.0
Tstar1 = 58.2
Tstar2 = 0.0654
Ttot = 426.7654
VI = 150000.0
VNI = 3570.0


#assignment rules
!F viruscopiesperml = VNI + VI