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

add unit conversions in qr-calculation

Browse source
This commit is contained in:
markus 2021-01-18 16:32:04 +01:00
parent 03cad9d505
commit ecc35ee03c
1 changed files with 5 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

5
cara/montecarlo.py
View file

@ -15,6 +15,11 @@ def calculate_qr(viral_load: float, emission: float, diameter: float, mask_effic
"""
Calculates the quantum generation rate given a set of parameters.
"""
# Unit conversions
diameter *= 1e-4
viral_load = 10 ** viral_load
emission = (emission * 60) if breathing_rate is None else (emission * 1e6)
volume = 4 * np.pi * (diameter / 2) ** 3 / 3
if breathing_rate is None:
breathing_rate = 1