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

fix compare_concentration_curves

Browse source
This commit is contained in:
markus 2021-02-24 13:18:20 +01:00
parent 04dfa59fce
commit 646b005baa
1 changed files with 6 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

9
cara/montecarlo.py
View file

@ -1036,16 +1036,19 @@ def compare_concentration_curves(exp_models: typing.List[MCExposureModel], label
times = np.arange(start, stop, TIME_STEP) times = np.arange(start, stop, TIME_STEP)
concentrations = [[np.mean(model.concentration_model.concentration(t)) for t in times] for model in exp_models] concentrations = [[np.median(model.concentration_model.concentration(t)) for t in times] for model in exp_models]
for c, label, color in zip(concentrations, labels, colors): for c, label, color in zip(concentrations, labels, colors):
plt.plot(times, c, label=label, color=color) plt.plot(times, c, label=label, color=color)
plt.legend() plt.legend()
factors = [0.6 * model.exposed.activity.inhalation_rate * (1 - model.exposed.mask.η_inhale) for model in exp_models]
qds = [[np.trapz(c[:i + 1], times[:i + 1]) * factor for i in range(len(times))]
for c, factor in zip(concentrations, factors)]
if show_qd: if show_qd:
plt.twinx() plt.twinx()
for c, label, color in zip(concentrations, labels, colors): for qd, label, color in zip(qds, labels, colors):
plt.plot(times, np.cumsum(c), label='qD - ' + label, color=color, linestyle='dotted') plt.plot(times, qd, label='qD - ' + label, color=color, linestyle='dotted')
plt.legend() plt.legend()
plt.show() plt.show()