lodging 0.0.12 → 0.0.13

data/features/lodging_impacts.feature

@@ -0,0 +1,74 @@
+Feature: Lodging Emissions Calculations
+  The lodging model should generate correct emission calculations
+
+  Background:
+    Given a Lodging
+
+  Scenario: Calculations starting from nothing
+    Given a lodging has nothing
+    When impacts are calculated
+    Then the amount of "carbon" should be within "0.005" of "28.62"
+
+  Scenario Outline: Calculations starting from date
+    Given it has "date" of "<date>"
+    And it has "timeframe" of "<timeframe>"
+    When impacts are calculated
+    Then the amount of "carbon" should be within "0.005" of "<carbon>"
+    Examples:
+      | date       | timeframe             | carbon |
+      | 2011-01-15 | 2011-01-01/2012-01-01 |  28.62 |
+      | 2012-01-15 | 2011-01-01/2012-01-01 |   0.0  |
+
+  Scenario: Calculations starting from rooms and duration
+    Given it has "rooms" of "2"
+    And it has "duration" of "172800"
+    When impacts are calculated
+    Then the amount of "carbon" should be within "0.005" of "114.47"
+
+  Scenario Outline: Calculations from fuzzy inference based on country degree days
+    Given it has "rooms" of "2"
+    And it has "duration" of "172800"
+    And it has "country.iso_3166_code" of "<country>"
+    And it has "state.postal_abbreviation" of "<state>"
+    When impacts are calculated
+    Then the amount of "carbon" should be within "0.005" of "<carbon>"
+    Examples:
+      | country | state | carbon |
+      | VI      |       | 114.47 |
+      | GB      |       | 189.56 |
+      | US      |       | 131.28 |
+      |         | CA    | 131.28 |
+
+  Scenario Outline: Calculations from fuzzy inference
+    Given it has "rooms" of "2"
+    And it has "duration" of "172800"
+    And it has "property.northstar_id" of "<id>"
+    And it has "zip_code.name" of "<zip>"
+    And it has "city" of "<city>"
+    And it has "state.postal_abbreviation" of "<state>"
+    When impacts are calculated
+    Then the amount of "carbon" should be within "0.005" of "<carbon>"
+    Examples:
+      | id | zip   | city          | state | carbon | notes |
+      | 1  | 94122 |               |       |  86.85 | dd from climate divizion; fuzzy from property attributes |
+      | 1  |       | San Francisco | CA    | 122.33 | dd from country; fuzzy from property attributes |
+      | 2  | 94133 |               |       |  94.24 | dd from country; fuzzy from property attributes |
+      | 2  |       | San Francisco | CA    |  94.24 | dd from country; fuzzy from property attributes |
+      | 3  | 94014 |               |       | 128.13 | dd from country; fuzzy from property attributes |
+
+  Scenario: Calculations from all user inputs
+    Given it has "rooms" of "2"
+    And it has "duration" of "172800"
+    And it has "country.iso_3166_code" of "GB"
+    And it has "heating_degree_days" of "2700"
+    And it has "cooling_degree_days" of "100"
+    And it has "property_rooms" of "100"
+    And it has "floors" of "3"
+    And it has "construction_year" of "1993"
+    And it has "ac_coverage" of "0.5"
+    And it has "refrigerator_coverage" of "0.6"
+    And it has "hot_tubs" of "6"
+    And it has "indoor_pools" of "5"
+    And it has "outdoor_pools" of "5"
+    When impacts are calculated
+    Then the amount of "carbon" should be within "0.005" of "208.83"

data/features/support/db/fixtures/census_divisions.csv

@@ -0,0 +1,2 @@
+number,census_region_number
+9,4

data/features/support/db/fixtures/climate_divisions.csv

@@ -0,0 +1,2 @@
+name,heating_degree_days,cooling_degree_days,state_postal_abbreviation
+CA4,1350,150,CA

data/features/support/db/fixtures/commercial_building_energy_consumption_survey_responses.csv

@@ -0,0 +1,11 @@
+id,heating_degree_days,cooling_degree_days,lodging_rooms,floors,construction_year,percent_cooled,natural_gas_per_room_night,fuel_oil_per_room_night,electricity_per_room_night,district_heat_per_room_night,weighting
+1,5500,400,10,1,1967,1,0,1.9,12,0,1700
+2,1500,1200,35,2,1988,0.6,0,0,14,0,420
+3,200,2700,100,4,2000,1,1.2,0,34,0,125
+4,2300,900,2,1,1992,1,0,0,23,0,1770
+5,2300,700,1000,25,1986,1,11,0.015,120,0,135
+6,3800,200,7,4,1881,0,0,0,9.2,0,1925
+7,800,1700,125,3,2002,1,1.9,0,27,0,240
+8,3300,400,10,2,1920,0,0,6.1,11,0,1925
+9,1000,1700,14,2,1950,1,9.7,0,83,180,365
+10,2400,800,450,25,1929,1,0.74,0.00092,41,97,120

data/features/support/db/fixtures/countries.csv

