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

changed linespace range

Browse source
This commit is contained in:
Luis Aleixo 2021-09-14 17:36:53 +02:00
parent c4d43c099c
commit cb7d8abeb8
1 changed files with 3 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

5
cara/results_paper.py
View file

@ -507,7 +507,7 @@ def calculate_deposition_factor():
br_heavy_exercise = breathing_heavy_exercise_exposure() br_heavy_exercise = breathing_heavy_exercise_exposure()
br_heavy_exercise_model = br_heavy_exercise.build_model(size=SAMPLE_SIZE) br_heavy_exercise_model = br_heavy_exercise.build_model(size=SAMPLE_SIZE)
rho_p = 1.2 rho_p = 1000
mu_air = 1.8*10**-5 mu_air = 1.8*10**-5
FRC = 0.003 FRC = 0.003
Vt = 0.0004 Vt = 0.0004
@ -515,7 +515,8 @@ def calculate_deposition_factor():
k = 1.38*10**-23 k = 1.38*10**-23
T = 300 T = 300
diameters = np.linspace(0.001, 100, 400) #particle diameter (multiply later by 10**(-6)) diameters = np.linspace(0.001, 0.01, 1000) #particle diameter (multiply later by 10**(-6))
diameters = np.append(diameters, np.linspace(0.01, 100, 400))
activity_fractions = [] activity_fractions = []
for scenario in (br_seated_model, br_light_activity_model, br_heavy_exercise_model): for scenario in (br_seated_model, br_light_activity_model, br_heavy_exercise_model):
BRk = scenario.exposed.activity.inhalation_rate BRk = scenario.exposed.activity.inhalation_rate