#essunger1



#reactions
CD4TproductionV:
 $pool = {1.0}V
  muastarstar * NN * Tastarstar
VinfectionCD4T:
 {1.0}V + {1.0}Ta = {1.0}Tastarstar
  ki * Ta * V
activationCD4TrestingMemoryInfected:
 {1.0}Tmstarstar = {1.0}Tastarstar
  kb * Tmstarstar
activationCD4TrestingMemoryPlain:
 {1.0}Tm = {1.0}Ta
  kb * Tm
activationCD4TrestingNaivePlain:
 {1.0}Tv = {1.0}Ta
  ka * Tv
clearanceV:
 {1.0}V = $pool
  muV * V
deathCD4TactiveInfected:
 {1.0}Tastarstar + {1.0}Ttot = $pool
  muastarstar * Tastarstar
deathCD4TactivePlain:
 {1.0}Ttot + {1.0}Ta = $pool
  mua * Ta
deathCD4TrestingMemoryInfected:
 {1.0}Tmstarstar + {1.0}Ttot = $pool
  mum * Tmstarstar
deathCD4TrestingMemoryPlain:
 {1.0}Tm + {1.0}Ttot = $pool
  mum * Tm
deathCD4TrestingNaivePlain:
 {1.0}Ttot + {1.0}Tv = $pool
  muv * Tv
fractionVinfectionCD4T:
 {1.0}Tastarstar = {1.0}Tmstarstar
  ki * phi * Ta * V
sourceCD4T:
 $pool = {1.0}Tv + {1.0}Ttot
  s
v10:
 {1.0}Ttot + {1.0}Tm = $pool
  2 * r * Ta * Ttot / Tmax
v15:
 {1.0}Ttot + {1.0}Tastarstar = $pool
  rstar * Tastarstar
v16:
 $pool = {1.0}Ttot + {1.0}Tastarstar
  rstar * Tastarstar * Ttot / Tmax
v17:
 $pool = {1.0}Ttot + {1.0}Tmstarstar
  2 * rstar * Tastarstar
v18:
 {1.0}Ttot + {1.0}Tmstarstar = $pool
  2 * rstar * Tastarstar * Ttot / Tmax
v5:
 {1.0}Ttot + {1.0}Ta = $pool
  r * Ta
v6:
 $pool = {1.0}Ttot + {1.0}Ta
  r * Ta * Ttot / Tmax
v9:
 $pool = {1.0}Ttot + {1.0}Tm
  2 * r * Ta


#parameters
NN = 600.0
Tmax = 1500.0
ka = 0.00495
kb = 0.0055
ki = 0.0007
muV = 0.75
mua = 0.02
muastarstar = 0.5
mum = 0.013
muv = 0.013
phi = 1.0
r = 1.5
rstar = 0.0
s = 8.076
default = 1.0


#initial values
Ta = 10.0
Tastarstar = 0.0
Tm = 540.0
Tmstarstar = 0.0
Ttot = 1000.0
Tv = 450.0
V = 5e-05


#assignment rules