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

Fixing bug pointed in issue #262

Browse source
This commit is contained in:
Nicolas Mounet 2022-04-25 17:25:59 +02:00
parent ac49e102d4
commit 2d1f7fa823
1 changed files with 4 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

7
cara/models.py
View file

@ -1356,7 +1356,7 @@ class ExposureModel:
# to perform properly the Monte-Carlo integration over
# particle diameters (doing things in another order would
# lead to wrong results for the probability of infection).
deposited_exposure += (np.array(short_range_jet_exposure
this_deposited_exposure = (np.array(short_range_jet_exposure
* fdep).mean()/dilution
- np.array(short_range_lr_exposure * fdep).mean()
* self.concentration_model.infected.activity.exhalation_rate
@ -1364,14 +1364,15 @@ class ExposureModel:
else:
# in the case of a single diameter or no diameter defined,
# one should not take any mean at this stage.
deposited_exposure += (short_range_jet_exposure
this_deposited_exposure = (short_range_jet_exposure
* fdep/dilution
- short_range_lr_exposure * fdep
* self.concentration_model.infected.activity.exhalation_rate
/dilution)
# multiply by the (diameter-independent) inhalation rate
deposited_exposure *= interaction.activity.inhalation_rate
deposited_exposure += (this_deposited_exposure *
interaction.activity.inhalation_rate)
# then we multiply by diameter-independent quantities: viral load
# and fraction of infected virions