diff --git a/conf/grid_description/griddes_ukmo602.txt b/conf/grid_description/griddes_ukmo602.txt
new file mode 100644
index 0000000000000000000000000000000000000000..1a79ef46a19adbc37e4a28f588fc5692f474087a
--- /dev/null
+++ b/conf/grid_description/griddes_ukmo602.txt
@@ -0,0 +1,17 @@
+#
+# gridID 2
+#
+gridtype  = lonlat
+gridsize  = 64800
+xsize     = 360
+ysize     = 180
+xname     = lon
+xlongname = "longitude"
+xunits    = "degrees_east"
+yname     = lat
+ylongname = "latitude"
+yunits    = "degrees_north"
+xfirst    = 0.5
+xinc      = 1
+yfirst    = 89.5
+yinc      = -1
diff --git a/modules/Loading/R/dates2load.R b/modules/Loading/R/dates2load.R
index 1d3fcbaa4b82d2693deba8c10b76353a276a4034..a21795865d00baae7e14e9c8b378a1341ea5a4bb 100644
--- a/modules/Loading/R/dates2load.R
+++ b/modules/Loading/R/dates2load.R
@@ -21,11 +21,22 @@ dates2load <- function(recipe, logger) {
 # temporal frequency of the variable requested
 # the details in $Time request: hcst period, forecast times ...
   temp_freq <- recipe$Analysis$Variables$freq
+  system_name <- recipe$Analysis$Datasets$System$name
   recipe <- recipe$Analysis$Time
   if (temp_freq == "monthly_mean") {
-   # hcst dates
+    # hcst dates
     file_dates <- paste0(strtoi(recipe$hcst_start):strtoi(recipe$hcst_end),
                          recipe$sdate)
+    # Note: UKMO models are always missing the start date of January 1993.
+    # This piece of code removes that date from the start dates.
+    UKMO_MODELS <- c("UK-MetOffice-Glosea600", "UK-MetOffice-Glosea601",
+                     "UK-MetOffice-Glosea602", "UK-MetOffice-Glosea603")
+    if ((system_name %in% UKMO_MODELS) && (file_dates[1] == "19930101")) {
+      file_dates <- file_dates[-1]
+      warn(logger,
+           paste("January 1993 start date is not available for", system_name,
+                 "and has been removed from the list of start dates."))
+    }
     ## TODO: Add dims?
     file_dates <- .add_dims(file_dates)
   } else if (temp_freq == "weekly_mean") {
@@ -83,7 +94,6 @@ dates2load <- function(recipe, logger) {
     file_dates.fcst <- NULL
   }
   return(list(hcst = file_dates, fcst = file_dates.fcst))
-  ## TODO: document header of fun
 }
 
 # adds the correspondent dims to each sdate array
diff --git a/recipes/atomic_recipes/recipe_ukmo-tas.yml b/recipes/atomic_recipes/recipe_ukmo-tas.yml
new file mode 100644
index 0000000000000000000000000000000000000000..bea2bda58447bd706e6e1405fc8caa2605042611
--- /dev/null
+++ b/recipes/atomic_recipes/recipe_ukmo-tas.yml
@@ -0,0 +1,57 @@
+Description:
+  Author: V. Agudetse
+  Description: Analysis of MF System 7 with temperature
+Analysis:
+  Horizon: Seasonal
+  Variables:
+    name: tas
+    freq: monthly_mean
+  Datasets:
+    System:
+      name: UK-MetOffice-Glosea600
+    Multimodel: False
+    Reference:
+      name: ERA5
+  Time:
+    sdate: '0101'
+    fcst_year: '2020'
+    hcst_start: '1993'
+    hcst_end: '2010'
+    ftime_min: 1
+    ftime_max: 2
+  Region:
+    latmin: 30
+    latmax: 70
+    lonmin: -20
+    lonmax: 40
+  Regrid:
+    method: bilinear
+    type: to_system
+  Workflow:
+    Anomalies:
+      compute: yes # yes/no, default yes
+      cross_validation: yes # yes/no, default yes
+      save: 'all' # 'all'/'none'/'exp_only'/'fcst_only'
+    Calibration:
+      method: mse_min
+      save: 'none' # 'all'/'none'/'exp_only'/'fcst_only'
+    Skill:
+      metric: BSS10 BSS90
+      save: 'all' # 'all'/'none'
+    Probabilities:
+      percentiles: [[1/3, 2/3], [1/10, 9/10]]
+      save: 'percentiles_only' # 'all'/'none'/'bins_only'/'percentiles_only'
+    Visualization:
+      plots: skill_metrics, forecast_ensemble_mean, most_likely_terciles
+      multi_panel: no
+      projection: lambert_europe
+    Indicators:
+      index: no
+  ncores: 10
+  remove_NAs: yes
+  Output_format: S2S4E
+Run:
+  Loglevel: INFO
+  Terminal: yes
+  output_dir: /esarchive/scratch/vagudets/auto-s2s-outputs/
+  code_dir: /esarchive/scratch/vagudets/repos/auto-s2s/