Add LPI #8

Author: mammatus95

.github/workflows/greetings.yml

@@ -12,5 +12,5 @@ jobs:
     - uses: actions/first-interaction@v1
       with:
         repo-token: ${{ secrets.GITHUB_TOKEN }}
-        issue-message: "Message that will be displayed on users' first issue"
-        pr-message: "Message that will be displayed on users' first pull request"
+        issue-message: "Thanks for you suggestions to improve the project further."
+        pr-message: "Thanks for interest in the project. We will review the pull request as quickly as possible."

.gitignore

@@ -6,18 +6,13 @@
 
 **log.txt
 
-*/*.png
-*/*/*.png
+**.png
+**.svg
+**.gif
 
-*/*.svg
-*/*/*.svg
+**.log
+**.pyc
 
-*.log
-*/*.log
-
-*.pyc
-*/*.pyc
-*/*/*.pyc
 
 src/modeldata
 src/run.yml

README.md

@@ -18,7 +18,7 @@ You need to addg Conda initialization to your etc/profile, as well.
 
 ![Map of Hodographs IFS](images/hodographmap_IFS.png)
 
-![Map of Hodographs GFS](images/hodographmap_GFS.png)
+![Map of Hodographs GFS](images/GFS.gif)
 
 
 ## How to run

images/GFS.gif

@@ -14,7 +14,7 @@ fi
 
 
 FP=${1}
-R=0 # select run 0 or 12z
+R=00 # select run 0 or 12z
 D=$(date +"%Y%m%d") # date in format YYYYMMDD
 #######################################################################
 
@@ -66,7 +66,7 @@ echo "Start downloading leadtime ${T}h"
 
 # ICON
 # single level
-for N in CAPE_ML CAPE_CON PS
+for N in CAPE_ML CAPE_CON PS LPI_CON_MAX
 do
   typeset -l nvar
   nvar=${N}
@@ -107,13 +107,16 @@ echo default_date: \"$(date +%Y-%m-%d)\" >> run.yml
 echo "Plot Hodograph Maps"
 # run python script
 #which python3
+python3 main.py ICON >> log.txt 2>&1
 python3 main.py IFS >> log.txt 2>&1
 python3 main.py GFS >> log.txt 2>&1
-python3 main.py ICON >> log.txt 2>&1
-
-echo "done with leadtime ${T}h on $(date)"
-ls -lh ./images/*${FP}.png
 
 rm run.yml
 # remove nwp files
 rm -r ${store_path}
+
+echo "done with leadtime ${T}h on $(date)"
+cd ..
+mkdir -p ./results_img
+mv ./src/images/*${FP}.png ./results_img
+ls -lh ./results_img/*${FP}.png

images/ICON.gif

@@ -6,6 +6,8 @@ hodomap : True
 
 threshold: 0
 
+lpi: True
+
 fontsize: 14
 titlesize: 18
 

run_script.sh

@@ -47,7 +47,7 @@ do
   echo "Start downloading leadtime ${T}h" 
   # ICON
   # single level
-  for N in CAPE_ML CAPE_CON PS
+  for N in CAPE_ML CAPE_CON PS LPI_CON_MAX
   do
     typeset -l nvar
     nvar=${N}

src/config.yml

@@ -60,7 +60,14 @@ def run_plots(model_obj, fp):
                 print(f"u_mean by lvl: {np.nanmean(u_fld, axis=(1, 2))}")
             plotlib.basic_plot(model_obj, cape_fld, u_fld, v_fld, lats, lons, fp,
                                threshold=config["threshold"])
-            plotlib.basic_plot_custarea(model_obj, cape_fld, u_fld, v_fld, lats, lons, fp,
+            
+            # Create a 2D array filled with NaNs
+            lpi_fld = np.full(cape_fld.shape, np.nan)
+            if config["lpi"] is True:
+                lpi_fld, _, _ = model_obj.open_icon_gribfile_single("LPI_CON_MAX", fp, path="./modeldata/")
+
+
+            plotlib.basic_plot_custarea(model_obj, cape_fld, lpi_fld, u_fld, v_fld, lats, lons, fp,
                                         threshold=config["threshold"])
         else:
             print("Wrong command line argument")

src/down_script.sh

@@ -260,7 +260,7 @@ def basic_plot(model_obj, cape_fld, u, v, lats, lons, hour, threshold=10., imfmt
     plt.close()
 
 
-def basic_plot_custarea(model_obj, cape_fld, u, v, lats, lons, hour, threshold=10., imfmt="png"):
+def basic_plot_custarea(model_obj, cape_fld, lpi_fld, u, v, lats, lons, hour, threshold=10., imfmt="png"):
     pres_levels = model_obj.getlevels()
     model_name = model_obj.getname()
     start = model_obj.getrun()
@@ -275,6 +275,12 @@ def basic_plot_custarea(model_obj, cape_fld, u, v, lats, lons, hour, threshold=1
     wx = ax.contourf(lons, lats, cape_fld[:, :], levels=clevs, transform=crs.PlateCarree(),
                      cmap=cmap, extend='max', alpha=0.4, antialiased=True)
 
+    # LPI contours
+    lpi_levs = np.array([50, 500])
+    cs = plt.contour(lons, lats, lpi_fld[:,:], levels=lpi_levs, transform=crs.PlateCarree(), colors='darkred', linewidths=1.4)
+    fmt = '%.f'
+    plt.clabel(cs, fontsize=9, inline=1, fmt=fmt)
+
     for i in range(340, 400, 5):
         for j in range(555, 665, 5):
             if np.mean(cape_fld[i-1:i+1, j-1:j+1]) > threshold:
@@ -285,8 +291,9 @@ def basic_plot_custarea(model_obj, cape_fld, u, v, lats, lons, hour, threshold=1
     cax = fig.add_axes([0.44, 0.03, 0.35, 0.05])
     fig.colorbar(wx, cax=cax, orientation='horizontal')
 
-    ax.annotate("CAPE ML(contour plot)", xy=(0.8, -0.07), xycoords='axes fraction', fontsize=14)
-    ax.annotate(r'in $J/kg$', xy=(0.8, -0.1), xycoords='axes fraction', fontsize=14)
+    ax.annotate("CAPE ML(contour plot)", xy=(0.8, -0.03), xycoords='axes fraction', fontsize=14)
+    ax.annotate(r'in $J/kg$', xy=(0.8, -0.06), xycoords='axes fraction', fontsize=14)
+    ax.annotate("LPI (darkred contours)", xy=(0.8, -0.1), xycoords='axes fraction', fontsize=14)
 
     ax.annotate('1000-850 hPa in magenta', xy=(0.02, -0.025), xycoords='axes fraction', fontsize=13)
     ax.annotate(' 850-600 hPa in blue', xy=(0.02, -0.05), xycoords='axes fraction', fontsize=13)