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

print values for paper

Browse source
This commit is contained in:
Andrejh 2021-03-22 17:41:51 +01:00
parent 95931f6c53
commit 6ff3b90ceb
1 changed files with 32 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

36
cara/model_scenarios_publication.py
View file

@ -2,7 +2,7 @@ from cara import models
from cara.montecarlo import *
######### Standard exposure models ###########
exposure_models = [MCExposureModel(
qR_models = [MCExposureModel(
concentration_model=MCConcentrationModel(
room=models.Room(volume=45),
ventilation=models.SlidingWindow(
@ -25,7 +25,35 @@ exposure_models = [MCExposureModel(
exposed=models.Population(
number=2,
presence=models.SpecificInterval(((0, 4), (5, 9))),
activity=models.Activity.types['Seated'],
mask=models.Mask.types['Type I']
activity=models.Activity.types['Light activity'],
mask=models.Mask.types['No mask']
)
) for a in (1, 2, 3)]
) for a in (1, 2, 3)]
qR_models_shout_light = [MCExposureModel(
concentration_model=MCConcentrationModel(
room=models.Room(volume=45),
ventilation=models.SlidingWindow(
active=models.PeriodicInterval(period=120, duration=120),
inside_temp=models.PiecewiseConstant((0, 24), (293,)),
outside_temp=models.PiecewiseConstant((0, 24), (283,)),
window_height=1.6, opening_length=0.2,
),
infected=MCInfectedPopulation(
number=1,
presence=models.SpecificInterval(((0, 4), (5, 9))),
masked=False,
virus=MCVirus(halflife=1.1, qID=qid),
expiratory_activity=3,
samples=2000000,
breathing_category=3,
expiratory_activity_weights=(0.7, 0.3, 0)
)
),
exposed=models.Population(
number=2,
presence=models.SpecificInterval(((0, 4), (5, 9))),
activity=models.Activity.types['Light activity'],
mask=models.Mask.types['No mask']
)
) for qid in (100, 500, 1000)]