Add fictitious load to account for land transport ICE emissions

Recently ICE vehicles were added, but the emissions were not accounted
for. Now, like aviation emissions, they are added as a negative load
to the "co2 atmosphere" bus.

config.default.yaml

@@ -370,6 +370,7 @@ plotting:
     "process emissions to atmosphere" : "#888888"
     "process emissions" : "#222222"
     "oil emissions" : "#666666"
+    "land transport oil emissions" : "#666666"
     "land transport fuel cell" : "#AAAAAA"
     "biogas" : "#800000"
     "solid biomass" : "#DAA520"

scripts/prepare_sector_network.py

@@ -1157,6 +1157,13 @@ def add_land_transport(network):
                      carrier="land transport oil",
                      p_set=ice_share/options['transport_internal_combustion_efficiency']*transport[nodes])
 
+        co2 = ice_share/options['transport_internal_combustion_efficiency']*transport[nodes].sum().sum()/8760.*costs.at["oil",'CO2 intensity']
+
+        network.madd("Load",
+                     ["land transport oil emissions"],
+                     bus="co2 atmosphere",
+                     carrier="land transport oil emissions",
+                     p_set=-co2)
 
 
 def add_heat(network):
@@ -1906,7 +1913,7 @@ if __name__ == "__main__":
                         floor_area = "resources/floor_area_{network}_s{simpl}_{clusters}.csv"
             ),
             output=['results/version-cb48be3/prenetworks/{network}_s{simpl}_{clusters}_lv{lv}__{sector_opts}_{planning_horizons}.nc']
-    
+
         )
         import yaml
         with open('config.yaml', encoding='utf8') as f: