def test_run_idf_5vertex():
"""Test the run_idf methods with 5-vertex apertures."""
# Get input Model
input_hb_model = './tests/json/model_colinear_verts.json'
with open(input_hb_model) as json_file:
data = json.load(json_file)
model = Model.from_dict(data)
# Get the input SimulationParameter
sim_par = SimulationParameter()
sim_par.output.add_zone_energy_use()
ddy_file = './tests/ddy/chicago.ddy'
sim_par.sizing_parameter.add_from_ddy_996_004(ddy_file)
sim_par.run_period.end_date = Date(1, 7)
# create the IDF string for simulation paramters and model
idf_str = '\n\n'.join((sim_par.to_idf(), model.to.idf(model)))
# write the final string into an IDF
idf = os.path.join(folders.default_simulation_folder, 'test_file_5verts',
'in.idf')
write_to_file(idf, idf_str, True)
# prepare the IDF for simulation
epw_file = './tests/simulation/chicago.epw'
prepare_idf_for_simulation(idf, epw_file)
# run the IDF through EnergyPlus
sql, zsz, rdd, html, err = run_idf(idf, epw_file)
assert os.path.isfile(sql)
assert os.path.isfile(err)
err_obj = Err(err)
assert 'EnergyPlus Completed Successfully' in err_obj.file_contents
def run_idf_and_report_errors(i):
"""Run an IDF file through EnergyPlus and report errors/warnings on this component."""
# process the additional strings
idf_i = idfs[i]
if add_str_ != [] and add_str_[0] is not None:
a_str = '/n'.join(add_str_)
with open(idf_i, "a") as idf_file:
idf_file.write(a_str)
sql_i, zsz_i, rdd_i, html_i, err_i = run_idf(idf_i,
_epw_file,
silent=silent)
# report any errors on this component
if err_i is not None:
err_obj = Err(err_i)
err_objs[i] = err_obj
for warn in err_obj.severe_errors:
give_warning(ghenv.Component, warn)
for error in err_obj.fatal_errors:
print err_obj.file_contents # print before raising the error
raise Exception(error)
# append everything to the global lists
sql[i] = sql_i
zsz[i] = zsz_i
rdd[i] = rdd_i
html[i] = html_i
err[i] = err_i
def simulate_idf(idf_file, epw_file, folder, log_file):
"""Simulate an IDF file in EnergyPlus.
\b
Args:
idf_file: Full path to a simulate-able .idf file.
epw_file: Full path to an .epw file.
"""
try:
# set the default folder to the default if it's not specified and copy the IDF
if folder is None:
proj_name = os.path.basename(idf_file).replace('.idf', '')
folder = os.path.join(folders.default_simulation_folder, proj_name)
preparedir(folder, remove_content=False)
idf = os.path.join(folder, 'in.idf')
shutil.copy(idf_file, idf)
# run the file through EnergyPlus
gen_files = [idf]
sql, eio, rdd, html, err = run_idf(idf, epw_file)
if err is not None and os.path.isfile(err):
gen_files.extend([sql, eio, rdd, html, err])
err_obj = Err(err)
for error in err_obj.fatal_errors:
log_file.write(
err_obj.file_contents) # log before raising the error
raise Exception(error)
else:
raise Exception('Running EnergyPlus failed.')
log_file.write(json.dumps(gen_files))
except Exception as e:
_logger.exception('IDF simulation failed.\n{}'.format(e))
sys.exit(1)
else:
sys.exit(0)
def test_run_idf():
"""Test the prepare_idf_for_simulation and run_idf methods."""
# Get input Model
room = Room.from_box('TinyHouseZone', 5, 10, 3)
room.properties.energy.program_type = office_program
room.properties.energy.add_default_ideal_air()
model = Model('TinyHouse', [room])
# Get the input SimulationParameter
sim_par = SimulationParameter()
sim_par.output.add_zone_energy_use()
ddy_file = './tests/ddy/chicago.ddy'
sim_par.sizing_parameter.add_from_ddy_996_004(ddy_file)
sim_par.run_period.end_date = Date(1, 7)
# create the IDF string for simulation paramters and model
idf_str = '\n\n'.join((sim_par.to_idf(), model.to.idf(model)))
# write the final string into an IDF
idf = os.path.join(folders.default_simulation_folder, 'test_file',
'in.idf')
write_to_file(idf, idf_str, True)
# prepare the IDF for simulation
epw_file = './tests/simulation/chicago.epw'
prepare_idf_for_simulation(idf, epw_file)
# run the IDF through EnergyPlus
sql, zsz, rdd, html, err = run_idf(idf, epw_file)
assert os.path.isfile(sql)
assert os.path.isfile(err)
err_obj = Err(err)
assert 'EnergyPlus Completed Successfully' in err_obj.file_contents
def run_osw_and_report_errors(i):
"""Run an OSW through OpenStudio CLI."""
osw = _osw[i]
osm_i, idf_i = run_osw(osw, silent=silent)
# process the additional strings
if add_str_ != [] and add_str_[0] is not None and idf is not None:
add_str = '/n'.join(add_str_)
with open(idf, "a") as idf_file:
idf_file.write(add_str)
osm.append(osm_i)
idf.append(idf_i)
# run the IDF through EnergyPlus
if run_:
sql_i, zsz_i, rdd_i, html_i, err_i = run_idf(idf_i, _epw_file, silent=silent)
# report any errors on this component
if err_i is not None:
err_obj = Err(err_i)
err_objs.append(err_obj)
for warn in err_obj.severe_errors:
give_warning(ghenv.Component, warn)
for error in err_obj.fatal_errors:
print err_obj.file_contents # print before raising the error
raise Exception(error)
# append everything to the global lists
sql.append(sql_i)
zsz.append(zsz_i)
rdd.append(rdd_i)
html.append(html_i)
err.append(err_i)
def run_osm_and_report_errors(i):
"""Run an OSW through OpenStudio CLI."""
# create a blank osw for the translation
osw_dict = {'seed_file': _osm[i], 'weather_file': _epw_file}
osw_directory = os.path.dirname(_osm[i])
sch_directory1 = os.path.join(os.path.dirname(osw_directory), 'schedules')
sch_directory2 = os.path.join(osw_directory, 'schedules')
if os.path.isdir(sch_directory1):
osw_dict['file_paths'] = [sch_directory1]
elif os.path.isdir(sch_directory2):
osw_dict['file_paths'] = [sch_directory2]
osw = os.path.join(osw_directory, 'workflow.osw')
with open(osw, 'w') as fp:
json.dump(osw_dict, fp, indent=4)
# get an IDF from the OSM using the OpenStudio CLI
osm_i, idf_i = run_osw(osw, silent=silent)
if idf_i is None:
log_osw = OSW(os.path.join(osw_directory, 'out.osw'))
errors = []
print(log_osw.stdout)
for error, tb in zip(log_osw.errors, log_osw.error_tracebacks):
print(tb)
errors.append(error)
raise Exception('Failed to run OpenStudio CLI:\n{}'.format(
'\n'.join(errors)))
# process the additional strings
if add_str_ != [] and add_str_[0] is not None and idf is not None:
a_str = '/n'.join(add_str_)
with open(idf_i, "a") as idf_file:
idf_file.write(a_str)
osm[i] = osm_i
idf[i] = idf_i
# run the IDF through EnergyPlus
if run_:
sql_i, zsz_i, rdd_i, html_i, err_i = run_idf(idf_i,
_epw_file,
silent=silent)
# report any errors on this component
if err_i is not None:
err_obj = Err(err_i)
err_objs[i] = err_obj
for warn in err_obj.severe_errors:
give_warning(ghenv.Component, warn)
for error in err_obj.fatal_errors:
print err_obj.file_contents # print before raising the error
raise Exception(error)
# append everything to the global lists
sql[i] = sql_i
zsz[i] = zsz_i
rdd[i] = rdd_i
html[i] = html_i
err[i] = err_i
def test_run_idf_window_ventilation():
"""Test the Model.to.idf and run_idf method with a model possessing operable windows."""
room = Room.from_box('TinyHouseZone', 5, 10, 3)
room.properties.energy.add_default_ideal_air()
south_face = room[3]
north_face = room[1]
south_face.apertures_by_ratio(0.4, 0.01)
north_face.apertures_by_ratio(0.4, 0.01)
south_face.apertures[0].is_operable = True
north_face.apertures[0].is_operable = True
heat_setpt = ScheduleRuleset.from_constant_value(
'House Heating', 20, schedule_types.temperature)
cool_setpt = ScheduleRuleset.from_constant_value(
'House Cooling', 28, schedule_types.temperature)
setpoint = Setpoint('House Setpoint', heat_setpt, cool_setpt)
room.properties.energy.setpoint = setpoint
vent_control = VentilationControl(22, 27, 12, 30)
room.properties.energy.window_vent_control = vent_control
ventilation = VentilationOpening(wind_cross_vent=True)
op_aps = room.properties.energy.assign_ventilation_opening(ventilation)
assert len(op_aps) == 2
model = Model('TinyHouse', [room])
# Get the input SimulationParameter
sim_par = SimulationParameter()
sim_par.output.add_zone_energy_use()
ddy_file = './tests/ddy/chicago.ddy'
sim_par.sizing_parameter.add_from_ddy_996_004(ddy_file)
sim_par.run_period.end_date = Date(6, 7)
sim_par.run_period.start_date = Date(6, 1)
# create the IDF string for simulation paramters and model
idf_str = '\n\n'.join((sim_par.to_idf(), model.to.idf(model)))
# write the final string into an IDF
idf = os.path.join(folders.default_simulation_folder,
'test_file_window_vent', 'in.idf')
write_to_file(idf, idf_str, True)
# prepare the IDF for simulation
epw_file = './tests/simulation/chicago.epw'
prepare_idf_for_simulation(idf, epw_file)
# run the IDF through EnergyPlus
sql, zsz, rdd, html, err = run_idf(idf, epw_file)
assert os.path.isfile(sql)
assert os.path.isfile(err)
err_obj = Err(err)
assert 'EnergyPlus Completed Successfully' in err_obj.file_contents
def test_run_idf_window_shade():
"""Test the Model.to.idf and run_idf method with a model that has a window shade."""
# Get input Model
room = Room.from_box('TinyHouseZone', 5, 10, 3)
room.properties.energy.program_type = office_program
room.properties.energy.add_default_ideal_air()
double_clear = WindowConstruction('Double Pane Clear',
[clear_glass, air_gap, clear_glass])
shade_mat = EnergyWindowMaterialShade('Low-e Diffusing Shade', 0.005, 0.15,
0.5, 0.25, 0.5, 0, 0.4, 0.2, 0.1,
0.75, 0.25)
sched = ScheduleRuleset.from_daily_values('NighSched', [
1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1
])
double_ec = WindowConstructionShade('Double Low-E Inside EC', double_clear,
shade_mat, 'Interior',
'OnIfHighSolarOnWindow', 200, sched)
south_face = room[3]
south_face.apertures_by_ratio(0.4, 0.01)
south_face.apertures[0].properties.energy.construction = double_ec
north_face = room[1]
north_face.apertures_by_ratio(0.4, 0.01)
model = Model('TinyHouse', [room])
# Get the input SimulationParameter
sim_par = SimulationParameter()
sim_par.output.add_zone_energy_use()
ddy_file = './tests/ddy/chicago.ddy'
sim_par.sizing_parameter.add_from_ddy_996_004(ddy_file)
sim_par.run_period.end_date = Date(1, 7)
# create the IDF string for simulation paramters and model
idf_str = '\n\n'.join((sim_par.to_idf(), model.to.idf(model)))
# write the final string into an IDF
idf = os.path.join(folders.default_simulation_folder, 'test_file_shd',
'in.idf')
write_to_file(idf, idf_str, True)
# prepare the IDF for simulation
epw_file = './tests/simulation/chicago.epw'
prepare_idf_for_simulation(idf, epw_file)
# run the IDF through EnergyPlus
sql, zsz, rdd, html, err = run_idf(idf, epw_file)
assert os.path.isfile(sql)
assert os.path.isfile(err)
err_obj = Err(err)
assert 'EnergyPlus Completed Successfully' in err_obj.file_contents
def simulate_osm(osm_file, epw_file, folder, log_file):
"""Simulate an OSM file in EnergyPlus.
\b
Args:
osm_file: Full path to a simulate-able .osm file.
epw_file: Full path to an .epw file.
"""
try:
# set the default folder to the default if it's not specified and copy the IDF
if folder is None:
proj_name = os.path.basename(osm_file).replace('.osm', '')
folder = os.path.join(folders.default_simulation_folder, proj_name)
preparedir(folder, remove_content=False)
base_osm = os.path.join(folder, 'in.osm')
shutil.copy(osm_file, base_osm)
# create a blank osw for the translation
osw_dict = {'seed_file': osm_file, 'weather_file': epw_file}
osw = os.path.join(folder, 'workflow.osw')
with open(osw, 'w') as fp:
json.dump(osw_dict, fp, indent=4)
# run the OSW through OpenStudio CLI
osm, idf = run_osw(osw)
# run the file through EnergyPlus
if idf is not None and os.path.isfile(idf):
gen_files = [osw, osm, idf]
sql, eio, rdd, html, err = run_idf(idf, epw_file)
if err is not None and os.path.isfile(err):
gen_files.extend([sql, eio, rdd, html, err])
err_obj = Err(err)
for error in err_obj.fatal_errors:
log_file.write(
err_obj.file_contents) # log before raising the error
raise Exception(error)
else:
raise Exception('Running EnergyPlus failed.')
else:
raise Exception('Running OpenStudio CLI failed.')
log_file.write(json.dumps(gen_files))
except Exception as e:
_logger.exception('OSM simulation failed.\n{}'.format(e))
sys.exit(1)
else:
sys.exit(0)
def test_run_idf_daylight_control():
"""Test Model.to.idf and run_idf with a model possessing daylight controls."""
room = Room.from_box('TinyHouseZone', 5, 10, 3)
room.properties.energy.program_type = office_program
room.properties.energy.add_default_ideal_air()
south_face = room[3]
south_face.apertures_by_ratio(0.4, 0.01)
room.properties.energy.add_daylight_control_to_center(0.8, 500, 0.5)
sens_pt = room.properties.energy.daylighting_control.sensor_position
assert sens_pt.x == pytest.approx(2.5, rel=1e-3)
assert sens_pt.y == pytest.approx(5, rel=1e-3)
assert sens_pt.z == pytest.approx(0.8, rel=1e-3)
room.move(Vector3D(5, 5, 0))
assert room.properties.energy.daylighting_control.sensor_position == \
sens_pt.move(Vector3D(5, 5, 0))
model = Model('TinyHouse', [room])
# Get the input SimulationParameter
sim_par = SimulationParameter()
sim_par.output.add_zone_energy_use()
ddy_file = './tests/ddy/chicago.ddy'
sim_par.sizing_parameter.add_from_ddy_996_004(ddy_file)
sim_par.run_period.end_date = Date(6, 7)
sim_par.run_period.start_date = Date(6, 1)
# create the IDF string for simulation paramters and model
idf_str = '\n\n'.join((sim_par.to_idf(), model.to.idf(model)))
# write the final string into an IDF
idf = os.path.join(folders.default_simulation_folder, 'test_file_daylight',
'in.idf')
write_to_file(idf, idf_str, True)
# prepare the IDF for simulation
epw_file = './tests/simulation/chicago.epw'
prepare_idf_for_simulation(idf, epw_file)
# run the IDF through EnergyPlus
sql, zsz, rdd, html, err = run_idf(idf, epw_file)
assert os.path.isfile(sql)
assert os.path.isfile(err)
err_obj = Err(err)
assert 'EnergyPlus Completed Successfully' in err_obj.file_contents
def run_idf_and_report_errors(i):
"""Run an IDF file through EnergyPlus and report errors/warnings on this component."""
idf = _idf[i]
sql_i, zsz_i, rdd_i, html_i, err_i = run_idf(idf, _epw_file)
# report any errors on this component
if err_i is not None:
err_obj = Err(err_i)
err_objs.append(err_obj)
for warn in err_obj.severe_errors:
give_warning(ghenv.Component, warn)
for error in err_obj.fatal_errors:
raise Exception(error)
# append everything to the global lists
sql.append(sql_i)
zsz.append(zsz_i)
rdd.append(rdd_i)
html.append(html_i)
err.append(err_i)
def test_run_idf_hot_water():
"""Test the Model.to.idf and run_idf method with a model possessing hot water."""
# Get input Model
room = Room.from_box('TinyHouseZone', 5, 10, 3)
room.properties.energy.program_type = office_program
simple_office = ScheduleDay(
'Simple Weekday', [0, 1, 0],
[Time(0, 0), Time(9, 0), Time(17, 0)])
schedule = ScheduleRuleset('Office Water Use', simple_office, None,
schedule_types.fractional)
shw = ServiceHotWater('Office Hot Water', 0.1, schedule)
room.properties.energy.service_hot_water = shw
room.properties.energy.add_default_ideal_air()
model = Model('TinyHouse', [room])
# Get the input SimulationParameter
sim_par = SimulationParameter()
sim_par.output.add_zone_energy_use()
ddy_file = './tests/ddy/chicago.ddy'
sim_par.sizing_parameter.add_from_ddy_996_004(ddy_file)
sim_par.run_period.end_date = Date(1, 7)
# create the IDF string for simulation paramters and model
idf_str = '\n\n'.join((sim_par.to_idf(), model.to.idf(model)))
# write the final string into an IDF
idf = os.path.join(folders.default_simulation_folder, 'test_file_shw',
'in.idf')
write_to_file(idf, idf_str, True)
# prepare the IDF for simulation
epw_file = './tests/simulation/chicago.epw'
prepare_idf_for_simulation(idf, epw_file)
# run the IDF through EnergyPlus
sql, zsz, rdd, html, err = run_idf(idf, epw_file)
assert os.path.isfile(sql)
assert os.path.isfile(err)
err_obj = Err(err)
assert 'EnergyPlus Completed Successfully' in err_obj.file_contents
def find_max_cooling_des_day(des_days, sim_par, base_strs):
"""Find the cooling design day with the highest coincident peak load."""
# create sizing parameters with all of the design days
sim_par_dup = sim_par.duplicate()
sim_par_dup.output.outputs = None
for dy in des_days:
sim_par_dup.sizing_parameter.add_design_day(dy)
# write the IDF and run the sizing calculation
idf_str_init = '\n\n'.join([sim_par_dup.to_idf()] + base_strs)
idf = os.path.join(directory, 'in.idf')
write_to_file_by_name(directory, 'in.idf', idf_str_init, True)
sql, zsz, rdd, html, err = run_idf(idf, silent=True)
# determine the design day with the highest peak using the sizing results
sql_obj = SQLiteResult(sql)
d_day_dict = {d_day.name.upper(): [0, d_day] for d_day in des_days}
peak_cool_dict = {}
for zs in sql_obj.zone_cooling_sizes:
d_day_dict[zs.design_day_name][0] += zs.calculated_design_load
peak_cool_dict[zs.zone_name] = zs.calculated_design_load
day_loads = list(d_day_dict.values())
day_loads.sort(key=lambda y: y[0])
return [day_loads[-1][1]], peak_cool_dict
# create the strings for simulation paramters and model
ver_str = energyplus_idf_version() if energy_folders.energyplus_version \
is not None else energyplus_idf_version(compatibe_ep_version)
sim_par_str = _sim_par_.to_idf()
model_str = _model.to.idf(_model,
schedule_directory=sch_directory,
patch_missing_adjacencies=True)
idf_str = '\n\n'.join([ver_str, sim_par_str, model_str])
# write the final string into an IDF
idf = os.path.join(directory, 'in.idf')
write_to_file_by_name(directory, 'in.idf', idf_str, True)
# run the IDF through EnergyPlus
silent = True if _run == 1 else False
sql, zsz, rdd, html, err = run_idf(idf, _epw_file, silent=silent)
if html is None and err is not None: # something went wrong; parse the errors
err_obj = Err(err)
print(err_obj.file_contents)
for error in err_obj.fatal_errors:
raise Exception(error)
# parse the result sql and get the monthly data collections
if os.name == 'nt': # we are on windows; use IronPython like usual
sql_obj = SQLiteResult(sql)
cool_init = sql_obj.data_collections_by_output_name(cool_out)
heat_init = sql_obj.data_collections_by_output_name(heat_out)
light_init = sql_obj.data_collections_by_output_name(light_out)
elec_equip_init = sql_obj.data_collections_by_output_name(el_equip_out)
gas_equip_init = sql_obj.data_collections_by_output_name(gas_equip_out)
process1_init = sql_obj.data_collections_by_output_name(process1_out)
def simulate_model(model_json, epw_file, sim_par_json, obj_per_model,
multiplier, no_plenum, no_cap, shade_dist, base_osw, folder,
log_file):
"""Simulate a Dragonfly Model JSON file in EnergyPlus.
\b
Args:
model_json: Full path to a Dragonfly Model JSON file.
epw_file: Full path to an .epw file.
"""
try:
# get a ddy variable that might get used later
epw_folder, epw_file_name = os.path.split(epw_file)
ddy_file = os.path.join(epw_folder,
epw_file_name.replace('.epw', '.ddy'))
# set the default folder to the default if it's not specified
if folder is None:
proj_name = \
os.path.basename(model_json).replace('.json', '').replace('.dfjson', '')
folder = os.path.join(folders.default_simulation_folder, proj_name,
'OpenStudio')
preparedir(folder, remove_content=False)
# process the simulation parameters and write new ones if necessary
def ddy_from_epw(epw_file, sim_par):
"""Produce a DDY from an EPW file."""
epw_obj = EPW(epw_file)
des_days = [
epw_obj.approximate_design_day('WinterDesignDay'),
epw_obj.approximate_design_day('SummerDesignDay')
]
sim_par.sizing_parameter.design_days = des_days
def write_sim_par(sim_par):
"""Write simulation parameter object to a JSON."""
sim_par_dict = sim_par.to_dict()
sp_json = os.path.abspath(
os.path.join(folder, 'simulation_parameter.json'))
with open(sp_json, 'w') as fp:
json.dump(sim_par_dict, fp)
return sp_json
if sim_par_json is None: # generate some default simulation parameters
sim_par = SimulationParameter()
sim_par.output.add_zone_energy_use()
sim_par.output.add_hvac_energy_use()
else:
with open(sim_par_json) as json_file:
data = json.load(json_file)
sim_par = SimulationParameter.from_dict(data)
if len(sim_par.sizing_parameter.design_days) == 0 and os.path.isfile(
ddy_file):
try:
sim_par.sizing_parameter.add_from_ddy_996_004(ddy_file)
except AssertionError: # no design days within the DDY file
ddy_from_epw(epw_file, sim_par)
elif len(sim_par.sizing_parameter.design_days) == 0:
ddy_from_epw(epw_file, sim_par)
sim_par_json = write_sim_par(sim_par)
# re-serialize the Dragonfly Model
with open(model_json) as json_file:
data = json.load(json_file)
model = Model.from_dict(data)
model.convert_to_units('Meters')
# convert Dragonfly Model to Honeybee
add_plenum = not no_plenum
cap = not no_cap
hb_models = model.to_honeybee(obj_per_model, shade_dist, multiplier,
add_plenum, cap)
# write out the honeybee JSONs
osms = []
idfs = []
sqls = []
for hb_model in hb_models:
model_dict = hb_model.to_dict(triangulate_sub_faces=True)
directory = os.path.join(folder, hb_model.identifier)
file_path = os.path.join(directory,
'{}.json'.format(hb_model.identifier))
preparedir(directory, remove_content=False) # create the directory
with open(file_path, 'w') as fp:
json.dump(model_dict, fp, indent=4)
# Write the osw file to translate the model to osm
osw = to_openstudio_osw(directory,
file_path,
sim_par_json,
base_osw=base_osw,
epw_file=epw_file)
# run the measure to translate the model JSON to an openstudio measure
if osw is not None and os.path.isfile(osw):
if base_osw is None: # separate the OS CLI run from the E+ run
osm, idf = run_osw(osw)
if idf is not None and os.path.isfile(idf):
sql, eio, rdd, html, err = run_idf(idf, epw_file)
osms.append(osm)
idfs.append(idf)
sqls.append(sql)
if err is None or not os.path.isfile(err):
raise Exception('Running EnergyPlus failed.')
else:
raise Exception('Running OpenStudio CLI failed.')
else: # run the whole simulation with the OpenStudio CLI
osm, idf = run_osw(osw, measures_only=False)
if idf is None or not os.path.isfile(idf):
raise Exception('Running OpenStudio CLI failed.')
sql, eio, rdd, html, err = \
output_energyplus_files(os.path.dirname(idf))
if err is None or not os.path.isfile(err):
raise Exception('Running EnergyPlus failed.')
osms.append(osm)
idfs.append(idf)
sqls.append(sql)
else:
raise Exception('Writing OSW file failed.')
log_file.write(json.dumps({'osm': osms, 'idf': idfs, 'sql': sqls}))
except Exception as e:
_logger.exception('Model translation failed.\n{}'.format(e))
sys.exit(1)
else:
sys.exit(0)
# create the strings for simulation paramters and model
ver_str = energyplus_idf_version() if energy_folders.energyplus_version \
is not None else energyplus_idf_version(compatibe_ep_version)
model_str = _model.to.idf(_model, schedule_directory=sch_directory)
idf_str = '\n\n'.join(
[ver_str, location.to_idf(),
_sim_par_.to_idf(), model_str])
# write the final string into an IDF
idf = os.path.join(directory, 'in.idf')
write_to_file_by_name(directory, 'in.idf', idf_str, True)
# run the IDF through EnergyPlus
silent = True if _run == 1 else False
sql, zsz, rdd, html, err = run_idf(idf, silent=silent)
if html is None and err is not None: # something went wrong; parse the errors
err_obj = Err(err)
print(err_obj.file_contents)
for error in err_obj.fatal_errors:
raise Exception(error)
# parse the result ZSZ and get the timestep data collections
zsz_obj = ZSZ(zsz)
cool_init = zsz_obj.cooling_load_data
heat_init = zsz_obj.heating_load_data
cooling = data_to_load(cool_init, 'Cooling', cool_ap)
heating = data_to_load(heat_init, 'Heating', heat_ap)
# parse the result sql and get the timestep data collections
if os.name == 'nt': # we are on windows; use IronPython like usual
# check the model to be sure there are no orphaned faces, apertures, or doors
assert len(_model.orphaned_faces) == 0, orphaned_warning('Face')
assert len(_model.orphaned_apertures) == 0, orphaned_warning('Aperture')
assert len(_model.orphaned_doors) == 0, orphaned_warning('Door')
# create the strings for simulation paramters and model
ver_str = energyplus_idf_version() if energy_folders.energyplus_version \
is not None else energyplus_idf_version((9, 2, 0))
sim_par_str = _sim_par_.to_idf()
model_str = _model.to.idf(_model, schedule_directory=sch_directory,
solar_distribution=_sim_par_.shadow_calculation.solar_distribution)
idf_str = '\n\n'.join([ver_str, sim_par_str, model_str, add_str])
# delete any existing files in the directory
nukedir(directory)
# write the final string into an IDF.
idf = os.path.join(directory, 'in.idf')
write_to_file_by_name(directory, 'in.idf', idf_str, True)
if run_:
sql, zsz, rdd, html, err = run_idf(idf, _epw_file)
if err is not None:
err_obj = Err(err)
print(err_obj.file_contents)
for warn in err_obj.severe_errors:
give_warning(ghenv.Component, warn)
for error in err_obj.fatal_errors:
raise Exception(error)
def simulate_model(model_json, epw_file, sim_par_json, base_osw, folder,
check_model, log_file):
"""Simulate a Model JSON file in EnergyPlus.
\b
Args:
model_json: Full path to a Model JSON file.
epw_file: Full path to an .epw file.
"""
try:
# get a ddy variable that might get used later
epw_folder, epw_file_name = os.path.split(epw_file)
ddy_file = os.path.join(epw_folder,
epw_file_name.replace('.epw', '.ddy'))
# set the default folder to the default if it's not specified
if folder is None:
proj_name = \
os.path.basename(model_json).replace('.json', '').replace('.hbjson', '')
folder = os.path.join(folders.default_simulation_folder, proj_name,
'OpenStudio')
preparedir(folder, remove_content=False)
# process the simulation parameters and write new ones if necessary
def ddy_from_epw(epw_file, sim_par):
"""Produce a DDY from an EPW file."""
epw_obj = EPW(epw_file)
des_days = [
epw_obj.approximate_design_day('WinterDesignDay'),
epw_obj.approximate_design_day('SummerDesignDay')
]
sim_par.sizing_parameter.design_days = des_days
def write_sim_par(sim_par):
"""Write simulation parameter object to a JSON."""
sim_par_dict = sim_par.to_dict()
sp_json = os.path.abspath(
os.path.join(folder, 'simulation_parameter.json'))
with open(sp_json, 'w') as fp:
json.dump(sim_par_dict, fp)
return sp_json
if sim_par_json is None: # generate some default simulation parameters
sim_par = SimulationParameter()
sim_par.output.add_zone_energy_use()
sim_par.output.add_hvac_energy_use()
else:
with open(sim_par_json) as json_file:
data = json.load(json_file)
sim_par = SimulationParameter.from_dict(data)
if len(sim_par.sizing_parameter.design_days) == 0 and os.path.isfile(
ddy_file):
try:
sim_par.sizing_parameter.add_from_ddy_996_004(ddy_file)
except AssertionError: # no design days within the DDY file
ddy_from_epw(epw_file, sim_par)
elif len(sim_par.sizing_parameter.design_days) == 0:
ddy_from_epw(epw_file, sim_par)
sim_par_json = write_sim_par(sim_par)
# run the Model re-serialization and check if specified
if check_model:
model_json = measure_compatible_model_json(model_json, folder)
# Write the osw file to translate the model to osm
osw = to_openstudio_osw(folder,
model_json,
sim_par_json,
base_osw=base_osw,
epw_file=epw_file)
# run the measure to translate the model JSON to an openstudio measure
if osw is not None and os.path.isfile(osw):
gen_files = [osw]
if base_osw is None: # separate the OS CLI run from the E+ run
osm, idf = run_osw(osw)
# run the resulting idf through EnergyPlus
if idf is not None and os.path.isfile(idf):
gen_files.extend([osm, idf])
sql, eio, rdd, html, err = run_idf(idf, epw_file)
if err is not None and os.path.isfile(err):
gen_files.extend([sql, eio, rdd, html, err])
else:
raise Exception('Running EnergyPlus failed.')
else:
raise Exception('Running OpenStudio CLI failed.')
else: # run the whole simulation with the OpenStudio CLI
osm, idf = run_osw(osw, measures_only=False)
if idf is not None and os.path.isfile(idf):
gen_files.extend([osm, idf])
else:
raise Exception('Running OpenStudio CLI failed.')
sql, eio, rdd, html, err = output_energyplus_files(
os.path.dirname(idf))
if os.path.isfile(err):
gen_files.extend([sql, eio, rdd, html, err])
else:
raise Exception('Running EnergyPlus failed.')
log_file.write(json.dumps(gen_files))
else:
raise Exception('Writing OSW file failed.')
except Exception as e:
_logger.exception('Model simulation failed.\n{}'.format(e))
sys.exit(1)
else:
sys.exit(0)
def simulate_model(model_json, epw_file, sim_par_json, base_osw, folder,
check_model, log_file):
"""Simulate a Model JSON file in EnergyPlus.
\n
Args:
model_json: Full path to a Model JSON file.\n
epw_file: Full path to an .epw file.
"""
try:
# check that the model JSON and the EPW file is there
assert os.path.isfile(model_json), \
'No Model JSON file found at {}.'.format(model_json)
assert os.path.isfile(epw_file), \
'No EPW file found at {}.'.format(epw_file)
# ddy variable that might get used later
epw_folder, epw_file_name = os.path.split(epw_file)
ddy_file = os.path.join(epw_folder,
epw_file_name.replace('.epw', '.ddy'))
# set the default folder to the default if it's not specified
if folder is None:
proj_name = os.path.basename(model_json).replace('.json', '')
folder = os.path.join(folders.default_simulation_folder, proj_name,
'OpenStudio')
preparedir(folder, remove_content=False)
# process the simulation parameters and write new ones if necessary
def write_sim_par(sim_par):
"""Write simulation parameter object to a JSON."""
sim_par_dict = sim_par.to_dict()
sp_json = os.path.abspath(
os.path.join(folder, 'simulation_parameter.json'))
with open(sp_json, 'w') as fp:
json.dump(sim_par_dict, fp)
return sp_json
if sim_par_json is None: # generate some default simulation parameters
sim_par = SimulationParameter()
sim_par.output.add_zone_energy_use()
sim_par.output.add_hvac_energy_use()
if os.path.isfile(ddy_file):
sim_par.sizing_parameter.add_from_ddy_996_004(ddy_file)
else:
raise ValueError(
'No sim-par-json was input and there is no .ddy file next to '
'the .epw.\nAt least one of these two cirtieria must be satisfied '
'for a successful simulation.')
sim_par_json = write_sim_par(sim_par)
else:
assert os.path.isfile(sim_par_json), \
'No simulation parameter file found at {}.'.format(sim_par_json)
with open(sim_par_json) as json_file:
data = json.load(json_file)
sim_par = SimulationParameter.from_dict(data)
if len(sim_par.sizing_parameter.design_days
) == 0 and os.path.isfile(ddy_file):
sim_par.sizing_parameter.add_from_ddy_996_004(ddy_file)
sim_par_json = write_sim_par(sim_par)
elif len(sim_par.sizing_parameter.design_days) == 0:
raise ValueError(
'No design days were found in the input sim-par-json and there is '
'no .ddy file next to the .epw.\nAt least one of these two cirtieria '
'must be satisfied for a successful simulation.')
# run the Model re-serialization and check if specified
if check_model:
model_json = measure_compatible_model_json(model_json, folder)
# Write the osw file to translate the model to osm
osw = to_openstudio_osw(folder,
model_json,
sim_par_json,
base_osw=base_osw,
epw_file=epw_file)
# run the measure to translate the model JSON to an openstudio measure
if osw is not None and os.path.isfile(osw):
gen_files = [osw]
if base_osw is None: # separate the OS CLI run from the E+ run
osm, idf = run_osw(osw)
# run the resulting idf through EnergyPlus
if idf is not None and os.path.isfile(idf):
gen_files.extend([osm, idf])
sql, eio, rdd, html, err = run_idf(idf, epw_file)
if err is not None and os.path.isfile(err):
gen_files.extend([sql, eio, rdd, html, err])
else:
raise Exception('Running EnergyPlus failed.')
else:
raise Exception('Running OpenStudio CLI failed.')
else: # run the whole simulation with the OpenStudio CLI
osm, idf = run_osw(osw, measures_only=False)
if idf is not None and os.path.isfile(idf):
gen_files.extend([osm, idf])
else:
raise Exception('Running OpenStudio CLI failed.')
sql, eio, rdd, html, err = output_energyplus_files(
os.path.dirname(idf))
if os.path.isfile(err):
gen_files.extend([sql, eio, rdd, html, err])
else:
raise Exception('Running EnergyPlus failed.')
log_file.write(json.dumps(gen_files))
else:
raise Exception('Writing OSW file failed.')
except Exception as e:
_logger.exception('Model simulation failed.\n{}'.format(e))
sys.exit(1)
else:
sys.exit(0)
