#lou1



#reactions
v_1:
 $pool = {1.0}St
  rt
v_10:
 {1.0}Sb = {1.0}Ib
  lambdaB * Sb
v_11:
 {1.0}Sb = $pool
  dM * Sb
v_12:
 {1.0}Ib = $pool
  dI * Ib
v_2:
 {1.0}St = {1.0}It
  lambdaT * St
v_3:
 {1.0}St = $pool
  dM * St
v_4:
 {1.0}It = $pool
  dI * It
v_5:
 $pool = {1.0}Sv
  rv
v_6:
 {1.0}Sv = {1.0}Iv
  lambdaV * Sv
v_7:
 {1.0}Sv = $pool
  dM * Sv
v_8:
 {1.0}Iv = $pool
  dI * Iv
v_9:
 $pool = {1.0}Sb
  rb


#parameters
betaBT = 0.2563
betaBV = 0.2826
betaTB = 0.5786
betaTV = 0.5653
betaVB = 0.5786
betaVT = 0.2563
betaVV = 0.4239
dI = 0.105
dM = 0.022
rb = 70817.0
rt = 96667.0
rv = 603950.0
default = 1.0


#initial values
Ib = 39455.0
It = 36931.0
Iv = 234800.0
Sb = 2590800.0
St = 2425100.0
Sv = 15418000.0


#assignment rules
!F lambdaT = (betaBT * Ib + betaVT * Iv) / (Sb + Ib + Sv + Iv)
!F lambdaV = (betaTV * It + betaBV * Ib + betaVV * Iv) / (St + It + Sv + Iv + Sb + Ib)
!F lambdaB = (betaTB * It + betaVB * Iv) / (St + It + Sv + Iv)
!F Total_infected_MSM = It + Iv + Ib