Plotting script radiosounding

.github/workflows/test.yml

@@ -15,7 +15,7 @@ jobs:
       - name: Setup Python
         uses: actions/setup-python@v1
         with:
-          python-version: 3.9
+          python-version: 3.11
 
       - name: Install dependencies with pip
         run: |
@@ -26,3 +26,23 @@ jobs:
       - run: flake8
       - run: mypy
       - run: pytest --cov --cov-fail-under=50
+  test_plotting:
+    runs-on: ubuntu-latest
+    steps:
+      - uses: actions/checkout@v2
+      - name: Setup Python
+        uses: actions/setup-python@v1
+        with:
+          python-version: 3.11
+      - name: Install dependencies with pip
+        run: |
+          python -m pip install -e .[dev]
+      - name: Create plots
+        run: |
+          python scripts/plot_radiosonde_data.py -i examples/level1/EUREC4A_BCO_Vaisala-RS_L1-ascent_20200126T2244_v3.0.0.nc -o artifacts
+      - name: Upload plots
+        uses: actions/upload-artifact@v3
+        with:
+          name: figures
+          path: |
+            artifacts

CHANGELOG.rst

@@ -2,6 +2,11 @@
 CHANGELOG
 =========
 
+unreleased
+------------------
+
+* Add scripts to plot radiosonde data incl. trajectory, skewT and profiles of measured quantities (:pr:`21`) `Laura Köhler`_.
+
 0.0.4 (2022-11-02)
 ------------------
 

requirements.txt

@@ -1,4 +1,5 @@
 bottleneck
+cartopy
 coverage==5.3.1
 dask
 metpy

scripts/README.md

@@ -0,0 +1,17 @@
+# Visualization of level 1 radiosonde data
+
+This folder contains a plotting script for converted level 1 sounding data
+
+Requirements: The script needs the `metpy` package (`conda install metpy`).
+
+Running the script:
+```
+python plot_radiosonde_data.py -i infile -o outdir
+```
+* `infile`: input file. It should include the direction (ascent or descent) in its name. If no direction is given in the name of the input file, no SkewT plot will be created.
+* `outdir`: output directory. Three subfolders `Quantities`, `SkewT`, and `Trajectories` will be created in the output directory if not already existing.
+
+The script produces 3 plots:
+ * Trajectory of the radiosonde in `Trajectories`
+ * Temperature T, dew point tau, relative humidity rh, water vapor mixing ratio mr, wind speed, and wind direction in `Quantities`
+ * SkewT diagram in `SkewT` (only if the direction, i.e. ascent or descent, is given in the input file name)