urbanopt-cli 0.8.0 → 0.8.1

data/example_files/mappers/Baseline.rb

@@ -37,11 +37,13 @@ require 'openstudio/load_flexibility_measures'
 
 require 'json'
 require 'rexml/document'
+require 'logger'
 
 module URBANopt
   module Scenario
     class BaselineMapper < SimulationMapperBase
       # class level variables
+      @@logger = Logger.new($stdout)
       @@instance_lock = Mutex.new
       @@osw = nil
       @@geometry = nil
@@ -383,222 +385,215 @@ module URBANopt
       end
 
       # epw_state to subregions mapping methods
-      #REK: Maybe we can move these method to the geojson gem
+      # REK: Maybe we can move these method to the geojson gem
       def get_future_emissions_region(feature)
         # Options are: AZNMc, CAMXc, ERCTc, FRCCc, MROEc, MROWc, NEWEc, NWPPc, NYSTc, RFCEc, RFCMc, RFCWc, RMPAc, SPNOc, SPSOc, SRMVc, SRMWc, SRSOc, SRTVc, and SRVCc
         # egrid subregions can map directly to zipcodes but not to states. Some state might include multiple egrid subregions. the default mapper prioritize the egrid subregion that is most common in the state (covers the biggest number of zipcodes)
-        future_emissions_mapping_hash =   
-        {'FL': 'FRCCc', #['FRCCc', 'SRSOc']
-        'MS': 'SRMVc', #['SRMVc', 'SRTVc']
-        'NE': 'MROWc', #['MROWc', 'RMPAc']
-        'OR': 'NWPPc',
-        'CA': 'CAMXc', #['CAMXc', 'NWPPc']
-        'VA': 'SRVCc', #['SRVCc', 'RFCWc', 'RFCEc'],
-        'AR': 'SRMVc', #['SRMVc', 'SPSOc']
-        'TX': 'ERCTc', #['ERCTc', 'SRMVc', 'SPSOc', 'AZNMc']
-        'OH': 'RFCWc', 
-        'UT': 'NWPPc',
-        'MT': 'NWPPc', #['NWPPc', 'MROWc']
-        'TN': 'SRTVc',
-        'ID': 'NWPPc',
-        'WI': 'MROEc', #['RFCWc', 'MROEc', 'MROWc']
-        'WV': 'RFCWc',
-        'NC': 'SRVCc',
-        'LA': 'SRMVc',
-        'IL': 'SRMWc', #['RFCWc', 'SRMWc']
-        'OK': 'SPSOc',
-        'IA': 'MROWc',
-        'WA': 'NWPPc',
-        'SD': 'MROWc', #['MROWc', 'RMPAc']
-        'MN': 'MROWc',
-        'KY': 'SRTVc', #['SRTVc', 'RFCWc']
-        'MI': 'RFCMc', #['RFCMc', 'MROEc']
-        'KS': 'SPNOc',
-        'NJ': 'RFCEc',
-        'NY': 'NYSTc',
-        'IN': 'RFCWc',
-        'VT': 'NEWEc',
-        'NM': 'AZNMc', #['AZNMc', 'SPSOc']
-        'WY': 'RMPAc', #['RMPAc', 'NWPPc']
-        'GA': 'SRSOc',
-        'MO': 'SRMWc', #['SRMWc', 'SPNOc']
-        'DC': 'RFCEc',
-        'SC': 'SRVCc',
-        'PA': 'RFCEc', #['RFCEc', 'RFCWc']
-        'CO': 'RMPAc',
-        'AZ': 'AZNMc',
-        'ME': 'NEWEc',
-        'AL': 'SRSOc',
-        'MD': 'RFCEc', #['RFCEc', 'RFCWc']
-        'NH': 'NEWEc',
-        'MA': 'NEWEc',
-        'ND': 'MROWc',
-        'NV': 'NWPPc', #['NWPPc', 'AZNMc']
-        'CT': 'NEWEc',
-        'DE': 'RFCEc',
-        'RI': 'NEWEc'}
-
-        #get the state from weather file
+        future_emissions_mapping_hash =
+          { 'FL': 'FRCCc', # ['FRCCc', 'SRSOc']
+            'MS': 'SRMVc', # ['SRMVc', 'SRTVc']
+            'NE': 'MROWc', # ['MROWc', 'RMPAc']
+            'OR': 'NWPPc',
+            'CA': 'CAMXc', # ['CAMXc', 'NWPPc']
+            'VA': 'SRVCc', # ['SRVCc', 'RFCWc', 'RFCEc'],
+            'AR': 'SRMVc', # ['SRMVc', 'SPSOc']
+            'TX': 'ERCTc', # ['ERCTc', 'SRMVc', 'SPSOc', 'AZNMc']
+            'OH': 'RFCWc',
+            'UT': 'NWPPc',
+            'MT': 'NWPPc', # ['NWPPc', 'MROWc']
+            'TN': 'SRTVc',
+            'ID': 'NWPPc',
+            'WI': 'MROEc', # ['RFCWc', 'MROEc', 'MROWc']
+            'WV': 'RFCWc',
+            'NC': 'SRVCc',
+            'LA': 'SRMVc',
+            'IL': 'SRMWc', # ['RFCWc', 'SRMWc']
+            'OK': 'SPSOc',
+            'IA': 'MROWc',
+            'WA': 'NWPPc',
+            'SD': 'MROWc', # ['MROWc', 'RMPAc']
+            'MN': 'MROWc',
+            'KY': 'SRTVc', # ['SRTVc', 'RFCWc']
+            'MI': 'RFCMc', # ['RFCMc', 'MROEc']
+            'KS': 'SPNOc',
+            'NJ': 'RFCEc',
+            'NY': 'NYSTc',
+            'IN': 'RFCWc',
+            'VT': 'NEWEc',
+            'NM': 'AZNMc', # ['AZNMc', 'SPSOc']
+            'WY': 'RMPAc', # ['RMPAc', 'NWPPc']
+            'GA': 'SRSOc',
+            'MO': 'SRMWc', # ['SRMWc', 'SPNOc']
+            'DC': 'RFCEc',
+            'SC': 'SRVCc',
+            'PA': 'RFCEc', # ['RFCEc', 'RFCWc']
+            'CO': 'RMPAc',
+            'AZ': 'AZNMc',
+            'ME': 'NEWEc',
+            'AL': 'SRSOc',
+            'MD': 'RFCEc', # ['RFCEc', 'RFCWc']
+            'NH': 'NEWEc',
+            'MA': 'NEWEc',
+            'ND': 'MROWc',
+            'NV': 'NWPPc', # ['NWPPc', 'AZNMc']
+            'CT': 'NEWEc',
+            'DE': 'RFCEc',
+            'RI': 'NEWEc' }
+
+        # get the state from weather file
         state = feature.weather_filename.split('_', -1)[1]
