Ano failling when doesn't using a climatology with the dimension member
Hi,
The Ano subroutine is failing when it is using a climatology without the dimension member. I do not know whether it is normal in the last s2dverification version, but it is contradictory with what is said in the documentation:
In the documentation: The Argument clim is the Climatologies from clim with the dimension c(nmod/nexp/nobs, nltime) up to c(nmod/nexp/nobs, nltime, nlevel, nlat, nlon)
Here is an example failing:
/esnas/scratch/mmenegoz/decadal_forecast/WMO_application/global/test.R
The following command fails:
clim_MO <- Clim(full_exp,obsforclim,memb="FALSE") # it means that I am using the same climatology for all the members. ano_exp <- Ano(full_exp,clim_MO$clim_exp)
whereas this command works fine:
ano_exp <- Ano(full_exp,InsertDim(clim_MO$clim_exp,2,5))
When working with lon/lat datasets, the command does not fail, but I wonder if we get correct anomalies when we do not specify that the climatology has a member dimension (even in the case that the climatology is the same for all the members).
@cprodhomme is also interested