Changes between Version 1 and Version 2 of Users/Reference Configurations/ORCA2_LIM3_PISCES


Ignore:
Timestamp:
2017-03-02T08:23:21+01:00 (4 years ago)
Author:
cbricaud
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • Users/Reference Configurations/ORCA2_LIM3_PISCES

    v1 v2  
    1 Configuration Description 
     1'''Configuration Description''' 
    22 
    33ORCA 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. 
    44 
    5 =Space-time domain= 
     5Space-time domain 
    66 
    7   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.[[BR]] 
     7  - 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.[[BR]] 
    88 
    99 
    10  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 [https://www.drakkar-ocean.eu DRAKKAR ] project). [[BR]] 
     10  - 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 [https://www.drakkar-ocean.eu DRAKKAR ] project). [[BR]] 
    1111 
    1212 
    13  The time step depends on the resolution. It is 1h36' for ORCA2 so that there is 15 time steps in one day.[[BR]] 
     13 - The time step depends on the resolution. It is 1h36' for ORCA2 so that there is 15 time steps in one day.[[BR]] 
    1414 
    1515 
    16 =Ocean Physics (for ORCA2)= 
     16Ocean Physics (for ORCA2) 
    1717 
    18     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.[[BR]] 
    19     isopycnal diffusion on tracers: the diffusion acts along the isopycnal surfaces (neutral surface) with a eddy diffusivity coefficient of 2000. m2/s.[[BR]] 
    20     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).[[BR]] 
    21     lateral boundary conditions : zero fluxes of heat and salt and no-slip conditions are applied through lateral solid boundaries.[[BR]] 
    22     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.[[BR]] 
    23     convection: the vertical eddy viscosity and diffusivity coefficients are increased to 1 m2/s in case of static instability.[[BR]] 
    24     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)[[BR]] 
     18    -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.[[BR]] 
     19    -isopycnal diffusion on tracers: the diffusion acts along the isopycnal surfaces (neutral surface) with a eddy diffusivity coefficient of 2000. m2/s.[[BR]] 
     20    -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).[[BR]] 
     21    -lateral boundary conditions : zero fluxes of heat and salt and no-slip conditions are applied through lateral solid boundaries.[[BR]] 
     22    -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.[[BR]] 
     23    -convection: the vertical eddy viscosity and diffusivity coefficients are increased to 1 m2/s in case of static instability.[[BR]] 
     24    -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)[[BR]] 
    2525 
    2626 
    27 =ORCA2_LIM is a reference configuration with the following characteristics:= 
     27ORCA2_LIM is a reference configuration with the following characteristics: 
    2828 
    29     global ocean configuration,[[BR]] 
    30     based on a tri-polar ORCA grid, with a 2° horizontal resolution,[[BR]] 
    31     31 vertical levels,[[BR]] 
    32     forced with climatological surface fields,[[BR]] 
    33     coupled to the sea-ice model LIM[[BR]] 
     29    -global ocean configuration,[[BR]] 
     30    -based on a tri-polar ORCA grid, with a 2° horizontal resolution,[[BR]] 
     31    -31 vertical levels,[[BR]] 
     32    -forced with climatological surface fields,[[BR]] 
     33    -coupled to the sea-ice model [ http://www.geosci-model-dev.net/8/2991/2015/ LIM ] [[BR]] 
     34    -coupled to TOP passive tracer transport module and [ http://www.geosci-model-dev.net/8/2465/2015/ PISCES biogeochemichal model ] .[[BR]]  
    3435 
    3536 
    36 =AGRIF demonstrator:= 
     37AGRIF demonstrator: 
    3738 
    3839From 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. 
    39 To set up this configuration, after extracting NEMO (see http://www.nemo-ocean.eu/Using-NEMO/User-Guides/Basics/NEMO-Quick-Start-Guide):[[BR]] 
     40To set up this configuration, after extracting NEMO :[[BR]] 
    4041 
    41     Build your AGRIF configuration directory from ORCA2_LIM, with the key_agrif cpp key activated:[[BR]] 
     42    -Build your AGRIF configuration directory from ORCA2_LIM, with the key_agrif cpp key activated:[[BR]] 
    4243 
    4344 
     
    4748 
    4849 
    49     Using the ORCA2_LIM input files and namelist, AGRIF test configuration is ready to run[[BR]] 
     50    -Using the ORCA2_LIM input files and namelist, AGRIF test configuration is ready to run[[BR]] 
     51 
     52 
     53'''Configuration Inputs'''  
     54 
     55 Input files can be downloaded from the DODS server of IDRIS, (CAUTION different versions are available):[[BR]] 
     56 
     57    ORCA2_LIM_nemo_v3.7.tar and INPUTS_PISCES_v3.6st.tar are dedicated to NEMO trunk version[[BR]] 
     58    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[[BR]] 
     59    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[[BR]] 
     60    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[[BR]] 
     61   
     62     
     63 Brief description of files in ORCA2_LIM_nemo_v3.2.tar file :[[BR]] 
     64 
     65           ahmcoef and ahmcoef.nc: 2D integer used to increase the horizontal eddy viscosity along western boundary in the equatorial strip 20°S-20°N[[BR]] 
     66           bathy_level.nc: the bathymetry in levels (Z-coordinates)[[BR]] 
     67           bathy_meter.nc and bathy_updated.nc : the bathymetry in meters (Partial-cells)[[BR]] 
     68           chlorophyll.nc : used for solar penetration scheme (traqsr) [[BR]] 
     69           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-[[BR]] 
     70           data_1m_potential_temperature_nomask.nc: Levitus climatological potential temperature[[BR]] 
     71           data_1m_salinity_nomask.nc: Levitus climatological salinity[[BR]] 
     72           EMPave_old.dat: fresshwater correction flux[[BR]] 
     73           geothermal_heating.nc: constant geothermal heat flux[[BR]] 
     74           K1rowdrg.nc, M2rowdrg.nc, mask_itf.nc: parametrization of vertical mixing due to tidal mixing (key zdftmx)  [[BR]] 
     75           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[[BR]] 
     76           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[[BR]] 
     77           resto.nc: damping zones for T&S( tradmp) [[BR]] 
     78 
     79  
     80           More informations on input and configuration files in NEMO [ ?????? reference manual ] 
     81  
     82 
     83ON-THE-FLY Interpolation 
     84 
     85NEMO 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.[[BR]] 
     86You 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.[[BR]] 
     87The 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[[BR]] 
    5088 
    5189 
     
    5391 
    5492 
    55  
    56