diff --git a/conf/archive_seasonal.yml b/conf/archive_seasonal.yml
index d4e2d56e6d25c29f937077713b605abf8e012dd6..03ed1636bc3ff5f6d03d7701afd5d83d75c4e3e8 100644
--- a/conf/archive_seasonal.yml
+++ b/conf/archive_seasonal.yml
@@ -25,10 +25,13 @@ gpfs:
       institution: "European Centre for Medium-Range Weather Forecasts"
       src: "exp/ecmwf/system51c3s/"
       monthly_mean: {"tas":"monthly_mean/tas_f6h/",
+		     "tasmax":"monthly_mean/tasmax_f24h/",
+                     "tasmin":"monthly_mean/tasmin_f24h/",
                      "tos":"monthly_mean/tos_f6h/",
                      "prlr":"monthly_mean/prlr_s0-24h/",
                      "sfcWind":"monthly_mean/sfcWind_f6h/",
-                     "psl":"monthly_mean/psl_f6h/"}
+                     "psl":"monthly_mean/psl_f6h/",
+                     "rsds":"monthly_mean/rsds_s0-24h/"}
       nmember:
         fcst: 51
         hcst: 25
diff --git a/operationals/bigpredidata/recipe_bigpredidata_oper.yml b/operationals/bigpredidata/recipe_bigpredidata_oper.yml
index 1f3f0d542efc4feaa7a5383102f52e13d9e0f160..4e45d3d0f1ec5c05891de3070e11fac47cef743b 100644
--- a/operationals/bigpredidata/recipe_bigpredidata_oper.yml
+++ b/operationals/bigpredidata/recipe_bigpredidata_oper.yml
@@ -24,8 +24,8 @@ Analysis:
 #      - '0301'
 #      - '0401'
 #      - '0501'
-      - '0601'
-#      - '0701'
+#      - '0601'
+      - '0701'
 #      - '0801'
 #      - '0901'
 #      - '1001'
@@ -37,8 +37,8 @@ Analysis:
     ftime_min: 1 # Mandatory, int: First leadtime time step in months 
     ftime_max: 6 # Mandatory, int: Last leadtime time step in months
   Region:
-    - {name: "Iberia", latmin: 36, latmax: 44, lonmin: -10, lonmax: 5}
-    - {name: "EU", latmin: 20, latmax: 80, lonmin: -20, lonmax: 40}
+  - {name: "Iberia", latmin: 36, latmax: 44, lonmin: -10, lonmax: 5}
+  - {name: "EU", latmin: 20, latmax: 80, lonmin: -20, lonmax: 40}
   Regrid:
     method: bilinear  # Mandatory, str: Interpolation method. See docu.
     type: "to_reference"
@@ -51,23 +51,28 @@ Analysis:
     Calibration:
       method: bias # Mandatory, str: bias, evmos, mse_min, crps_min, rpc_based
       cross_validation: yes
-      save: none 
+      save: none
     Skill:
       metric: rpss
       save: all 
       cross_validation: yes
     Probabilities:
-      percentiles: [[1/3, 2/3]] # frac: Quantile thresholds.
+      percentiles: #[[1/3, 2/3]]
+        Terciles: [1/3, 2/3]
+        P10: [1/10]
+        P90: [9/10]
       save:  all
     Indicators:
       index: no
     Visualization:
-      plots: most_likely_terciles
+      plots: most_likely_terciles, extreme_probabilities
+      brks_min: 5
+      brks_max: 95
       NA_color: white 
       multi_panel: no
       dots: no
       mask_terciles: yes
-      shapefile: recintos_provinciales_inspire_peninbal_etrs89.shp
+      shapefile: shapefiles/bigpredidata/recintos_provinciales_inspire_peninbal_etrs89.shp
       file_format: PNG # Final file format of the plots. Formats available: PNG, JPG, JPEG, EPS. Defaults to PDF.
   ncores: 12 # Optional, int: number of cores, defaults to 1
   remove_NAs: TRUE # Optional, bool: Whether NAs are removed, defaults to FALSE
@@ -82,7 +87,7 @@ Run:
   autosubmit: yes
   # fill only if using autosubmit
   auto_conf:
