Version 2 (modified by cbricaud, 4 years ago) (diff)

Configuration Description

ORCA is the generic name given to global ocean configurations. Its specificity lies on the horizontal curvilinear mesh used to overcome the North Pole singularity found for geographical meshes. LIM (Louvain-la-Neuve sea-ice model, multi-category model LIM3 is used) is a thermodynamic-dynamic sea ice model specifically designed for climate studies. A brief description of the model is given here.

Space-time domain

  • The horizontal resolution available through the standard configuration is ORCA2. it is based on a 2 degrees Mercator mesh, (i.e. variation of meridian scale factor as cosinus of the latitude). In the northern hemisphere the mesh has two poles so that the ratio of anisotropy is nearly one everywhere. The mean grid spacing is about 2/3 of the nominal value: for example it is 1.3 degrees for ORCA2. Other resolutions (ORCA4, ORCA05 and ORCA025) are running or under development within specific projects. In the coarse resolution version (i.e. ORCA2 and ORCA4) the meridional grid spacing is increased near the equator to improve the equatorial dynamics. Figures in pdf format of mesh and bathymetry can be found and downloaded here. The sea-ice model runs on the same grid.
  • The vertical domain spreads from the surface to a depth of 5000m. There are 31 levels, with 10 levels in the top 100m. The vertical mesh is deduced from a mathematical function of z (Madec and Imbard 1996). The ocean surface corresponds to the w-level k=1, and the ocean bottom to the w-level k=31. The last T-level (k=31) is thus always in the ground.The depths of the vertical levels and the associated scale factors can be viewed .Higher vertical resolution is used in ORCA025 and ORCA12 ( see DRAKKAR project).
  • The time step depends on the resolution. It is 1h36' for ORCA2 so that there is 15 time steps in one day.

Ocean Physics (for ORCA2)

-horizontal diffusion on momentum: the eddy viscosity coefficient depends on the geographical position. It is taken as 40000. m2/s, reduced in the equator regions (2000 m2/s) excepted near the western boundaries.
-isopycnal diffusion on tracers: the diffusion acts along the isopycnal surfaces (neutral surface) with a eddy diffusivity coefficient of 2000. m2/s.
-Eddy induced velocity parametrization with a coefficient that depends on the growth rate of baroclinic instabilities (it usually varies from 15 m2/s to 3000m2/s).
-lateral boundary conditions : zero fluxes of heat and salt and no-slip conditions are applied through lateral solid boundaries.
-bottom boundary condition : zero fluxes of heat and salt are applied through the ocean bottom. The Beckmann [19XX] simple bottom boundary layer parameterization is applied along continental slopes. A linear friction is applied on momentum.
-convection: the vertical eddy viscosity and diffusivity coefficients are increased to 1 m2/s in case of static instability.
-forcings: the ocean receives heat, freshwater, and momentum fluxes from the atmosphere and/or the sea-ice. The solar radiation penetrates the top meters of the ocean. The downward irradiance I(z) is formulated with two extinction coefficients [Paulson and Simpson, 1977], whose values correspond to a Type I water in Jerlov's classification (i.e the most transparent water)

ORCA2_LIM is a reference configuration with the following characteristics:

-global ocean configuration,
-based on a tri-polar ORCA grid, with a 2° horizontal resolution,
-31 vertical levels,
-forced with climatological surface fields,
-coupled to the sea-ice model [ http://www.geosci-model-dev.net/8/2991/2015/ LIM ]
-coupled to TOP passive tracer transport module and [ http://www.geosci-model-dev.net/8/2465/2015/ PISCES biogeochemichal model ] .

AGRIF demonstrator:

From the ORCA2_LIM configuration, a demonstrator using AGRIF nesting can be activated. It includes the global ORCA2_LIM configuration and a nested grid in the Agulhas region. To set up this configuration, after extracting NEMO :

-Build your AGRIF configuration directory from ORCA2_LIM, with the key_agrif cpp key activated:

cd NEMOGCM/CONFIG ; ./makenemo -r ORCA2_LIM -n AGRIF add_key "key_agrif"

-Using the ORCA2_LIM input files and namelist, AGRIF test configuration is ready to run

Configuration Inputs

Input files can be downloaded from the DODS server of IDRIS, (CAUTION different versions are available):

ORCA2_LIM_nemo_v3.7.tar and INPUTS_PISCES_v3.6st.tar are dedicated to NEMO trunk version
ORCA2_LIM_nemo_v3.6.tar and INPUTS_PISCES_v3.6st.tar are dedicated to NEMO version nemo_v_3.6 It includes a README file with files description
ORCA2_LIM_nemo_v3.4.tar and INPUTS_PISCES_v3.4.tar are dedicated to NEMO version tagged nemo_v_3.5 and nemo_v3.4 It includes a README file with files description
ORCA2_LIM_nemo_v3.3.tar and INPUTS_PISCES_v3.tar are dedicated to NEMO version tagged nemo_v3.3 It includes a README file with files description

Brief description of files in ORCA2_LIM_nemo_v3.2.tar file :

ahmcoef and ahmcoef.nc: 2D integer used to increase the horizontal eddy viscosity along western boundary in the equatorial strip 20°S-20°N
bathy_level.nc: the bathymetry in levels (Z-coordinates)
bathy_meter.nc and bathy_updated.nc : the bathymetry in meters (Partial-cells)
chlorophyll.nc : used for solar penetration scheme (traqsr)
coordinates.nc: scale factors e[123][tuvwf] and geographical glam[tuvwf] and gphi[tuvwf] location of each grid points T- U- V- W- and F-
data_1m_potential_temperature_nomask.nc: Levitus climatological potential temperature
data_1m_salinity_nomask.nc: Levitus climatological salinity
EMPave_old.dat: fresshwater correction flux
geothermal_heating.nc: constant geothermal heat flux
K1rowdrg.nc, M2rowdrg.nc, mask_itf.nc: parametrization of vertical mixing due to tidal mixing (key zdftmx)
ncar_precip.15JUNE2009_fill.nc, ncar_rad.15JUNE2009_fill.nc, q_10.15JUNE2009_fill.nc, t_10.15JUNE2009_fill.nc, u_10.15JUNE2009_fill.nc, v_10.15JUNE2009_fill.nc; forcing data set for CORE bulk formulae
weights_core_orca2_bicubic_noc.nc and weights_core_orca2_bilinear_noc.nc: weight ( used in interpolation-on-the-fly) to interpolate forcing fields on the model grid
resto.nc: damping zones for T&S( tradmp)

More informations on input and configuration files in NEMO [ ?????? reference manual ]

ON-THE-FLY Interpolation

NEMO allows to use the interpolation on the fly option allowing to interpolate input data during the run. If you want to use this option you need files giving informations on weights, which have been created.
You can find here 2 weights files bil_weights for scalar field (bilinear interpolation) and bic_weights for vector field (bicubic interpolation) [ http://prodn.idris.fr/thredds/catalog/ipsl_public/reee512/ORCA2_ONTHEFLY/WEIGHTS/catalog.html Weights ] files.
The data files used are COREII forcing (downloaded from http://data1.gdfl.noaa.gov/nomads/forms/mom4/COREv2) extrapolated on continents, ready to be used for on the fly option: [ http://prodn.idris.fr/thredds/catalog/ipsl_public/reee512/ORCA2_ONTHEFLY/FILLED_FILES/catalog.html COREII ] forcing files extrapolated on continents

Attachments (1)