solve conflicts
This commit is contained in:
Andre Henriques
commit
cfe14d7eeb
3 changed files with 760 additions and 620 deletions
|
|
@ -1,179 +1,223 @@
|
|||
from cara.montecarlo import *
|
||||
from cara.model_scenarios_IGH_paper import *
|
||||
|
||||
##Fig 1a
|
||||
#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_windows_open_breaks[0],
|
||||
# #Fig 1c
|
||||
# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_windows_open_breaks[0],
|
||||
# classroom_model_IGH_no_mask_windows_open_alltimes[0]],
|
||||
# labels=['Windows closed', 'Window open during breaks', 'Window open at all times'],
|
||||
# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'],
|
||||
# title='No mask - Spring/Summer period'
|
||||
# title=''
|
||||
# )
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_open_breaks[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_open_alltimes[0])
|
||||
# print("\n############ Fig 1c ###########\n")
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_open_breaks[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_open_alltimes[0])
|
||||
#
|
||||
##Fig 1b
|
||||
#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_2windows_open_breaks[0],
|
||||
# classroom_model_IGH_no_mask_2windows_open_alltimes[0]],
|
||||
# labels=['Windows closed', '2 windows open during breaks', '2 windows open at all times'],
|
||||
# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'],
|
||||
# title='No mask - Spring/Summer period'
|
||||
# )
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_2windows_open_breaks[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_2windows_open_alltimes[0])
|
||||
|
||||
##Fig 1c
|
||||
#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_windows_open_breaks_winter[0],
|
||||
# #Fig 1d
|
||||
# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed_winter[0],classroom_model_IGH_no_mask_windows_open_breaks_winter[0],
|
||||
# classroom_model_IGH_no_mask_windows_open_alltimes_winter[0], classroom_model_IGH_no_mask_windows_fully_open_alltimes_winter[0]],
|
||||
# labels=['Windows closed', 'Window fully open during breaks', 'Window slightly open at all times', 'Window fully open at all times'],
|
||||
# colors=['tomato', 'lightskyblue', '#1f77b4', 'limegreen', 'seagreen', 'lightskyblue', 'deepskyblue'],
|
||||
# title='No mask - Winter period'
|
||||
# title=''
|
||||
# )
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_open_breaks_winter[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_open_alltimes_winter[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_fully_open_alltimes_winter[0])
|
||||
# print("\n############ Fig 1d ###########\n")
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_closed_winter[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_open_breaks_winter[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_open_alltimes_winter[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_fully_open_alltimes_winter[0])
|
||||
#
|
||||
|
||||
##Fig 1d
|
||||
#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_2windows_open_breaks_winter[0],
|
||||
# #Fig 1e
|
||||
# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_2windows_open_breaks[0],
|
||||
# classroom_model_IGH_no_mask_2windows_open_alltimes[0]],
|
||||
# labels=['Windows closed', '2 windows open during breaks', '2 windows open at all times'],
|
||||
# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'],
|
||||
# title=''
|
||||
# )
|
||||
# print("\n############ Fig 1e ###########\n")
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_2windows_open_breaks[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_2windows_open_alltimes[0])
|
||||
#
|
||||
# #Fig 1f
|
||||
# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed_winter[0],classroom_model_IGH_no_mask_2windows_open_breaks_winter[0],
|
||||
# classroom_model_IGH_no_mask_2windows_open_alltimes_winter[0], classroom_model_IGH_no_mask_2windows_fully_open_alltimes_winter[0]],
|
||||
# labels=['Windows closed', '2 windows fully open during breaks', '2 windows slighty open at all times', '2 windows fully open at all times'],
|
||||
# colors=['tomato', 'lightskyblue', '#1f77b4', 'limegreen', 'seagreen', 'lightskyblue', 'deepskyblue'],
|
||||
# title='No mask - Winter period'
|
||||
# title=''
|
||||
# )
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_2windows_open_breaks_winter[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_2windows_open_alltimes_winter[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_2windows_fully_open_alltimes_winter[0])
|
||||
# print("\n############ Fig 1f ###########\n")
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_closed_winter[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_2windows_open_breaks_winter[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_2windows_open_alltimes_winter[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_2windows_fully_open_alltimes_winter[0])
|
||||
#
|
||||
##Fig 1e
|
||||
#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_6windows_open_breaks[0],
|
||||
# #Fig 1g
|
||||
# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_6windows_open_breaks[0],
|
||||
# classroom_model_IGH_no_mask_6windows_open_alltimes[0]],
|
||||
# labels=['Windows closed', '6 windows open during breaks', '6 windows open at all times'],
|
||||
# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'],
|
||||
# title='No mask - Spring/Summer period'
|
||||
# title=''
|
||||
# )
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_6windows_open_alltimes[0])
|
||||
# print("\n############ Fig 1g ###########\n")
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_6windows_open_alltimes[0])
|
||||
#
|
||||
##Fig 1f
|
||||
#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_6windows_open_breaks_winter[0],
|
||||
# #Fig 1h
|
||||
# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed_winter[0],classroom_model_IGH_no_mask_6windows_open_breaks_winter[0],
|
||||
# classroom_model_IGH_no_mask_6windows_open_alltimes_winter[0], classroom_model_IGH_no_mask_6windows_fully_open_alltimes_winter[0]],
|
||||
# labels=['Windows closed', '6 windows fully open during breaks', '6 windows slighty open at all times', '6 windows fully open at all times'],
|
||||
# colors=['tomato', 'lightskyblue', '#1f77b4', 'limegreen', 'seagreen', 'lightskyblue', 'deepskyblue'],
|
||||
# title='No mask - Winter period'
|
||||
# title=''
|
||||
# )
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_winter[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_6windows_open_alltimes_winter[0])
|
||||
#rint_qd_info(classroom_model_IGH_no_mask_6windows_fully_open_alltimes_winter[0])
|
||||
# print("\n############ Fig 1h ###########\n")
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_closed_winter[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_winter[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_6windows_open_alltimes_winter[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_6windows_fully_open_alltimes_winter[0])
|
||||
#
|
||||
##Fig 2
|
||||
#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_windows_closed_1HEPA[0],
|
||||
# #Fig 1i
|
||||
# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],
|
||||
# classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass[0]],
|
||||
# labels=['Windows closed',
|
||||
# '6 windows open during breaks + end of classes'],
|
||||
# colors=['tomato', '#1f77b4'],
|
||||
# title=''
|
||||
# )
|
||||
# print("\n############ Fig 1i ###########\n")
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass[0])
|
||||
#
|
||||
# #Fig 1j
|
||||
# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed_winter[0],
|
||||
# classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass_winter[0]],
|
||||
# labels=['Windows closed',
|
||||
# '6 windows open during breaks + end of classes'],
|
||||
# colors=['tomato', '#1f77b4'],
|
||||
# title=''
|
||||
# )
|
||||
# print("\n############ Fig 1j ###########\n")
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_closed_winter[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass_winter[0])
|
||||
#
|
||||
#
|
||||
#
|
||||
# #Fig 2a
|
||||
# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_windows_closed_1HEPA[0],
|
||||
# classroom_model_IGH_no_mask_windows_closed_2HEPA[0]],
|
||||
# labels=['Windows closed', 'Windows closed + 1 HEPA filter (2.5 ACH)', 'Windows closed + 2 HEPA filter (5 ACH)'],
|
||||
# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'],
|
||||
# title='No mask'
|
||||
# title=''
|
||||
# )
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_closed_1HEPA[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_closed_2HEPA[0])
|
||||
# print("\n############ Fig 2a ###########\n")
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_closed_1HEPA[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_closed_2HEPA[0])
|
||||
#
|
||||
##Fig 3
|
||||
#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_with_mask_windows_closed[0]],
|
||||
# #Fig 2b
|
||||
# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed_winter[0],classroom_model_IGH_no_mask_windows_closed_1HEPA_winter[0],
|
||||
# classroom_model_IGH_no_mask_windows_closed_2HEPA_winter[0]],
|
||||
# labels=['Windows closed', 'Windows closed + 1 HEPA filter (2.5 ACH)', 'Windows closed + 2 HEPA filter (5 ACH)'],
|
||||
# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'],
|
||||
# title=''
|
||||
# )
|
||||
# print("\n############ Fig 2b ###########\n")
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_closed_winter[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_closed_1HEPA_winter[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_closed_2HEPA_winter[0])
|
||||
#
|
||||
# #Fig 2c
|
||||
# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_with_mask_windows_closed[0]],
|
||||
# labels=['Windows closed', 'Windows closed + surgical type masks'],
|
||||
# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'],
|
||||
# title='Use of mask'
|
||||
# title=''
|
||||
# )
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
#print_qd_info(classroom_model_IGH_with_mask_windows_closed[0])
|
||||
# print("\n############ Fig 2c ###########\n")
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
# print_qd_info(classroom_model_IGH_with_mask_windows_closed[0])
|
||||
#
|
||||
##Fig 4a
|
||||
#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_2windows_open_breaks_winter_bis[0],
|
||||
# #Fig 2d
|
||||
# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed_winter[0],classroom_model_IGH_with_mask_windows_closed_winter[0]],
|
||||
# labels=['Windows closed', 'Windows closed + surgical type masks'],
|
||||
# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'],
|
||||
# title=''
|
||||
# )
|
||||
# print("\n############ Fig 2d ###########\n")
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_closed_winter[0])
|
||||
# print_qd_info(classroom_model_IGH_with_mask_windows_closed_winter[0])
|
||||
#
|
||||
#
|
||||
# #Fig 3b
|
||||
# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed_winter[0],classroom_model_IGH_no_mask_2windows_open_alltimes_winter_bis[0],
|
||||
# classroom_model_IGH_with_mask_2windows_open_alltimes_winter[0]],
|
||||
# labels=['Windows closed', '2 windows open at all times', '2 windows open at all times + surgical masks'],
|
||||
# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'],
|
||||
# title=''
|
||||
# )
|
||||
# print("\n############ Fig 3b ###########\n")
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_closed_winter[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_2windows_open_alltimes_winter_bis[0])
|
||||
# print_qd_info(classroom_model_IGH_with_mask_2windows_open_alltimes_winter[0])
|
||||
#
|
||||
# #Fig 3c
|
||||
# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed_winter[0],classroom_model_IGH_with_mask_windows_closed_1HEPA[0],
|
||||
# classroom_model_IGH_with_mask_windows_closed_2HEPA[0]],
|
||||
# labels=['No masks', '1 HEPA filter (2.5 ACH) + surgical masks', '2 HEPA filter (5 ACH) + surgical masks'],
|
||||
# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'],
|
||||
# title=''
|
||||
# )
|
||||
# print("\n############ Fig 3c ###########\n")
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
# print_qd_info(classroom_model_IGH_with_mask_windows_closed_1HEPA[0])
|
||||
# print_qd_info(classroom_model_IGH_with_mask_windows_closed_2HEPA[0])
|
||||
#
|
||||
# #Fig 3d
|
||||
# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed_winter[0],classroom_model_IGH_with_mask_2windows_open_alltimes_winter_1HEPA[0],
|
||||
# classroom_model_IGH_with_mask_2windows_open_alltimes_winter_2HEPA[0]],
|
||||
# labels=['Windows closed', '2 windows open at all times + 1 HEPA filter (2.5 ACH) + surgical masks', '2 windows open at all times + 2 HEPA filter (5 ACH) + surgical masks'],
|
||||
# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'],
|
||||
# title=''
|
||||
# )
|
||||
# print("\n############ Fig 3d ###########\n")
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_closed_winter[0])
|
||||
# print_qd_info(classroom_model_IGH_with_mask_2windows_open_alltimes_winter_1HEPA[0])
|
||||
# print_qd_info(classroom_model_IGH_with_mask_2windows_open_alltimes_winter_2HEPA[0])
|
||||
#
|
||||
#
|
||||
# #Hourly breaks
|
||||
# print("\n############ Hourly breaks - summer ###########\n")
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_open_breaks_endOfClass[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_2windows_open_breaks_endOfClass[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass[0])
|
||||
#
|
||||
# print("\n############ Hourly breaks - winter ###########\n")
|
||||
# print_qd_info(classroom_model_IGH_no_mask_windows_open_breaks_endOfClass_winter[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_2windows_open_breaks_endOfClass_winter[0])
|
||||
# print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass_winter[0])
|
||||
|
||||
|
||||
|
||||
#>>>>>>>>>>> Aux >>>>>>>>
|
||||
|
||||
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_2windows_open_breaks_endOfClass_winter[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_open_breaks_endOfClass_winter[0])
|
||||
|
||||
# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_2windows_open_breaks_winter_bis[0],
|
||||
# classroom_model_IGH_with_mask_2windows_open_breaks_winter[0]],
|
||||
# labels=['Windows closed', '2 windows open during breaks', '2 windows open during breaks + surgical masks'],
|
||||
# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'],
|
||||
# title='Combination of measures (winter period) '
|
||||
# )
|
||||
#
|
||||
##Fig 4b
|
||||
#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_with_mask_2windows_open_breaks_winter_1HEPA[0],
|
||||
|
||||
# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_with_mask_2windows_open_breaks_winter_1HEPA[0],
|
||||
# classroom_model_IGH_with_mask_2windows_open_breaks_winter_2HEPA[0]],
|
||||
# labels=['Windows closed', '2 windows open during breaks + 1 HEPA filter (2.5 ACH) + surgical masks', '2 windows open during breaks + 2 HEPA filter (5 ACH) + surgical masks'],
|
||||
# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'],
|
||||
# title='Combination of measures (winter period) '
|
||||
# )
|
||||
#
|
||||
##Fig 4c
|
||||
#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_2windows_open_alltimes_winter_bis[0],
|
||||
# classroom_model_IGH_with_mask_2windows_open_alltimes_winter[0]],
|
||||
# labels=['Windows closed', '2 windows open at all times', '2 windows open at all times + surgical masks'],
|
||||
# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'],
|
||||
# title='Combination of measures (winter period) '
|
||||
# )
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_2windows_open_alltimes_winter_bis[0])
|
||||
#print_qd_info(classroom_model_IGH_with_mask_2windows_open_alltimes_winter[0])
|
||||
#
|
||||
##Fig 4d
|
||||
#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_with_mask_2windows_open_alltimes_winter_1HEPA[0],
|
||||
# classroom_model_IGH_with_mask_2windows_open_alltimes_winter_2HEPA[0]],
|
||||
# labels=['Windows closed', '2 windows open at all times + 1 HEPA filter (2.5 ACH) + surgical masks', '2 windows open at all times + 2 HEPA filter (5 ACH) + surgical masks'],
|
||||
# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'],
|
||||
# title='Combination of measures (winter period) '
|
||||
# )
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
#print_qd_info(classroom_model_IGH_with_mask_2windows_open_alltimes_winter_1HEPA[0])
|
||||
#print_qd_info(classroom_model_IGH_with_mask_2windows_open_alltimes_winter_2HEPA[0])
|
||||
#
|
||||
|
||||
#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],
|
||||
# classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass[0]],
|
||||
# labels=['Windows closed',
|
||||
# '6 windows open during breaks + end of classes'],
|
||||
# colors=['tomato', '#1f77b4'],
|
||||
# title='Without mask - Spring/Summer period'
|
||||
# )
|
||||
#
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
##print_qd_info(classroom_model_IGH_no_mask_6windows_open_alltimes[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass[0])
|
||||
|
||||
#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],
|
||||
# classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass_winter[0]],
|
||||
# labels=['Windows closed',
|
||||
# '6 windows open during breaks + end of classes'],
|
||||
# colors=['tomato', '#1f77b4'],
|
||||
# title='Without mask - Winter period'
|
||||
# )
|
||||
#
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_6windows_open_alltimes_winter[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass_winter[0])
|
||||
|
||||
#print_qd_info(classroom_model_IGH_no_mask_2windows_open_breaks_endOfClass_winter[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_open_breaks_endOfClass_winter[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_2windows_open_breaks_endOfClass_winter[0])
|
||||
#print_qd_info(classroom_model_IGH_no_mask_windows_open_breaks_endOfClass_winter[0])
|
||||
|
||||
compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_with_mask_windows_closed_1HEPA[0],
|
||||
classroom_model_IGH_with_mask_windows_closed_2HEPA[0]],
|
||||
labels=['No masks', '1 HEPA filter (2.5 ACH) + surgical masks', '2 HEPA filter (5 ACH) + surgical masks'],
|
||||
colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'],
|
||||
title='Combination of measures (windows closed) '
|
||||
)
|
||||
|
||||
print_qd_info(classroom_model_IGH_no_mask_windows_closed[0])
|
||||
print_qd_info(classroom_model_IGH_with_mask_windows_closed_1HEPA[0])
|
||||
print_qd_info(classroom_model_IGH_with_mask_windows_closed_2HEPA[0])
|
||||
|
||||
|
||||
#compare_concentration_curves([classroom_model_IGH_no_mask_windows_closed[1],classroom_model_IGH_no_mask_2windows_open_breaks[1],classroom_model_IGH_no_mask_2windows_open_breaks_endOfClass[1],
|
||||
# classroom_model_IGH_no_mask_2windows_open_alltimes[1], classroom_model_IGH_no_mask_windows_closed_1HEPA[1],
|
||||
|
|
|
|||
File diff suppressed because it is too large
Load diff
|
|
@ -1202,12 +1202,13 @@ def compare_concentration_curves_virus_IGH_paper(exp_models: typing.List[MCExpos
|
|||
times = np.arange(start, stop, TIME_STEP)
|
||||
|
||||
concentrations = [[np.mean(model.concentration_model.concentration(t))*exp_models[0].concentration_model.virus.qID for t in times] for model in exp_models]
|
||||
fig, ax = plt.subplots()
|
||||
fig, ax = plt.subplots(figsize=(8.5,6))
|
||||
for c, label, color in zip(concentrations, labels, colors):
|
||||
ax.plot(times, c, label=label, color=color)
|
||||
|
||||
ax.legend(loc='upper left')
|
||||
ax.set_ylim(ax.get_ylim()[0], ax.get_ylim()[1] * 1.2)
|
||||
#ax.set_ylim(ax.get_ylim()[0], ax.get_ylim()[1] * 1.2)
|
||||
ax.set_ylim(ax.get_ylim()[0], 230)
|
||||
ax.spines["right"].set_visible(False)
|
||||
|
||||
factors = [0.6 * model.exposed.activity.inhalation_rate * (1 - model.exposed.mask.η_inhale) for model in exp_models]
|
||||
|
|
@ -1230,7 +1231,8 @@ def compare_concentration_curves_virus_IGH_paper(exp_models: typing.List[MCExpos
|
|||
ax1.spines["right"].set_linestyle("--")
|
||||
ax1.spines["right"].set_linestyle((0,(1,5)))
|
||||
ax1.set_ylabel('Mean cumulative dose\n(virion)', fontsize=14)
|
||||
ax1.set_ylim(ax1.get_ylim()[0], ax1.get_ylim()[1] * 1.2)
|
||||
#ax1.set_ylim(ax1.get_ylim()[0], ax1.get_ylim()[1] * 1.2)
|
||||
ax1.set_ylim(ax.get_ylim()[0], 300)
|
||||
|
||||
#ax2 = ax.twinx()
|
||||
#ax2.spines["right"].set_position(("axes", 1.15))
|
||||
|
|
|
|||
Loading…
Reference in a new issue