# MAIN for scenario in scenarios: # scenario name, read and modify data for scenario sce = scenario.__name__ data = vicus.read_excel(filename) scenario(data) # create model, solve it, read results model = vicus.create_model(data, timesteps) prob = model.create() optim = SolverFactory('glpk') # cplex, glpk, gurobi, ... result = optim.solve(prob, tee=True) prob.load(result) # write report to spreadsheet vicus.report(prob, os.path.join('results', '{}.xlsx').format(sce), ['ElecAC', 'ElecDC']) # example: customise color of commodity vicus.COLOURS['ElecDC'] = (0, 121, 239) # timeseries plot for each commodity for co in ['ElecAC', 'ElecDC']: fig = vicus.plot(prob, co=co) for ext in ['png', 'pdf']: fig_filename = os.path.join('results', '{}-{}.{}').format(sce, co, ext) fig.savefig(fig_filename, bbox_inches='tight')
# MAIN for scenario in scenarios: # scenario name, read and modify data for scenario sce = scenario.__name__ data = vicus.read_excel(filename) scenario(data) # create model, solve it, read results model = vicus.create_model(data, timesteps) prob = model.create() optim = SolverFactory('glpk') # cplex, glpk, gurobi, ... result = optim.solve(prob, tee=True) prob.load(result) # write report to spreadsheet vicus.report(prob, '{}.xlsx'.format(sce), ['ElecAC', 'ElecDC']) # example: customise color of commodity vicus.COLOURS['ElecDC'] = (0, 121, 239) # timeseries plot for each commodity for co in ['ElecAC', 'ElecDC']: fig = vicus.plot(prob, co=co) for ext in ['png', 'pdf']: fig_filename = '{}-{}.{}'.format(sce, co, ext) fig.savefig(fig_filename, bbox_inches='tight')
# MAIN for scenario in scenarios: # scenario name, read and modify data for scenario sce = scenario.__name__ data = vicus.read_excel(filename) scenario(data) # create model, solve it, read results model = vicus.create_model(data, timesteps) prob = model.create() optim = SolverFactory('glpk') # cplex, glpk, gurobi, ... result = optim.solve(prob, tee=True) prob.load(result) # write report to spreadsheet vicus.report(prob, os.path.join('results', '{}.xlsx').format(sce), ['ElecAC', 'ElecDC']) # example: customise color of commodity vicus.COLOURS['ElecDC'] = (0, 121, 239) # timeseries plot for each commodity for co in ['ElecAC', 'ElecDC']: fig = vicus.plot(prob, co=co) for ext in ['png', 'pdf']: fig_filename = os.path.join( 'results', '{}-{}.{}').format(sce, co, ext) fig.savefig(fig_filename, bbox_inches='tight')