Homogenization of basins files used by Earthdiagnostics
Hi Earthdiagnostics potential users, @vlapin, @rfernand, @vsicardi , @jacosta , @yruprich , @eexarchou , @emoreno, @etourign
Continuing in the line of trying to homogenize our diagnoses, I raise this issue about the basin masks used in Earthdiagnostics, which is now implemented in the workflow of auto-ecearth.
As you may know, you can ask for diagnostics (regmean
but also moc
) to be computed in specific regions.
Earthdiagnostics uses for that the basins.Ec${version}_${RES}.nc
in /esarchive/autosubmit/con_files/
.
Right now, there is:
- basins.Ec3.2_O1L75.nc
- basins.Ec3.2_O25L75.nc
- basins.Ec3.6_12L75.nc (that @emoreno recently produced).
BUT the different regions are not consistent between grids!
For instance, the first region, named AMV_North_Atlantic
looks like, depending on the grid (from ORCA1 to ORCA012):
These basins
files are produced by this script and uses the basins.yml for the definition of the regions.
Three things here:
-
the regions should be consistent between grids except that this
basins.yml
file uses for the main oceans (actually onlyAtlantic_Ocean
is called in thebasins.yml
file) the basins from Nemo directly, that is from this file: `/esarchive/autosubmit/con_files/new_maskglo.Ec3.Screenshot_from_2022-02-23_14-56-20 -
some regions such as
AMV_North_Atlantic
may have been created in a different way for ORCA1 and ORCA025, that's why it's different from the one in ORCA12 (which uses the correct definition inbasins.yml
). -
some regions already included in ORCA1 and ORCA025 are not present in the new ORCA12 file (
basins.Ec3.6_12L75.nc
), because these regions haven't been updated in ourbasins.yml
in cpg_shared_tool project. This file shows the different regions present or not: Earthdiagnostics_regions.ods
In any case, it can be the moment to redefine/clean the regions we want. Also because for instance, when calling for the 'moc' function in Earthdiagnostics, if you call it for the North_Atlantic_Ocean
region, you compute it over the North Atlantic + the whole Arctic in ORCA1 and ORCA12, over only the North Atlantic (not the highest latitudes) in ORCA025 when we may want the North_Atlantic + the extension of the Arctic only in the Atlantic part, not including the other Pacific part of the Arctic.
It can be done easily:
- we have to clean or update the basins.yml file
- we have to check the
new_maskglo.{version}_{res}.nc
file for the latest 3.6 version of NEMO and see if they are coherent between resolution (I think the issue is only for ORCA025 here). - we have to produce the basins.Ec{version}_{res}.nc for this latest 3.6 version with the
create_basins_mask.py
Who's with me?
PYI @portega