Indices: Error when trying to compute and plot Niño1+2 and Niño3

Hi @aho,

As discussed offline, this issue is about the error that appears when trying to use the Indices module to compute some of the El Niño indices using tos as the variable.

Current Behavior

When trying to compute Niño1+2 and Niño3, the following error arises. The weird part is that it seems to happen even if I set plot_ts and plot_sp to FALSE, so it may not be coming from the plotting part of the code? But don't take my word on it. Niño3.4 and Niño4 don't raise this error.

> nino_indices <- Indices(recipe = recipe, data = data)
Error in plot.window(...) : need finite 'ylim' values
In addition: Warning messages:
1: In CST_SaveExp(data = nino$hcst, destination = file_dest, startdates = as.vector(file_dates),  :
  Metadata is not found for any variable.
2: In CST_SaveExp(data = nino$obs, destination = file_dest, startdates = as.vector(file_dates),  :
  Metadata is not found for any variable.
3: In max(abs(fcst) + 0.1, abs(obs) + 0.1, na.rm = TRUE) :
  no non-missing arguments to max; returning -Inf
4: In max(abs(fcst) + 0.1, abs(obs) + 0.1, na.rm = TRUE) :
  no non-missing arguments to max; returning -Inf

Steps To Reproduce

Recipe:

/esarchive/scratch/vagudets/auto-s2s-tests/recipes/elnino_bug.yml

Script:

source("modules/Loading/Loading.R")
source("modules/Anomalies/Anomalies.R")
source("modules/Indices/Indices.R")
recipe_file <- "/esarchive/scratch/vagudets/auto-s2s-tests/recipes/elnino_bug.yml"
recipe <- prepare_outputs(recipe_file)

data <- Loading(recipe = recipe)
data <- Anomalies(recipe = recipe, data = data)
nino_indices <- Indices(recipe = recipe, data = data)

Branch/SUNSET Version: master
Environment: WS and Nord3v2

If you need anything, please let me know.

Thanks,

Victòria

FYI @nperez

Edited Nov 23, 2023 by vagudets