# -*- coding: utf-8 -*-
"""
Created on Thu Mar 24 10:39:28 2016
@author: Jonatan
Compare results for different accumulator sizes
"""
from __future__ import division
from sim_p2h import sim_p2h, power_scen_to_res, temp_to_heat
#scen = [1, 2, 3, 4]
scen = ['Conservative', 'High Wind', 'High Wind & Solar']
ac_sizes = [1, 10]
P_res = {}
cc, dhd = temp_to_heat(exclude_waste=False)
for s in scen:
P_res[s] = power_scen_to_res(s)
for a in ac_sizes:
sim_p2h(P_res[s], cc, dhd, s, ac_size=a, ac_transfer=a)
# -*- coding: utf-8 -*-
"""
Created on Wed Apr 06 09:13:25 2016
@author: Jonatan
"""
from __future__ import division
from sim_p2h import temp_to_heat, power_scen_to_res, sim_p2h
#scen = [1, 2, 3, 4]
scen = ['Conservative', 'High Wind', 'High Wind & Solar']
min_hydro = [0, 1.682, 4]
# first calculate heat curves
county_curves, design_heat_dict = temp_to_heat()
for s in scen:
for h in min_hydro:
P_res = power_scen_to_res(s, hydro_min = h)
sim_p2h(P_res, county_curves, design_heat_dict, s, hydro_min = h)
# -*- coding: utf-8 -*-
"""
Created on Wed Mar 23 14:34:08 2016
@author: Jonatan
"""
from __future__ import division
from sim_p2h import temp_to_heat, power_scen_to_res, sim_p2h
#scen = [0, 1, 2, 3, 4]
#scen = ['Current', 'Conservative', 'High RES', 'High RES & Surplus',
# 'Extreme RES']
scen = ['Current', 'Conservative', 'High Wind',
'High Wind & Solar']
cc, dhd = temp_to_heat()
for s in scen:
P_res = power_scen_to_res(s)
result = sim_p2h(P_res, cc, dhd, s)
my_scenario["Other"] = 15
my_scenario["Wind"] = 60
my_scenario["Solar"] = 25
my_scenario["Consumption"] = 170
# give a name for your scenario
my_scenario["Name"] = "Jonatan's scenario"
# choose scenario to compare to
comp_scen = 3
# read heat curve
heat_curves, design_heat_loads = temp_to_heat()
# calculate power residual
P_res = power_scen_to_res(scenario=my_scenario)
# simulate P2H, results are written to file
sim_p2h(P_res, heat_curves, design_heat_loads, my_scenario)
# store results in a dictionary
results = {}
# read my scenario
results[my_scenario['Name']] = read_results(my_scenario["Name"])
# read other scenario
results['Scenario_{}'.format(comp_scen)] = read_results(comp_scen)
# compare scenarios
plot_compare_scenarios(results, my_scenario['Name'])