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 }, 'aw': { 'dens': 6.964, 'mu': 52.5, 'sd': 27.3, 'min': 11.1, 'max': 112.5 },