@@ -0,0 +1,4 @@
+iso_3166_code,iso_3166_alpha_3_code,heating_degree_days,cooling_degree_days,electricity_emission_factor,electricity_loss_factor,lodging_occupancy_rate,lodging_natural_gas_intensity,lodging_fuel_oil_intensity,lodging_electricity_intensity,lodging_district_heat_intensity
+US,USA,2200,880,0.59000,0.06025,0.6,2.0,0.4,33.9,1.8
+GB,GBR,2800,60,0.9,0.1,0.5,3.0,0.5,60.0,0.0
+VI,VIR

data/features/support/db/fixtures/egrid_countries.csv

@@ -0,0 +1,2 @@
+name,generation,imports,consumption
+US,400000,3100,379000

data/features/support/db/fixtures/egrid_regions.csv

@@ -0,0 +1,2 @@
+name,loss_factor
+W,0.05

data/features/support/db/fixtures/egrid_subregions.csv

@@ -0,0 +1,3 @@
+abbreviation,egrid_region_name,net_generation,electricity_emission_factor
+CAMX,W,1000,0.31
+RMPA,W,1000,0.87

data/features/support/db/fixtures/fuel_years.csv

@@ -0,0 +1,2 @@
+name,fuel_name,year,energy_content,energy_content_units,co2_emission_factor,co2_emission_factor_units
+Pipeline Natural Gas 2008,Pipeline Natural Gas,2008,38,megajoules_per_cubic_metre,1.9,kilograms_per_cubic_metre

data/features/support/db/fixtures/fuels.csv

@@ -0,0 +1,4 @@
+name,energy_content,co2_emission_factor
+Residual Fuel Oil No. 6,42,3.0
+District Heat,1,0.077
+Pipeline Natural Gas

data/features/support/db/fixtures/lodging_classes.csv

@@ -0,0 +1,4 @@
+name
+Hotel
+Inn
+Motel

data/features/support/db/fixtures/lodging_properties.csv

@@ -0,0 +1,5 @@
+northstar_id,name,city,locality,postcode,country_iso_3166_alpha_3_code,lodging_rooms,floors,construction_year,ac_coverage,lodging_class_name,fridge_coverage,mini_bar_coverage,hot_tubs,pools_indoor,pools_outdoor
+1,Courtyard by Marriott,San Francisco,California,94122,USA,100,3,1993,0.5,Motel,0.5,0.6,6,6,6
+2,Hilton San Francisco,San Francisco,California,94133,USA,50,2,2002,1,Hotel,0.5,,0,0,0
+3,Pacific Inn,Daly City,California,94014,USA,25,,,,Inn,,0.5,1,1,1
+4,,,,,,,,,,,0.6,0.5,5,5,5

data/features/support/db/fixtures/states.csv

@@ -0,0 +1,2 @@
+postal_abbreviation,name,census_division_number
+CA,California,9

data/features/support/db/fixtures/zip_codes.csv

@@ -0,0 +1,4 @@
+name,state_postal_abbreviation,description,egrid_subregion_abbreviation,climate_division_name
+94122,CA,San Francisco,CAMX,CA4
+94133,CA,San Francisco
+94014,CA,Daly City

data/features/support/env.rb

@@ -3,6 +3,23 @@ Bundler.setup
 
 require 'cucumber'
 require 'cucumber/formatter/unicode'
+require 'rspec'
+require 'rspec/expectations'
+require 'cucumber/rspec/doubles'
+
+require 'data_miner'
+DataMiner.logger = Logger.new(nil)
 
 require 'sniff'
-Sniff.init File.join(File.dirname(__FILE__), '..', '..'), :earth => [:hospitality, :fuel, :locality], :cucumber => true
+
+# need to require this before initializing Sniff for LodgingProperty fixture to be loaded
+require File.expand_path('../lodging_property', __FILE__)
+
+Sniff.init File.join(File.dirname(__FILE__), '..', '..'),
+  :adapter => 'mysql2',
+  :database => 'test_lodging',
+  :username => 'root',
+  :password => 'password',
+  :earth => [:hospitality, :fuel, :locality],
+  :cucumber => true,
+  :logger => 'log/test_log.txt'

data/features/support/lodging_property.rb

@@ -0,0 +1,26 @@
+# LodgingProperty is defined in CM1 rather than Earth, so we need to define it here for testing
+class LodgingProperty < ActiveRecord::Base
+  col :northstar_id
+  col :name
+  col :city
+  col :locality # state / province / etc.
+  col :postcode # zip code / postal code / etc.
+  col :country_iso_3166_alpha_3_code
+  col :chain_name
+  col :lodging_rooms,     :type => :integer
+  col :floors,            :type => :integer
+  col :construction_year, :type => :integer
+  col :lodging_class_name
+  col :ac_coverage,       :type => :float
+  col :mini_bar_coverage, :type => :float
+  col :fridge_coverage,   :type => :float
+  col :hot_tubs,          :type => :float # float b/c fallback needs to be a float
+  col :pools_indoor,      :type => :float # float b/c fallback needs to be a float
+  col :pools_outdoor,     :type => :float # float b/c fallback needs to be a float
+  
+  # based on LodgingProperty as of 2012-02-21
+  falls_back_on :fridge_coverage => 0.6,
+                :hot_tubs        => 0.3,
+                :pools_indoor    => 0.3,
+                :pools_outdoor   => 0.6
+end

