update paper code, new release

Author: junzis

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/plots.ipynb

@@ -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
+```