def test_import_stat():
"""Test import standard stat."""
relative_path = './tests/fixtures/stat/chicago.stat'
abs_path = os.path.abspath(relative_path)
stat_rel = STAT(relative_path)
stat = STAT(abs_path)
# Test imports don't break
assert stat.file_path == abs_path
assert stat_rel.file_path == os.path.normpath(relative_path)
assert isinstance(stat, STAT)
def test_import_STAT(self):
"""Test import standard stat."""
relative_path = './tests/stat/chicago.stat'
abs_path = os.path.abspath(relative_path)
stat_rel = STAT(relative_path)
stat = STAT(abs_path)
# Test imports don't break
assert stat.file_path == abs_path
assert stat_rel.file_path == os.path.normpath(relative_path)
assert hasattr(stat, 'isStat')
def test_annual_cooling_design_days(self):
relative_path = './tests/stat/tokyo.stat'
stat = STAT(relative_path)
ann_cdd_04 = stat.annual_cooling_design_day_004
ann_cdd_10 = stat.annual_cooling_design_day_010
assert ann_cdd_04.day_type == 'SummerDesignDay'
assert ann_cdd_04.name == '0.4% Cooling Design Day for TOKYO HYAKURI'
assert ann_cdd_04.analysis_period.st_month == \
ann_cdd_04.analysis_period.end_month == \
ann_cdd_04.sky_condition.month == 8
assert ann_cdd_04.analysis_period.st_day == \
ann_cdd_04.analysis_period.end_day == 21
assert ann_cdd_04.humidity_condition.hum_type == 'Wetbulb'
assert ann_cdd_04.dry_bulb_condition.dry_bulb_max == 32.1
assert ann_cdd_04.humidity_condition.hum_value == 26.0
assert ann_cdd_04.dry_bulb_condition.dry_bulb_range == 7.7
assert ann_cdd_04.wind_condition.wind_speed == 4.8
assert ann_cdd_04.wind_condition.wind_direction == 210
assert ann_cdd_04.sky_condition.solar_model == 'ASHRAETau'
assert ann_cdd_04.sky_condition.tau_b == 0.546
assert ann_cdd_04.sky_condition.tau_d == 1.827
assert ann_cdd_10.name == '1.0% Cooling Design Day for TOKYO HYAKURI'
assert ann_cdd_10.dry_bulb_condition.dry_bulb_max == 30.9
assert ann_cdd_10.humidity_condition.hum_value == 25.8
def test_annual_heating_design_days(self):
relative_path = './tests/stat/chicago.stat'
stat = STAT(relative_path)
ann_hdd_96 = stat.annual_heating_design_day_996
ann_hdd_90 = stat.annual_heating_design_day_990
assert ann_hdd_96.day_type == 'WinterDesignDay'
assert ann_hdd_96.name == '99.6% Heating Design Day for Chicago Ohare Intl Ap'
assert ann_hdd_96.analysis_period.st_month == \
ann_hdd_96.analysis_period.end_month == \
ann_hdd_96.sky_condition.month == 1
assert ann_hdd_96.analysis_period.st_day == \
ann_hdd_96.analysis_period.end_day == 21
assert ann_hdd_96.humidity_condition.hum_type == 'Wetbulb'
assert ann_hdd_96.dry_bulb_condition.dry_bulb_max == \
ann_hdd_96.humidity_condition.hum_value == -20
assert ann_hdd_96.dry_bulb_condition.dry_bulb_range == 0
assert ann_hdd_96.wind_condition.wind_speed == 4.9
assert ann_hdd_96.wind_condition.wind_direction == 270
assert ann_hdd_96.sky_condition.solar_model == 'ASHRAEClearSky'
assert ann_hdd_96.sky_condition.clearness == 0
assert ann_hdd_90.name == '99.0% Heating Design Day for Chicago Ohare Intl Ap'
assert ann_hdd_90.dry_bulb_condition.dry_bulb_max == \
ann_hdd_90.humidity_condition.hum_value == -16.6
def test_annual_heating_design_days():
"""Test the annual heating design days within the stat object."""
relative_path = './tests/fixtures/stat/chicago.stat'
stat = STAT(relative_path)
ann_hdd_96 = stat.annual_heating_design_day_996
ann_hdd_90 = stat.annual_heating_design_day_990
assert ann_hdd_96.day_type == 'WinterDesignDay'
assert ann_hdd_96.name == '99.6% Heating Design Day for Chicago Ohare Intl Ap'
assert ann_hdd_96.analysis_period.st_month == \
ann_hdd_96.analysis_period.end_month == \
ann_hdd_96.sky_condition.date.month == 1
assert ann_hdd_96.analysis_period.st_day == \
ann_hdd_96.analysis_period.end_day == 21
assert ann_hdd_96.humidity_condition.humidity_type == 'Wetbulb'
assert ann_hdd_96.dry_bulb_condition.dry_bulb_max == \
ann_hdd_96.humidity_condition.humidity_value == -20
assert ann_hdd_96.dry_bulb_condition.dry_bulb_range == 0
assert ann_hdd_96.wind_condition.wind_speed == 4.9
assert ann_hdd_96.wind_condition.wind_direction == 270
assert isinstance(ann_hdd_96.sky_condition, ASHRAEClearSky)
assert ann_hdd_96.sky_condition.clearness == 0
assert ann_hdd_90.name == '99.0% Heating Design Day for Chicago Ohare Intl Ap'
assert ann_hdd_90.dry_bulb_condition.dry_bulb_max == \
ann_hdd_90.humidity_condition.humidity_value == -16.6
def test_annual_cooling_design_days():
"""Test the annual cooling design days within the stat object."""
relative_path = './tests/fixtures/stat/tokyo.stat'
stat = STAT(relative_path)
ann_cdd_04 = stat.annual_cooling_design_day_004
ann_cdd_10 = stat.annual_cooling_design_day_010
assert ann_cdd_04.day_type == 'SummerDesignDay'
assert ann_cdd_04.name == '0.4% Cooling Design Day for TOKYO HYAKURI'
assert ann_cdd_04.analysis_period.st_month == \
ann_cdd_04.analysis_period.end_month == \
ann_cdd_04.sky_condition.date.month == 8
assert ann_cdd_04.analysis_period.st_day == \
ann_cdd_04.analysis_period.end_day == 21
assert ann_cdd_04.humidity_condition.humidity_type == 'Wetbulb'
assert ann_cdd_04.dry_bulb_condition.dry_bulb_max == 32.1
assert ann_cdd_04.humidity_condition.humidity_value == 26.0
assert ann_cdd_04.dry_bulb_condition.dry_bulb_range == 7.7
assert ann_cdd_04.wind_condition.wind_speed == 4.8
assert ann_cdd_04.wind_condition.wind_direction == 210
assert isinstance(ann_cdd_04.sky_condition, ASHRAETau)
assert ann_cdd_04.sky_condition.tau_b == 0.546
assert ann_cdd_04.sky_condition.tau_d == 1.827
assert ann_cdd_10.name == '1.0% Cooling Design Day for TOKYO HYAKURI'
assert ann_cdd_10.dry_bulb_condition.dry_bulb_max == 30.9
assert ann_cdd_10.humidity_condition.humidity_value == 25.8
def test_dict_methods():
"""Test dictionary serialization methods."""
relative_path = './tests/fixtures/stat/chicago.stat'
stat_obj = STAT(relative_path)
stat_dict = stat_obj.to_dict()
rebuilt_stat = STAT.from_dict(stat_dict)
assert stat_dict == rebuilt_stat.to_dict()
def test_json_methods(self):
"""Test JSON serialization methods"""
relative_path = './tests/stat/chicago.stat'
stat_obj = STAT(relative_path)
stat_json = stat_obj.to_json()
rebuilt_stat = STAT.from_json(stat_json)
assert stat_json == rebuilt_stat.to_json()
def test_monthly_cooling_design_days(self):
relative_path = './tests/stat/chicago.stat'
stat = STAT(relative_path)
m_ddy_050 = stat.monthly_cooling_design_days_050
m_ddy_100 = stat.monthly_cooling_design_days_100
m_ddy_020 = stat.monthly_cooling_design_days_020
m_ddy_004 = stat.monthly_cooling_design_days_004
assert len(m_ddy_050) == len(m_ddy_100) == len(m_ddy_020) == \
len(m_ddy_004) == 12
ddy_path = './tests/ddy/chicago_monthly.ddy'
monthly_ddy = DDY(stat.location, m_ddy_050)
monthly_ddy.save(ddy_path)
def test_monthly_cooling_design_days():
"""Test the monthly cooling design days within the stat object."""
relative_path = './tests/fixtures/stat/chicago.stat'
stat = STAT(relative_path)
m_ddy_050 = stat.monthly_cooling_design_days_050
m_ddy_100 = stat.monthly_cooling_design_days_100
m_ddy_020 = stat.monthly_cooling_design_days_020
m_ddy_004 = stat.monthly_cooling_design_days_004
assert len(m_ddy_050) == len(m_ddy_100) == len(m_ddy_020) == \
len(m_ddy_004) == 12
ddy_path = './tests/fixtures/ddy/chicago_monthly.ddy'
monthly_ddy = DDY(stat.location, m_ddy_050)
monthly_ddy.save(ddy_path)
def test_typical_extreme_weeks(self):
relative_path = './tests/stat/chicago.stat'
stat = STAT(relative_path)
extreme_cold = stat.extreme_cold_week
extreme_hot = stat.extreme_hot_week
typical_winter = stat.typical_winter_week
typical_spring = stat.typical_spring_week
typical_summer = stat.typical_summer_week
typical_autumn = stat.typical_autumn_week
assert str(extreme_cold) == '1/27 to 2/2 between 0 to 23 @1'
assert str(extreme_hot) == '7/13 to 7/19 between 0 to 23 @1'
assert str(typical_winter) == '12/22 to 12/28 between 0 to 23 @1'
assert str(typical_spring) == '4/26 to 5/2 between 0 to 23 @1'
assert str(typical_summer) == '8/24 to 8/30 between 0 to 23 @1'
assert str(typical_autumn) == '10/27 to 11/2 between 0 to 23 @1'
def test_typical_extreme_weeks():
"""Test the typical and extreme weeks within the stat object."""
relative_path = './tests/fixtures/stat/chicago.stat'
stat = STAT(relative_path)
extreme_cold = stat.extreme_cold_week
extreme_hot = stat.extreme_hot_week
typical_winter = stat.typical_winter_week
typical_spring = stat.typical_spring_week
typical_summer = stat.typical_summer_week
typical_autumn = stat.typical_autumn_week
assert str(extreme_cold) == '1/27 to 2/2 between 0 and 23 @1'
assert str(extreme_hot) == '7/13 to 7/19 between 0 and 23 @1'
assert str(typical_winter) == '12/22 to 12/28 between 0 and 23 @1'
assert str(typical_spring) == '4/26 to 5/2 between 0 and 23 @1'
assert str(typical_summer) == '8/24 to 8/30 between 0 and 23 @1'
assert str(typical_autumn) == '10/27 to 11/2 between 0 and 23 @1'
def test_stat_location(self):
relative_path = './tests/stat/tokyo.stat'
stat = STAT(relative_path)
# Test accuracy of import
assert len(stat.header) == 10
assert stat.location.city == 'TOKYO HYAKURI'
assert stat.location.country == 'JPN'
assert stat.location.source == 'IWEC Data'
assert stat.location.station_id == '477150'
assert stat.location.latitude == approx(36.18, rel=1e-1)
assert stat.location.longitude == approx(140.42, rel=1e-1)
assert stat.location.time_zone == 9
assert stat.location.elevation == 35
assert stat.monthly_tau_beam == \
[0.341, 0.392, 0.48, 0.601, 0.654, 0.632, 0.553, 0.546, 0.473,
0.423, 0.388, 0.338]
assert stat.monthly_tau_diffuse == \
[2.214, 1.999, 1.765, 1.546, 1.501, 1.591, 1.811, 1.827, 2.024,
2.083, 2.105, 2.267]
assert stat.ashrae_climate_zone == '4A'
assert stat.koppen_climate_zone == 'Cfa'
try:
from ladybug.stat import STAT
from ladybug.wea import Wea
except ImportError as e:
raise ImportError('\nFailed to import ladybug:\n\t{}'.format(e))
try:
from ladybug_rhino.grasshopper import all_required_inputs
except ImportError as e:
raise ImportError('\nFailed to import ladybug_rhino:\n\t{}'.format(e))
if all_required_inputs(ghenv.Component):
stat_obj = STAT(_stat_file)
# output location and climate zone
location = stat_obj.location
ashrae_zone = stat_obj.ashrae_climate_zone
koppen_zone = stat_obj.koppen_climate_zone
# output clear sky radiation
try: # first see if we can get the values from monthly optical depths
wea = Wea.from_stat_file(_stat_file)
except: # no optical data was found; use the original clear sky
wea = Wea.from_ashrae_clear_sky(location)
clear_dir_norm_rad = wea.direct_normal_irradiance
clear_diff_horiz_rad = wea.diffuse_horizontal_irradiance
# output design day objects
except ImportError as e:
raise ImportError('\nFailed to import ladybug_rhino:\n\t{}'.format(e))
if all_required_inputs(ghenv.Component) and _write:
# set default values
_percentile_ = 0.4 if _percentile_ is None else _percentile_
_folder_ = os.path.join(folders.default_epw_folder, 'ddy') if _folder_ \
is None else _folder_
if _weather_file.lower().endswith('.epw'):
f_name = os.path.basename(_weather_file.lower()).replace(
'.epw', '.ddy')
epw_data = True
elif _weather_file.lower().endswith('.stat'):
f_name = os.path.basename(_weather_file.lower()).replace(
'.stat', '.ddy')
epw_data = False
else:
raise ValueError('Failed to recognize the file type of "{}".\n'
'Must end in .epw or .stat.'.format(_weather_file))
f_path = os.path.join(_folder_, f_name)
# create the DDY file
if epw_data:
epw = EPW(_weather_file)
ddy_file = epw.to_ddy_monthly_cooling(f_path, _percentile_, 2) \
if monthly_cool_ else epw.to_ddy(f_path, _percentile_)
else:
stat = STAT(_weather_file)
ddy_file = stat.to_ddy_monthly_cooling(f_path, _percentile_, 2) \
if monthly_cool_ else stat.to_ddy(f_path, _percentile_)
