update paper code, new release

paper/contrail.py

@@ -18,19 +18,29 @@
 
 
 def saturation_pressure_over_water(temperature):
- return pressure_steam * 10 ** (
- -7.90298 * (temperature_steam / temperature - 1)
- + 5.02808 * np.log10(temperature_steam / temperature)
- - 1.3816e-7 * (10 ** (11.344 * (1 - temperature / temperature_steam)) - 1)
- + 8.1328e-3 * (10 ** (-3.49149 * (temperature_steam / temperature - 1)) - 1)
+ # Murphy and Koop 2005
+ return np.exp(
+ 54.842763
+ - 6763.22 / temperature
+ - 4.210 * np.log(temperature)
+ + 0.000367 * temperature
+ + np.tanh(0.0415 * (temperature - 218.8))
+ * (
+ 53.878
+ - 1331.22 / temperature
+ - 9.44523 * np.log(temperature)
+ + 0.014025 * temperature
+ )
  )
 
 
 def saturation_pressure_over_ice(temperature):
- return pressure_ice_point * 10 ** (
- -9.09718 * (temperature_ice_point / temperature - 1)
- - 3.56654 * np.log10(temperature_ice_point / temperature)
- + 0.876793 * (1 - temperature / temperature_ice_point)
+ # Murphy and Koop 2005
+ return np.exp(
+ 9.550426
+ - 5723.265 / temperature
+ + 3.53068 * np.log(temperature)
+ - 0.00728332 * temperature
  )
 
 

paper/readme.md

@@ -2,6 +2,8 @@
 
 This repository accompanies the paper "Fast Contrail Estimation with OpenSky Data", providing Python notebooks, scripts, and datasets integral to this paper and its findings.
 
+The paper can be accessed at: https://journals.open.tudelft.nl/joas/article/view/7264
+
 ## Python Notebooks and Scripts
 
 - `plot.ipynb`: This Jupyter notebook contains the code to replicate the figures presented in the paper.
@@ -17,3 +19,14 @@ This repository accompanies the paper "Fast Contrail Estimation with OpenSky Dat
 ## Installation Note
 
 To reproduce the results presented, it is necessary to install the `fastmeteo` Python library. Installation instructions and additional information are available at: https://github.com/junzis/fastmeteo
+
+You will also need `traffic` library for trajectory visualization. More details at: https://traffic-viz.github.io/
+
+Sample installation:
+
+```
+mamba create -n fastcontrail python=3.11 -c conda-forge
+mamba activate fastcontrail
+mamba install traffic
+pip install git+https://github.com/junzis/fastmeteo
+```