logger = logging.getLogger()
info = logging.StreamHandler()
info.setLevel(logging.INFO)
logger.addHandler(info)
logging.info("START: %s", time.asctime())
###############################################################################
# initial conditions
initial_populations = []
# instantiate an SW kernel to grab an 'x' array...
sw = ABBSW()
sw.L, sw.U = L, U
sw.setup(MidPoint(), N)
# dens: # / 10 m^2
scenarios = {
'C': {
'sw': { 'dens': 8.846, 'mu': 25.7, 'sd': 17.1, 'min': 3.8, 'max': 53.5 },
'aw': { 'dens': 3.308, 'mu': 30.5, 'sd': 22.9, 'min': 9.4, 'max': 76.1 },
},
'CD': {
'sw': { 'dens': 3.778, 'mu': 25.7, 'sd': 17.1, 'min': 3.8, 'max': 53.5 },
'aw': { 'dens': 7.278, 'mu': 63.8, 'sd': 28.5, 'min': 23.7, 'max': 109.2 },
},
'D': {
'sw': { 'dens': 0.036, 'mu': 25.7, 'sd': 17.1, 'min': 3.8, 'max': 53.5 },
