Skip to content

Homogenization of basins files used by Earthdiagnostics

Hi Earthdiagnostics potential users, @vlapin, @rfernand, @vsicardi , @jacosta , @yruprich , @eexarchou , @emoreno, @etourign

hi @mcastril , @eferre1,

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): Screenshot_from_2022-02-23_13-39-06

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 only Atlantic_Ocean is called in the basins.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 in basins.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 our basins.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