#essunger3



#reactions
CD4TproductionV:
 $pool = {1.0}V
  muAA * N0 * (1 + gamma * time**2 / (time**2 + Tc**2)) * Tastarstar
VinfectionCD4T:
 {1.0}Ta + {1.0}V = {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
  muVV * V
deathCD4TactiveInfected:
 {1.0}Ttot + {1.0}Tastarstar = $pool
  muAA * Tastarstar
deathCD4TactivePlain:
 {1.0}Ta + {1.0}Ttot = $pool
  muA * Ta
deathCD4TrestingMemoryInfected:
 {1.0}Tmstarstar + {1.0}Ttot = $pool
  muM * Tmstarstar
deathCD4TrestingMemoryPlain:
 {1.0}Ttot + {1.0}Tm = $pool
  muM * Tm
deathCD4TrestingNaivePlain:
 {1.0}Ttot + {1.0}Tv = $pool
  muV * Tv
fractionVinfectionCD4T:
 {1.0}Tastarstar = {1.0}Tmstarstar
  kI * phi * Ta * V
growthCD4T:
 {1.0}Ta = {1.0}Ttot
  r * Ta * (1 - Ttot / Tmax)
sourceCD4T:
 $pool = {1.0}Ttot + {1.0}Tv
  s0 * s1 / (s1 + V**(1/3))
v15:
 $pool = {1.0}Tm
  2*r*Ta*(1 - Ttot/Tmax)
v16:
 $pool = {1.0}Tmstarstar
  2*rStar*Tastarstar*(1 - Ttot/Tmax)
v9:
 {1.0}Tastarstar = {1.0}Ttot
  rStar*Tastarstar*(1 - Ttot/Tmax)


#parameters
N0 = 600.0
Tc = 15000.0
Tmax = 1500.0
gamma = 130.0
kA = 0.00495
kB = 0.0055
kI = 0.0007
muA = 0.02
muAA = 0.5
muM = 0.013
muV = 0.013
muVV = 0.75
phi = 1.0
r = 1.5
rStar = 0.0
s0 = 8.076
s1 = 2.0
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