-      
-        #find region input based on the state
+
+        # find region input based on the state
         region = future_emissions_mapping_hash[state.to_sym]
 
-        puts "emissions_future_subregion for #{state} is assigned to: #{region}"
-        puts "You can overwrite this assigned input by specifiying the emissions_future_subregion input in the FeatureFile"
+        @@logger.warn("emissions_future_subregion for #{state} is assigned to: #{region}. Note: Not all states have a 1 to 1 mapping with a subregion. Some states('ND','IN', 'MN', 'SD', 'IA', 'WV', 'OH', 'NE' ) include 2 subregions.
+        The default mapper maps to the subregion that includes the most zipcodes in the corresponding state. You can overwrite this assigned input by specifiying the emissions_future_subregion input in the FeatureFile.")
 
         return region
-
       end
 
       def get_hourly_historical_emissions_region(feature)
-      
         # Options are: California, Carolinas, Central, Florida, Mid-Atlantic, Midwest, New England, New York, Northwest, Rocky Mountains, Southeast, Southwest, Tennessee, and Texas
-        # There is no "correct" mapping of eGrid to AVERT regions as they are both large geographical areas that partially overlap. 
-        # Mapping is done using mapping tools from eGrid and AVERT (ZipCode for eGrid and fraction of state for AVERT). 
+        # There is no "correct" mapping of eGrid to AVERT regions as they are both large geographical areas that partially overlap.
+        # Mapping is done using mapping tools from eGrid and AVERT (ZipCode for eGrid and fraction of state for AVERT).
         # Mapped based on the maps of each set of regions:
-        hourly_historical_mapping_hash = 
-        {'FL': 'Florida',
-        'MS': 'Midwest',
-        'NE': 'Midwest',#MRWO could be Midwest / Central
-        'OR': 'Northwest',
-        'CA': 'California',
-        'VA': 'Carolinas',
-        'AR': 'Midwest',
-        'TX': 'Texas',
-        'OH': 'Midwest',#RFCW could be Midwest / Mid Atlantic
-        'UT': 'Northwest',
-        'MT': 'Northwest',
-        'TN': 'Tennessee',
-        'ID': 'Northwest',
-        'WI': 'Midwest',
-        'WV': 'Midwest', #RFCW could be Midwest / Mid Atlantic
-        'NC': 'Carolinas',
-        'LA': 'Midwest',
-        'IL': 'Midwest',
-        'OK': 'Central',
-        'IA': 'Midwest', #MRWO could be Midwest / Central
-        'WA': 'Northwest',
-        'SD': 'Midwest',#MRWO could be Midwest / Central
-        'MN': 'Midwest',#MRWO could be Midwest / Central
-        'KY': 'Tennessee',
-        'MI': 'Midwest',
-        'KS': 'Central',
-        'NJ': 'Mid-Atlantic',
-        'NY': 'New York',
-        'IN': 'Midwest', #RFCW could be Midwest / Mid Atlantic
-        'VT': 'New England',
-        'NM': 'Southwest',
-        'WY': 'Rocky Mountains',
-        'GA': 'SRSO',
-        'MO': 'Midwest',
-        'DC': 'Mid-Atlantic',
-        'SC': 'Carolinas',
-        'PA': 'Mid-Atlantic',
-        'CO': 'Rocky Mountains',
-        'AZ': 'Southwest',
-        'ME': 'New England',
-        'AL': 'Southeast',
-        'MD': 'Mid-Atlantic',
-        'NH': 'New England',
-        'MA': 'New England',
-        'ND': 'Midwest',#MRWO could be Midwest / Central
-        'NV': 'Northwest',
-        'CT': 'New England',
-        'DE': 'Mid-Atlantic',
-        'RI': 'New England'}
-
-        #get the state from weather file
+        hourly_historical_mapping_hash =
+          { 'FL': 'Florida',
+            'MS': 'Midwest',
+            'NE': 'Midwest',  # MRWO could be Midwest / Central
+            'OR': 'Northwest',
+            'CA': 'California',
+            'VA': 'Carolinas',
+            'AR': 'Midwest',
+            'TX': 'Texas',
+            'OH': 'Midwest',  # RFCW could be Midwest / Mid Atlantic
+            'UT': 'Northwest',
+            'MT': 'Northwest',
+            'TN': 'Tennessee',
+            'ID': 'Northwest',
+            'WI': 'Midwest',
+            'WV': 'Midwest', # RFCW could be Midwest / Mid Atlantic
+            'NC': 'Carolinas',
+            'LA': 'Midwest',
+            'IL': 'Midwest',
+            'OK': 'Central',
+            'IA': 'Midwest', # MRWO could be Midwest / Central
+            'WA': 'Northwest',
+            'SD': 'Midwest',  # MRWO could be Midwest / Central
+            'MN': 'Midwest',  # MRWO could be Midwest / Central
+            'KY': 'Tennessee',
+            'MI': 'Midwest',
+            'KS': 'Central',
+            'NJ': 'Mid-Atlantic',
+            'NY': 'New York',
+            'IN': 'Midwest', # RFCW could be Midwest / Mid Atlantic
+            'VT': 'New England',
+            'NM': 'Southwest',
+            'WY': 'Rocky Mountains',
+            'GA': 'SRSO',
+            'MO': 'Midwest',
+            'DC': 'Mid-Atlantic',
+            'SC': 'Carolinas',
+            'PA': 'Mid-Atlantic',
+            'CO': 'Rocky Mountains',
+            'AZ': 'Southwest',
+            'ME': 'New England',
+            'AL': 'Southeast',
+            'MD': 'Mid-Atlantic',
+            'NH': 'New England',
+            'MA': 'New England',
+            'ND': 'Midwest', # MRWO could be Midwest / Central
+            'NV': 'Northwest',
+            'CT': 'New England',
+            'DE': 'Mid-Atlantic',
+            'RI': 'New England' }
+
+        # get the state from weather file
         state = feature.weather_filename.split('_', -1)[1]
 
-        #find region input based on the state
-        region = hourly_historical_mapping_hash[state.to_sym]       
-        puts "emissions_hourly_historical_subregion for #{state} is assigned to: #{region}"
-        puts "You can overwrite this assigned input by specifiying the emissions_hourly_historical_subregion input in the FeatureFile"
-
-        return region   
+        # find region input based on the state
+        region = hourly_historical_mapping_hash[state.to_sym]
+        @@logger.warn("emissions_hourly_historical_subregion for #{state} is assigned to: #{region}. Note: Not all states have a 1 to 1 mapping with a subregion. Some states('ND','IN', 'MN', 'SD', 'IA', 'WV', 'OH', 'NE' ) include 2 subregions.
+        The default mapper maps to the subregion that includes the most zipcodes in the corresponding state. You can overwrite this assigned input by specifiying the emissions_hourly_historical_subregion input in the FeatureFile.")
 
+        return region
       end
 
-
       def get_annual_historical_emissions_region(feature)
-      
         # Options are: AKGD, AKMS, AZNM, CAMX, ERCT, FRCC, HIMS, HIOA, MROE, MROW, NEWE, NWPP, NYCW, NYLI, NYUP, RFCE, RFCM, RFCW, RMPA, SPNO, SPSO, SRMV, SRMW, SRSO, SRTV, and SRVC
         # egrid subregions can map directly to zipcodes but not to states. Some state might include multiple egrid subregions. the default mapper prioritize the egrid subregion that is most common in the state (covers the biggest number of zipcodes)
-        annual_historical_mapping_hash = 
-        {'FL': 'FRCC',
-        'MS': 'SRMV',
-        'NE': 'MROW',
-        'OR': 'NWPP',
-        'CA': 'CAMX',
-        'VA': 'SRVC',
-        'AR': 'SRMV',
-        'TX': 'ERCT',
-        'OH': 'RFCW',
-        'UT': 'NWPP',
-        'MT': 'NWPP',
-        'TN': 'SRTV',
-        'ID': 'NWPP',
-        'WI': 'MROE',
-        'WV': 'RFCW',
-        'NC': 'SRVC',
-        'LA': 'SRMV',
-        'IL': 'SRMW',
-        'OK': 'SPSO',
-        'IA': 'MROW',
-        'WA': 'NWPP',
-        'SD': 'MROW',
-        'MN': 'MROW',
-        'KY': 'SRTV',
-        'MI': 'RFCM',
-        'KS': 'SPNO',
-        'NJ': 'RFCE',
-        'NY': 'NYCW',
-        'IN': 'RFCW',
-        'VT': 'NEWE',
-        'NM': 'AZNM',
-        'WY': 'RMPA',
-        'GA': 'SRSO',
-        'MO': 'SRMW',
-        'DC': 'RFCE',
-        'SC': 'SRVC',
-        'PA': 'RFCE',
-        'CO': 'RMPA',
-        'AZ': 'AZNM',
-        'ME': 'NEWE',
-        'AL': 'SRSO',
-        'MD': 'RFCE',
-        'NH': 'NEWE',
-        'MA': 'NEWE',
-        'ND': 'MROW',
-        'NV': 'NWPP',
-        'CT': 'NEWE',
-        'DE': 'RFCE',
-        'RI': 'NEWE'}
-        #get the state from weather file
+        annual_historical_mapping_hash =
+          { 'FL': 'FRCC',
+            'MS': 'SRMV',
+            'NE': 'MROW',
+            'OR': 'NWPP',
+            'CA': 'CAMX',
+            'VA': 'SRVC',
+            'AR': 'SRMV',
+            'TX': 'ERCT',
+            'OH': 'RFCW',
+            'UT': 'NWPP',
+            'MT': 'NWPP',
+            'TN': 'SRTV',
+            'ID': 'NWPP',
+            'WI': 'MROE',
+            'WV': 'RFCW',
+            'NC': 'SRVC',
+            'LA': 'SRMV',
+            'IL': 'SRMW',
+            'OK': 'SPSO',
+            'IA': 'MROW',
+            'WA': 'NWPP',
+            'SD': 'MROW',
+            'MN': 'MROW',
+            'KY': 'SRTV',
+            'MI': 'RFCM',
+            'KS': 'SPNO',
+            'NJ': 'RFCE',
+            'NY': 'NYCW',
+            'IN': 'RFCW',
+            'VT': 'NEWE',
+            'NM': 'AZNM',
+            'WY': 'RMPA',
+            'GA': 'SRSO',
+            'MO': 'SRMW',
+            'DC': 'RFCE',
+            'SC': 'SRVC',
+            'PA': 'RFCE',
+            'CO': 'RMPA',
+            'AZ': 'AZNM',
+            'ME': 'NEWE',
+            'AL': 'SRSO',
+            'MD': 'RFCE',
+            'NH': 'NEWE',
+            'MA': 'NEWE',
+            'ND': 'MROW',
+            'NV': 'NWPP',
+            'CT': 'NEWE',
+            'DE': 'RFCE',
+            'RI': 'NEWE' }
+        # get the state from weather file
         state = feature.weather_filename.split('_', -1)[1]
-      
-        #finf region input based on the state
+
+        # finf region input based on the state
         region = annual_historical_mapping_hash[state.to_sym]
 
-        puts "emissions_annual_historical_subregion for #{state} is assigned to: #{region}"
-        puts "You can overwrite this assigned input by specifiying the emissions_annual_historical_subregion input in the FeatureFile"
-        
+        @@logger.warn("electricity_emissions_annual_historical_subregion for #{state} is assigned to: #{region}. Note: Not all states have a 1 to 1 mapping with a subregion. Some states('ND','IN', 'MN', 'SD', 'IA', 'WV', 'OH', 'NE' ) include 2 subregions.
+        The default mapper maps to the subregion that includes the most zipcodes in the corresponding state. You can overwrite this assigned input by specifiying the electricity_emissions_annual_historical_subregion input in the FeatureFile.")
+
         return region
+      end
 
-      end 
-      
-      def is_defined(feature, method_name, raise_error=true)
-        begin
-          if feature.method_missing(method_name)
-            return true
-          end
-        rescue NoMethodError
-          if raise_error
-            raise "*** ERROR *** #{method_name} is not set on this feature"
-          end
-          return false
+      def is_defined(feature, method_name, raise_error = true)
+        if feature.method_missing(method_name)
+          return true
         end
+      rescue NoMethodError
+        if raise_error
+          raise "*** ERROR *** #{method_name} is not set on this feature"
+        end
+
+        return false
       end
 
       def create_osw(scenario, features, feature_names)
@@ -641,8 +636,8 @@ module URBANopt
             # Check for required residential fields
             is_defined(feature, :number_of_stories_above_ground)
             is_defined(feature, :foundation_type)
-            
-            if not is_defined(feature, :hpxml_directory, false)
+
+            if !is_defined(feature, :hpxml_directory, false)
               # check additional fields when HPXML dir is not given
               is_defined(feature, :attic_type)
               is_defined(feature, :number_of_bedrooms)
@@ -782,7 +777,7 @@ module URBANopt
 
             args[:geometry_unit_num_occupants] = 'auto'
             begin
-              args[:geometry_unit_num_occupants] = "#{feature.number_of_occupants / args[:geometry_building_num_units]}"
+              args[:geometry_unit_num_occupants] = (feature.number_of_occupants / args[:geometry_building_num_units]).to_s
             rescue StandardError
             end
 
@@ -1018,11 +1013,10 @@ module URBANopt
           elsif commercial_building_types.include? building_type
             # set_run_period
             OpenStudio::Extension.set_measure_argument(osw, 'set_run_period', '__SKIP__', false)
-            
             # can enable reporting (commercial building types only for now)
-            #OpenStudio::Extension.set_measure_argument(osw, 'openstudio_results', '__SKIP__', false)
-            #OpenStudio::Extension.set_measure_argument(osw, 'envelope_and_internal_load_breakdown', '__SKIP__', false)
-            #OpenStudio::Extension.set_measure_argument(osw, 'generic_qaqc', '__SKIP__', false)
+            # OpenStudio::Extension.set_measure_argument(osw, 'openstudio_results', '__SKIP__', false)
+            # OpenStudio::Extension.set_measure_argument(osw, 'envelope_and_internal_load_breakdown', '__SKIP__', false)
+            # OpenStudio::Extension.set_measure_argument(osw, 'generic_qaqc', '__SKIP__', false)
 
             begin
               timesteps_per_hour = feature.timesteps_per_hour
@@ -1353,111 +1347,144 @@ module URBANopt
           else
             raise "Building type #{building_type} not currently supported."
           end
+
         end
-        
-        ####### Emissions Adition
+
+        ######## Emissions Adition from add_ems_emissions_reporting
         if feature_type == 'Building'
 
+          # emissions options
+          future_regions = ['AZNMc', 'CAMXc', 'ERCTc', 'FRCCc', 'MROEc', 'MROWc', 'NEWEc', 'NWPPc', 'NYSTc', 'RFCEc', 'RFCMc', 'RFCWc', 'RMPAc', 'SPNOc', 'SPSOc', 'SRMVc', 'SRMWc', 'SRSOc', 'SRTVc', 'SRVCc']
+          hourly_historical_regions = ['California', 'Carolinas', 'Central', 'Florida', 'Mid-Atlantic', 'Midwest', 'New England', 'New York', 'Northwest', 'Rocky Mountains', 'Southeast', 'Southwest', 'Tennessee', 'Texas']
+          annual_historical_regions = ['AKGD', 'AKMS', 'AZNM', 'CAMX', 'ERCT', 'FRCC', 'HIMS', 'HIOA', 'MROE', 'MROW', 'NEWE', 'NWPP', 'NYCW', 'NYLI', 'NYUP', 'RFCE', 'RFCM', 'RFCW', 'RMPA', 'SPNO', 'SPSO', 'SRMV', 'SRMW', 'SRSO', 'SRTV', 'SRVC']
+          annual_historical_years = ['2007', '2009', '2010', '2012', '2014', '2016', '2018', '2019']
+          future_years = ['2020', '2022', '2024', '2026', '2028', '2030', '2032', '2034', '2036', '2038', '2040', '2042', '2044', '2046', '2048', '2050']
+          hourly_historical_years = ['2019']
+
           # add Emissions
           emissions = nil
 
           begin
             emissions = feature.emissions
-          rescue
+          rescue StandardError
           end
-          
+
           if emissions != true
-            puts "Emissions is not activated for this feature. Please set emissions to true in the the Feature properties in the GeoJSON file to add emissions results."
+            @@logger.info('Emissions is not activated for this feature. Please set emissions to true in the the Feature properties in the GeoJSON file to add emissions results.')
 
           elsif emissions == true
-            
-            #activate emissions measure
+
+            # activate emissions measure
             OpenStudio::Extension.set_measure_argument(osw, 'add_ems_emissions_reporting', '__SKIP__', false)
 
-            #get emissions inputs if they are available or get them from the mapping methods if the are not
+            # get emissions inputs if they are available or get them from the mapping methods if the are not
             begin
-              emissions_future_subregion = feature.emissions_future_subregion
-            rescue
-              puts "\nemission_future_subregion is not assigned for feature #{feature_id}. Defining subregion based on the State...."
-              emissions_future_subregion =  get_future_emissions_region(feature)
+              electricity_emissions_future_subregion = feature.electricity_emissions_future_subregion
+            rescue StandardError
+              @@logger.info("\nelectricity_emission_future_subregion is not assigned for feature #{feature_id}. Defining subregion based on the State....")
+              electricity_emissions_future_subregion = get_future_emissions_region(feature)
             end
 
             begin
-              emissions_hourly_historical_subregion = feature.emissions_hourly_historical_subregion
-            rescue
-              puts "\nemissions_hourly_historical_subregion is not assigned for feature #{feature_id}. Defining subregion based on the State...."
-              emissions_hourly_historical_subregion =  get_hourly_historical_emissions_region(feature)
+              electricity_emissions_hourly_historical_subregion = feature.electricity_emissions_hourly_historical_subregion
+            rescue StandardError
+              @@logger.info("\nelectricity_emissions_hourly_historical_subregion is not assigned for feature #{feature_id}. Defining subregion based on the State....")
+              electricity_emissions_hourly_historical_subregion = get_hourly_historical_emissions_region(feature)
             end
 
             begin
-              emissions_annual_historical_subregion = feature.emissions_annual_historical_subregion
-            rescue
-              puts "\nemissions_annual_historical_subregion is not assigned for feature #{feature_id}. Defining subregion based on the State...."
-              emissions_annual_historical_subregion =  get_annual_historical_emissions_region(feature)
+              electricity_emissions_annual_historical_subregion = feature.electricity_emissions_annual_historical_subregion
+            rescue StandardError
+              @@logger.info("\nelectricity_emissions_annual_historical_subregion is not assigned for feature #{feature_id}. Defining subregion based on the State....")
+              electricity_emissions_annual_historical_subregion = get_annual_historical_emissions_region(feature)
             end
 
             begin
-              emissions_future_year = feature.emissions_future_year
-            rescue
-              puts "emissions_future_year should be assigned !"
+              electricity_emissions_future_year = feature.electricity_emissions_future_year
+            rescue StandardError
+              @@logger.info("\nelectricity_emissions_future_year was not assigned by the user. The assigned default value is 2030")
+              electricity_emissions_future_year = '2030'
             end
 
             begin
-              emissions_hourly_historical_year = feature.emissions_hourly_historical_year
-            rescue
-              puts "emissions_hourly_historical_year should be assigned !"
+              electricity_emissions_hourly_historical_year = feature.electricity_emissions_hourly_historical_year
+            rescue StandardError
+              @@logger.info("\nelectricity_emissions_hourly_historical_year was not assigned by the user. The assigned default value is 2019")
+              electricity_emissions_hourly_historical_year = '2019'
             end
 
             begin
-              emissions_annual_historical_year = feature.emissions_annual_historical_year
-            rescue
-              puts "emissions_annual_historical_year should be assigned !"
+              electricity_emissions_annual_historical_year = feature.electricity_emissions_annual_historical_year
+            rescue StandardError
+              @@logger.info("\nelectricity_emissions_annual_historical_year was not assigned by the user. The assigned default value is 2019")
+              electricity_emissions_annual_historical_year = '2019'
             end
 
-            puts "\n building #{feature_id} emission inputs summarry: 
-              emissions_future_subregion = #{emissions_future_subregion}; 
-              emissions_hourly_historical_subregion = #{emissions_hourly_historical_subregion}; 
-              emissions_annual_historical_subregion = #{emissions_annual_historical_subregion}; 
-              emissions_future_year = #{emissions_future_year};
-              emissions_hourly_historical_year = #{emissions_hourly_historical_year};
-              emissions_annual_historical_year = #{emissions_annual_historical_year} \n "
+            # puts "\n building #{feature_id} emission inputs summarry:
+            # electricity_emissions_future_subregion = #{electricity_emissions_future_subregion};
+            #   electricity_emissions_hourly_historical_subregion = #{electricity_emissions_hourly_historical_subregion};
+            #   electricity_emissions_annual_historical_subregion = #{electricity_emissions_annual_historical_subregion};
+            #   electricity_emissions_future_year = #{electricity_emissions_future_year};
+            #   electricity_emissions_hourly_historical_year = #{electricity_emissions_hourly_historical_year};
+            #   electricity_emissions_annual_historical_year = #{electricity_emissions_annual_historical_year}\n "
 
             ## Assign the OS measure arguments
             begin
-              
-              #emissions_future_subregion
-              if !emissions_future_subregion.nil? && !emissions_future_subregion.empty?
-                OpenStudio::Extension.set_measure_argument(osw, 'add_ems_emissions_reporting', 'future_subregion', emissions_future_subregion)
+              # emissions_future_subregion
+              if !electricity_emissions_future_subregion.nil? && !electricity_emissions_future_subregion.empty?
+                if future_regions.include? electricity_emissions_future_subregion
+                  OpenStudio::Extension.set_measure_argument(osw, 'add_ems_emissions_reporting', 'future_subregion', electricity_emissions_future_subregion)
+                else
+                  @@logger.error(" '#{electricity_emissions_future_subregion}' is not valid option for electricity_emissions_future_subregion. Please choose an input from #{future_regions}")
+                end
               end
 
-              #hourly_historical_subregion
-              if !emissions_hourly_historical_subregion.nil? && !emissions_hourly_historical_subregion.empty?
-                OpenStudio::Extension.set_measure_argument(osw, 'add_ems_emissions_reporting', 'hourly_historical_subregion', emissions_hourly_historical_subregion)
+              # hourly_historical_subregion
+              if !electricity_emissions_hourly_historical_subregion.nil? && !electricity_emissions_hourly_historical_subregion.empty?
+                if hourly_historical_regions.include? electricity_emissions_hourly_historical_subregion
+                  OpenStudio::Extension.set_measure_argument(osw, 'add_ems_emissions_reporting', 'hourly_historical_subregion', electricity_emissions_hourly_historical_subregion)
+                else
+                  @@logger.error(" '#{electricity_emissions_hourly_historical_subregion}' is not valid option for electricity_emissions_hourly_historical_subregion. Please choose an input from #{hourly_historical_regions}")
+                end
               end
 
-              #annual_historical_subregion
-              if !emissions_annual_historical_subregion.nil? && !emissions_annual_historical_subregion.empty?
-                OpenStudio::Extension.set_measure_argument(osw, 'add_ems_emissions_reporting', 'annual_historical_subregion', emissions_annual_historical_subregion)
+              # annual_historical_subregion
+              if !electricity_emissions_annual_historical_subregion.nil? && !electricity_emissions_annual_historical_subregion.empty?
+                if annual_historical_regions.include? electricity_emissions_annual_historical_subregion
+                  OpenStudio::Extension.set_measure_argument(osw, 'add_ems_emissions_reporting', 'annual_historical_subregion', electricity_emissions_annual_historical_subregion)
+                else
+                  @@logger.error(" '#{electricity_emissions_annual_historical_subregion}' is not valid option for electricity_emissions_annual_historical_subregion. Please choose an input from #{annual_historical_regions}")
+                end
               end
 
-              #future_year
-              if !emissions_future_year.nil? && !emissions_future_year.empty?
-                OpenStudio::Extension.set_measure_argument(osw, 'add_ems_emissions_reporting', 'future_year', emissions_future_year)              
+              # future_year
+              if !electricity_emissions_future_year.nil? && !electricity_emissions_future_year.empty?
+
+                if future_years.include? electricity_emissions_future_year
+                  OpenStudio::Extension.set_measure_argument(osw, 'add_ems_emissions_reporting', 'future_year', electricity_emissions_future_year)
+                else
+                  @@logger.error(" '#{electricity_emissions_future_year}' is not valid option for electricity_emissions_future_year. Please choose an input from #{future_years}")
+                end
               end
 
-              #hourly_historical_year
-              if !emissions_hourly_historical_year.nil? && !emissions_hourly_historical_year.empty?
-                OpenStudio::Extension.set_measure_argument(osw, 'add_ems_emissions_reporting', 'hourly_historical_year', emissions_hourly_historical_year)
-              else
-                
+              # hourly_historical_year
+              if !electricity_emissions_hourly_historical_year.nil? && !electricity_emissions_hourly_historical_year.empty?
+                if hourly_historical_years.include? electricity_emissions_hourly_historical_year
+                  OpenStudio::Extension.set_measure_argument(osw, 'add_ems_emissions_reporting', 'hourly_historical_year', electricity_emissions_hourly_historical_year)
+                else
+                  @@logger.error(" '#{electricity_emissions_hourly_historical_year}' is not valid option for electricity_emissions_hourly_historical_year. Please choose an input from #{hourly_historical_years}")
+                end
               end
-              
-              #annual_historical_year'
-              if !emissions_annual_historical_year.nil? && !emissions_annual_historical_year.empty?
-                OpenStudio::Extension.set_measure_argument(osw, 'add_ems_emissions_reporting', 'annual_historical_year', emissions_annual_historical_year)                
-              end     
 
-            rescue
+              # annual_historical_year
+              if !electricity_emissions_annual_historical_year.nil? && !electricity_emissions_annual_historical_year.empty?
+                if annual_historical_years.include? electricity_emissions_annual_historical_year
+                  OpenStudio::Extension.set_measure_argument(osw, 'add_ems_emissions_reporting', 'annual_historical_year', electricity_emissions_annual_historical_year)
+                else
+                  @@logger.error("'#{electricity_emissions_annual_historical_year}' is not valid option for electricity_emissions_annual_historical_year. Please choose an input from #{annual_historical_years}")
+                end
+              end
+            rescue StandardError
             end
 
           end