use skin temp instead of soil layer 1 plus 1K as previously (#168)

Update dephy and kpt conversion routines to use skin temperature instead of soil layer 1 plus 1K as previously for surface temperature.

Co-authored-by: Steven Boeing <stevenboeing@localhost.localdomain>
Co-authored-by: Leif Denby <leif@denby.eu>

lagtraj/domain/sources/era5/download.py

@@ -335,6 +335,7 @@ def _build_single_level_an_query(date, bbox, latlon_sampling):
     # 187 Medium cloud cover  [(0 - 1)],  mcc
     # 188 High cloud cover  [(0 - 1)],  hcc
     # 198 Skin reservoir content [m of water equivalent], src
+    # 235 Skin temperature [K], skt
     # 236 Soil temperature level 41 [K], stl4
     # 238 Temperature of snow layer [K], tsn
     # 243 Forecast albedo [(0 - 1)], fal
@@ -351,8 +352,8 @@ def _build_single_level_an_query(date, bbox, latlon_sampling):
             "/35.128/39.128/40.128/41.128/42.128/66.128/67.128"
             "/129.128/134.128/136.128/139.128/141.128/151.128/159.128"
             "/160.128/161.128/162.128/163.128/164.128/170.128/172.128"
-            "/183.128/186.128/187.128/188.128/198.128/236.128/238.128"
-            "/243.128/244.128/245.128"
+            "/183.128/186.128/187.128/188.128/198.128/235.128/236.128"
+            "/238.128/243.128/244.128/245.128"
         ),
         "stream": "oper",
         "time": (

lagtraj/forcings/conversion/targets/dephy.py

@@ -428,8 +428,7 @@ def unit_check(unit_guess, variable):
     # Heat roughness, derive from "flsr" variable
     z0th_traj = np.exp(ds_era5["flsr_mean"].values)
     ds_dephy["z0th_traj"] = forcing_field_dephy(z0th_traj, "z0th_traj")
-    # Include same t_skin correction used for DALES, may need further work
-    ts = ds_era5["stl1_mean"].values + 1.0
+    ts = ds_era5["skt_mean"].values
     ds_dephy["ts"] = forcing_field_dephy(ts, "ts")
     # Surface fluxes: obtain from time mean in ERA data, change sign for dephy!
     sfc_sens_flx = -central_estimate(ds_era5["msshf_mean"].values)

lagtraj/forcings/conversion/targets/kpt.py

@@ -490,8 +490,7 @@ def from_era5(ds_era5, da_levels, parameters, metadata):
     # Heat roughness, derive from "flsr" variable
     ds_kpt["heat_rough"] = np.exp(ds_era5["flsr_mean"])
     ds_kpt["heat_rough"].attrs.update(**kpt_attributes["heat_rough"])
-    # Apply correction to t_skin (see output files)
-    ds_kpt["t_skin"] = ds_era5["stl1_mean"] + 1.0
+    ds_kpt["t_skin"] = ds_era5["skt_mean"]
     ds_kpt["t_skin"].attrs.update(**kpt_attributes["t_skin"])
     # Surface fluxes: obtain from time mean in ERA data, do not change sign
     sfc_sens_flx = central_estimate(ds_era5["msshf_mean"].values)