-    script: operational/script_bigpredidata_oper.R # replace with the path to your script
+    script: operationals/bigpredidata/script_bigpredidata_oper.R # replace with the path to your script
     expid: a8kg # replace with your EXPID
     hpc_user: bsc032413 # replace with your hpc username
     wallclock: 03:00 # hh:mm
diff --git a/operationals/bigpredidata/recipe_bigpredidata_oper_subseasonal.yml b/operationals/bigpredidata/recipe_bigpredidata_oper_subseasonal.yml
index 264321fe16777664cfb1bd33f51299d0fde776d6..83e9f86c4449f3fb8d3caba156a1a8160b4cf002 100755
--- a/operationals/bigpredidata/recipe_bigpredidata_oper_subseasonal.yml
+++ b/operationals/bigpredidata/recipe_bigpredidata_oper_subseasonal.yml
@@ -50,17 +50,22 @@ Analysis:
       save: all 
       cross_validation: yes
     Probabilities:
-      percentiles: [[1/3, 2/3]]
+      percentiles: #[[1/3, 2/3]]
+        Terciles: [1/3, 2/3]
+        P10: [1/10]
+        P90: [9/10]
       save:  all
     Indicators:
       index: no
     Visualization:
-      plots: most_likely_terciles
+      plots: most_likely_terciles, extreme_probabilities
+      brks_min: 5
+      brks_max: 95
       NA_color: white 
       multi_panel: no
       dots: no
       mask_terciles: yes
-      shapefile: recintos_provinciales_inspire_peninbal_etrs89.shp
+      shapefile: shapefiles/bigpredidata/recintos_provinciales_inspire_peninbal_etrs89.shp
       file_format: PNG # Final file format of the plots. Formats available: PNG, JPG, JPEG, EPS. Defaults to PDF.
   ncores: 12 # Optional, int: number of cores, defaults to 1
   remove_NAs: TRUE # Optional, bool: Whether NAs are removed, defaults to FALSE
@@ -75,7 +80,7 @@ Run:
   autosubmit: yes 
   # fill only if using autosubmit
   auto_conf:
-    script: operational/script_bigpredidata_oper.R # replace with the path to your script
+    script: operationals/bigpredidata/script_bigpredidata_oper.R # replace with the path to your script
     expid: a8kh # replace with your EXPID
     hpc_user: bsc032413 # replace with your hpc username
     wallclock: 05:00 # hh:mm
diff --git a/operationals/bigpredidata/script_bigpredidata_oper.R b/operationals/bigpredidata/script_bigpredidata_oper.R
index cec71a8b9ac3940cbd6f43f5cc08181069c48232..d34e7991f31f02bc84905188d43938a98254548a 100644
--- a/operationals/bigpredidata/script_bigpredidata_oper.R
+++ b/operationals/bigpredidata/script_bigpredidata_oper.R
@@ -10,10 +10,14 @@ recipe_file <- args[1]
 #recipe_file <- "/esarchive/scratch/ptrascas/R/dev-test_bigpredidata/sunset/sunset/recipes/recipe_bigpredidata_oper.yml"
 recipe <- read_atomic_recipe(recipe_file)
 
-if (recipe$Analysis$Region$name == "EU") {
-  recipe$Analysis$Regrid$type <- "to_system"
-    # "/gpfs/projects/bsc32/esarchive_cache/exp/ncep/cfs-v2/weekly_mean/s2s/tas_f24h/tas_20050109.nc"
+if (recipe$Analysis$Region$name == "EU"){
+	if (recipe$Analysis$Horizon == "subseasonal"){
+  recipe$Analysis$Regrid$type <- "/gpfs/projects/bsc32/esarchive_cache/exp/ncep/cfs-v2/weekly_mean/s2s/tas_f24h/tas_20050109.nc"
+} else if (recipe$Analysis$Horizon == "seasonal"){
+  recipe$Analysis$Regrid$type <- "/gpfs/projects/bsc32/esarchive_cache/exp/ecmwf/system51c3s/monthly_mean/tasmax_f24h/tasmax_20250701.nc"
 }
+}
+
 
 #recipe <- prepare_outputs(recipe_file)
 # Load datasets