class Standard # Helper method to set the weather file, import the design days, set # water mains temperature, and set ground temperature. # Based on ChangeBuildingLocation measure by Nicholas Long def model_add_design_days_and_weather_file(model, climate_zone, epw_file) success = true require_relative 'Weather.stat_file' # Remove any existing Design Day objects that are in the file model.getDesignDays.each(&:remove) OpenStudio.logFree(OpenStudio::Info, 'openstudio.weather.Model', "Started adding weather file for climate zone: #{climate_zone}.") # Define the weather file for each climate zone climate_zone_weather_file_map = { 'ASHRAE 169-2006-1A' => 'USA_FL_Miami.Intl.AP.722020_TMY3.epw', 'ASHRAE 169-2006-1B' => 'SAU_Riyadh.404380_IWEC.epw', 'ASHRAE 169-2006-2A' => 'USA_TX_Houston-Bush.Intercontinental.AP.722430_TMY3.epw', 'ASHRAE 169-2006-2B' => 'USA_AZ_Phoenix-Sky.Harbor.Intl.AP.722780_TMY3.epw', 'ASHRAE 169-2006-3A' => 'USA_TN_Memphis.Intl.AP.723340_TMY3.epw', 'ASHRAE 169-2006-3B' => 'USA_TX_El.Paso.Intl.AP.722700_TMY3.epw', 'ASHRAE 169-2006-3C' => 'USA_CA_San.Francisco.Intl.AP.724940_TMY3.epw', 'ASHRAE 169-2006-4A' => 'USA_MD_Baltimore-Washington.Intl.AP.724060_TMY3.epw', 'ASHRAE 169-2006-4B' => 'USA_NM_Albuquerque.Intl.AP.723650_TMY3.epw', 'ASHRAE 169-2006-4C' => 'USA_OR_Salem-McNary.Field.726940_TMY3.epw', 'ASHRAE 169-2006-5A' => 'USA_IL_Chicago-OHare.Intl.AP.725300_TMY3.epw', 'ASHRAE 169-2006-5B' => 'USA_ID_Boise.Air.Terminal.726810_TMY3.epw', 'ASHRAE 169-2006-5C' => 'CAN_BC_Vancouver.718920_CWEC.epw', 'ASHRAE 169-2006-6A' => 'USA_VT_Burlington.Intl.AP.726170_TMY3.epw', 'ASHRAE 169-2006-6B' => 'USA_MT_Helena.Rgnl.AP.727720_TMY3.epw', 'ASHRAE 169-2006-7A' => 'USA_MN_Duluth.Intl.AP.727450_TMY3.epw', 'ASHRAE 169-2006-7B' => 'USA_MN_Duluth.Intl.AP.727450_TMY3.epw', 'ASHRAE 169-2006-8A' => 'USA_AK_Fairbanks.Intl.AP.702610_TMY3.epw', 'ASHRAE 169-2006-8B' => 'USA_AK_Fairbanks.Intl.AP.702610_TMY3.epw', # For measure input 'NECB HDD Method' => epw_file.to_s, # For testing 'NECB-CNEB-5' => epw_file.to_s, 'NECB-CNEB-6' => epw_file.to_s, 'NECB-CNEB-7a' => epw_file.to_s, 'NECB-CNEB-7b' => epw_file.to_s, 'NECB-CNEB-8' => epw_file.to_s, # For DEER 'CEC T24-CEC1' => 'ARCATA_725945_CZ2010.epw', 'CEC T24-CEC2' => 'SANTA-ROSA_724957_CZ2010.epw', 'CEC T24-CEC3' => 'OAKLAND_724930_CZ2010.epw', 'CEC T24-CEC4' => 'SAN-JOSE-REID_724946_CZ2010.epw', 'CEC T24-CEC5' => 'SANTA-MARIA_723940_CZ2010.epw', 'CEC T24-CEC6' => 'TORRANCE_722955_CZ2010.epw', 'CEC T24-CEC7' => 'SAN-DIEGO-LINDBERGH_722900_CZ2010.epw', 'CEC T24-CEC8' => 'FULLERTON_722976_CZ2010.epw', 'CEC T24-CEC9' => 'BURBANK-GLENDALE_722880_CZ2010.epw', 'CEC T24-CEC10' => 'RIVERSIDE_722869_CZ2010.epw', 'CEC T24-CEC11' => 'RED-BLUFF_725910_CZ2010.epw', 'CEC T24-CEC12' => 'SACRAMENTO-EXECUTIVE_724830_CZ2010.epw', 'CEC T24-CEC13' => 'FRESNO_723890_CZ2010.epw', 'CEC T24-CEC14' => 'PALMDALE_723820_CZ2010.epw', 'CEC T24-CEC15' => 'PALM-SPRINGS-INTL_722868_CZ2010.epw', 'CEC T24-CEC16' => 'BLUE-CANYON_725845_CZ2010.epw' } # Get the weather file name from the hash weather_file_name = if epw_file.nil? || (epw_file.to_s.strip == '') climate_zone_weather_file_map[climate_zone] else epw_file.to_s end if weather_file_name.nil? OpenStudio.logFree(OpenStudio::Error, 'openstudio.weather.Model', "Could not determine the weather file for climate zone: #{climate_zone}.") success = false end # Define where the weather files lives weather_dir = nil if File.dirname(__FILE__)[0] == ':' # running embedded copy of the gem # load weather file from embedded files epw_string = load_resource_relative("../../../data/weather/#{weather_file_name}") ddy_string = load_resource_relative("../../../data/weather/#{weather_file_name.gsub('.epw', '.ddy')}") stat_string = load_resource_relative("../../../data/weather/#{weather_file_name.gsub('.epw', '.stat')}") # extract to local weather dir weather_dir = File.expand_path(File.join(Dir.pwd, 'extracted_files/weather/')) puts "Extracting weather files to #{weather_dir}" FileUtils.mkdir_p(weather_dir) File.open("#{weather_dir}/#{weather_file_name}", 'wb') { |f| f << epw_string; f.flush } File.open("#{weather_dir}/#{weather_file_name.gsub('.epw', '.ddy')}", 'wb') { |f| f << ddy_string; f.flush } File.open("#{weather_dir}/#{weather_file_name.gsub('.epw', '.stat')}", 'wb') { |f| f << stat_string; f.flush } else # loaded gem from system path top_dir = File.expand_path('../../..', File.dirname(__FILE__)) weather_dir = File.expand_path("#{top_dir}/data/weather") end # Add Weather File unless (Pathname.new weather_dir).absolute? weather_dir = File.expand_path(File.join(File.dirname(__FILE__), weather_dir)) end weather_file = File.join(weather_dir, weather_file_name) epw_file = OpenStudio::EpwFile.new(weather_file) OpenStudio::Model::WeatherFile.setWeatherFile(model, epw_file).get weather_name = "#{epw_file.city}_#{epw_file.stateProvinceRegion}_#{epw_file.country}" weather_lat = epw_file.latitude weather_lon = epw_file.longitude weather_time = epw_file.timeZone weather_elev = epw_file.elevation # Add or update site data site = model.getSite site.setName(weather_name) site.setLatitude(weather_lat) site.setLongitude(weather_lon) site.setTimeZone(weather_time) site.setElevation(weather_elev) # Add SiteWaterMainsTemperature -- via parsing of STAT file. stat_filename = "#{File.join(File.dirname(weather_file), File.basename(weather_file, '.*'))}.stat" if File.exist? stat_filename stat_file = EnergyPlus::StatFile.new(stat_filename) water_temp = model.getSiteWaterMainsTemperature water_temp.setAnnualAverageOutdoorAirTemperature(stat_file.mean_dry_bulb) water_temp.setMaximumDifferenceInMonthlyAverageOutdoorAirTemperatures(stat_file.delta_dry_bulb) # OpenStudio::logFree(OpenStudio::Info, "openstudio.weather.Model", "Mean dry bulb is #{stat_file.mean_dry_bulb}") # OpenStudio::logFree(OpenStudio::Info, "openstudio.weather.Model", "Delta dry bulb is #{stat_file.delta_dry_bulb}") else OpenStudio.logFree(OpenStudio::Warn, 'openstudio.weather.Model', "Could not find .stat file for weather, will use default water mains temperatures which may be inaccurate for the location.") success = false end # Load in the ddy file based on convention that it is in # the same directory and has the same basename as the epw file. ddy_file = "#{File.join(File.dirname(weather_file), File.basename(weather_file, '.*'))}.ddy" if File.exist? ddy_file ddy_model = OpenStudio::EnergyPlus.loadAndTranslateIdf(ddy_file).get ddy_model.getObjectsByType('OS:SizingPeriod:DesignDay'.to_IddObjectType).each do |d| # Import the 99.6% Heating and 0.4% Cooling design days ddy_list = /(Htg 99.6. Condns DB)|(Clg .4% Condns DB=>MWB)/ if d.name.get =~ ddy_list model.addObject(d.clone) # OpenStudio::logFree(OpenStudio::Info, 'openstudio.weather.Model', "Added #{d.name} design day.") end end # Check to ensure that some design days were added if model.getDesignDays.size.zero? OpenStudio.logFree(OpenStudio::Error, 'openstudio.weather.Model', "No design days were loaded, check syntax of .ddy file: #{ddy_file}.") end else OpenStudio.logFree(OpenStudio::Error, 'openstudio.weather.Model', "Could not find .ddy file for: #{ddy_file}.") puts "Could not find .ddy file for: #{ddy_file}." success = false end return success end def model_add_ground_temperatures(model, building_type, climate_zone) ground_temp_vals = model_find_object(standards_data['ground_temperatures'], 'template' => template, 'climate_zone' => climate_zone, 'building_type' => building_type) if ground_temp_vals && ground_temp_vals['jan'] ground_temp = model.getSiteGroundTemperatureBuildingSurface ground_temp.setJanuaryGroundTemperature(ground_temp_vals['jan']) ground_temp.setFebruaryGroundTemperature(ground_temp_vals['feb']) ground_temp.setMarchGroundTemperature(ground_temp_vals['mar']) ground_temp.setAprilGroundTemperature(ground_temp_vals['apr']) ground_temp.setMayGroundTemperature(ground_temp_vals['may']) ground_temp.setJuneGroundTemperature(ground_temp_vals['jun']) ground_temp.setJulyGroundTemperature(ground_temp_vals['jul']) ground_temp.setAugustGroundTemperature(ground_temp_vals['aug']) ground_temp.setSeptemberGroundTemperature(ground_temp_vals['sep']) ground_temp.setOctoberGroundTemperature(ground_temp_vals['oct']) ground_temp.setNovemberGroundTemperature(ground_temp_vals['nov']) ground_temp.setDecemberGroundTemperature(ground_temp_vals['dec']) else OpenStudio.logFree(OpenStudio::Warn, 'openstudio.weather.Model', 'Could not find ground temperatures; will use generic temperatures, which will skew results.') ground_temp = model.getSiteGroundTemperatureBuildingSurface ground_temp.setJanuaryGroundTemperature(19.527) ground_temp.setFebruaryGroundTemperature(19.502) ground_temp.setMarchGroundTemperature(19.536) ground_temp.setAprilGroundTemperature(19.598) ground_temp.setMayGroundTemperature(20.002) ground_temp.setJuneGroundTemperature(21.640) ground_temp.setJulyGroundTemperature(22.225) ground_temp.setAugustGroundTemperature(22.375) ground_temp.setSeptemberGroundTemperature(21.449) ground_temp.setOctoberGroundTemperature(20.121) ground_temp.setNovemberGroundTemperature(19.802) ground_temp.setDecemberGroundTemperature(19.633) end end # Gets the maximum OA dry bulb temperatures # for all WinterDesignDays in the model. # # @return [Array] an array of OA temperatures in C def heating_design_outdoor_temperatures heating_design_outdoor_temps = [] getDesignDays.each do |dd| next unless dd.dayType == 'WinterDesignDay' heating_design_outdoor_temps << dd.maximumDryBulbTemperature end return heating_design_outdoor_temps end end # ********************************************************************* # * Copyright (c) 2008-2015, Natural Resources Canada # * All rights reserved. # * # * This library is free software; you can redistribute it and/or # * modify it under the terms of the GNU Lesser General Public # * License as published by the Free Software Foundation; either # * version 2.1 of the License, or (at your option) any later version. # * # * This library is distributed in the hope that it will be useful, # * but WITHOUT ANY WARRANTY; without even the implied warranty of # * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU # * Lesser General Public License for more details. # * # * You should have received a copy of the GNU Lesser General Public # * License along with this library; if not, write to the Free Software # * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA # **********************************************************************/ # This module has been created to make it easier to manipulate weather files can contains region specific data. module BTAP module Environment require_relative 'Weather.stat_file' # rubocop:disable Style/MutableConstant # rubocop:enable Style/MutableConstant # this method is used to populate user interfaces if needed from the hash above. def self.get_canadian_weather_file_names() canadian_file_names = [] Dir.glob("#{File.dirname(__FILE__)}/../../../**/*.epw").each do |file| canadian_file_names << File.basename(file).to_s end return canadian_file_names end # This method will create a climate index file. # @author phylroy.lopez@nrcan.gc.ca # @param folder [String] # @param output_file [String] def self.create_climate_index_file(folder = "#{File.dirname(__FILE__)}/../../../weather", output_file = 'C:/test/phylroy.csv') data = '' counter = 0 File.open(output_file, 'w') do |file| puts "outpus #{output_file}" data << "file,location_name,energy_plus_location_name,country,state_province_region,city,hdd10,hdd18,cdd10,cdd18,latitude,longitude,elevation, deltaDB, climate_zone, cz_standard, summer_wet_months, winter_dry_months,autumn_months, spring_months, typical_summer_wet_week, typical_winter_dry_week, typical_autumn_week, typical_spring_week, heating_design_info[1],cooling_design_info[1],extremes_design_info[1],db990\n" BTAP::FileIO.get_find_files_from_folder_by_extension(folder, 'epw').sort.each do |wfile| wf = BTAP::Environment::WeatherFile.new(wfile) data << "#{File.basename(wfile)}, #{wf.location_name}\,#{wf.energy_plus_location_name},#{wf.country}, #{wf.state_province_region}, #{wf.city}, #{wf.hdd10}, #{wf.hdd18},#{wf.cdd10},#{wf.cdd18},#{wf.latitude}, #{wf.longitude}, #{wf.elevation}, #{wf.delta_dry_bulb} ,#{wf.climate_zone},#{wf.standard},#{wf.summer_wet_months}, #{wf.winter_dry_months},#{wf.autumn_months}, #{wf.spring_months}, #{wf.typical_summer_wet_week}, #{wf.typical_winter_dry_week}, #{wf.typical_autumn_week}, #{wf.typical_spring_week},#{wf.heating_design_info[1]},#{wf.cooling_design_info[1]},#{wf.extremes_design_info[1]},#{wf.db990}\n" counter += 1 end file.write(data) end puts "parsed #{counter} weather files." end # This method will create a climate index file. # @author phylroy.lopez@nrcan.gc.ca # @param folder [String] # @param output_file [String] def self.create_climate_json_file(folder = "#{File.dirname(__FILE__)}/../../../weather", output_file = 'C:/test/phylroy.csv') data_array = [] File.open(output_file, 'w') do |file| BTAP::FileIO.get_find_files_from_folder_by_extension(folder, 'epw').sort.each do |wfile| wf = BTAP::Environment::WeatherFile.new(wfile) data = {} data_array << data data['file'] = File.basename(wfile).encode('UTF-8') data['location_name'] = wf.location_name.force_encoding('ISO-8859-1').encode('UTF-8') data['energy_plus_location_name'] = wf.energy_plus_location_name.force_encoding('ISO-8859-1').encode('UTF-8') data['country'] = wf.country.force_encoding('ISO-8859-1').encode('UTF-8') data['state_province_region'] = wf.state_province_region.force_encoding('ISO-8859-1').encode('UTF-8') data['city'] = wf.city.force_encoding('ISO-8859-1').encode('UTF-8') data['hdd10'] = wf.hdd10 data['hdd18'] = wf.hdd18 data['cdd10'] = wf.cdd10 data['cdd18'] = wf.cdd18 data['latitude'] = wf.latitude data['longitude'] = wf.longitude data['elevation'] = wf.delta_dry_bulb data['climate_zone'] = wf.climate_zone.force_encoding('ISO-8859-1').encode('UTF-8') data['standard'] = wf.standard data['summer_wet_months'] = wf.summer_wet_months.force_encoding('ISO-8859-1').encode('UTF-8') data['winter_dry_months'] = wf.autumn_months.force_encoding('ISO-8859-1').encode('UTF-8') data['spring_months'] = wf.spring_months.force_encoding('ISO-8859-1').encode('UTF-8') data['typical_summer_wet_week'] = wf.typical_summer_wet_week data['typical_winter_dry_week'] = wf.typical_winter_dry_week data['typical_autumn_week'] = wf.typical_autumn_week data['typical_spring_week'] = wf.typical_spring_week data['wf.heating_design_info[1]'] = wf.heating_design_info[1] data['cooling_design_info[1]'] = wf.cooling_design_info[1] data['extremes_design_info[1]'] = wf.extremes_design_info[1] data['db990'] = wf.db990 end File.write(output_file,JSON.pretty_generate(data_array)) end end class WeatherFile attr_accessor :location_name, :energy_plus_location_name, :latitude, :longitude, :elevation, :city, :state_province_region, :country, :hdd18, :cdd18, :hdd10, :cdd10, :heating_design_info, :cooling_design_info, :extremes_design_info, :monthly_dry_bulb, :delta_dry_bulb, :climate_zone, :standard, :summer_wet_months, :winter_dry_months, :autumn_months, :spring_months, :typical_summer_wet_week, :typical_winter_dry_week, :typical_autumn_week, :typical_spring_week, :epw_filepath, :ddy_filepath, :stat_filepath, :db990 YEAR = 0 MONTH = 1 DAY = 2 HOUR = 3 MINUTE = 4 DATA_SOURCE = 5 DRY_BULB_TEMPERATURE = 6 DEW_POINT_TEMPERATURE = 7 RELATIVE_HUMIDITY = 8 ATMOSPHERIC_STATION_PRESSURE = 9 EXTRATERRESTRIAL_HORIZONTAL_RADIATION = 10 # not used EXTRATERRESTRIAL_DIRECT_NORMAL_RADIATION = 11 # not used HORIZONTAL_INFRARED_RADIATION_INTENSITY = 12 GLOBAL_HORIZONTAL_RADIATION = 13 # not used DIRECT_NORMAL_RADIATION = 14 DIFFUSE_HORIZONTAL_RADIATION = 15 GLOBAL_HORIZONTAL_ILLUMINANCE = 16 # not used DIRECT_NORMAL_ILLUMINANCE = 17 # not used DIFFUSE_HORIZONTAL_ILLUMINANCE = 18 # not used ZENITH_LUMINANCE = 19 # not used WIND_DIRECTION = 20 WIND_SPEED = 21 TOTAL_SKY_COVER = 22 # not used OPAQUE_SKY_COVER = 23 # not used VISIBILITY = 24 # not used CEILING_HEIGHT = 25 # not used PRESENT_WEATHER_OBSERVATION = 26 PRESENT_WEATHER_CODES = 27 PRECIPITABLE_WATER = 28 # not used AEROSOL_OPTICAL_DEPTH = 29 # not used SNOW_DEPTH = 30 DAYS_SINCE_LAST_SNOWFALL = 31 # not used ALBEDO = 32 # not used LIQUID_PRECIPITATION_DEPTH = 33 LIQUID_PRECIPITATION_QUANTITY = 34 # This method initializes and returns self. # @author phylroy.lopez@nrcan.gc.ca # @param weather_file [String] # @return [String] self def initialize(weather_file) # Define the openstudio-standards weather location top_dir = File.expand_path('../../..', File.dirname(__FILE__)) weather_dir = "#{top_dir}/data/weather" # First check if the epw file exists at a full path. If not found there, # check for the file in the openstudio-standards/data/weather directory. weather_file = weather_file.to_s @epw_filepath = nil @ddy_filepath = nil @stat_filepath = nil if File.exist?(weather_file) @epw_filepath = weather_file.to_s @ddy_filepath = weather_file.sub('epw', 'ddy').to_s @stat_filepath = weather_file.sub('epw', 'stat').to_s elsif File.exist?("#{weather_dir}/#{weather_file}") @epw_filepath = "#{weather_dir}/#{weather_file}" @ddy_filepath = "#{weather_dir}/#{weather_file.sub('epw', 'ddy')}" @stat_filepath = "#{weather_dir}/#{weather_file.sub('epw', 'stat')}" else raise("Could not find weather file #{weather_file}. Make sure file path is correct.") end # Ensure that epw, ddy, and stat file all exist raise("Weather file #{@epw_filepath} not found.") unless File.exist?(@epw_filepath) && @epw_filepath.downcase.include?('.epw') raise("Weather file ddy #{@ddy_filepath} not found.") unless File.exist?(@ddy_filepath) && @ddy_filepath.downcase.include?('.ddy') raise("Weather file stat #{@stat_filepath} not found.") unless File.exist?(@stat_filepath) && @stat_filepath.downcase.include?('.stat') # load file objects. @epw_file = OpenStudio::EpwFile.new(OpenStudio::Path.new(@epw_filepath)) if OpenStudio::EnergyPlus.loadAndTranslateIdf(@ddy_filepath).empty? raise "Unable to load ddy idf file#{@ddy_filepath}." else @ddy_file = OpenStudio::EnergyPlus.loadAndTranslateIdf(@ddy_filepath).get end @stat_file = EnergyPlus::StatFile.new(@stat_filepath) # assign variables. @latitude = @epw_file.latitude @longitude = @epw_file.longitude @elevation = @epw_file.elevation @city = @epw_file.city @state_province_region = @epw_file.stateProvinceRegion @country = @epw_file.country @hdd18 = @stat_file.hdd18 @cdd18 = @stat_file.cdd18 @hdd10 = @stat_file.hdd10 @cdd10 = @stat_file.cdd10 @heating_design_info = @stat_file.heating_design_info @cooling_design_info = @stat_file.cooling_design_info @extremes_design_info = @stat_file.extremes_design_info @monthly_dry_bulb = @stat_file.monthly_dry_bulb @mean_dry_bulb = @stat_file.mean_dry_bulb @delta_dry_bulb = @stat_file.delta_dry_bulb @location_name = "#{@country}-#{@state_province_region}-#{@city}" @energy_plus_location_name = "#{@city}_#{@state_province_region}_#{@country}" @climate_zone = @stat_file.climate_zone @standard = @stat_file.standard @summer_wet_months = @stat_file.summer_wet_months @winter_dry_months = @stat_file.winter_dry_months @autumn_months = @stat_file.autumn_months @spring_months = @stat_file.spring_months @typical_summer_wet_week = @stat_file.typical_summer_wet_week @typical_winter_dry_week = @stat_file.typical_winter_dry_week @typical_autumn_week = @stat_file.typical_autumn_week @typical_spring_week = @stat_file.typical_spring_week @db990 = @heating_design_info[2] return self end # This method returns the Thermal Zone based on cdd10 and hdd18 # @author padmassun.rajakareyar@canada.ca # @return [String] thermal_zone def a169_2006_climate_zone cdd10 = self.cdd10.to_f hdd18 = self.hdd18.to_f if cdd10 > 6000 # Extremely Hot Humid (0A), Dry (0B) return 'ASHRAE 169-2006-0A' elsif (cdd10 > 5000) && (cdd10 <= 6000) # Very Hot Humid (1A), Dry (1B) return 'ASHRAE 169-2006-1A' elsif (cdd10 > 3500) && (cdd10 <= 5000) # Hot Humid (2A), Dry (2B) return 'ASHRAE 169-2006-2A' elsif ((cdd10 > 2500) && (cdd10 < 3500)) && (hdd18 <= 2000) # Warm Humid (3A), Dry (3B) return 'ASHRAE 169-2006-3A' # and 'ASHRAE 169-2006-3B' elsif (cdd10 <= 2500) && (hdd18 <= 2000) # Warm Marine (3C) return 'ASHRAE 169-2006-3C' elsif ((cdd10 > 1500) && (cdd10 < 3500)) && ((hdd18 > 2000) && (hdd18 <= 3000)) # Mixed Humid (4A), Dry (4B) return 'ASHRAE 169-2006-4A' # and 'ASHRAE 169-2006-4B' elsif (cdd10 <= 1500) && ((hdd18 > 2000) && (hdd18 <= 3000)) # Mixed Marine return 'ASHRAE 169-2006-4C' elsif ((cdd10 > 1000) && (cdd10 <= 3500)) && ((hdd18 > 3000) && (hdd18 <= 4000)) # Cool Humid (5A), Dry (5B) return 'ASHRAE 169-2006-5A' # and 'ASHRAE 169-2006-5B' elsif (cdd10 <= 1000) && ((hdd18 > 3000) && (hdd18 <= 4000)) # Cool Marine (5C) return 'ASHRAE 169-2006-5C' elsif (hdd18 > 4000) && (hdd18 <= 5000) # Cold Humid (6A), Dry (6B) return 'ASHRAE 169-2006-6A' # and 'ASHRAE 169-2006-6B' elsif (hdd18 > 5000) && (hdd18 <= 7000) # Very Cold (7) return 'ASHRAE 169-2006-7A' elsif hdd18 > 7000 # Subarctic/Arctic (8) return 'ASHRAE 169-2006-8A' else # raise ("invalid cdd10 of #{cdd10} or hdd18 of #{hdd18}") return '[INVALID]' end end # This method will set the weather file and returns a log string. # @author phylroy.lopez@nrcan.gc.ca # @param model [OpenStudio::model::Model] A model object # @return [String] log def set_weather_file(model, runner = nil) BTAP.runner_register('Info', 'BTAP::Environment::WeatherFile::set_weather', runner) OpenStudio::Model::WeatherFile.setWeatherFile(model, @epw_file) building_name = model.building.get.name weather_file_path = model.weatherFile.get.path.get BTAP.runner_register('Info', "Set model \"#{building_name}\" to weather file #{weather_file_path}.\n", runner) # Add or update site data site = model.getSite site.setName("#{@epw_file.city}_#{@epw_file.stateProvinceRegion}_#{@epw_file.country}") site.setLatitude(@epw_file.latitude) site.setLongitude(@epw_file.longitude) site.setTimeZone(@epw_file.timeZone) site.setElevation(@epw_file.elevation) BTAP.runner_register('Info', 'Setting water main temperatures via parsing of STAT file.', runner) water_temp = model.getSiteWaterMainsTemperature water_temp.setAnnualAverageOutdoorAirTemperature(@stat_file.mean_dry_bulb) water_temp.setMaximumDifferenceInMonthlyAverageOutdoorAirTemperatures(@stat_file.delta_dry_bulb) BTAP.runner_register('Info', "SiteWaterMainsTemperature.AnnualAverageOutdoorAirTemperature = #{@stat_file.mean_dry_bulb}.", runner) BTAP.runner_register('Info', "SiteWaterMainsTemperature.MaximumDifferenceInMonthlyAverageOutdoorAirTemperatures = #{@stat_file.delta_dry_bulb}.", runner) # Remove all the Design Day objects that are in the file model.getObjectsByType('OS:SizingPeriod:DesignDay'.to_IddObjectType).each(&:remove) # Load in the ddy file based on convention that it is in the same directory and has the same basename as the weather @ddy_file.getObjectsByType('OS:SizingPeriod:DesignDay'.to_IddObjectType).each do |d| # grab only the ones that matter ddy_list = /(Htg 99.6. Condns DB)|(Clg .4. Condns WB=>MDB)|(Clg .4% Condns DB=>MWB)/ if d.name.get =~ ddy_list BTAP.runner_register('Info', "Adding design day '#{d.name}'.", runner) # add the object to the existing model model.addObject(d.clone) end end return true end # This method scans the epw file into memory. # @author phylroy.lopez@nrcan.gc.ca def scan @filearray = [] file = File.new(@epw_filepath, 'r') while (line = file.gets) @filearray.push(line.split(',')) end file.close end # This method will sets column to a value. # @author phylroy.lopez@nrcan.gc.ca # @param column [String] # @param value [Fixnum] def setcolumntovalue(column, value) @filearray.each do |line| unless line.first =~ /\D(.*)/ line[column] = value end end end # This method will eliminate all radiation from the weather and returns self. # @author phylroy.lopez@nrcan.gc.ca # @return [String] self def eliminate_all_radiation scan if @filearray.nil? setcolumntovalue(EXTRATERRESTRIAL_HORIZONTAL_RADIATION, '0') # not used setcolumntovalue(EXTRATERRESTRIAL_DIRECT_NORMAL_RADIATION, '0') # not used setcolumntovalue(HORIZONTAL_INFRARED_RADIATION_INTENSITY, '315') setcolumntovalue(GLOBAL_HORIZONTAL_RADIATION, '0') # not used setcolumntovalue(DIRECT_NORMAL_RADIATION, '0') setcolumntovalue(DIFFUSE_HORIZONTAL_RADIATION, '0') setcolumntovalue(TOTAL_SKY_COVER, '10') # not used setcolumntovalue(OPAQUE_SKY_COVER, '10') # not used setcolumntovalue(VISIBILITY, '0') # not used setcolumntovalue(CEILING_HEIGHT, '0') # not used # lux values setcolumntovalue(GLOBAL_HORIZONTAL_ILLUMINANCE, '0') # not used setcolumntovalue(DIRECT_NORMAL_ILLUMINANCE, '0') # not used setcolumntovalue(DIFFUSE_HORIZONTAL_ILLUMINANCE, '0') # not used setcolumntovalue(ZENITH_LUMINANCE, '0') # not used return self end # This method will eliminate solar radiation and returns self. # @author phylroy.lopez@nrcan.gc.ca # @return [String] self def eliminate_only_solar_radiation scan if @filearray.nil? setcolumntovalue(GLOBAL_HORIZONTAL_RADIATION, '0') # not used setcolumntovalue(DIRECT_NORMAL_RADIATION, '0') setcolumntovalue(DIFFUSE_HORIZONTAL_RADIATION, '0') return self end # This method will eliminate all radiation except solar and returns self. # @author phylroy.lopez@nrcan.gc.ca # @return [String] self def eliminate_all_radiation_except_solar scan if @filearray.nil? setcolumntovalue(EXTRATERRESTRIAL_HORIZONTAL_RADIATION, '0') # not used setcolumntovalue(EXTRATERRESTRIAL_DIRECT_NORMAL_RADIATION, '0') # not used setcolumntovalue(HORIZONTAL_INFRARED_RADIATION_INTENSITY, '315') setcolumntovalue(TOTAL_SKY_COVER, '10') # not used setcolumntovalue(OPAQUE_SKY_COVER, '10') # not used setcolumntovalue(VISIBILITY, '0') # not used setcolumntovalue(CEILING_HEIGHT, '0') # not used # lux values setcolumntovalue(GLOBAL_HORIZONTAL_ILLUMINANCE, '0') # not used setcolumntovalue(DIRECT_NORMAL_ILLUMINANCE, '0') # not used setcolumntovalue(DIFFUSE_HORIZONTAL_ILLUMINANCE, '0') # not used setcolumntovalue(ZENITH_LUMINANCE, '0') # not used return self end # This method will eliminate percipitation and returns self. # @author phylroy.lopez@nrcan.gc.ca # @return [String] self def eliminate_percipitation scan if @filearray.nil? setcolumntovalue(PRESENT_WEATHER_OBSERVATION, '0') setcolumntovalue(PRESENT_WEATHER_CODES, '999999999') # no weather. Clear day. setcolumntovalue(SNOW_DEPTH, '0') setcolumntovalue(LIQUID_PRECIPITATION_DEPTH, '0') setcolumntovalue(LIQUID_PRECIPITATION_QUANTITY, '0') return self end # This method eliminates wind and returns self. # @author phylroy.lopez@nrcan.gc.ca # @return [String] self def eliminate_wind scan if @filearray.nil? setcolumntovalue(WIND_DIRECTION, '0') setcolumntovalue(WIND_SPEED, '0') return self end # This method sets Constant Dry and Dew Point Temperature Humidity And Pressure and returns self. # @author phylroy.lopez@nrcan.gc.ca # @param dbt [Float] dry bulb temperature # @param dpt [Float] dew point temperature # @param hum [Fixnum] humidity # @param press [Fixnum] pressure # @return [String] self def set_constant_dry_and_dewpoint_temperature_humidity_pressure(dbt = '0.0', dpt = '-1.1', hum = '92', press = '98500') scan if @filearray.nil? setcolumntovalue(DRY_BULB_TEMPERATURE, dbt) setcolumntovalue(DEW_POINT_TEMPERATURE, dpt) setcolumntovalue(RELATIVE_HUMIDITY, hum) setcolumntovalue(ATMOSPHERIC_STATION_PRESSURE, press) return self end # This method writes to a file. # @author phylroy.lopez@nrcan.gc.ca # @param filename [String] def writetofile(filename) scan if @filearray.nil? begin FileUtils.mkdir_p(File.dirname(filename)) file = File.open(filename, 'w') @filearray.each do |line| firstvalue = true newline = '' line.each do |value| if firstvalue == true firstvalue = false else newline += ',' end newline += value end file.puts(newline) end rescue IOError => e # some error occur, dir not writable etc. ensure file.close unless file.nil? end # copies original file FileUtils.cp(@ddy_filepath, "#{File.dirname(filename)}/#{File.basename(filename, '.epw')}.ddy") FileUtils.cp(@stat_filepath, "#{File.dirname(filename)}/#{File.basename(filename, '.epw')}.stat") end end # Environment end end