data/{lib/test_support → features/support}/lodging_record.rb

@@ -1,9 +1,6 @@
-require 'active_record'
-require 'falls_back_on'
 require 'lodging'
-require 'sniff'
 
 class LodgingRecord < ActiveRecord::Base
-  include Sniff::Emitter
+  include BrighterPlanet::Emitter
   include BrighterPlanet::Lodging
 end

data/lib/lodging.rb

@@ -3,5 +3,6 @@ require 'emitter'
 module BrighterPlanet
   module Lodging
     extend BrighterPlanet::Emitter
+    scope 'The lodging emission estimate is the anthropogenic emissions from lodging room energy use. It includes CO2 emissions from direct fuel combustion and indirect fuel combustion to generate purchased electricity.'
   end
 end

data/lib/lodging/cbecs.rb

@@ -0,0 +1,31 @@
+require 'fuzzy_infer'
+
+class CommercialBuildingEnergyConsumptionSurveyResponse < ActiveRecord::Base
+  fuzzy_infer :target     => [:natural_gas_per_room_night, :fuel_oil_per_room_night, :electricity_per_room_night, :district_heat_per_room_night], # list of columns that this model is designed to infer
+              :basis      => [:heating_degree_days, :cooling_degree_days, :lodging_rooms, :floors, :construction_year, :percent_cooled],          # list of columns that are believed to affect energy use (aka MU)
+              :sigma      => "(STDDEV_SAMP(:column)/5)+(ABS(AVG(:column)-:value)/3)",                                                                  # empirically determined formula (SQL!) that captures the desired sample size once all the weights are compiled, across the full range of possible mu values
+              :membership => :energy_use_membership,                                                                                              # name of instance method to be called on kernel
+              :weight     => :weighting                                                                                                           # (optional) a pre-existing row weighting, if any, provided by the dataset authors
+  
+  # empirically determined formula that minimizes variance between real and predicted energy use
+  # SQL! - :heating_degree_days_n_w will be replaced with, for example, `tmp_table_9301293.hdd_normalized_weight`
+  def energy_use_membership(basis)
+    keys = basis.keys
+    
+    unless keys.include?(:heating_degree_days) and keys.include?(:cooling_degree_days)
+      raise ArgumentError, "[lodging] Must provide at least :heating_degree_days and :cooling_degree_days"
+    end
+    
+    formula = ['(POW(:heating_degree_days_n_w, 0.8) + POW(:cooling_degree_days_n_w, 0.8))']
+    if keys.include?(:lodging_rooms) and keys.include?(:floors)
+      formula << '(POW(:lodging_rooms_n_w, 0.8) + POW(:floors_n_w, 0.8))'
+    elsif keys.include? :lodging_rooms
+      formula << 'POW(:lodging_rooms_n_w, 0.8)'
+    elsif keys.include? :floors
+      formula << 'POW(:floors_n_w, 0.8)'
+    end
+    formula << 'POW(:percent_cooled_n_w, 0.8)' if keys.include?(:percent_cooled)
+    formula << 'POW(:construction_year_n_w, 0.8)' if keys.include?(:construction_year)
+    formula.join(' * ')
+  end
+end

data/lib/lodging/characterization.rb

@@ -3,11 +3,27 @@ module BrighterPlanet
     module Characterization
       def self.included(base)
         base.characterize do
-          has :lodging_class
+          has :date
+          has :rooms
+          has :duration, :measures => :time
           has :zip_code
+          has :city
           has :state
-          has :rooms
-          has :duration
+          has :country
+          has :lodging_class
+          has :heating_degree_days
+          has :cooling_degree_days
+          has :property do |prop|
+            prop.name
+          end
+          has :property_rooms
+          has :floors
+          has :construction_year
+          has :ac_coverage
+          has :refrigerator_coverage
+          has :hot_tubs
+          has :outdoor_pools
+          has :indoor_pools
         end
       end
     end

data/lib/lodging/data.rb

@@ -2,16 +2,30 @@ module BrighterPlanet
   module Lodging
     module Data
       def self.included(base)
+        base.col :date,     :type => :date
+        base.col :rooms,    :type => :integer
+        base.col :duration, :type => :integer
+        base.col :zip_code_name
+        base.col :city
+        base.col :state_postal_abbreviation
+        base.col :country_iso_3166_code
+        base.col :lodging_class_name
+        base.col :heating_degree_days, :type => :float
+        base.col :cooling_degree_days, :type => :float
+        base.col :property_northstar_id
+        base.col :property_rooms,        :type => :integer
+        base.col :floors,                :type => :integer
+        base.col :construction_year,     :type => :integer
+        base.col :ac_coverage,           :type => :float
+        base.col :refrigerator_coverage, :type => :float
+        base.col :hot_tubs,              :type => :integer
+        base.col :outdoor_pools,         :type => :integer
+        base.col :indoor_pools,          :type => :integer
+        
         base.data_miner do
