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

fix time conversion factor

Browse source
This commit is contained in:
markus 2021-01-18 17:00:33 +01:00
parent b3327506fe
commit 11561c10e8
1 changed files with 3 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
cara/montecarlo.py
View file

@ -20,7 +20,7 @@ def calculate_qr(viral_load: float, emission: float, diameter: float, mask_effic
# Unit conversions
diameter *= 1e-4
viral_load = 10 ** viral_load
emission = (emission * 60) if breathing_rate is None else (emission * 1e6)
emission = (emission * 3600) if breathing_rate is None else (emission * 1e6)
volume = 4 * np.pi * (diameter / 2) ** 3 / 3
if breathing_rate is None:
@ -49,3 +49,5 @@ def generate_qr_values(samples: int, expiratory_activity: int, qid: int = 100) -
breathing_rate = 1
return qr_func(viral_load, emissions, diameters, mask_efficiency, qid)
print(generate_qr_values(samples=5, expiratory_activity=1))