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

Corrected handling of min-1/hr-1 half life

Browse source 
Changed half life 'limit' value to 6.43 (as paper model, used when formula gives negative values).
This commit is contained in:
jdevine 2022-05-09 14:04:27 +02:00
parent 645bdd0197
commit f518a893d5
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/models.py
View file

@ -459,11 +459,14 @@ class SARSCoV2(Virus):
CERN-OPEN-2021-004, DOI: 10.17181/CERN.1GDQ.5Y75)
"""
# Updated to use the formula from Dabish et al. with correction https://doi.org/10.1080/02786826.2020.1829536
# with a minimum at hl = 1.1. Note that humidity is in percentage and inside_temp in °C.
return np.maximum(1.1, (0.693/((0.16030 + 0.04018*(((inside_temp-273.15)-20.615)/10.585)
# with a maximum at hl = 6.43. Note that humidity is in percentage and inside_temp in °C.
hl_calc = ((0.693/((0.16030 + 0.04018*(((inside_temp-273.15)-20.615)/10.585)
+0.02176*(((humidity*100)-45.235)/28.665)
-0.14369
-0.02636*((inside_temp-273.15)-20.615)/10.585))))
-0.02636*((inside_temp-273.15)-20.615)/10.585)))/60)
if (hl_calc <= 0):
hl_calc = 6.43
return hl_calc
Virus.types = {