From 11561c10e8e94df15652c3937b3673767935bd96 Mon Sep 17 00:00:00 2001 From: markus Date: Mon, 18 Jan 2021 17:00:33 +0100 Subject: [PATCH] fix time conversion factor --- cara/montecarlo.py | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/cara/montecarlo.py b/cara/montecarlo.py index 5ef26711..02d16026 100644 --- a/cara/montecarlo.py +++ b/cara/montecarlo.py @@ -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))