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

move buonanno model into function

Browse source
This commit is contained in:
markus 2021-02-02 13:31:24 +01:00
parent 5f95c410e9
commit ee17e9f937
1 changed files with 18 additions and 17 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

35
cara/montecarlo.py
View file

@ -338,6 +338,24 @@ def print_qr_info(qr_values: np.ndarray) -> None:
print(f"qR_{quantile} = {np.quantile(qr_values, quantile)}")
def buaonanno_exposure_model():
return MCExposureModel(
concentration_model=MCConcentrationModel(
room=models.Room(volume=800),
ventilation=models.AirChange(active=models.PeriodicInterval(period=120, duration=120),
air_exch=0.5),
infected=BuonannoSpecificInfectedPopulation(virus=MCVirus(halflife=1.1),
samples=1)
),
exposed=models.Population(
number=1,
presence=models.SpecificInterval(((0, 2),)),
activity=models.Activity.types['Light activity'],
mask=models.Mask.types['No mask']
)
)
baseline_mc_exposure_model = MCExposureModel(
concentration_model=MCConcentrationModel(
room=models.Room(volume=75),
@ -366,23 +384,6 @@ baseline_mc_exposure_model = MCExposureModel(
)
)
buonanno_exposure_model = MCExposureModel(
concentration_model=MCConcentrationModel(
room=models.Room(volume=800),
ventilation=models.AirChange(active=models.PeriodicInterval(period=120, duration=120),
air_exch=0.5),
infected=BuonannoSpecificInfectedPopulation(virus=MCVirus(halflife=1.1),
samples=1)
),
exposed=models.Population(
number=1,
presence=models.SpecificInterval(((0, 2),)),
activity=models.Activity.types['Light activity'],
mask=models.Mask.types['No mask']
)
)
durations = tuple(np.linspace(0, 120, 10))
models = [