Use ifs name 'si' for wind speed instead of 'sp'

Author: dnerini

src/plotting/colormap_defaults.py

@@ -17,8 +17,8 @@ def _fallback():
     "2t": load_ncl_colormap("t2m_29lev.ct") | {"units": "degC"},
     "10v": load_ncl_colormap("modified_uv_17lev.ct") | {"units": "m/s"},
     "10u": load_ncl_colormap("modified_uv_17lev.ct") | {"units": "m/s"},
-    "10sp": load_ncl_colormap("modified_uv_17lev.ct") | {"units": "m/s"},
-    "10si": {"cmap": plt.get_cmap("GnBu", 11), "vmin": 0, "vmax": 25},
+    "10si": load_ncl_colormap("modified_uv_17lev.ct") | {"units": "m/s"},
+    # "10si": {"cmap": plt.get_cmap("GnBu", 11), "vmin": 0, "vmax": 25},
     "t_850": {"cmap": plt.get_cmap("inferno", 11), "vmin": 220, "vmax": 310},
     "z_850": {"cmap": plt.get_cmap("coolwarm", 11), "vmin": 8000, "vmax": 17000},
     "q_925": load_ncl_colormap("RH_6lev.ct"),

workflow/Snakefile

@@ -52,7 +52,7 @@ rule showcase_all:
             / "showcases/{run_id}/{init_time}/{init_time}_{param}_{region}.gif",
             init_time=[t.strftime("%Y%m%d%H%M") for t in REFTIMES],
             run_id=collect_all_candidates(),
-            param=["tp", "2t", "10sp"],
+            param=["tp", "2t", "10si"],
             region=["globe", "europe", "switzerland"],
         ),
 

workflow/scripts/plot_forecast_frame.mo.py

@@ -76,9 +76,9 @@ def _(ArgumentParser, Path):
 def _(grib_dir, init_time, lead_time, load_state_from_grib, param):
     # load grib file
     grib_file = grib_dir / f"{init_time}_{lead_time}.grib"
-    if param == "10sp":
+    if param == "10si":
         paramlist = ["10u", "10v"]
-    elif param == "sp":
+    elif param == "si":
         paramlist = ["u", "v"]
     else:
         paramlist = [param]
@@ -159,21 +159,21 @@ def _m_to_mm(arr):
     def preprocess_field(param: str, state: dict):
         """
         - Temperatures (2t, 2d, t, d): K -> °C
-        - Wind speed at 10m (10sp): m/s -> kn, sqrt(10u^2 + 10v^2)
-        - Wind speed (sp): m/s -> kn, sqrt(u^2 + v^2)
+        - Wind speed at 10m (10si): m/s -> kn, sqrt(10u^2 + 10v^2)
+        - Wind speed (si): m/s -> kn, sqrt(u^2 + v^2)
         Returns: (field_array, units_override or None)
         """
         fields = state["fields"]
         # temperature variables
         if param in ("2t", "2d", "t", "d"):
             return _k_to_c(fields[param]), "°C"
         # 10m wind speed (allow legacy 'uv' alias)
-        if param == "10sp":
+        if param == "10si":
             u = fields["10u"]
             v = fields["10v"]
             return np.sqrt(u**2 + v**2), "m/s"
         # wind speed from standard-level components
-        if param == "sp":
+        if param == "si":
             u = fields["u"]
             v = fields["v"]
             return np.sqrt(u**2 + v**2), "m/s"