Merge pull request #18 from dlr-pa/feature/cont

Integration of AirClim 2.1 contrail module

.gitignore

@@ -26,6 +26,7 @@ example/repository/*.nc
 !example/repository/ch4_bg.nc
 !example/repository/n2o_bg.nc
 !example/repository/resp_RF.nc
+!example/repository/resp_cont.nc
 tests/repository/*.*
 tests/repository/cache
 # exceptions

README.md

@@ -155,7 +155,9 @@ Major planned software releases and milestones for the project planning are:
 ## References
 - Grewe, V., & Stenke, A. (2008). AirClim: an efficient tool for climate evaluation of aircraft technology. Atmospheric Chemistry and Physics, 8(16), 4621-4639.
 - Dahlmann, K. (2011). A method for the efficient evaluation of climate optimisation measures for air transport [Eine Methode zur effizienten Bewertung von Maßnahmen zur Klimaoptimierung des Luftverkehrs] (Doctoral dissertation, Ph. D. Thesis, Ludwig-Maximilians-Universität München, Munich).
+- Hüttenhofer, L. (2013). Parametrisierung von Kondensstreifenzirren für AirClim 2.0 (Bachelor Thesis, Ludwig-Maximilians-Universität München, Munich). 
 - Dahlmann, K., Grewe, V., Frömming, C., & Burkhardt, U. (2016). Can we reliably assess climate mitigation options for air traffic scenarios despite large uncertainties in atmospheric processes?. Transportation Research Part D: Transport and Environment, 46, 40-55.
+- Grewe, V., Bock, L., Burkhardt, U., et al. (2017). Assessing the climate impact of the AHEAD multi-fuel blended wing body. Meteorologische Zeitschrift, 26(6), 711-725.
 - Leipold, A. et al. (2021) DEPA 2050 – Development Pathways for Aviation up to 2050 (Final Report). https://elib.dlr.de/142185/
 
 

example/example.toml

@@ -3,23 +3,23 @@
 # Species considered
 [species]
 # Species defined in emission inventories
-# possible values: "CO2", "H2O", "NOx"
-inv = ["CO2", "H2O", "NOx"]
+# possible values: "CO2", "H2O", "NOx", "distance"
+inv = ["CO2", "H2O", "NOx", "distance"]
 # Assumed NOx species in emission inventory
 # possible values: "NO", "NO2"
 nox = "NO"
 # Output / response species
 # possible values: "CO2", "H2O"
-out = ["CO2", "H2O", "CH4"]
+out = ["CO2", "H2O", "cont"]
 
 # Emission inventories                                                    
 [inventories]
 dir = "repository/"
 files = [
-    "rnd_inv_2020.nc",
-    "rnd_inv_2030.nc",
-    "rnd_inv_2040.nc",
-    "rnd_inv_2050.nc",
+    "emi_inv_2020.nc",
+#     "emi_inv_2030.nc",
+#     "emi_inv_2040.nc",
+    "emi_inv_2050.nc",
     #    "rnd_inv_2060.nc",
     #    "rnd_inv_2070.nc",
     #    "rnd_inv_2080.nc",
@@ -77,6 +77,12 @@ CH4.response_grid = "2D"
 CH4.tau.file = "repository/resp_ch4.nc"
 CH4.rf.method = "Etminan_2016"
 
+cont.response_grid = "cont"
+cont.resp.file = "repository/resp_cont.nc"
+cont.G_comp = 0.04  # conventional: 0.04; hydrogen: 0.12 
+cont.eff_fac = 1.0  # efficiency factor compared to kerosene (1.0)
+cont.PMrel = 1.0  # relative PM emissions compared to kerosene (1.0)
+
 # Temperature options
 [temperature]
 # valid methods: "Boucher&Reddy"
@@ -88,6 +94,7 @@ CO2.lambda = 0.73
 H2O.efficacy = 1.14
 O3.efficacy = 1.37
 CH4.efficacy = 1.14
+cont.efficacy = 0.59
 
 # Climate metrics options
 [metrics]

openairclim/__init__.py

@@ -14,6 +14,7 @@
 from openairclim.calc_response import *  # noqa: F401, F403
 from openairclim.calc_co2 import *  # noqa: F401, F403
 from openairclim.calc_ch4 import *  # noqa: F401, F403
+from openairclim.calc_cont import *  # noqa: F401, F403
 from openairclim.calc_dt import *  # noqa: F401, F403
 from openairclim.calc_metric import *  # noqa: F401, F403
 from openairclim.uncertainties import *  # noqa: F401, F403