pingu_98/cara
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Using concentration_limit method in concentration (ConcentrationModel)

Browse source
This commit is contained in:
Nicolas Mounet 2021-05-04 06:43:27 +02:00
parent 68893f905f
commit 94d01104d0
1 changed files with 1 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

3
cara/models.py
View file

@ -694,14 +694,13 @@ class ConcentrationModel:
if time == 0:
return 0.0
IVRR = self.infectious_virus_removal_rate(time)
V = self.room.volume
t_last_state_change = self.last_state_change(time)
concentration_at_last_state_change = self.concentration(t_last_state_change)
delta_time = time - t_last_state_change
fac = np.exp(-IVRR * delta_time)
concentration_limit = (self.infected.emission_rate(time)) / (IVRR * V)
concentration_limit = self.concentration_limit(time)
return concentration_limit * (1 - fac) + concentration_at_last_state_change * fac