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

Merge remote-tracking branch 'origin/feature/mc' into feature/mc

Browse source
This commit is contained in:
Andrejh 2021-02-26 16:03:07 +01:00
commit 762c147956
1 changed files with 6 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

7
cara/montecarlo.py
View file

@ -1063,8 +1063,13 @@ def compare_concentration_curves(exp_models: typing.List[MCExposureModel], label
ax.set_ylim(ax.get_ylim()[0], ax.get_ylim()[1] * 1.2)
factors = [0.6 * model.exposed.activity.inhalation_rate * (1 - model.exposed.mask.η_inhale) for model in exp_models]
present_indexes = np.array([exp_models[0].exposed.person_present(t) for t in times])
modified_concentrations = [np.array(c) for c in concentrations]
for mc in modified_concentrations:
mc[~present_indexes] = 0
qds = [[np.trapz(c[:i + 1], times[:i + 1]) * factor for i in range(len(times))]
for c, factor in zip(concentrations, factors)]
for c, factor in zip(modified_concentrations, factors)]
plt.suptitle("")
plt.xlabel("Time ($h$)")