Just skip

.codespell-exclude

@@ -19,7 +19,7 @@
     # allvars = Series(["TEMP", "Q", "PRES"])
     #         "p": dset["PRES"][:].compute().values
         var_list.append("pres")  # 'pressure'
-        if var == "pres":  # Insert special versions.
+        if var == "pres":
         "SIZ",
                 df.loc[con, "variable"] = "Caf"
                 df.loc[con, "variable"] = "Laf"

monetio/models/_wrfchem_mm.py

@@ -95,15 +95,12 @@ def open_mfdataset(
         # Additional defaults for satellite analysis
         var_list.append("zstag")  # 'height'
 
-    # Make sure we get 'temperature_k' if requested
-    if "temperature_k" in var_list:
-        var_list.remove("temperature_k")
-        if "tk" not in var_list:
-            var_list.append("tk")
-
     var_wrf_list = []
     for var in var_list:
-        if var == "pres":  # Insert special versions.
+        if not surf_only_nc and var in {"temperature_k", "pres_pa_mid"}:
+            # These will already be included in the final result
+            continue
+        if var == "pres":
             var_wrf = getvar(
                 wrflist, var, timeidx=ALL_TIMES, method="cat", squeeze=False, units="Pa"
             )