-          schema do
-            string  'lodging_class_name'
-            string  'zip_code_name'
-            string  'state_postal_abbreviation'
-            integer 'rooms'
-            integer 'duration'
+          process 'pull orphans' do
+            Fuel.run_data_miner!
           end
-
-          process :run_data_miner_on_belongs_to_associations
         end
       end
     end

data/lib/lodging/impact_model.rb

@@ -0,0 +1,544 @@
+# Copyright © 2010 Brighter Planet.
+# See LICENSE for details.
+# Contact Brighter Planet for dual-license arrangements.
+
+## Lodging impact model
+# This model is used by the [Brighter Planet](http://brighterplanet.com) [CM1 web service](http://carbon.brighterplanet.com) to calculate the impacts of a lodging (e.g. a hotel stay), such as energy use, greenhouse gas emissions, and water use.
+
+##### Timeframe
+# The model calculates impacts that occured during a particular time period (`timeframe`).
+# For example if the `timeframe` is February 2010, a lodging that occurred (`date`) on February 15, 2010 will have impacts, but a lodging that occurred on January 31, 2010 will have zero impacts.
+#
+# The default `timeframe` is the current calendar year.
+
+##### Calculations
+# The final impacts are the result of the calculations below. These are performed in reverse order, starting with the last calculation listed and finishing with the greenhouse gas emissions calculation.
+#
+# Each calculation listing shows:
+#
+# * value returned (*units of measurement*)
+# * description of the value
+# * calculation methods, listed from most to least preferred
+#
+# Some methods use `values` returned by prior calculations. If any of these `values` are unknown the method is skipped.
+# If all the methods for a calculation are skipped, the value the calculation would return is unknown.
+
+##### Standard compliance
+# When compliance with a particular standard is requested, all methods that do not comply with that standard are ignored.
+# Thus any `values` a method needs will have been calculated using a compliant method or will be unknown.
+# To see which standards a method complies with, look at the `:complies =>` section of the code in the right column.
+#
+# Client input complies with all standards.
+
+##### Collaboration
+# Contributions to this impact model are actively encouraged and warmly welcomed. This library includes a comprehensive test suite to ensure that your changes do not cause regressions. All changes should include test coverage for new functionality. Please see [sniff](https://github.com/brighterplanet/sniff#readme), our emitter testing framework, for more information.
+require 'lodging/cbecs'
+
+module BrighterPlanet
+  module Lodging
+    module ImpactModel
+      def self.included(base)
+        base.decide :impact, :with => :characteristics do
+          # * * *
+          
+          #### Carbon (*kg CO<sub>2</sub>e*)
+          # *The lodging's total anthropogenic greenhouse gas emissions during `timeframe`.*
+          committee :carbon do
+            # Multiply `natural gas use` (*m<sup>3</sup>*) by natural gas' emission factor (*kg CO<sub>2</sub> / m<sup>3</sup>*) to give *kg CO<sub>2</sub>.*
+            # Multiply `fuel oil use` (*l*) by fuel oil's emission factor (*kg CO<sub>2</sub> / l*) to give *kg CO<sub>2</sub>.*
+            # Multiply `electricity use` (*kWh*) by `electricity emission factor` (*kg CO<sub>2</sub> / kWh*) to give *kg CO<sub>2</sub>.*
+            # Multiply `district heat use` (*MJ*) by district heat's emission factor` (*kg CO<sub>2</sub> / MJ*) to give *kg CO<sub>2</sub>.*
+            # Sum to give *kg CO<sub>2</sub>e.*
+            quorum 'from natural gas use, fuel oil use, electricity use, district heat use, and electricity emission factor', :needs => [:natural_gas_use, :fuel_oil_use, :electricity_use, :district_heat_use, :electricity_emission_factor],
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:natural_gas_use] * Fuel.find_by_name('Pipeline Natural Gas').co2_emission_factor +
+                characteristics[:fuel_oil_use] * Fuel.find_by_name('Residual Fuel Oil No. 6').co2_emission_factor +
+                characteristics[:electricity_use] * characteristics[:electricity_emission_factor] +
+                characteristics[:district_heat_use] * Fuel.find_by_name('District Heat').co2_emission_factor
+            end
+          end
+          
+          #### Electricity emission factor (*kg CO<sub>2</sub>e / kWh*)
+          # *A greenhouse gas emission factor for electricity used by the lodging.*
+          committee :electricity_emission_factor do
+            # Multiply the `eGRID subregion` electricity emission factors by 1 minus the the subregion's eGRID region loss factor (to account for transmission and distribution losses) to give *kg CO<sub>2</sub>e / kWh*.
+            quorum 'from eGRID subregion', :needs => :egrid_subregion,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:egrid_subregion].electricity_emission_factor / (1 - characteristics[:egrid_subregion].egrid_region.loss_factor)
+            end
+            
+            # Otherwise look up the `country` electricity emission factor (*kg CO<sub>2</sub>e / kWh*).
+            quorum 'from country', :needs => :country,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                if (emission_factor = characteristics[:country].electricity_emission_factor).present? and (loss_factor = characteristics[:country].electricity_loss_factor).present?
+                  emission_factor / (1 - loss_factor)
+                end
+            end
+            
+            # Otherwise use a global average electricity emission factor (*kg CO<sub>2</sub>e / kWh*).
+            quorum 'default',
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do
+                Country.fallback.electricity_emission_factor / (1 - Country.fallback.electricity_loss_factor)
+            end
+          end
+          
+          #### Natural gas use (*m<sup>3</sup>*)
+          # *The lodging's natural gas use during `timeframe`.*
+          committee :natural_gas_use do
+            # Multiply `room nights` (*room-nights*) by `natural gas intensity` (*m<sup>3</sup> / occupied room-night*) to give *m<sup>3</sup>*.
+            quorum 'from adjusted fuel intensities and room nights', :needs => [:adjusted_fuel_intensities, :room_nights],
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:room_nights] * characteristics[:adjusted_fuel_intensities][:natural_gas]
+            end
+          end
+          
+          #### Fuel oil use (*l*)
+          # *The lodging's fuel oil use during `timeframe`.*
+          committee :fuel_oil_use do
+            # Multiply `room nights` (*room-nights*) by `fuel oil intensity` (*l / occupied room-night*) to give *l*.
+            quorum 'from adjusted fuel intensities and room nights', :needs => [:adjusted_fuel_intensities, :room_nights],
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:room_nights] * characteristics[:adjusted_fuel_intensities][:fuel_oil]
+            end
+          end
+          
+          #### Electricity use (*kWh*)
+          # *The lodging's electricity use during `timeframe`.*
+          committee :electricity_use do
+            # Multiply `room nights` (*room-nights*) by `electricity intensity` (*kWh / occupied room-night*) to give *kWh*.
+            quorum 'from adjusted fuel intensities and room nights', :needs => [:adjusted_fuel_intensities, :room_nights],
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:room_nights] * characteristics[:adjusted_fuel_intensities][:electricity]
+            end
+          end
+          
+          #### District heat use (*MJ*)
+          # *The lodging's district heat use during `timeframe`.*
+          committee :district_heat_use do
+            # Multiply `room nights` (*room-nights*) by `district heat intensity` (*MJ / occupied room-night*) to give *MJ*.
+            quorum 'from adjusted fuel intensities and room nights', :needs => [:adjusted_fuel_intensities, :room_nights],
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:room_nights] * characteristics[:adjusted_fuel_intensities][:district_heat]
+            end
+          end
+          
+          #### Adjusted fuel intensities (*various*)
+          # *The lodging's use per occupied room night of a variety of fuels, adjusted by number of pools, mini-fridges, etc.*
+          #
+          # - Natural gas intensity: *m<sup>3</sup> / occupied room-night*
+          # - Fuel oil intensity: *l / occupied room-night*
+          # - Electricity intensity: *kWh / occupied room-night*
+          # - District heat intensity: *MJ / occupied room-night*
+          committee :adjusted_fuel_intensities do
+            # Adjust `fuel intensities` based on any amenity adjustments:
+            quorum 'from fuel intensities and amenity adjustments',
+              :needs => :fuel_intensities, :appreciates => [:indoor_pool_adjustment, :outdoor_pool_adjustment, :refrigerator_adjustment, :hot_tub_adjustment],
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                gas = Fuel.find_by_name('Pipeline Natural Gas')
+                oil = Fuel.find_by_name('Residual Fuel Oil No. 6')
+                
+                # Duplicate the current intensities so we don't overwrite them.
+                intensities = characteristics[:fuel_intensities].dup
+                
+                # TODO combine pool adjustments - so if we're below average in indoor but above in outdoor we don't accidentally subtract, get less than zero, discard some energy use, and then add energy back for outdoor pools
+                
+                # Cycle through each adjustment...
+                characteristics.each do |characteristic, value|
+                  unless characteristic == :fuel_intensities
+                    value.each do |fuel, adjustment|
+                      if fuel == :pool_energy
+                        # If the adjustment relates to pool energy, convert natural gas and fuel oil inensities from physical units to energy units.
+                        gas_energy = intensities[:natural_gas] * gas.energy_content
+                        oil_energy = intensities[:fuel_oil]    * oil.energy_content
+                        
+                        # Apply the adjustment to natural gas intensity.
+                        # If adjusted natural gas intensity reaches zero, apply any remaining adjustment to fuel oil intensity.
+                        # If adjusted fuel oil intensity also reaches zero, discard any remaining adjustment.
+                        if (gas_energy += adjustment) < 0.0
+                          if (oil_energy += gas_energy) < 0.0
+                            oil_energy = 0.0
+                          end
+                          gas_energy = 0
+                        end
+                        
+                        # Convert adjusted natural gas and fuel oil intensities from energy units back to physical units.
+                        intensities[:natural_gas] = gas_energy / gas.energy_content
+                        intensities[:fuel_oil]    = oil_energy / oil.energy_content
+                      # Otherwise apply the adjustment to the appropriate fuel intensity. If adjusted fuel intensity reaches zero, discard any remaining adjustment.
+                      elsif (intensities[fuel] += adjustment) < 0.0
+                        intensities[fuel] = 0.0
+                      end
+                    end
+                  end
+                end
+                intensities
+            end
+          end
+          
+          #### Outdoor pools adjustment (*MJ / occupied room-night*)
+          # *Adjusts the natural gas intensity based on the number of outdoor pools.*
+          committee :outdoor_pool_adjustment do
+            # Assume outdoor pool energy intensity of 329,917 *BTU / night* per [Energy Star](http://www.energystar.gov/ia/business/evaluate_performance/swimming_pool_tech_desc.pdf).
+            # Calculate the difference between `outdoor pools` and average outdoor pools.
+            # Multiply the difference by outdoor pool energy intensity (*MJ / night*) and divide by `property rooms` and `occupancy rate` to give *MJ / occupied room-night*.
+            quorum 'from outdoor pools, property rooms, and occupancy rate', :needs => [:outdoor_pools, :property_rooms, :occupancy_rate],
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                difference = characteristics[:outdoor_pools] - LodgingProperty.fallback.pools_outdoor
+                { :pool_energy => difference * 329_917.btus.to(:megajoules) / characteristics[:property_rooms] / characteristics[:occupancy_rate] }
+            end
+          end
+          
+          #### Indoor pools adjustment (*MJ / occupied room-night*)
+          # *Adjusts the natural gas intensity based on the number of indoor pools.*
+          committee :indoor_pool_adjustment do
+            # Assume indoor pool energy intensity of 2,770,942 *BTU / night* per [Energy Star](http://www.energystar.gov/ia/business/evaluate_performance/swimming_pool_tech_desc.pdf).
+            # Calculate the difference between `indoor pools` and average indoor pools.
+            # Multiply the difference by indoor pool energy intensity (*MJ / night*) and divide by `property rooms` and `occupancy rate` to give *MJ / occupied room-night*.
+            quorum 'from indoor pools, property rooms, and occupancy rate', :needs => [:indoor_pools, :property_rooms, :occupancy_rate],
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                difference = characteristics[:indoor_pools] - LodgingProperty.fallback.pools_indoor
+                { :pool_energy => difference * 2_770_942.btus.to(:megajoules) / characteristics[:property_rooms] / characteristics[:occupancy_rate] }
+            end
+          end
+          
+          #### Hot tub adjustment (*kWh / occupied room-night*)
+          # *Adjusts the electricity intensity based on the number of hot tubs.*
+          committee :hot_tub_adjustment do
+            # Calculate the difference between `hot tubs` and average hot tubs.
+            # Assume hot tub electricity intensity of 6.3 *kWh / night* per [LBL residential energy data sourcebook, p128](http://enduse.lbl.gov/info/LBNL-40297.pdf).
+            # Multiply the difference by the hot tub electricity intensity (*kWh / night*) and divide by `property rooms` and `occupancy rate` to give *kWh / occupied room-night*.
+            quorum 'from hot tubs, property rooms, and occupancy rate', :needs => [:hot_tubs, :property_rooms, :occupancy_rate],
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                difference = characteristics[:hot_tubs] - LodgingProperty.fallback.hot_tubs
+                { :electricity => difference * 6.3 / characteristics[:property_rooms] / characteristics[:occupancy_rate] }
+            end
+          end
+          
+          #### Refrigerator adjustment (*kWh / occupied room-night*)
+          # *Adjusts the electricity intensity based on refrigerator coverage.*
+          committee :refrigerator_adjustment do
+            # Calculate the difference between `refrigerator coverage` (*refrigerators / room*) and average refrigerator coverage (*refrigerators / room*).
+            # Assume an auto-defrost compact fridge electricity intensity of 1.18 *kWh / refrigerator night* per [Energy Star](http://www.energystar.gov/ia/business/bulk_purchasing/bpsavings_calc/Bulk_Purchasing_CompactRefrig_Sav_Calc.xls).
+            # Multiply the difference (*refrigerators / room*) by the refrigerator electricity intensity (*kWh / refrigerator night*) and divide by `occupancy rate` to give *kWh / occupied room-night*.
+            quorum 'from refrigerator coverage and occupancy rate', :needs => [:refrigerator_coverage, :occupancy_rate],
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                difference = characteristics[:refrigerator_coverage] - LodgingProperty.fallback.fridge_coverage
+                { :electricity => (difference * 1.18 / characteristics[:occupancy_rate]) }
+            end
+          end
+          
+          #### Fuel intensities (*various*)
+          # *The lodging's use per occupied room night of a variety of fuels.*
+          #
+          # - Natural gas intensity: *m<sup>3</sup> / occupied room-night*
+          # - Fuel oil intensity: *l / occupied room-night*
+          # - Electricity intensity: *kWh / occupied room-night*
+          # - District heat intensity: *MJ / occupied room-night*
+          committee :fuel_intensities do
+            # If we know `heating degree days` and `cooling degree days`, calculate fuel intensities from [CBECS 2003](http://data.brighterplanet.com/commercial_building_energy_consumption_survey_responses) data using fuzzy inference and adjust the inferred values based on `occupancy rate`.
+            quorum 'from degree days, occupancy rate, and user inputs', :needs => [:heating_degree_days, :cooling_degree_days, :occupancy_rate], :appreciates => [:property_rooms, :floors, :construction_year, :ac_coverage],
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                inputs = characteristics.to_hash.inject({}) do |memo, (characteristic, value)|
+                  case characteristic
+                  when :property_rooms
+                    memo[:lodging_rooms] = value
+                  when :ac_coverage
+                    memo[:percent_cooled] = value
+                  when :occupancy_rate
+                    # Don't include `occupancy rate` in inputs to fuzzy inference
+                  else
+                    memo[characteristic] = value
+                  end
+                  memo
+                end
+                
+                kernel = CommercialBuildingEnergyConsumptionSurveyResponse.new(inputs)
+                n, f, e, d = kernel.fuzzy_infer(:natural_gas_per_room_night, :fuel_oil_per_room_night, :electricity_per_room_night, :district_heat_per_room_night)
+                {
+                  :natural_gas   => n / characteristics[:occupancy_rate],
+                  :fuel_oil      => f / characteristics[:occupancy_rate],
+                  :electricity   => e / characteristics[:occupancy_rate],
+                  :district_heat => d / characteristics[:occupancy_rate]
+                }
+            end
+            
+            # Otherwise use global averages.
+            quorum 'default',
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                intensities = {
+                  :natural_gas   => Country.fallback.lodging_natural_gas_intensity,
+                  :fuel_oil      => Country.fallback.lodging_fuel_oil_intensity,
+                  :electricity   => Country.fallback.lodging_electricity_intensity,
+                  :district_heat => Country.fallback.lodging_district_heat_intensity,
+                }
+            end
+          end
+          
+          #### Occupancy rate
+          # *The percent of the proprety's rooms that are occupied on an average night.*
+          committee :occupancy_rate do
+            # Look up the `country` average lodging occupancy rate.
+            quorum 'from country', :needs => :country,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:country].lodging_occupancy_rate
+            end
+            
+            # Otherwise use a global average.
+            quorum 'default', :complies => [:ghg_protocol_scope_3, :iso, :tcr] do
+              Country.fallback.lodging_occupancy_rate
+            end
+          end
+          
+          #### Outdoor pools
+          # *The number of outdoor pools.*
+          committee :outdoor_pools do
+            # Use client input, if available.
+            
+            # Otherwise look up the `property` number of outdoor pools, but set a ceiling of 5 outdoor pools.
+            quorum 'from property', :needs => :property,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                unless characteristics[:property].pools_outdoor.nil?
+                  [characteristics[:property].pools_outdoor.to_f, 5].min
+                end
+            end
+          end
+          
+          #### Indoor pools
+          # *The number of indoor pools.*
+          committee :indoor_pools do
+            # Use client input, if available.
+            
+            # Otherwise look up the `property` number of indoor pools, but set a ceiling of 5 indoor pools.
+            quorum 'from property', :needs => :property,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                unless characteristics[:property].pools_indoor.nil?
+                  [characteristics[:property].pools_indoor.to_f, 5].min
+                end
+            end
+          end
+          
+          #### Hot tubs
+          # *The number of hot tubs.*
+          committee :hot_tubs do
+            # Use client input, if available.
+            
+            # Otherwise look up the `property` number of hot tubs.
+            quorum 'from property', :needs => :property,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:property].hot_tubs
+            end
+          end
+          
+          #### Refrigerator coverage
+          # *The percentage of property rooms that have refrigerators.*
+          committee :refrigerator_coverage do
+            # Use client input, if available.
+            # Otherwise take whichever is greater of the `property` refrigerator coverage and mini bar coverage.
+            quorum 'from property', :needs => :property,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                if characteristics[:property].mini_bar_coverage || characteristics[:property].fridge_coverage
+                  [characteristics[:property].mini_bar_coverage.to_f, characteristics[:property].fridge_coverage.to_f].max
+                end
+            end
+          end
+          
+          #### A/C coverage
+          # *The percentage of property rooms that are air conditioned.*
+          committee :ac_coverage do
+            # Use client input, if available.
+            
+            # Otherwise look up the `property` A/C coverage.
+            quorum 'from property', :needs => :property,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:property].ac_coverage
+            end
+          end
+          
+          #### Construction year
+          # *The year the property was built.*
+          committee :construction_year do
+            # Use client input, if available.
+            
+            # Otherwise look up the year the `property` was built.
+            quorum 'from property', :needs => :property,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:property].construction_year
+            end
+          end
+          
+          #### Floors
+          # *The number of floors.*
+          committee :floors do
+            # Use client input, if available.
+            
+            # Otherwise look up the `property` number of floors.
+            quorum 'from property', :needs => :property,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:property].floors
+            end
+          end
+          
+          #### Property rooms
+          # *The number of guest rooms in the property.*
+          committee :property_rooms do
+            # Use client input, if available.
+            
+            # Otherwise look up the `property` number of rooms.
+            quorum 'from property', :needs => :property,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:property].lodging_rooms
+            end
+          end
+          
+          #### Property
+          # *The property where the stay occurred.*
+          #
+          # Use client input, if available
+          
+          #### Heating degree days
+          # *The average number of annual heating degree days (base 18°C) at the lodging's location.*
+          committee :heating_degree_days do
+            # Use client input, if available
+            
+            # Otherwise look up the `climate division` heating degree days.
+            quorum 'from climate division', :needs => :climate_division,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:climate_division].heating_degree_days
+            end
+            
+            # Otherwise look up the `country` heating degree days.
+            quorum 'from country', :needs => :country,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:country].heating_degree_days
+            end
+          end
+          
+          #### Cooling degree days
+          # *The average number of annual cooling degree days (base 18°C) at the lodging's location.*
+          committee :cooling_degree_days do
+            # Use client input, if available
+            
+            # Otherwise look up the `climate division` cooling degree days.
+            quorum 'from climate division', :needs => :climate_division,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:climate_division].cooling_degree_days
+            end
+            
+            # Otherwise look up the `country` cooling degree days.
+            quorum 'from country', :needs => :country,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:country].cooling_degree_days
+            end
+          end
+          
+          #### eGRID subregion
+          # *The [eGRID subregion](http://data.brighterplanet.com/egrid_subregions).*
+          committee :egrid_subregion do
+            # Look up the `zip code` eGRID subregion.
+            quorum 'from zip code', :needs => :zip_code,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:zip_code].egrid_subregion
+            end
+          end
+          
+          #### Country
+          # *The [country](http://data.brighterplanet.com/countries).*
+          committee :country do
+            # Use client input, if available.
+            
+            # Otherwise if state is defined then the country is the United States.
+            quorum 'from state', :needs => :state,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                Country.united_states
+            end
+          end
+          
+          #### State
+          # *The [US state](http://data.brighterplanet.com/states).*
+          committee :state do
+            # Use client input, if available.
+            
+            # Otherwise look up the `zip code` state.
+            quorum 'from zip code', :needs => :zip_code,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:zip_code].state
+            end
+          end
+          
+          #### City
+          # *The city.*
+          committee :city do
+            # Use client input, if available.
+            
+            # Otherwise look up the `zip code` description.
+            quorum 'from zip code', :needs => :zip_code,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:zip_code].description
+            end
+          end
+          
+          #### Climate division
+          # *The [US climate division](http://data.brighterplanet.com/climate_divisions).*
+          committee :climate_division do
+            # Look up the `zip code` climate division.
+            quorum 'from zip code', :needs => :zip_code,
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics|
+                characteristics[:zip_code].climate_division
+            end
+          end
+          
+          #### Zip code
+          # *The [US zip code](http://data.brighterplanet.com/zip_codes).*
+          #
+          # Use client input, if available.
+          
+          #### Room nights (*room-nights*)
+          # *The stay's room-nights that occurred during `timeframe`.*
+          committee :room_nights do
+            # If `date` falls within `timeframe`, divide `duration` (*seconds*) by 86,400 (*seconds / night*) and multiply by `rooms` to give *room-nights*.
+            # Otherwise `room nights` is zero.
+            quorum 'from rooms, duration, date, and timeframe', :needs => [:rooms, :duration, :date],
+              :complies => [:ghg_protocol_scope_3, :iso, :tcr] do |characteristics, timeframe|
+                date = characteristics[:date].is_a?(Date) ? characteristics[:date] : Date.parse(characteristics[:date].to_s)
+                timeframe.include?(date) ? characteristics[:duration] / 86400.0 * characteristics[:rooms] : 0
+            end
+          end
+          
+          #### Duration (*seconds*)
+          # *The stay's duration.*
+          committee :duration do
+            # Use client input, if available.
+            
+            # Otherwise use 86400 *seconds* (1 night).
+            quorum 'default' do
+              86400.0
+            end
+          end
+          
+          #### Rooms
+          # *The number of rooms used.*
+          committee :rooms do
+            # Use client input, if available.
+            
+            # Otherwise use 1 *room*.
+            quorum 'default' do
+              1
+            end
+          end
+          
+          #### Date (*date*)
+          # *The day the stay occurred.*
+          committee :date do
+            # Use client input, if available.
+            
+            # Otherwise use the first day of `timeframe`.
+            quorum 'from timeframe',
+              :complies => [:ghg_protocol_scope_3, :iso] do |characteristics, timeframe|
+                timeframe.from
+            end
+          end
+        end
+      end
+    end
+  end
+end