New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 12928 for NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/cfgs/SHARED – NEMO

Ignore:
Timestamp:
2020-05-14T21:46:00+02:00 (4 years ago)
Author:
smueller
Message:

Synchronizing with /NEMO/trunk@12925 (ticket #2170)

Location:
NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser
Files:
11 edited
1 copied

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser

    • Property svn:externals
      •  

        old new  
        66^/vendors/FCM@HEAD            ext/FCM 
        77^/vendors/IOIPSL@HEAD         ext/IOIPSL 
         8 
         9# SETTE 
         10^/utils/CI/sette@HEAD         sette 
  • NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/cfgs/SHARED/README.rst

    r10598 r12928  
    33*********** 
    44 
     5.. todo:: 
     6 
     7 
     8 
    59.. contents:: 
    6            :local: 
     10   :local: 
    711 
    812Output of diagnostics in NEMO is usually done using XIOS. 
    9 This is an efficient way of writing diagnostics because the time averaging, file writing and even some simple arithmetic or regridding is carried out in parallel to the NEMO model run. 
     13This is an efficient way of writing diagnostics because 
     14the time averaging, file writing and even some simple arithmetic or regridding is carried out in 
     15parallel to the NEMO model run. 
    1016This page gives a basic introduction to using XIOS with NEMO. 
    11 Much more information is available from the XIOS homepage above and from the NEMO manual. 
     17Much more information is available from the :xios:`XIOS homepage<>` above and from the NEMO manual. 
    1218 
    13 Use of XIOS for diagnostics is activated using the pre-compiler key ``key_iomput``.  
     19Use of XIOS for diagnostics is activated using the pre-compiler key ``key_iomput``. 
    1420 
    1521Extracting and installing XIOS 
    16 ------------------------------ 
     22============================== 
    1723 
    18241. Install the NetCDF4 library. 
    19    If you want to use single file output you will need to compile the HDF & NetCDF libraries to allow parallel IO. 
    20 2. Download the version of XIOS that you wish to use. The recommended version is now XIOS 2.5: 
    21     
    22 .. code-block:: console 
     25   If you want to use single file output you will need to compile the HDF & NetCDF libraries to 
     26   allow parallel IO. 
     272. Download the version of XIOS that you wish to use. 
     28   The recommended version is now XIOS 2.5: 
    2329 
    24    $ svn co http://forge.ipsl.jussieu.fr/ioserver/svn/XIOS/branchs/xios-2.5 xios-2.5 
     30   .. code-block:: console 
    2531 
    26 and follow the instructions in `XIOS documentation <http://forge.ipsl.jussieu.fr/ioserver/wiki/documentation>`_ to compile it. 
    27    If you find problems at this stage, support can be found by subscribing to the `XIOS mailing list <http://forge.ipsl.jussieu.fr/mailman/listinfo.cgi/xios-users>`_ and sending a mail message to it.  
     32      $ svn co http://forge.ipsl.jussieu.fr/ioserver/svn/XIOS/branchs/xios-2.5 
     33 
     34and follow the instructions in :xios:`XIOS documentation <wiki/documentation>` to compile it. 
     35If you find problems at this stage, support can be found by subscribing to 
     36the :xios:`XIOS mailing list <../mailman/listinfo.cgi/xios-users>` and sending a mail message to it. 
    2837 
    2938XIOS Configuration files 
    3039------------------------ 
    3140 
    32 XIOS is controlled using xml input files that should be copied to your model run directory before running the model. 
    33 Examples of these files can be found in the reference configurations (``cfgs``). The XIOS executable expects to find a file called ``iodef.xml`` in the model run directory. 
    34 In NEMO we have made the decision to use include statements in the ``iodef.xml`` file to include ``field_def_nemo-oce.xml`` (for physics), ``field_def_nemo-ice.xml`` (for ice), ``field_def_nemo-pisces.xml`` (for biogeochemistry) and ``domain_def.xml`` from the /cfgs/SHARED directory. 
    35 Most users will not need to modify ``domain_def.xml`` or ``field_def_nemo-???.xml`` unless they want to add new diagnostics to the NEMO code. 
    36 The definition of the output files is organized into separate ``file_definition.xml`` files which are included in the ``iodef.xml`` file. 
     41XIOS is controlled using XML input files that should be copied to 
     42your model run directory before running the model. 
     43Examples of these files can be found in the reference configurations (:file:`./cfgs`). 
     44The XIOS executable expects to find a file called :file:`iodef.xml` in the model run directory. 
     45In NEMO we have made the decision to use include statements in the :file:`iodef.xml` file to include: 
     46 
     47- :file:`field_def_nemo-oce.xml` (for physics), 
     48- :file:`field_def_nemo-ice.xml` (for ice), 
     49- :file:`field_def_nemo-pisces.xml` (for biogeochemistry) and 
     50- :file:`domain_def.xml` from the :file:`./cfgs/SHARED` directory. 
     51 
     52Most users will not need to modify :file:`domain_def.xml` or :file:`field_def_nemo-???.xml` unless 
     53they want to add new diagnostics to the NEMO code. 
     54The definition of the output files is organized into separate :file:`file_definition.xml` files which 
     55are included in the :file:`iodef.xml` file. 
    3756 
    3857Modes 
    39 ----- 
     58===== 
    4059 
    4160Detached Mode 
     
    4463In detached mode the XIOS executable is executed on separate cores from the NEMO model. 
    4564This is the recommended method for using XIOS for realistic model runs. 
    46 To use this mode set ``using_server`` to ``true`` at the bottom of the ``iodef.xml`` file: 
     65To use this mode set ``using_server`` to ``true`` at the bottom of the :file:`iodef.xml` file: 
    4766 
    4867.. code-block:: xml 
    4968 
    50    <variable id="using_server" type="boolean">true</variable> 
     69   <variable id="using_server" type="boolean">true</variable> 
    5170 
    52 Make sure there is a copy (or link to) your XIOS executable in the working directory and in your job submission script allocate processors to XIOS. 
     71Make sure there is a copy (or link to) your XIOS executable in the working directory and 
     72in your job submission script allocate processors to XIOS. 
    5373 
    5474Attached Mode 
     
    5676 
    5777In attached mode XIOS runs on each of the cores used by NEMO. 
    58 This method is less efficient than the detached mode but can be more convenient for testing or with small configurations. 
    59 To activate this mode simply set ``using_server`` to false in the ``iodef.xml`` file 
     78This method is less efficient than the detached mode but can be more convenient for testing or 
     79with small configurations. 
     80To activate this mode simply set ``using_server`` to false in the :file:`iodef.xml` file 
    6081 
    6182.. code-block:: xml 
    6283 
    63    <variable id="using_server" type="boolean">false</variable> 
     84   <variable id="using_server" type="boolean">false</variable> 
    6485 
    6586and don't allocate any cores to XIOS. 
    66 Note that due to the different domain decompositions between XIOS and NEMO if the total number of cores is larger than the number of grid points in the j direction then the model run will fail. 
     87 
     88.. note:: 
     89 
     90   Due to the different domain decompositions between XIOS and NEMO, 
     91   if the total number of cores is larger than the number of grid points in the ``j`` direction then 
     92   the model run will fail. 
    6793 
    6894Adding new diagnostics 
    69 ---------------------- 
     95====================== 
    7096 
    7197If you want to add a NEMO diagnostic to the NEMO code you will need to do the following: 
    7298 
    73991. Add any necessary code to calculate you new diagnostic in NEMO 
    74 2. Send the field to XIOS using ``CALL iom_put( 'field_id', variable )`` where ``field_id`` is a unique id for your new diagnostics and variable is the fortran variable containing the data. 
    75    This should be called at every model timestep regardless of how often you want to output the field. No time averaging should be done in the model code.  
    76 3. If it is computationally expensive to calculate your new diagnostic you should also use "iom_use" to determine if it is requested in the current model run. For example, 
    77     
    78 .. code-block:: fortran 
     1002. Send the field to XIOS using ``CALL iom_put( 'field_id', variable )`` where 
     101   ``field_id`` is a unique id for your new diagnostics and 
     102   variable is the fortran variable containing the data. 
     103   This should be called at every model timestep regardless of how often you want to output the field. 
     104   No time averaging should be done in the model code. 
     1053. If it is computationally expensive to calculate your new diagnostic 
     106   you should also use "iom_use" to determine if it is requested in the current model run. 
     107   For example, 
    79108 
    80       IF iom_use('field_id') THEN 
    81          !Some expensive computation 
    82          !... 
    83          !... 
    84          iom_put('field_id', variable) 
    85       ENDIF 
     109   .. code-block:: fortran 
    86110 
    87 4. Add a variable definition to the ``field_def_nemo-???.xml`` file. 
    88 5. Add the variable to the ``iodef.xml`` or ``file_definition.xml`` file. 
     111      IF iom_use('field_id') THEN 
     112         !Some expensive computation 
     113         !... 
     114         !... 
     115    iom_put('field_id', variable) 
     116      ENDIF 
     117 
     1184. Add a variable definition to the :file:`field_def_nemo-???.xml` file. 
     1195. Add the variable to the :file:`iodef.xml` or :file:`file_definition.xml` file. 
  • NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/cfgs/SHARED/domain_def_nemo.xml

    r12225 r12928  
    1111 
    1212     <domain id="1point" domain_ref="grid_T" > 
    13        <zoom_domain ibegin="1" jbegin="1" ni="1" nj="1"/> 
     13       <zoom_domain ibegin="139" jbegin="119" ni="1" nj="1"/> 
    1414     </domain> 
    1515     <domain id="danger_T" domain_ref="grid_T" > 
     
    191191     <domain id="EqW" domain_ref="grid_W" > <zoom_domain id="EqW"/> </domain> 
    192192 
    193       
    194      <!--   zonal mean grid   --> 
    195      <domain id="gznl" long_name="gznl"/> 
    196      <domain id="znl_T" domain_ref="gznl" > <zoom_domain id="znl_T"/> </domain> 
    197      <domain id="znl_W" domain_ref="gznl" > <zoom_domain id="znl_W"/> </domain> 
     193              <!--   zonal mean grid   --> 
     194     <domain_group id="gznl"> 
     195        <domain id="gznl" long_name="gznl"/> 
     196        <domain id="ptr" domain_ref="gznl" >  
     197            <zoom_domain id="ptr" ibegin="0000" jbegin="0" ni="1" nj="0000" />  
     198        </domain> 
     199      </domain_group> 
    198200 
    199201      
  • NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/cfgs/SHARED/field_def_nemo-ice.xml

    r12178 r12928  
    9292          <field id="qemp_ice"     long_name="Downward Heat Flux from E-P over ice"                                                                              unit="W/m2" /> 
    9393          <field id="albedo"       long_name="Mean albedo over sea ice and ocean"                                                                                unit=""     /> 
     94          <field id="Cd_ice"       long_name="Momentum turbulent exchange coefficient"                                                                           unit=""     /> 
     95          <field id="Ch_ice"       long_name="Heat turbulent exchange coefficient"                                                                               unit=""     /> 
    9496      
    9597     <!-- trends --> 
     
    228230 
    229231     <!-- momentum (advection) --> 
    230      <field id="xmtrpice"     long_name="X-component of ice mass transport"                       standard_name="ice_x_transport"                       unit="kg/s" /> 
    231           <field id="ymtrpice"     long_name="Y-component of ice mass transport"                       standard_name="ice_y_transport"                       unit="kg/s" /> 
    232           <field id="xmtrpsnw"     long_name="X-component of snw mass transport"                       standard_name="snw_x_transport"                       unit="kg/s" /> 
    233           <field id="ymtrpsnw"     long_name="Y-component of snw mass transport"                       standard_name="snw_y_transport"                       unit="kg/s" /> 
    234           <field id="xatrp"        long_name="X-component of ice area transport"                       standard_name="area_x_transport"                      unit="m2/s" /> 
    235           <field id="yatrp"        long_name="Y-component of ice area transport"                       standard_name="area_y_transport"                      unit="m2/s" /> 
    236           <field id="xmtrptot"     long_name="X-component of sea-ice mass transport"                   standard_name="sea_ice_x_transport"                   unit="kg/s" > xmtrpice + xmtrpsnw </field> 
    237           <field id="ymtrptot"     long_name="Y-component of sea-ice mass transport"                   standard_name="sea_ice_y_transport"                   unit="kg/s" > ymtrpice + ymtrpsnw </field> 
    238      
     232        <field id="xmtrpice"     long_name="X-component of ice mass transport"                       standard_name="ice_x_transport"          grid_ref="grid_U_2D"              unit="kg/s" /> 
     233          <field id="ymtrpice"     long_name="Y-component of ice mass transport"                       standard_name="ice_y_transport"          grid_ref="grid_V_2D"              unit="kg/s" /> 
     234          <field id="xmtrpsnw"     long_name="X-component of snw mass transport"                       standard_name="snw_x_transport"          grid_ref="grid_U_2D"              unit="kg/s" /> 
     235          <field id="ymtrpsnw"     long_name="Y-component of snw mass transport"                       standard_name="snw_y_transport"          grid_ref="grid_V_2D"              unit="kg/s" /> 
     236          <field id="xatrp"        long_name="X-component of ice area transport"                       standard_name="area_x_transport"         grid_ref="grid_U_2D"              unit="m2/s" /> 
     237          <field id="yatrp"        long_name="Y-component of ice area transport"                       standard_name="area_y_transport"         grid_ref="grid_V_2D"              unit="m2/s" /> 
     238          <field id="xmtrptot"     long_name="X-component of sea-ice mass transport"                   standard_name="sea_ice_x_transport"      grid_ref="grid_U_2D"              unit="kg/s" > xmtrpice + xmtrpsnw </field> 
     239        <field id="ymtrptot"     long_name="Y-component of sea-ice mass transport"                   standard_name="sea_ice_y_transport"      grid_ref="grid_V_2D"              unit="kg/s" > ymtrpice + ymtrpsnw </field> 
     240 
     241     <!-- Ice transport through straits" --> 
     242         <field id="xmtrpice_ave"     long_name="Monthly average of x-ice mass transport"   field_ref="xmtrpice"          grid_ref="grid_U_2D"        freq_op="1mo" freq_offset="_reset_"  > @xmtrpice </field> 
     243         <field id="xmtrpice_section"                                                                                     grid_ref="grid_U_scalar"  > xmtrpice_ave </field> 
     244         <field id="xmtrpice_strait"                                                        field_ref="xmtrpice_section"  grid_ref="grid_U_4strait_ice"  /> 
     245         <field id="xstrait_mifl"                                                           field_ref="xmtrpice_strait"   grid_ref="grid_U_4strait_ice_hsum" unit="kg/s" detect_missing_value="true" > this * maskMFO_u_ice </field> 
     246 
     247         <field id="ymtrpice_ave"     long_name="Monthly average of y-ice mass transport"   field_ref="ymtrpice"          grid_ref="grid_V_2D"        freq_op="1mo" freq_offset="_reset_"  > @ymtrpice </field> 
     248         <field id="ymtrpice_section"                                                                                     grid_ref="grid_V_scalar"  > ymtrpice_ave </field> 
     249         <field id="ymtrpice_strait"                                                        field_ref="ymtrpice_section"  grid_ref="grid_V_4strait_ice"  /> 
     250    <field id="ystrait_mifl"                                                           field_ref="ymtrpice_strait"   grid_ref="grid_V_4strait_ice_hsum" unit="kg/s" detect_missing_value="true"  > this * maskMFO_v_ice </field> 
     251 
     252         <field id="xmtrpsnw_ave"     long_name="Monthly average of x-snow mass transport"  field_ref="xmtrpsnw"          grid_ref="grid_U_2D"        freq_op="1mo" freq_offset="_reset_"  > @xmtrpsnw </field> 
     253         <field id="xmtrpsnw_section"                                                                                     grid_ref="grid_U_scalar"  > xmtrpsnw_ave </field> 
     254         <field id="xmtrpsnw_strait"                                                        field_ref="xmtrpsnw_section"  grid_ref="grid_U_4strait_ice"  /> 
     255         <field id="xstrait_msfl"                                                           field_ref="xmtrpsnw_strait"   grid_ref="grid_U_4strait_ice_hsum" unit="kg/s" detect_missing_value="true" > this * maskMFO_u_ice </field> 
     256 
     257         <field id="ymtrpsnw_ave"     long_name="Monthly average of y-snow mass transport"  field_ref="ymtrpsnw"          grid_ref="grid_V_2D"        freq_op="1mo" freq_offset="_reset_"  > @ymtrpsnw </field> 
     258         <field id="ymtrpsnw_section"                                                                                     grid_ref="grid_V_scalar"  > ymtrpsnw_ave </field> 
     259         <field id="ymtrpsnw_strait"                                                        field_ref="ymtrpsnw_section"  grid_ref="grid_V_4strait_ice"  /> 
     260    <field id="ystrait_msfl"                                                           field_ref="ymtrpsnw_strait"   grid_ref="grid_V_4strait_ice_hsum" unit="kg/s" detect_missing_value="true"  > this * maskMFO_v_ice </field> 
     261 
     262         <field id="xatrp_ave"        long_name="Monthly average of x-ice area transport"   field_ref="xatrp"             grid_ref="grid_U_2D"        freq_op="1mo" freq_offset="_reset_"  > @xatrp </field> 
     263         <field id="xatrp_section"                                                                                        grid_ref="grid_U_scalar"  > xatrp_ave </field> 
     264         <field id="xatrp_strait"                                                           field_ref="xatrp_section"     grid_ref="grid_U_4strait_ice"  /> 
     265         <field id="xstrait_arfl"                                                           field_ref="xatrp_strait"      grid_ref="grid_U_4strait_ice_hsum" unit="kg/s" detect_missing_value="true" > this * maskMFO_u_ice </field> 
     266 
     267         <field id="yatrp_ave"        long_name="Monthly average of y-ice area transport"   field_ref="yatrp"             grid_ref="grid_V_2D"        freq_op="1mo" freq_offset="_reset_"  > @yatrp </field> 
     268         <field id="yatrp_section"                                                                                        grid_ref="grid_V_scalar"  > yatrp_ave </field> 
     269         <field id="yatrp_strait"                                                           field_ref="yatrp_section"     grid_ref="grid_V_4strait_ice"  /> 
     270    <field id="ystrait_arfl"                                                           field_ref="yatrp_strait"      grid_ref="grid_V_4strait_ice_hsum" unit="m2/s" detect_missing_value="true"  > this * maskMFO_v_ice </field> 
     271 
     272         <field id="strait_mifl"      long_name="Sea ice mass flux through straits"      standard_name="sea_ice_mass_transport_across_line"   unit="kg/s"  freq_op="1mo"  grid_ref="grid_4strait_ice" > xstrait_mifl + ystrait_mifl </field> 
     273         <field id="strait_msfl"      long_name="Snow mass flux through straits"         standard_name="snow_mass_transport_across_line"      unit="kg/s"  freq_op="1mo"  grid_ref="grid_4strait_ice" > xstrait_msfl + ystrait_msfl </field> 
     274    <field id="strait_arfl"      long_name="Sea ice area flux through straits"      standard_name="sea_area_mass_transport_across_line"  unit="m2/s"  freq_op="1mo"  grid_ref="grid_4strait_ice" > xstrait_arfl + ystrait_arfl </field> 
     275 
    239276   </field_group> <!-- SBC_2D --> 
    240277 
     
    296333          <field id="sbgheat_tot"     long_name="global mean snow heat content"                                                             unit="1e20J"    /> 
    297334 
    298      <!-- available later --> 
    299      <!-- 
    300           <field id="strait_mifl"     long_name="Sea ice mass flux through straits"      standard_name="sea_ice_mass_transport_across_line" unit="kg/s"  grid_ref="grid_4strait" /> 
    301           <field id="strait_arfl"     long_name="Sea ice area flux through straits"      standard_name="sea_ice_area_transport_across_line" unit="m2/s"  grid_ref="grid_4strait" />   
    302           <field id="strait_msfl"     long_name="Sea ice snow flux through straits"      standard_name="snow_mass_transport_across_line"    unit="kg/s"  grid_ref="grid_4strait" /> 
    303           --> 
    304335   </field_group> 
    305336   
  • NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/cfgs/SHARED/field_def_nemo-oce.xml

    r12316 r12928  
    11<?xml version="1.0"?>  
    22    <!-- $id$ --> 
     3 
     4    <field_definition level="1" prec="4" operation="average" enabled=".TRUE." default_value="1.e20" > <!-- time step automaticaly defined --> 
     5 
     6    <!-- 
     7       ===================================================================================================== 
     8       =                                   Configurable diagnostics                                        = 
     9       ===================================================================================================== 
     10      --> 
     11 
     12    <field_group id="diamlr_fields"> 
     13 
     14    <!-- 
     15       ===================================================================================================== 
     16                         Configuration of multiple-linear-regression analysis (diamlr) 
     17       ===================================================================================================== 
     18        
     19       This field group configures diamlr for tidal harmonic analysis of field 
     20       ssh: in addition to a regressor for fitting the mean value (diamlr_r101), 
     21       it includes the regressors for the analysis of the tidal constituents 
     22       that are available in the tidal-forcing implementation (see 
     23       ./src/OCE/SBC/tide.h90). 
     24 
     25      --> 
     26 
     27      <!-- Time --> 
     28      <field id="diamlr_time" grid_ref="diamlr_grid_T_2D" prec="8" /> 
     29 
     30      <!-- Regressors for tidal harmonic analysis --> 
     31      <field id="diamlr_r001"  field_ref="diamlr_time" expr="sin( __TDE_M2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:M2"   /> 
     32      <field id="diamlr_r002"  field_ref="diamlr_time" expr="cos( __TDE_M2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:M2"   /> 
     33      <field id="diamlr_r003"  field_ref="diamlr_time" expr="sin( __TDE_N2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:N2"   /> 
     34      <field id="diamlr_r004"  field_ref="diamlr_time" expr="cos( __TDE_N2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:N2"   /> 
     35      <field id="diamlr_r005"  field_ref="diamlr_time" expr="sin( __TDE_2N2_omega__  * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:2N2"  /> 
     36      <field id="diamlr_r006"  field_ref="diamlr_time" expr="cos( __TDE_2N2_omega__  * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:2N2"  /> 
     37      <field id="diamlr_r007"  field_ref="diamlr_time" expr="sin( __TDE_S2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:S2"   /> 
     38      <field id="diamlr_r008"  field_ref="diamlr_time" expr="cos( __TDE_S2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:S2"   /> 
     39      <field id="diamlr_r009"  field_ref="diamlr_time" expr="sin( __TDE_K2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:K2"   /> 
     40      <field id="diamlr_r010"  field_ref="diamlr_time" expr="cos( __TDE_K2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:K2"   /> 
     41      <field id="diamlr_r011"  field_ref="diamlr_time" expr="sin( __TDE_K1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:K1"   /> 
     42      <field id="diamlr_r012"  field_ref="diamlr_time" expr="cos( __TDE_K1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:K1"   /> 
     43      <field id="diamlr_r013"  field_ref="diamlr_time" expr="sin( __TDE_O1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:O1"   /> 
     44      <field id="diamlr_r014"  field_ref="diamlr_time" expr="cos( __TDE_O1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:O1"   /> 
     45      <field id="diamlr_r015"  field_ref="diamlr_time" expr="sin( __TDE_Q1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:Q1"   /> 
     46      <field id="diamlr_r016"  field_ref="diamlr_time" expr="cos( __TDE_Q1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:Q1"   /> 
     47      <field id="diamlr_r017"  field_ref="diamlr_time" expr="sin( __TDE_P1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:P1"   /> 
     48      <field id="diamlr_r018"  field_ref="diamlr_time" expr="cos( __TDE_P1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:P1"   /> 
     49      <field id="diamlr_r019"  field_ref="diamlr_time" expr="sin( __TDE_M4_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:M4"   /> 
     50      <field id="diamlr_r020"  field_ref="diamlr_time" expr="cos( __TDE_M4_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:M4"   /> 
     51      <field id="diamlr_r021"  field_ref="diamlr_time" expr="sin( __TDE_Mf_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:Mf"   /> 
     52      <field id="diamlr_r022"  field_ref="diamlr_time" expr="cos( __TDE_Mf_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:Mf"   /> 
     53      <field id="diamlr_r023"  field_ref="diamlr_time" expr="sin( __TDE_Mm_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:Mm"   /> 
     54      <field id="diamlr_r024"  field_ref="diamlr_time" expr="cos( __TDE_Mm_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:Mm"   /> 
     55      <field id="diamlr_r025"  field_ref="diamlr_time" expr="sin( __TDE_Msqm_omega__ * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:Msqm" /> 
     56      <field id="diamlr_r026"  field_ref="diamlr_time" expr="cos( __TDE_Msqm_omega__ * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:Msqm" /> 
     57      <field id="diamlr_r027"  field_ref="diamlr_time" expr="sin( __TDE_Mtm_omega__  * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:Mtm"  /> 
     58      <field id="diamlr_r028"  field_ref="diamlr_time" expr="cos( __TDE_Mtm_omega__  * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:Mtm"  /> 
     59      <field id="diamlr_r029"  field_ref="diamlr_time" expr="sin( __TDE_S1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:S1"   /> 
     60      <field id="diamlr_r030"  field_ref="diamlr_time" expr="cos( __TDE_S1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:S1"   /> 
     61      <field id="diamlr_r031"  field_ref="diamlr_time" expr="sin( __TDE_MU2_omega__  * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:MU2"  /> 
     62      <field id="diamlr_r032"  field_ref="diamlr_time" expr="cos( __TDE_MU2_omega__  * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:MU2"  /> 
     63      <field id="diamlr_r033"  field_ref="diamlr_time" expr="sin( __TDE_NU2_omega__  * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:NU2"  /> 
     64      <field id="diamlr_r034"  field_ref="diamlr_time" expr="cos( __TDE_NU2_omega__  * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:NU2"  /> 
     65      <field id="diamlr_r035"  field_ref="diamlr_time" expr="sin( __TDE_L2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:L2"   /> 
     66      <field id="diamlr_r036"  field_ref="diamlr_time" expr="cos( __TDE_L2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:L2"   /> 
     67      <field id="diamlr_r037"  field_ref="diamlr_time" expr="sin( __TDE_T2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:T2"   /> 
     68      <field id="diamlr_r038"  field_ref="diamlr_time" expr="cos( __TDE_T2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:T2"   /> 
     69      <field id="diamlr_r101"  field_ref="diamlr_time" expr="diamlr_time^0.0"                         enabled=".TRUE."  comment="mean"              /> 
     70 
     71      <!-- Fields selected for regression analysis --> 
     72      <field id="diamlr_f001" field_ref="ssh" enabled=".TRUE." /> 
     73 
     74    </field_group> 
    375     
    476    <!--  
     
    880============================================================================================================ 
    981    --> 
    10     <field_definition level="1" prec="4" operation="average" enabled=".TRUE." default_value="1.e20" > <!-- time step automaticaly defined --> 
     82 
     83    <field_group id="diadetide_fields"> 
     84 
     85    <!-- 
     86       ===================================================================================================== 
     87                 Weight fields for the computation of daily detided model diagnostics (diadetide) 
     88       ===================================================================================================== 
     89 
     90       --> 
     91 
     92      <field id="diadetide_weight" grid_ref="diadetide_grid_T_2D" enabled=".TRUE." /> 
     93      <field id="diadetide_weight_grid_T_2D" field_ref="diadetide_weight" grid_ref="diadetide_grid_T_2D" enabled=".TRUE." > this </field> 
     94      <field id="diadetide_weight_grid_U_2D" field_ref="diadetide_weight" grid_ref="diadetide_grid_U_2D" enabled=".TRUE." > this </field> 
     95      <field id="diadetide_weight_grid_V_2D" field_ref="diadetide_weight" grid_ref="diadetide_grid_V_2D" enabled=".TRUE." > this </field> 
     96      <field id="diadetide_weight_grid_T_3D" field_ref="diadetide_weight" grid_ref="diadetide_grid_2D_to_grid_T_3D" enabled=".TRUE." > this </field> 
     97      <field id="diadetide_weight_grid_U_3D" field_ref="diadetide_weight" grid_ref="diadetide_grid_2D_to_grid_U_3D" enabled=".TRUE." > this </field> 
     98      <field id="diadetide_weight_grid_V_3D" field_ref="diadetide_weight" grid_ref="diadetide_grid_2D_to_grid_V_3D" enabled=".TRUE." > this </field> 
     99      <field id="diadetide_weight_grid_W_3D" field_ref="diadetide_weight" grid_ref="diadetide_grid_2D_to_grid_W_3D" enabled=".TRUE." > this </field> 
     100 
     101    </field_group> 
    11102 
    12103    <!--  
     
    19110       
    20111      <field_group id="grid_T" grid_ref="grid_T_2D" > 
    21         <field id="e3t"          long_name="T-cell thickness"                    standard_name="cell_thickness"     unit="m"   grid_ref="grid_T_3D" /> 
    22         <field id="e3t_surf"     long_name="T-cell thickness"   field_ref="e3t"  standard_name="cell_thickness"     unit="m"   grid_ref="grid_T_SFC"/> 
    23         <field id="e3t_0"        long_name="Initial T-cell thickness"            standard_name="ref_cell_thickness" unit="m"   grid_ref="grid_T_3D" /> 
    24  
     112        <field id="e3t"          long_name="T-cell thickness"                    standard_name="cell_thickness"        unit="m"   grid_ref="grid_T_3D" /> 
     113        <field id="e3ts"         long_name="T-cell thickness"   field_ref="e3t"  standard_name="cell_thickness"        unit="m"   grid_ref="grid_T_SFC"/> 
     114        <field id="e3t_0"        long_name="Initial T-cell thickness"            standard_name="ref_cell_thickness"    unit="m"   grid_ref="grid_T_3D" /> 
     115        <field id="e3tb"         long_name="bottom T-cell thickness"             standard_name="bottom_cell_thickness" unit="m"   grid_ref="grid_T_2D"/>  
     116        <field id="e3t_300"      field_ref="e3t"                grid_ref="grid_T_zoom_300"       detect_missing_value="true" /> 
     117        <field id="e3t_vsum300"  field_ref="e3t_300"            grid_ref="grid_T_vsum"   detect_missing_value="true" /> 
     118   <field id="masscello"    long_name="Sea Water Mass per unit area"   standard_name="sea_water_mass_per_unit_area"   unit="kg/m2"   grid_ref="grid_T_3D"/> 
     119        <field id="volcello"     long_name="Ocean Volume"                   standard_name="ocean_volume"   unit="m3"       grid_ref="grid_T_3D"/>  
    25120        <field id="toce"         long_name="temperature"                         standard_name="sea_water_potential_temperature"   unit="degC"     grid_ref="grid_T_3D"/> 
    26121        <field id="toce_e3t"     long_name="temperature (thickness weighted)"                                                      unit="degC"     grid_ref="grid_T_3D" > toce * e3t </field > 
     
    28123        <field id="soce_e3t"     long_name="salinity    (thickness weighted)"                                                      unit="1e-3"     grid_ref="grid_T_3D" > soce * e3t </field > 
    29124 
     125        <field id="toce_e3t_300"      field_ref="toce_e3t"          unit="degree_C"     grid_ref="grid_T_zoom_300"      detect_missing_value="true" /> 
     126        <field id="toce_e3t_vsum300"  field_ref="toce_e3t_300"      unit="degress_C*m"  grid_ref="grid_T_vsum"  detect_missing_value="true" /> 
     127        <field id="toce_vmean300"     field_ref="toce_e3t_vsum300"  unit="degree_C"     grid_ref="grid_T_vsum"  detect_missing_value="true" > toce_e3t_vsum300/e3t_vsum300 </field> 
     128 
     129        <!-- AGRIF sponge --> 
     130        <field id="agrif_spt"         long_name=" AGRIF t-sponge coefficient"   unit=" " /> 
     131 
    30132        <!-- t-eddy viscosity coefficients (ldfdyn) --> 
    31133   <field id="ahmt_2d"      long_name=" surface t-eddy viscosity coefficient"   unit="m2/s or m4/s"                      /> 
    32134   <field id="ahmt_3d"      long_name=" 3D      t-eddy viscosity coefficient"   unit="m2/s or m4/s"  grid_ref="grid_T_3D"/> 
    33135 
    34         <field id="sst"          long_name="sea surface temperature"                            standard_name="sea_surface_temperature"             unit="degC"     /> 
     136        <field id="sst"          long_name="Bulk sea surface temperature"                       standard_name="bulk_sea_surface_temperature"        unit="degC"     /> 
     137        <field id="t_skin"       long_name="Skin temperature aka SSST"                          standard_name="skin_temperature"                    unit="degC"     /> 
    35138        <field id="sst2"         long_name="square of sea surface temperature"                  standard_name="square_of_sea_surface_temperature"   unit="degC2"     > sst * sst </field > 
    36139        <field id="sstmax"       long_name="max of sea surface temperature"   field_ref="sst"   operation="maximum"                                                 /> 
     
    53156        <field id="taubot"       long_name="bottom stress module"                                                                                   unit="N/m2"     />  
    54157 
     158         <!-- Case EOS = TEOS-10 : output potential temperature --> 
     159   <field id="toce_pot"     long_name="Sea Water Potential Temperature"              standard_name="sea_water_potential_temperature"   unit="degC"     grid_ref="grid_T_3D"/> 
     160        <field id="sst_pot"      long_name="potential sea surface temperature"             standard_name="sea_surface_temperature"             unit="degC"     /> 
     161        <field id="tosmint_pot"  long_name="vertical integral of potential temperature times density"   standard_name="integral_wrt_depth_of_product_of_density_and_potential_temperature"  unit="(kg m2) degree_C" /> 
     162 
     163 
    55164        <field id="ssh"          long_name="sea surface height"                                 standard_name="sea_surface_height_above_geoid"             unit="m" /> 
    56165        <field id="ssh2"         long_name="square of sea surface height"                       standard_name="square_of_sea_surface_height_above_geoid"   unit="m2" > ssh * ssh </field > 
     
    92201        <field id="topthdep"     long_name="Top of Thermocline Depth (dT = -0.2 wrt 10m)"   standard_name="ocean_mixed_layer_thickness_defined_by_temperature"                                unit="m"                         /> 
    93202        <field id="pycndep"      long_name="Pycnocline Depth (dsigma[dT=-0.2] wrt 10m)"     standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta"                                unit="m"                         /> 
    94         <field id="BLT"          long_name="Barrier Layer Thickness"                                                                                                                          unit="m"                          > topthdep - pycndep </field> 
     203        <field id="BLT"          long_name="Barrier Layer Thickness"                                                                                                                          unit="m"   > topthdep - pycndep </field> 
    95204        <field id="tinv"         long_name="Max of vertical invertion of temperature"                                                                                                         unit="degC"                      /> 
    96205        <field id="depti"        long_name="Depth of max. vert. inv. of temperature"                                                                                                          unit="m"                         /> 
    97         <field id="20d"          long_name="Depth of 20C isotherm"                          standard_name="depth_of_isosurface_of_sea_water_potential_temperature"                            unit="m"      axis_ref="iax_20C" /> 
    98         <field id="28d"          long_name="Depth of 28C isotherm"                          standard_name="depth_of_isosurface_of_sea_water_potential_temperature"                            unit="m"      axis_ref="iax_28C" /> 
     206   <field id="20d"          long_name="Depth of 20C isotherm"                          standard_name="depth_of_isosurface_of_sea_water_potential_temperature"                            unit="m"   axis_ref="iax_20C" /> 
     207        <field id="26d"          long_name="Depth of 26C isotherm"                          standard_name="depth_of_isosurface_of_sea_water_potential_temperature"                            unit="m"   axis_ref="iax_26C"   /> 
     208        <field id="28d"          long_name="Depth of 28C isotherm"                          standard_name="depth_of_isosurface_of_sea_water_potential_temperature"                            unit="m"   axis_ref="iax_28C" /> 
    99209        <field id="hc300"        long_name="Heat content 0-300m"                            standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content"   unit="J/m2"                      /> 
     210        <field id="hc700"        long_name="Heat content 0-700m"                            standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content"   unit="J/m2"                      /> 
     211        <field id="hc2000"       long_name="Heat content 0-2000m"                           standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content"   unit="J/m2"                      /> 
    100212 
    101213        <!-- variables available with diaar5 --> 
     
    119231 
    120232      <field_group id="Tides_T" grid_ref="grid_T_2D" operation="once" > 
    121         <!-- tidal composante --> 
    122         <field id="M2x"          long_name="M2 Elevation harmonic real part "                             unit="m"        /> 
    123         <field id="M2y"          long_name="M2 Elevation harmonic imaginary part"                         unit="m"        /> 
    124         <field id="S2x"          long_name="S2 Elevation harmonic real part "                             unit="m"        /> 
    125         <field id="S2y"          long_name="S2 Elevation harmonic imaginary part"                         unit="m"        /> 
    126         <field id="N2x"          long_name="N2 Elevation harmonic real part "                             unit="m"        /> 
    127         <field id="N2y"          long_name="N2 Elevation harmonic imaginary part"                         unit="m"        /> 
    128         <field id="K1x"          long_name="K1 Elevation harmonic real part "                             unit="m"        /> 
    129         <field id="K1y"          long_name="K1 Elevation harmonic imaginary part"                         unit="m"        /> 
    130         <field id="O1x"          long_name="O1 Elevation harmonic real part "                             unit="m"        /> 
    131         <field id="O1y"          long_name="O1 Elevation harmonic imaginary part"                         unit="m"        /> 
    132         <field id="Q1x"          long_name="Q1 Elevation harmonic real part "                             unit="m"        /> 
    133         <field id="Q1y"          long_name="Q1 Elevation harmonic imaginary part"                         unit="m"        /> 
    134         <field id="M4x"          long_name="M4 Elevation harmonic real part "                             unit="m"        /> 
    135         <field id="M4y"          long_name="M4 Elevation harmonic imaginary part"                         unit="m"        /> 
    136         <field id="K2x"          long_name="K2 Elevation harmonic real part "                             unit="m"        /> 
    137         <field id="K2y"          long_name="K2 Elevation harmonic imaginary part"                         unit="m"        /> 
    138         <field id="P1x"          long_name="P1 Elevation harmonic real part "                             unit="m"        /> 
    139         <field id="P1y"          long_name="P1 Elevation harmonic imaginary part"                         unit="m"        /> 
    140         <field id="Mfx"          long_name="Mf Elevation harmonic real part "                             unit="m"        /> 
    141         <field id="Mfy"          long_name="Mf Elevation harmonic imaginary part"                         unit="m"        /> 
    142         <field id="Mmx"          long_name="Mm Elevation harmonic real part "                             unit="m"        /> 
    143         <field id="Mmy"          long_name="Mm Elevation harmonic imaginary part"                         unit="m"        /> 
    144       </field_group> 
    145      
    146       <field_group id="Tides_U" grid_ref="grid_U_2D" operation="once" > 
    147         <field id="M2x_u"        long_name="M2 current barotrope along i-axis harmonic real part "        unit="m/s"      /> 
    148         <field id="M2y_u"        long_name="M2 current barotrope along i-axis harmonic imaginary part "   unit="m/s"      /> 
    149         <field id="S2x_u"        long_name="S2 current barotrope along i-axis harmonic real part "        unit="m/s"      /> 
    150         <field id="S2y_u"        long_name="S2 current barotrope along i-axis harmonic imaginary part "   unit="m/s"      /> 
    151         <field id="N2x_u"        long_name="N2 current barotrope along i-axis harmonic real part "        unit="m/s"      /> 
    152         <field id="N2y_u"        long_name="N2 current barotrope along i-axis harmonic imaginary part "   unit="m/s"      /> 
    153         <field id="K1x_u"        long_name="K1 current barotrope along i-axis harmonic real part "        unit="m/s"      /> 
    154         <field id="K1y_u"        long_name="K1 current barotrope along i-axis harmonic imaginary part "   unit="m/s"      /> 
    155         <field id="O1x_u"        long_name="O1 current barotrope along i-axis harmonic real part "        unit="m/s"      /> 
    156         <field id="O1y_u"        long_name="O1 current barotrope along i-axis harmonic imaginary part "   unit="m/s"      /> 
    157         <field id="Q1x_u"        long_name="Q1 current barotrope along i-axis harmonic real part "        unit="m/s"      /> 
    158         <field id="Q1y_u"        long_name="Q1 current barotrope along i-axis harmonic imaginary part "   unit="m/s"      /> 
    159         <field id="M4x_u"        long_name="M4 current barotrope along i-axis harmonic real part "        unit="m/s"      /> 
    160         <field id="M4y_u"        long_name="M4 current barotrope along i-axis harmonic imaginary part "   unit="m/s"      /> 
    161         <field id="K2x_u"        long_name="K2 current barotrope along i-axis harmonic real part "        unit="m/s"      /> 
    162         <field id="K2y_u"        long_name="K2 current barotrope along i-axis harmonic imaginary part "   unit="m/s"      /> 
    163         <field id="P1x_u"        long_name="P1 current barotrope along i-axis harmonic real part "        unit="m/s"      /> 
    164         <field id="P1y_u"        long_name="P1 current barotrope along i-axis harmonic imaginary part "   unit="m/s"      /> 
    165         <field id="Mfx_u"        long_name="Mf current barotrope along i-axis harmonic real part "        unit="m/s"      /> 
    166         <field id="Mfy_u"        long_name="Mf current barotrope along i-axis harmonic imaginary part "   unit="m/s"      /> 
    167         <field id="Mmx_u"        long_name="Mm current barotrope along i-axis harmonic real part "        unit="m/s"      /> 
    168         <field id="Mmy_u"        long_name="Mm current barotrope along i-axis harmonic imaginary part "   unit="m/s"      /> 
    169       </field_group> 
    170      
    171       <field_group id="Tides_V" grid_ref="grid_V_2D" operation="once" > 
    172         <field id="M2x_v"        long_name="M2 current barotrope along j-axis harmonic real part "        unit="m/s"      /> 
    173         <field id="M2y_v"        long_name="M2 current barotrope along j-axis harmonic imaginary part "   unit="m/s"      /> 
    174         <field id="S2x_v"        long_name="S2 current barotrope along j-axis harmonic real part "        unit="m/s"      /> 
    175         <field id="S2y_v"        long_name="S2 current barotrope along j-axis harmonic imaginary part "   unit="m/s"      /> 
    176         <field id="N2x_v"        long_name="N2 current barotrope along j-axis harmonic real part "        unit="m/s"      /> 
    177         <field id="N2y_v"        long_name="N2 current barotrope along j-axis harmonic imaginary part "   unit="m/s"      /> 
    178         <field id="K1x_v"        long_name="K1 current barotrope along j-axis harmonic real part "        unit="m/s"      /> 
    179         <field id="K1y_v"        long_name="K1 current barotrope along j-axis harmonic imaginary part "   unit="m/s"      /> 
    180         <field id="O1x_v"        long_name="O1 current barotrope along j-axis harmonic real part "        unit="m/s"      /> 
    181         <field id="O1y_v"        long_name="O1 current barotrope along j-axis harmonic imaginary part "   unit="m/s"      /> 
    182         <field id="Q1x_v"        long_name="Q1 current barotrope along j-axis harmonic real part "        unit="m/s"      /> 
    183         <field id="Q1y_v"        long_name="Q1 current barotrope along j-axis harmonic imaginary part "   unit="m/s"      /> 
    184         <field id="M4x_v"        long_name="M4 current barotrope along j-axis harmonic real part "        unit="m/s"      /> 
    185         <field id="M4y_v"        long_name="M4 current barotrope along j-axis harmonic imaginary part "   unit="m/s"      /> 
    186         <field id="K2x_v"        long_name="K2 current barotrope along j-axis harmonic real part "        unit="m/s"      /> 
    187         <field id="K2y_v"        long_name="K2 current barotrope along j-axis harmonic imaginary part "   unit="m/s"      /> 
    188         <field id="P1x_v"        long_name="P1 current barotrope along j-axis harmonic real part "        unit="m/s"      /> 
    189         <field id="P1y_v"        long_name="P1 current barotrope along j-axis harmonic imaginary part "   unit="m/s"      /> 
    190         <field id="Mfx_v"        long_name="Mf current barotrope along j-axis harmonic real part "        unit="m/s"      /> 
    191         <field id="Mfy_v"        long_name="Mf current barotrope along j-axis harmonic imaginary part "   unit="m/s"      /> 
    192         <field id="Mmx_v"        long_name="Mm current barotrope along j-axis harmonic real part "        unit="m/s"      /> 
    193         <field id="Mmy_v"        long_name="Mm current barotrope along j-axis harmonic imaginary part "   unit="m/s"      />    
     233        <!-- Tidal potential --> 
     234        <field id="tide_pot"      long_name="Total tidal potential" unit="m" /> 
     235        <field id="tide_pot_M2"   long_name="M2 tidal potential"    unit="m" /> 
     236        <field id="tide_pot_N2"   long_name="N2 tidal potential"    unit="m" /> 
     237        <field id="tide_pot_2N2"  long_name="2N2 tidal potential"   unit="m" /> 
     238        <field id="tide_pot_S2"   long_name="S2 tidal potential"    unit="m" /> 
     239        <field id="tide_pot_K2"   long_name="K2 tidal potential"    unit="m" /> 
     240        <field id="tide_pot_K1"   long_name="K1 tidal potential"    unit="m" /> 
     241        <field id="tide_pot_O1"   long_name="O1 tidal potential"    unit="m" /> 
     242        <field id="tide_pot_Q1"   long_name="Q1 tidal potential"    unit="m" /> 
     243        <field id="tide_pot_P1"   long_name="P1 tidal potential"    unit="m" /> 
     244        <field id="tide_pot_M4"   long_name="M4 tidal potential"    unit="m" /> 
     245        <field id="tide_pot_Mf"   long_name="Mf tidal potential"    unit="m" /> 
     246        <field id="tide_pot_Mm"   long_name="Mm tidal potential"    unit="m" /> 
     247        <field id="tide_pot_Msqm" long_name="Msqm tidal potential"  unit="m" /> 
     248        <field id="tide_pot_Mtm"  long_name="Mtm tidal potential"   unit="m" /> 
     249        <field id="tide_pot_S1"   long_name="S1 tidal potential"    unit="m" /> 
     250        <field id="tide_pot_MU2"  long_name="MU2 tidal potential"   unit="m" /> 
     251        <field id="tide_pot_NU2"  long_name="NU2 tidal potential"   unit="m" /> 
     252        <field id="tide_pot_L2"   long_name="L2 tidal potential"    unit="m" /> 
     253        <field id="tide_pot_T2"   long_name="T2 tidal potential"    unit="m" /> 
    194254      </field_group> 
    195255 
     
    305365          <field id="runoffs"      long_name="River Runoffs"                        standard_name="water_flux_into_sea_water_from_rivers"                                unit="kg/m2/s"   /> 
    306366          <field id="precip"       long_name="Total precipitation"                  standard_name="precipitation_flux"                                                   unit="kg/m2/s"   /> 
     367          <field id="wclosea"      long_name="closed sea empmr correction"          standard_name="closea_empmr"                                                         unit="kg/m2/s"   /> 
    307368      
    308369          <field id="qt"           long_name="Net Downward Heat Flux"                standard_name="surface_downward_heat_flux_in_sea_water"                              unit="W/m2"                           /> 
     
    311372          <field id="qsr3d"        long_name="Shortwave Radiation 3D distribution"   standard_name="downwelling_shortwave_flux_in_sea_water"                              unit="W/m2"      grid_ref="grid_T_3D" /> 
    312373          <field id="qrp"          long_name="Surface Heat Flux: Damping"            standard_name="heat_flux_into_sea_water_due_to_newtonian_relaxation"                 unit="W/m2"                           /> 
     374          <field id="qclosea"      long_name="closed sea heat content flux"          standard_name="closea_heat_content_downward_flux"                                    unit="W/m2"     /> 
    313375          <field id="erp"          long_name="Surface Water Flux: Damping"           standard_name="water_flux_out_of_sea_water_due_to_newtonian_relaxation"              unit="kg/m2/s"                        /> 
    314376          <field id="taum"         long_name="wind stress module"                    standard_name="magnitude_of_surface_downward_stress"                                 unit="N/m2"                           /> 
     
    319381 
    320382          <!-- * variable related to ice shelf forcing * --> 
    321           <field id="fwfisf"       long_name="Ice shelf melting"                             unit="kg/m2/s"  /> 
    322           <field id="fwfisf3d"     long_name="Ice shelf melting"                             unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
    323           <field id="qlatisf"      long_name="Ice shelf latent heat flux"                    unit="W/m2"     /> 
    324           <field id="qlatisf3d"    long_name="Ice shelf latent heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
    325           <field id="qhcisf"       long_name="Ice shelf heat content flux"                   unit="W/m2"     /> 
    326           <field id="qhcisf3d"     long_name="Ice shelf heat content flux"                   unit="W/m2"     grid_ref="grid_T_3D" /> 
    327           <field id="isfgammat"    long_name="transfert coefficient for isf (temperature) "  unit="m/s"      /> 
    328           <field id="isfgammas"    long_name="transfert coefficient for isf (salinity)    "  unit="m/s"      /> 
    329           <field id="stbl"         long_name="salinity in the Losh tbl                    "  unit="PSU"      /> 
    330           <field id="ttbl"         long_name="temperature in the Losh tbl                 "  unit="C"        /> 
    331           <field id="utbl"         long_name="zonal current in the Losh tbl at T point    "  unit="m/s"      /> 
    332           <field id="vtbl"         long_name="merid current in the Losh tbl at T point    "  unit="m/s"      /> 
    333           <field id="thermald"     long_name="thermal driving of ice shelf melting        "  unit="C"        /> 
    334           <field id="tfrz"         long_name="top freezing point (used to compute melt)   "  unit="C"        /> 
    335           <field id="tinsitu"      long_name="top insitu temperature (used to cmpt melt)  "  unit="C"        /> 
    336           <field id="ustar"        long_name="ustar at T point used in ice shelf melting  "  unit="m/s"      /> 
     383          <field id="isftfrz_cav"     long_name="freezing point temperature at ocean/isf interface"                unit="degC"     /> 
     384          <field id="isftfrz_par"     long_name="freezing point temperature in the parametrization boundary layer" unit="degC"     /> 
     385          <field id="fwfisf_cav"      long_name="Ice shelf melt rate"                           unit="kg/m2/s"  /> 
     386          <field id="fwfisf_par"      long_name="Ice shelf melt rate"                           unit="kg/m2/s"  /> 
     387          <field id="qoceisf_cav"     long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     /> 
     388          <field id="qoceisf_par"     long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     /> 
     389          <field id="qlatisf_cav"     long_name="Ice shelf latent heat flux"                    unit="W/m2"     /> 
     390          <field id="qlatisf_par"     long_name="Ice shelf latent heat flux"                    unit="W/m2"     /> 
     391          <field id="qhcisf_cav"      long_name="Ice shelf heat content flux of injected water" unit="W/m2"     /> 
     392          <field id="qhcisf_par"      long_name="Ice shelf heat content flux of injected water" unit="W/m2"     /> 
     393          <field id="fwfisf3d_cav"    long_name="Ice shelf melt rate"                           unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
     394          <field id="fwfisf3d_par"    long_name="Ice shelf melt rate"                           unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
     395          <field id="qoceisf3d_cav"   long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     396          <field id="qoceisf3d_par"   long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     397          <field id="qlatisf3d_cav"   long_name="Ice shelf latent heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     398          <field id="qlatisf3d_par"   long_name="Ice shelf latent heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     399          <field id="qhcisf3d_cav"    long_name="Ice shelf heat content flux of injected water" unit="W/m2"     grid_ref="grid_T_3D" /> 
     400          <field id="qhcisf3d_par"    long_name="Ice shelf heat content flux of injected water" unit="W/m2"     grid_ref="grid_T_3D" /> 
     401          <field id="ttbl_cav"        long_name="temperature in Losch tbl"                      unit="degC"     /> 
     402          <field id="ttbl_par"        long_name="temperature in the parametrisation boundary layer" unit="degC" /> 
     403          <field id="isfthermald_cav" long_name="thermal driving of ice shelf melting"          unit="degC"     /> 
     404          <field id="isfthermald_par" long_name="thermal driving of ice shelf melting"          unit="degC"     /> 
     405          <field id="isfgammat"       long_name="Ice shelf heat-transfert velocity"             unit="m/s"      /> 
     406          <field id="isfgammas"       long_name="Ice shelf salt-transfert velocity"             unit="m/s"      /> 
     407          <field id="stbl"            long_name="salinity in the Losh tbl"                      unit="1e-3"     /> 
     408          <field id="utbl"            long_name="zonal current in the Losh tbl at T point"      unit="m/s"      /> 
     409          <field id="vtbl"            long_name="merid current in the Losh tbl at T point"      unit="m/s"      /> 
     410          <field id="isfustar"        long_name="ustar at T point used in ice shelf melting"    unit="m/s"      /> 
     411          <field id="qconisf"         long_name="Conductive heat flux through the ice shelf"    unit="W/m2"     /> 
    337412 
    338413          <!-- *_oce variables available with ln_blk_clio or ln_blk_core --> 
     414          <field id="rho_air"      long_name="Air density at 10m above sea surface"         standard_name="rho_air_10m"                                        unit="kg/m3" /> 
     415          <field id="dt_skin"      long_name="SSST-SST temperature difference"              standard_name="SSST-SST"                                             unit="K"   /> 
    339416          <field id="qlw_oce"      long_name="Longwave Downward Heat Flux over open ocean"  standard_name="surface_net_downward_longwave_flux"                 unit="W/m2"  /> 
    340417          <field id="qsb_oce"      long_name="Sensible Downward Heat Flux over open ocean"  standard_name="surface_downward_sensible_heat_flux"                unit="W/m2"  /> 
    341418          <field id="qla_oce"      long_name="Latent Downward Heat Flux over open ocean"    standard_name="surface_downward_latent_heat_flux"                  unit="W/m2"  /> 
     419          <field id="evap_oce"     long_name="Evaporation over open ocean"                  standard_name="evaporation"                                        unit="kg/m2/s" /> 
    342420          <field id="qt_oce"       long_name="total flux at ocean surface"                  standard_name="surface_downward_heat_flux_in_sea_water"            unit="W/m2"  /> 
    343421          <field id="qsr_oce"      long_name="solar heat flux at ocean surface"             standard_name="net_downward_shortwave_flux_at_sea_water_surface"   unit="W/m2"  /> 
     
    357435 
    358436          <!-- available if key_oasis3 + conservative method --> 
    359           <field id="rain"          long_name="Liquid precipitation"                                     standard_name="rainfall_flux"                                                                 unit="kg/m2/s"  /> 
     437     <field id="rain"          long_name="Liquid precipitation"                                     standard_name="rainfall_flux"                                                                 unit="kg/m2/s"  /> 
     438          <field id="rain_ao_cea"   long_name="Liquid precipitation over ice-free ocean (cell average)"  standard_name="rainfall_flux"                                                                 unit="kg/m2/s"  /> 
    360439          <field id="evap_ao_cea"   long_name="Evaporation over ice-free ocean (cell average)"           standard_name="water_evaporation_flux"                                                        unit="kg/m2/s"  /> 
    361440          <field id="isnwmlt_cea"   long_name="Snow over Ice melting (cell average)"                     standard_name="surface_snow_melt_flux"                                                        unit="kg/m2/s"  /> 
     
    364443          <field id="hflx_rain_cea" long_name="heat flux due to rainfall"                                standard_name="temperature_flux_due_to_rainfall_expressed_as_heat_flux_into_sea_water"        unit="W/m2"     /> 
    365444          <field id="hflx_evap_cea" long_name="heat flux due to evaporation"                             standard_name="temperature_flux_due_to_evaporation_expressed_as_heat_flux_out_of_sea_water"   unit="W/m2"     /> 
     445          <field id="hflx_prec_cea" long_name="heat flux due to all precip"                              standard_name="temperature_flux_due_to_all_precip_expressed_as_heat_flux_into_sea_water"      unit="W/m2"     /> 
    366446          <field id="hflx_snow_cea" long_name="heat flux due to snow falling"                            standard_name="heat_flux_onto_ocean_and_ice_due_to_snow_thermodynamics"                       unit="W/m2"     /> 
    367447          <field id="hflx_snow_ai_cea" long_name="heat flux due to snow falling over ice"                standard_name="heat_flux_onto_ice_due_to_snow_thermodynamics"                                 unit="W/m2"     /> 
     
    400480 
    401481      </field_group> <!-- SBC --> 
    402  
     482       
     483      <!-- ABL --> 
     484      <field_group id="ABL" > <!-- time step automaticaly defined based on nn_fsbc --> 
     485 
     486   <!-- variables available with ABL on atmospheric T grid--> 
     487   <field_group id="grid_ABL3D" grid_ref="grid_TA_3D" > 
     488          <field id="u_abl"      long_name="ABL i-horizontal velocity"     standard_name="abl_x_velocity" unit="m/s"      /> 
     489          <field id="v_abl"      long_name="ABL j-horizontal velocity"     standard_name="abl_y_velocity" unit="m/s"      /> 
     490          <field id="t_abl"      long_name="ABL potential temperature"     standard_name="abl_theta"      unit="K"        /> 
     491          <field id="q_abl"      long_name="ABL specific humidity"         standard_name="abl_qspe"       unit="kg/kg"    /> 
     492          <!-- debug (to be removed) --> 
     493          <field id="u_dta"      long_name="DTA i-horizontal velocity"     standard_name="dta_x_velocity" unit="m/s"      /> 
     494          <field id="v_dta"      long_name="DTA j-horizontal velocity"     standard_name="dta_y_velocity" unit="m/s"      /> 
     495          <field id="t_dta"      long_name="DTA potential temperature"     standard_name="dta_theta"      unit="K"        /> 
     496          <field id="q_dta"      long_name="DTA specific humidity"         standard_name="dta_qspe"       unit="kg/kg"    /> 
     497          <field id="coeft"      long_name="ABL nudging coefficient"       standard_name="coeft"          unit=""         /> 
     498          <field id="tke_abl"    long_name="ABL turbulent kinetic energy"  standard_name="abl_tke"        unit="m2/s2"    /> 
     499          <field id="avm_abl"    long_name="ABL turbulent viscosity"       standard_name="abl_avm"        unit="m2/s"     /> 
     500          <field id="avt_abl"    long_name="ABL turbulent diffusivity"     standard_name="abl_avt"        unit="m2/s"     /> 
     501          <field id="mxl_abl"    long_name="ABL mixing length"             standard_name="abl_mxl"        unit="m"        /> 
     502   </field_group> 
     503 
     504   <field_group id="grid_ABL2D" grid_ref="grid_TA_2D" > 
     505          <field id="pblh"       long_name="ABL height"                    standard_name="abl_height"     unit="m"        /> 
     506          <field id="uz1_abl"    long_name="ABL i-horizontal velocity"     standard_name="abl_x_velocity" unit="m/s"      /> 
     507          <field id="vz1_abl"    long_name="ABL j-horizontal velocity"     standard_name="abl_y_velocity" unit="m/s"      /> 
     508          <field id="uvz1_abl"   long_name="ABL wind speed module"         standard_name="abl_wind_speed" unit="m/s"       > sqrt( uz1_abl^2 + vz1_abl^2 ) </field> 
     509          <field id="tz1_abl"    long_name="ABL potential temperature"     standard_name="abl_theta"      unit="K"        /> 
     510          <field id="qz1_abl"    long_name="ABL specific humidity"         standard_name="abl_qspe"       unit="kg/kg"    /> 
     511          <field id="uz1_dta"    long_name="DTA i-horizontal velocity"     standard_name="dta_x_velocity" unit="m/s"      /> 
     512          <field id="vz1_dta"    long_name="DTA j-horizontal velocity"     standard_name="dta_y_velocity" unit="m/s"      /> 
     513          <field id="uvz1_dta"   long_name="DTA wind speed module"         standard_name="dta_wind_speed" unit="m/s"       > sqrt( uz1_dta^2 + vz1_dta^2 ) </field>  
     514          <field id="tz1_dta"    long_name="DTA potential temperature"     standard_name="dta_theta"      unit="K"        /> 
     515          <field id="qz1_dta"    long_name="DTA specific humidity"         standard_name="dta_qspe"       unit="kg/kg"    /> 
     516          <!-- debug (to be removed) --> 
     517          <field id="uz1_geo"    long_name="GEO i-horizontal velocity"     standard_name="geo_x_velocity" unit="m/s"      /> 
     518          <field id="vz1_geo"    long_name="GEO j-horizontal velocity"     standard_name="geo_y_velocity" unit="m/s"      /> 
     519          <field id="uvz1_geo"   long_name="GEO wind speed module"         standard_name="geo_wind_speed" unit="m/s"       > sqrt( uz1_geo^2 + vz1_geo^2 ) </field>  
     520   </field_group> 
     521 
     522      </field_group> <!-- ABL --> 
     523 
     524       
    403525      <!-- U grid --> 
    404526       
    405527      <field_group id="grid_U"   grid_ref="grid_U_2D"> 
    406         <field id="e3u"          long_name="U-cell thickness"                                       standard_name="cell_thickness"              unit="m"          grid_ref="grid_U_3D" /> 
    407         <field id="e3u_0"        long_name="Initial U-cell thickness"                               standard_name="ref_cell_thickness"          unit="m"          grid_ref="grid_U_3D"/> 
    408         <field id="utau"         long_name="Wind Stress along i-axis"                               standard_name="surface_downward_x_stress"   unit="N/m2"                            /> 
    409         <field id="uoce"         long_name="ocean current along i-axis"                             standard_name="sea_water_x_velocity"        unit="m/s"        grid_ref="grid_U_3D" /> 
    410         <field id="uoce_e3u"     long_name="ocean current along i-axis  (thickness weighted)"                                                   unit="m/s"        grid_ref="grid_U_3D"  > uoce * e3u </field> 
     528        <field id="e2u"           long_name="U-cell width in meridional direction"                   standard_name="cell_width"                  unit="m"                               /> 
     529        <field id="e3u"           long_name="U-cell thickness"                                       standard_name="cell_thickness"              unit="m"          grid_ref="grid_U_3D" /> 
     530        <field id="e3u_0"         long_name="Initial U-cell thickness"                               standard_name="ref_cell_thickness"          unit="m"          grid_ref="grid_U_3D"/> 
     531        <field id="utau"          long_name="Wind Stress along i-axis"                               standard_name="surface_downward_x_stress"   unit="N/m2"                            /> 
     532        <field id="uoce"          long_name="ocean current along i-axis"                             standard_name="sea_water_x_velocity"        unit="m/s"        grid_ref="grid_U_3D" /> 
     533      <field id="uoce_e3u"      long_name="ocean current along i-axis  (thickness weighted)"                                                   unit="m/s"        grid_ref="grid_U_3D"  > uoce * e3u </field> 
     534      <field id="uoce_e3u_vsum" long_name="ocean current along i-axis * e3u summed on the vertical"  field_ref="uoce_e3u"    unit="m3/s"       grid_ref="grid_U_vsum"/> 
     535        <field id="uocetr_vsum"   long_name="ocean transport along i-axis  summed on the vertical"         field_ref="e2u"       unit="m3/s"> this * uoce_e3u_vsum  </field>  
     536 
     537        <field id="uocetr_vsum_op"    long_name="ocean current along i-axis * e3u * e2u summed on the vertical"  read_access="true"  freq_op="1mo"    field_ref="e2u"       unit="m3/s"> @uocetr_vsum </field> 
     538        <field id="uocetr_vsum_cumul" long_name="ocean current along i-axis * e3u * e2u cumulated from southwest point" freq_offset="_reset_" operation="instant" freq_op="1mo"  unit="m3/s" />  
     539        <field id="msftbarot"         long_name="ocean_barotropic_mass_streamfunction"   unit="kg s-1" > uocetr_vsum_cumul * $rau0 </field> 
     540 
     541 
    411542        <field id="ssu"          long_name="ocean surface current along i-axis"                                                                 unit="m/s"                             /> 
    412543        <field id="sbu"          long_name="ocean bottom current along i-axis"                                                                  unit="m/s"                             /> 
     
    415546        <field id="uocet"        long_name="ocean transport along i-axis times temperature (CRS)"                                               unit="degC*m/s"   grid_ref="grid_U_3D" /> 
    416547        <field id="uoces"        long_name="ocean transport along i-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_U_3D" /> 
    417  
     548        <field id="ssuww"        long_name="ocean surface wind work along i-axis"                   standard_name="surface_x_wind_work"         unit="N/m*s"                            > utau * ssu </field> 
     549        <!-- AGRIF sponge --> 
     550        <field id="agrif_spu"    long_name=" AGRIF u-sponge coefficient"   unit=" " /> 
    418551        <!-- u-eddy diffusivity coefficients (available if ln_traldf_OFF=F) --> 
    419552        <field id="ahtu_2d"      long_name=" surface u-eddy diffusivity coefficient"   unit="m2/s or m4/s" /> 
     
    427560 
    428561        <!-- uoce_eiv: available EIV (ln_ldfeiv=T and ln_ldfeiv_dia=T) --> 
    429         <field id="uoce_eiv"     long_name="EIV ocean current along i-axis"   standard_name="bolus_sea_water_x_velocity"   unit="m/s"   grid_ref="grid_U_3D" /> 
     562      <field id="uoce_eiv"      long_name="EIV ocean current along i-axis"                                  standard_name="bolus_sea_water_x_velocity"                     unit="m/s"   grid_ref="grid_U_3D" /> 
     563        <field id="ueiv_masstr"   long_name="EIV Ocean Mass X Transport"                                      standard_name="bolus_ocean_mass_x_transport"                   unit="kg/s"  grid_ref="grid_U_3D" /> 
     564        <field id="ueiv_heattr"   long_name="ocean bolus heat transport along i-axis"                         standard_name="ocean_heat_x_transport_due_to_bolus_advection"  unit="W"                         /> 
     565        <field id="ueiv_salttr"   long_name="ocean bolus salt transport along i-axis"                         standard_name="ocean_salt_x_transport_due_to_bolus_advection"  unit="Kg"                        /> 
     566        <field id="ueiv_heattr3d" long_name="ocean bolus heat transport along i-axis"                         standard_name="ocean_heat_x_transport_due_to_bolus_advection"  unit="W"     grid_ref="grid_U_3D" /> 
     567        <field id="ueiv_salttr3d" long_name="ocean bolus salt transport along i-axis"                         standard_name="ocean_salt_x_transport_due_to_bolus_advection"  unit="kg"    grid_ref="grid_U_3D" /> 
    430568 
    431569        <!-- uoce_bbl: available with ln_trabbl=T and nn_bbl_adv=1 --> 
     
    441579        <field id="utbl"         long_name="zonal current in the Losh tbl"     unit="m/s" /> 
    442580 
    443         <field id="u_masstr"      long_name="Ocean Mass X Transport"                                          standard_name="ocean_mass_x_transport"                         unit="kg/s" grid_ref="grid_U_3D" /> 
     581        <!-- variables available with diaar5 --> 
     582        <field id="u_masstr"      long_name="Ocean Mass X Transport"                                          standard_name="ocean_mass_x_transport"                         unit="kg/s"  grid_ref="grid_U_3D" /> 
    444583        <field id="u_masstr_vint" long_name="vertical integral of ocean eulerian mass transport along i-axis" standard_name="vertical_integral_of_ocean_mass_x_transport"    unit="kg/s"                      /> 
    445584        <field id="u_heattr"      long_name="ocean eulerian heat transport along i-axis"                      standard_name="ocean_heat_x_transport"                         unit="W"                         /> 
     
    447586        <field id="uadv_heattr"   long_name="ocean advective heat transport along i-axis"                     standard_name="advectice_ocean_heat_x_transport"               unit="W"                         /> 
    448587        <field id="uadv_salttr"   long_name="ocean advective salt transport along i-axis"                     standard_name="advectice_ocean_salt_x_transport"               unit="1e-3*kg/s"                 /> 
    449         <field id="ueiv_heattr"   long_name="ocean bolus heat transport along i-axis"                         standard_name="ocean_heat_x_transport_due_to_bolus_advection"  unit="W"                         /> 
    450         <field id="ueiv_salttr"   long_name="ocean bolus salt transport along i-axis"                         standard_name="ocean_salt_x_transport_due_to_bolus_advection"  unit="Kg"                        /> 
    451         <field id="ueiv_heattr3d" long_name="ocean bolus heat transport along i-axis"                         standard_name="ocean_heat_x_transport_due_to_bolus_advection"  unit="W"    grid_ref="grid_U_3D" /> 
    452         <field id="ueiv_salttr3d" long_name="ocean bolus salt transport along i-axis"                         standard_name="ocean_salt_x_transport_due_to_bolus_advection"  unit="kg"   grid_ref="grid_U_3D" /> 
    453588        <field id="udiff_heattr"  long_name="ocean diffusion heat transport along i-axis"                     standard_name="ocean_heat_x_transport_due_to_diffusion"        unit="W"                         /> 
    454589        <field id="udiff_salttr"  long_name="ocean diffusion salt transport along i-axis"                     standard_name="ocean_salt_x_transport_due_to_diffusion"        unit="1e-3*kg/s"                 /> 
     
    458593       
    459594      <field_group id="grid_V"   grid_ref="grid_V_2D"> 
     595        <field id="e1v"          long_name="V-cell width in longitudinal direction"                 standard_name="cell_width"                  unit="m"                              /> 
    460596        <field id="e3v"          long_name="V-cell thickness"                                       standard_name="cell_thickness"              unit="m"          grid_ref="grid_V_3D" /> 
    461597        <field id="e3v_0"        long_name="Initial V-cell thickness"                               standard_name="ref_cell_thickness"          unit="m"          grid_ref="grid_V_3D" /> 
     
    469605        <field id="vocet"        long_name="ocean transport along j-axis times temperature (CRS)"                                               unit="degC*m/s"   grid_ref="grid_V_3D" /> 
    470606        <field id="voces"        long_name="ocean transport along j-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_V_3D" /> 
    471  
     607        <field id="ssvww"        long_name="ocean surface wind work along j-axis"                   standard_name="surface_y_wind_work"         unit="N/m*s"                            > vtau * ssv </field> 
     608        <!-- AGRIF sponge --> 
     609        <field id="agrif_spv"    long_name=" AGRIF v-sponge coefficient"   unit=" " /> 
    472610        <!-- v-eddy diffusivity coefficients (available if ln_traldf_OFF=F) --> 
    473611        <field id="ahtv_2d"      long_name=" surface v-eddy diffusivity coefficient"     unit="m2/s or (m4/s)^1/2" /> 
     
    481619 
    482620        <!-- voce_eiv: available EIV (ln_ldfeiv=T and ln_ldfeiv_dia=T)  --> 
    483         <field id="voce_eiv"     long_name="EIV ocean current along j-axis"   standard_name="bolus_sea_water_y_velocity"   unit="m/s"   grid_ref="grid_V_3D" /> 
     621   <field id="voce_eiv"     long_name="EIV ocean current along j-axis"  standard_name="bolus_sea_water_y_velocity"     unit="m/s"   grid_ref="grid_V_3D" /> 
     622        <field id="veiv_masstr"  long_name="EIV Ocean Mass Y Transport"      standard_name="bolus_ocean_mass_y_transport"   unit="kg/s"  grid_ref="grid_V_3D" /> 
     623        <field id="veiv_heattr"   long_name="ocean bolus heat transport along j-axis"       standard_name="ocean_heat_y_transport_due_to_bolus_advection"   unit="W"                         /> 
     624        <field id="veiv_salttr"   long_name="ocean bolus salt transport along j-axis"       standard_name="ocean_salt_x_transport_due_to_bolus_advection"   unit="Kg"                        /> 
     625        <field id="veiv_heattr3d" long_name="ocean bolus heat transport along j-axis"       standard_name="ocean_heat_y_transport_due_to_bolus_advection"   unit="W"    grid_ref="grid_V_3D" /> 
     626        <field id="veiv_salttr3d" long_name="ocean bolus salt transport along j-axis"       standard_name="ocean_salt_y_transport_due_to_bolus_advection"   unit="kg"   grid_ref="grid_V_3D" /> 
     627 
    484628 
    485629        <!-- voce_bbl: available with ln_trabbl=T and nn_bbl_adv=1 --> 
     
    501645        <field id="vadv_heattr"   long_name="ocean advective heat transport along j-axis"   standard_name="advectice_ocean_heat_y_transport"                unit="W"                         /> 
    502646        <field id="vadv_salttr"   long_name="ocean advective salt transport along j-axis"   standard_name="advectice_ocean_salt_y_transport"                unit="1e-3*kg/s"                 /> 
    503         <field id="veiv_heattr"   long_name="ocean bolus heat transport along j-axis"       standard_name="ocean_heat_y_transport_due_to_bolus_advection"   unit="W"                         /> 
    504         <field id="veiv_salttr"   long_name="ocean bolus salt transport along j-axis"       standard_name="ocean_salt_x_transport_due_to_bolus_advection"   unit="Kg"                        /> 
    505         <field id="veiv_heattr3d" long_name="ocean bolus heat transport along j-axis"       standard_name="ocean_heat_y_transport_due_to_bolus_advection"   unit="W"    grid_ref="grid_V_3D" /> 
    506         <field id="veiv_salttr3d" long_name="ocean bolus salt transport along j-axis"       standard_name="ocean_salt_y_transport_due_to_bolus_advection"   unit="kg"   grid_ref="grid_V_3D" /> 
    507647        <field id="vdiff_heattr"  long_name="ocean diffusion heat transport along j-axis"   standard_name="ocean_heat_y_transport_due_to_diffusion"         unit="W"                         /> 
    508648        <field id="vdiff_salttr"  long_name="ocean diffusion salt transport along j-axis"   standard_name="ocean_salt_y_transport_due_to_diffusion"         unit="1e-3*kg/s"                 /> 
     
    513653      <field_group id="grid_W" grid_ref="grid_W_3D"> 
    514654        <field id="e3w"          long_name="W-cell thickness"                              standard_name="cell_thickness"                         unit="m"    /> 
    515         <field id="woce"         long_name="ocean vertical velocity"                       standard_name="upward_sea_water_velocity"              unit="m/s"  /> 
     655   <field id="woce"         long_name="ocean vertical velocity"                       standard_name="upward_sea_water_velocity"              unit="m/s"  /> 
     656   <field id="woce_e3w"     long_name="ocean vertical velocity * e3w"                                                                        unit="m2/s"  > woce * e3w </field>   
    516657        <field id="wocetr_eff"   long_name="effective ocean vertical transport"                                                                   unit="m3/s" /> 
    517658 
    518659        <!-- woce_eiv: available with EIV  (ln_ldfeiv=T and ln_ldfeiv_dia=T)  --> 
    519660        <field id="woce_eiv"     long_name="EIV ocean vertical velocity"                    standard_name="bolus_upward_sea_water_velocity"       unit="m/s"  /> 
    520  
    521         <field id="avt"          long_name="vertical eddy diffusivity"                      standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
     661        <field id="weiv_masstr"  long_name="EIV Upward Ocean Mass Transport"  standard_name="bolus_upward_ocean_mass_transport"             unit="kg/s"   /> 
     662        <field id="weiv_heattr3d" long_name="ocean bolus heat transport"    standard_name="ocean_heat_z_transport_due_to_bolus_advection"   unit="W"    /> 
     663        <field id="weiv_salttr3d" long_name="ocean bolus salt transport"    standard_name="ocean_salt_z_transport_due_to_bolus_advection"   unit="kg"   /> 
     664 
     665   <field id="avt"          long_name="vertical eddy diffusivity"                      standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
     666        <field id="avt_e3w"      long_name="vertical heat diffusivity * e3w"                unit="m3/s" > avt * e3w </field>      
    522667        <field id="logavt"       long_name="logarithm of vertical eddy diffusivity"         standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
    523668        <field id="avm"          long_name="vertical eddy viscosity"                        standard_name="ocean_vertical_momentum_diffusivity"   unit="m2/s" /> 
     669        <field id="avm_e3w"      long_name="vertical eddy viscosity * e3w"   unit="m3/s" > avm * e3w </field>  
    524670 
    525671        <!-- avs: /= avt with ln_zdfddm=T --> 
    526672        <field id="avs"          long_name="salt vertical eddy diffusivity"                 standard_name="ocean_vertical_salt_diffusivity"       unit="m2/s" /> 
    527         <field id="logavs"       long_name="logarithm of salt vertical eddy diffusivity"    standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
     673        <field id="avs_e3w"      long_name="vertical salt diffusivity * e3w"   unit="m3/s" > avs * e3w </field>  
     674   <field id="logavs"       long_name="logarithm of salt vertical eddy diffusivity"    standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
    528675 
    529676        <!-- avt_evd and avm_evd: available with ln_zdfevd --> 
    530677        <field id="avt_evd"      long_name="convective enhancement of vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_convection"     unit="m2/s" /> 
    531         <field id="avm_evd"      long_name="convective enhancement of vertical viscosity"   standard_name="ocean_vertical_momentum_diffusivity_due_to_convection"   unit="m2/s" /> 
     678        <field id="avt_evd_e3w"  long_name="convective enhancement to vertical diffusivity * e3w "    unit="m3/s" > avt_evd * e3w </field>  
     679   <field id="avm_evd"      long_name="convective enhancement of vertical viscosity"   standard_name="ocean_vertical_momentum_diffusivity_due_to_convection"   unit="m2/s" /> 
    532680 
    533681        <!-- avt_tide: available with ln_zdfiwm=T --> 
     
    548696       
    549697      <!-- F grid --> 
     698      <!-- AGRIF sponge --> 
     699      <field id="agrif_spf"    long_name=" AGRIF f-sponge coefficient"   unit=" " /> 
    550700      <!-- f-eddy viscosity coefficients (ldfdyn) --> 
    551701      <field id="ahmf_2d"      long_name=" surface f-eddy viscosity coefficient"   unit="m2/s or m4/s" /> 
     
    560710        <field id="masstot"    long_name="global total mass"                            standard_name="sea_water_mass"                                 unit="kg"   /> 
    561711        <field id="temptot"    long_name="global mean temperature"                      standard_name="sea_water_potential_temperature"                unit="degC" /> 
    562         <field id="saltot"     long_name="global mean salinity"                         standard_name="sea_water_salinity"                             unit="1e-3" /> 
    563         <field id="fram_trans" long_name="Sea Ice Mass Transport Through Fram Strait"   standard_name="sea_ice_transport_across_line"                  unit="kg/s" /> 
     712   <field id="saltot"     long_name="global mean salinity"                         standard_name="sea_water_salinity"                             unit="1e-3" /> 
     713        <field id="ssttot"     long_name="global mean sea surface temperature"          standard_name="sea_water_conservative_temperature"             unit="degC" /> 
     714         <!-- EOS = TEOS-10 --> 
     715        <field id="temptot_pot" long_name="global mean potential temperature"            standard_name="sea_water_potential_temperature"               unit="degC" /> 
    564716 
    565717         <!-- available with ln_diahsb --> 
     
    578730        <field id="bgmissal"     long_name="global mean salinity error due to free surface (linssh true)"                                                                          unit="1e-3"     /> 
    579731      </field_group> 
     732 
     733 
     734      <!-- transects --> 
     735      <field_group id="oce_straits"> 
     736         <field id="uoce_e3u_ave"         long_name="Monthly average of u*e3u"                        field_ref="uoce_e3u"                    freq_op="1mo"   freq_offset="_reset_" > @uoce_e3u </field> 
     737         <field id="uoce_e3u_ave_vsum"    long_name="Vertical sum of u*e3u"                           field_ref="uoce_e3u_ave"         grid_ref="grid_U_vsum"     /> 
     738         <field id="uocetr_vsum_section"  long_name="Total 2D transport in i-direction"               field_ref="uoce_e3u_ave_vsum"    grid_ref="grid_U_scalar"  detect_missing_value="true"> this * e2u </field> 
     739         <field id="uocetr_strait"        long_name="Total transport across lines in i-direction"     field_ref="uocetr_vsum_section"  grid_ref="grid_U_4strait" /> 
     740         <field id="u_masstr_strait"      long_name="Sea water transport across line in i-direction"  field_ref="uocetr_strait"        grid_ref="grid_U_4strait_hsum" unit="kg/s"> this * maskMFO_u * $rau0 </field> 
     741 
     742         <field id="voce_e3v_ave"         long_name="Monthly average of v*e3v"                        field_ref="voce_e3v"                    freq_op="1mo"   freq_offset="_reset_" > @voce_e3v </field> 
     743         <field id="voce_e3v_ave_vsum"    long_name="Vertical sum of v*e3v"                           field_ref="voce_e3v_ave"         grid_ref="grid_V_vsum"      /> 
     744         <field id="vocetr_vsum_section"  long_name="Total 2D transport of in j-direction"            field_ref="voce_e3v_ave_vsum"    grid_ref="grid_V_scalar"  detect_missing_value="true"> this * e1v </field> 
     745         <field id="vocetr_strait"        long_name="Total transport across lines in j-direction"     field_ref="vocetr_vsum_section"  grid_ref="grid_V_4strait"  /> 
     746         <field id="v_masstr_strait"      long_name="Sea water transport across line in j-direction"  field_ref="vocetr_strait"        grid_ref="grid_V_4strait_hsum" unit="kg/s"> this * maskMFO_v * $rau0 </field> 
     747 
     748         <field id="masstr_strait"        long_name="Sea water transport across line"                                                  grid_ref="grid_4strait"  > u_masstr_strait + v_masstr_strait </field> 
     749      </field_group> 
     750 
    580751       
    581752      <!-- variables available with ln_floats --> 
     
    613784      <!-- Poleward transport : ptr -->      
    614785      <field_group id="diaptr" >   
    615         <field id="zomsfglo"          long_name="Meridional Stream-Function: Global"           unit="Sv"        grid_ref="gznl_W_3D" /> 
    616         <field id="zomsfatl"          long_name="Meridional Stream-Function: Atlantic"         unit="Sv"        grid_ref="gznl_W_3D" /> 
    617         <field id="zomsfpac"          long_name="Meridional Stream-Function: Pacific"          unit="Sv"        grid_ref="gznl_W_3D" /> 
    618         <field id="zomsfind"          long_name="Meridional Stream-Function: Indian"           unit="Sv"        grid_ref="gznl_W_3D" /> 
    619         <field id="zomsfipc"          long_name="Meridional Stream-Function: Pacific+Indian"   unit="Sv"        grid_ref="gznl_W_3D" /> 
    620         <field id="zotemglo"          long_name="Zonal Mean Temperature : Global"              unit="degree_C"  grid_ref="gznl_T_3D" /> 
    621         <field id="zotematl"          long_name="Zonal Mean Temperature : Atlantic"            unit="degree_C"  grid_ref="gznl_T_3D" /> 
    622         <field id="zotempac"          long_name="Zonal Mean Temperature : Pacific"             unit="degree_C"  grid_ref="gznl_T_3D" /> 
    623         <field id="zotemind"          long_name="Zonal Mean Temperature : Indian"              unit="degree_C"  grid_ref="gznl_T_3D" /> 
    624         <field id="zotemipc"          long_name="Zonal Mean Temperature : Pacific+Indian"      unit="degree_C"  grid_ref="gznl_T_3D" /> 
    625         <field id="zosalglo"          long_name="Zonal Mean Salinity : Global"                 unit="0.001"     grid_ref="gznl_T_3D" /> 
    626         <field id="zosalatl"          long_name="Zonal Mean Salinity : Atlantic"               unit="0.001"     grid_ref="gznl_T_3D" /> 
    627         <field id="zosalpac"          long_name="Zonal Mean Salinity : Pacific"                unit="0.001"     grid_ref="gznl_T_3D" /> 
    628         <field id="zosalind"          long_name="Zonal Mean Salinity : Indian"                 unit="0.001"     grid_ref="gznl_T_3D" /> 
    629         <field id="zosalipc"          long_name="Zonal Mean Salinity : Pacific+Indian"         unit="0.001"     grid_ref="gznl_T_3D" /> 
    630         <field id="zosrfglo"          long_name="Zonal Mean Surface"                           unit="m2"        grid_ref="gznl_T_3D" /> 
    631         <field id="zosrfatl"          long_name="Zonal Mean Surface : Atlantic"                unit="m2"        grid_ref="gznl_T_3D" /> 
    632         <field id="zosrfpac"          long_name="Zonal Mean Surface : Pacific"                 unit="m2"        grid_ref="gznl_T_3D" /> 
    633         <field id="zosrfind"          long_name="Zonal Mean Surface : Indian"                  unit="m2"        grid_ref="gznl_T_3D" /> 
    634         <field id="zosrfipc"          long_name="Zonal Mean Surface : Pacific+Indian"          unit="m2"        grid_ref="gznl_T_3D" /> 
    635         <field id="sophtadv"          long_name="Advective Heat Transport"                     unit="PW"        grid_ref="gznl_T_2D" /> 
    636         <field id="sophtadv_atl"      long_name="Advective Heat Transport: Atlantic"           unit="PW"        grid_ref="gznl_T_2D" /> 
    637         <field id="sophtadv_pac"      long_name="Advective Heat Transport: Pacific"            unit="PW"        grid_ref="gznl_T_2D" /> 
    638         <field id="sophtadv_ind"      long_name="Advective Heat Transport: Indian"             unit="PW"        grid_ref="gznl_T_2D" /> 
    639         <field id="sophtadv_ipc"      long_name="Advective Heat Transport: Pacific+Indian"     unit="PW"        grid_ref="gznl_T_2D" /> 
    640         <field id="sophtldf"          long_name="Diffusive Heat Transport"                     unit="PW"        grid_ref="gznl_T_2D" /> 
    641         <field id="sophtldf_atl"      long_name="Diffusive Heat Transport: Atlantic"           unit="PW"        grid_ref="gznl_T_2D" /> 
    642         <field id="sophtldf_pac"      long_name="Diffusive Heat Transport: Pacific"            unit="PW"        grid_ref="gznl_T_2D" /> 
    643         <field id="sophtldf_ind"      long_name="Diffusive Heat Transport: Indian"             unit="PW"        grid_ref="gznl_T_2D" /> 
    644         <field id="sophtldf_ipc"      long_name="Diffusive Heat Transport: Pacific+Indian"     unit="PW"        grid_ref="gznl_T_2D" /> 
    645         <field id="sophtove"          long_name="Overturning Heat Transport"                   unit="PW"        grid_ref="gznl_T_2D" /> 
    646         <field id="sophtove_atl"      long_name="Overturning Heat Transport: Atlantic"         unit="PW"        grid_ref="gznl_T_2D" /> 
    647         <field id="sophtove_pac"      long_name="Overturning Heat Transport: Pacific"          unit="PW"        grid_ref="gznl_T_2D" /> 
    648         <field id="sophtove_ind"      long_name="Overturning Heat Transport: Indian"           unit="PW"        grid_ref="gznl_T_2D" /> 
    649         <field id="sophtove_ipc"      long_name="Overturning Heat Transport: Pacific+Indian"   unit="PW"        grid_ref="gznl_T_2D" /> 
    650         <field id="sophtbtr"          long_name="Barotropic Heat Transport"                    unit="PW"        grid_ref="gznl_T_2D" /> 
    651         <field id="sophtbtr_atl"      long_name="Barotropic Heat Transport: Atlantic"          unit="PW"        grid_ref="gznl_T_2D" /> 
    652         <field id="sophtbtr_pac"      long_name="Barotropic Heat Transport: Pacific"           unit="PW"        grid_ref="gznl_T_2D" /> 
    653         <field id="sophtbtr_ind"      long_name="Barotropic Heat Transport: Indian"            unit="PW"        grid_ref="gznl_T_2D" /> 
    654         <field id="sophtbtr_ipc"      long_name="Barotropic Heat Transport: Pacific+Indian"    unit="PW"        grid_ref="gznl_T_2D" /> 
    655         <field id="sophteiv"          long_name="Heat Transport from mesoscale eddy advection"                     unit="PW"       grid_ref="gznl_T_2D" /> 
    656         <field id="sophteiv_atl"      long_name="Heat Transport from mesoscale eddy advection: Atlantic"           unit="PW"       grid_ref="gznl_T_2D" /> 
    657         <field id="sophteiv_pac"      long_name="Heat Transport from mesoscale eddy advection: Pacific"            unit="PW"       grid_ref="gznl_T_2D" /> 
    658         <field id="sophteiv_ind"      long_name="Heat Transport from mesoscale eddy advection: Indian"             unit="PW"       grid_ref="gznl_T_2D" /> 
    659         <field id="sophteiv_ipc"      long_name="Heat Transport from mesoscale eddy advection: Pacific+Indian"     unit="PW"       grid_ref="gznl_T_2D" /> 
    660         <field id="sopstadv"          long_name="Advective Salt Transport"                     unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    661         <field id="sopstadv_atl"      long_name="Advective Salt Transport: Atlantic"           unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    662         <field id="sopstadv_pac"      long_name="Advective Salt Transport: Pacific"            unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    663         <field id="sopstadv_ind"      long_name="Advective Salt Transport: Indian"             unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    664         <field id="sopstadv_ipc"      long_name="Advective Salt Transport: Pacific+Indian"     unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    665         <field id="sopstove"          long_name="Overturning Salt Transport"                   unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    666         <field id="sopstove_atl"      long_name="Overturning Salt Transport: Atlantic"         unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    667         <field id="sopstove_pac"      long_name="Overturning Salt Transport: Pacific"          unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    668         <field id="sopstove_ind"      long_name="Overturning Salt Transport: Indian"           unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    669         <field id="sopstove_ipc"      long_name="Overturning Salt Transport: Pacific+Indian"   unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    670         <field id="sopstbtr"          long_name="Barotropic Salt Transport"                    unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    671         <field id="sopstbtr_atl"      long_name="Barotropic Salt Transport: Atlantic"          unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    672         <field id="sopstbtr_pac"      long_name="Barotropic Salt Transport: Pacific"           unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    673         <field id="sopstbtr_ind"      long_name="Barotropic Salt Transport: Indian"            unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    674         <field id="sopstbtr_ipc"      long_name="Barotropic Salt Transport: Pacific+Indian"    unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    675         <field id="sopstldf"          long_name="Diffusive Salt Transport"                     unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    676         <field id="sopstldf_atl"      long_name="Diffusive Salt Transport: Atlantic"           unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    677         <field id="sopstldf_pac"      long_name="Diffusive Salt Transport: Pacific"            unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    678         <field id="sopstldf_ind"      long_name="Diffusive Salt Transport: Indian"             unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    679         <field id="sopstldf_ipc"      long_name="Diffusive Salt Transport: Pacific+Indian"     unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    680         <field id="sopsteiv"          long_name="Salt Transport from mesoscale eddy advection"                     unit="Giga g/s"       grid_ref="gznl_T_2D" /> 
    681         <field id="sopsteiv_atl"      long_name="Salt Transport from mesoscale eddy advection: Atlantic"           unit="Giga g/s"       grid_ref="gznl_T_2D" /> 
    682         <field id="sopsteiv_pac"      long_name="Salt Transport from mesoscale eddy advection: Pacific"            unit="Giga g/s"       grid_ref="gznl_T_2D" /> 
    683         <field id="sopsteiv_ind"      long_name="Salt Transport from mesoscale eddy advection: Indian"             unit="Giga g/s"       grid_ref="gznl_T_2D" /> 
    684         <field id="sopsteiv_ipc"      long_name="Salt Transport from mesoscale eddy advection: Pacific+Indian"     unit="Giga g/s"       grid_ref="gznl_T_2D" />        
    685       </field_group> 
     786        <field id="zomsf"         long_name="Overturning Stream-Function : All basins"                     unit="Sv"         grid_ref="grid_znl_W_3D" /> 
     787        <field id="zotem"         long_name="Zonal Mean Temperature : All basins"                          unit="degree_C"   grid_ref="grid_znl_T_3D" /> 
     788        <field id="zosal"         long_name="Zonal Mean Salinity : All basins"                             unit="0.001"      grid_ref="grid_znl_T_3D" /> 
     789        <field id="zosrf"         long_name="Zonal Mean Surface : All basins"                              unit="m2"         grid_ref="grid_znl_T_3D" /> 
     790        <field id="sophtove"      long_name="Overturning Heat Transport: All basins"                       unit="PW"         grid_ref="grid_znl_T_2D" /> 
     791        <field id="sopstove"      long_name="Overturning Salt Transport: All basins"                       unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
     792        <field id="sophtbtr"      long_name="Barotropic Heat Transport: All basins"                        unit="PW"         grid_ref="grid_znl_T_2D" /> 
     793        <field id="sopstbtr"      long_name="Barotropic Salt Transport: All basins"                        unit="Giga g/s"   grid_ref="grid_znl_T_2D" />  
     794        <field id="sophtadv"      long_name="Advective Heat Transport: All basins"                         unit="PW"         grid_ref="grid_znl_T_2D" /> 
     795        <field id="sopstadv"      long_name="Advective Salt Transport: All basins"                         unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
     796        <field id="sophtldf"      long_name="Diffusive Heat Transport: All basins"                         unit="PW"         grid_ref="grid_znl_T_2D" /> 
     797        <field id="sopstldf"      long_name="Diffusive Salt Transport: All basins"                         unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
     798        <field id="sophtvtr"      long_name="Heat Transport : All basins"                                  unit="PW"         grid_ref="grid_znl_T_2D" /> 
     799        <field id="sopstvtr"      long_name="Salt Transport : All basins"                                  unit="Giga g/s"   grid_ref="grid_znl_T_2D" />   
     800        <field id="sophteiv"      long_name="Heat Transport from mesoscale eddy advection: All basins"     unit="PW"         grid_ref="grid_znl_T_2D" /> 
     801        <field id="sopsteiv"      long_name="Salt Transport from mesoscale eddy advection : All basins"    unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
     802   <field id="sopstadv"      long_name="Advective Salt Transport"                                     unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
     803        <field id="sophtgyre"     long_name="Overturning heat transport due to gyre" field_ref="sophtove"  unit="W"          grid_ref="grid_znl_T_2D" > sophtvtr - sophtove  </field> 
     804   <field id="sopstgyre"     long_name="Overturning salt transport due to gyre" field_ref="sopstove"  unit="kg/s"       grid_ref="grid_znl_T_2D" > sophtvtr - sopstove  </field> 
     805       </field_group> 
     806 
     807       <field_group id="constant_fields"  grid_ref="grid_T_2D"  operation="once" > 
     808         <field id="bathy"       long_name="Sea floor depth below geoid"               standard_name="sea_floor_depth_below_geoid"               unit="m"/> 
     809         <field id="areacello"   long_name="Horizontal area of ocean grid cells"       standard_name="cell_area"                                 unit="m2" /> 
     810         <field id="hfgeou"      long_name="Upward geothermal heat flux at sea floor"  standard_name="upward_geothermal_heat_flux_at_sea_floor"  unit="W/m2"/> 
     811         <field id="basins"      long_name="ocean tracer region masks"                 standard_name="ocean_tracer_region_masks"   unit="none" grid_ref="grid_basin"  /> 
     812       </field_group> 
     813 
    686814 
    687815    <!--  
     
    716844     <field id="ttrd_totad"    long_name="temperature-trend: total advection"         standard_name="tendency_of_sea_water_salinity_due_to_advection"              unit="degC/s"                        /> 
    717845     <field id="strd_totad"    long_name="salinity   -trend: total advection"         standard_name="tendency_of_sea_water_salinity_due_to_advection"              unit="1e-3/s"                        /> 
    718      <field id="ttrd_sad"      long_name="temperature-trend: surface adv. (linssh true)"                                                                                unit="degC/s"   grid_ref="grid_T_2D" /> 
    719      <field id="strd_sad"      long_name="salinity   -trend: surface adv. (linssh true)"                                                                                unit="1e-3/s"   grid_ref="grid_T_2D" /> 
     846     <field id="ttrd_sad"      long_name="temperature-trend: surface adv. (linssh true)"                                                                           unit="degC/s"   grid_ref="grid_T_2D" /> 
     847     <field id="strd_sad"      long_name="salinity   -trend: surface adv. (linssh true)"                                                                           unit="1e-3/s"   grid_ref="grid_T_2D" /> 
    720848     <field id="ttrd_ldf"      long_name="temperature-trend: lateral  diffusion"      standard_name="tendency_of_sea_water_temperature_due_to_horizontal_mixing"   unit="degC/s"                        /> 
    721849     <field id="strd_ldf"      long_name="salinity   -trend: lateral  diffusion"      standard_name="tendency_of_sea_water_salinity_due_to_horizontal_mixing"      unit="1e-3/s"                        /> 
     
    783911     <field id="ttrd_npc_e3t"      unit="degC/s * m"  >  ttrd_npc * e3t </field> 
    784912     <field id="strd_npc_e3t"      unit="1e-3/s * m"  >  strd_npc * e3t </field> 
    785      <field id="ttrd_qns_e3t"      unit="degC/s * m"  >  ttrd_qns * e3t_surf </field> 
    786      <!-- <field id="strd_cdt_e3t"      unit="degC/s * m"  >  strd_cdt * e3t_surf </field> --> 
     913     <field id="ttrd_qns_e3t"      unit="degC/s * m"  >  ttrd_qns * e3ts </field> 
     914     <!-- <field id="strd_cdt_e3t"      unit="degC/s * m"  >  strd_cdt * e3ts </field> --> 
    787915     <field id="ttrd_qsr_e3t"      unit="degC/s * m"  >  ttrd_qsr * e3t </field> 
    788916     <field id="ttrd_bbc_e3t"      unit="degC/s * m"  >  ttrd_bbc * e3t </field> 
     
    9941122    <field_group id="groupU" > 
    9951123      <field field_ref="uoce"         name="uo"      long_name="sea_water_x_velocity"      /> 
    996       <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"    /> 
    9971124      <field field_ref="utau"         name="tauuo"   long_name="surface_downward_x_stress" /> 
    9981125    </field_group> 
     
    10001127    <field_group id="groupV" > 
    10011128      <field field_ref="voce"         name="vo"      long_name="sea_water_y_velocity"      /> 
    1002       <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"    /> 
    10031129      <field field_ref="vtau"         name="tauvo"   long_name="surface_downward_y_stress" /> 
    10041130    </field_group> 
     
    10061132    <field_group id="groupW" > 
    10071133      <field field_ref="woce"         name="wo"       long_name="ocean vertical velocity"  /> 
    1008     </field_group> 
    1009  
    1010     <!-- TMB diagnostic output --> 
    1011     <field_group  id="1h_grid_T_tmb" grid_ref="grid_T_2D" operation="instant"> 
    1012       <field id="top_temp"           name="votemper_top"  unit="degC"  /> 
    1013       <field id="mid_temp"           name="votemper_mid"  unit="degC"  /> 
    1014       <field id="bot_temp"           name="votemper_bot"  unit="degC"  /> 
    1015       <field id="top_sal"            name="vosaline_top"  unit="psu"   /> 
    1016       <field id="mid_sal"            name="vosaline_mid"  unit="psu"   /> 
    1017       <field id="bot_sal"            name="vosaline_bot"  unit="psu"   /> 
    1018       <field id="sshnmasked"         name="sossheig"      unit="m"     />  
    10191134    </field_group> 
    10201135 
  • NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/cfgs/SHARED/field_def_nemo-pisces.xml

    r12178 r12928  
    6060       <field id="NFe"          long_name="Nano iron Concentration"                  unit="mmol/m3" /> 
    6161       <field id="NFe_e3t"      long_name="NFe * e3t"                                unit="mmol/m2"  > NFe * e3t </field > 
    62        <field id="NCHL"         long_name="Nano chlorophyl Concentration"            unit="mg/m3"   /> 
     62       <field id="NCHL"         long_name="Nano chlorophyl Concentration"            unit="gChl/m3"   /> 
    6363       <field id="NCHL_e3t"     long_name="NCHL * e3t"                               unit="mmol/m2"  > NCHL * e3t </field > 
    64        <field id="DCHL"         long_name="Diatoms chlorophyl Concentration"         unit="mg/m3"   /> 
     64       <field id="DCHL"         long_name="Diatoms chlorophyl Concentration"         unit="gChl/m3"   /> 
    6565       <field id="DCHL_e3t"     long_name="DCHL * e3t"                               unit="mmol/m2"  > DCHL * e3t </field > 
    6666       <field id="NO3"          long_name="Nitrate Concentration"                    unit="mmol/m3" /> 
     
    9999       <field id="PFe"          long_name="Picophytoplankton Fe biomass"             unit="mmol/m3" /> 
    100100       <field id="PFe_e3t"      long_name="PFe * e3t"                                unit="mmol/m2"  > PFe * e3t </field > 
    101        <field id="PCHL"         long_name="Picophytoplankton Chl biomass"            unit="mg/m3" /> 
     101       <field id="PCHL"         long_name="Picophytoplankton Chl biomass"            unit="gChl/m3" /> 
    102102       <field id="PCHL_e3t"     long_name="PCHL * e3t"                               unit="mmol/m2"  > PCHL * e3t </field > 
    103103 
     
    230230       <field id="Kg"          long_name="Gas transfer"                            unit="mol/m2/s/uatm"                   /> 
    231231       <field id="Dpco2"       long_name="Delta CO2"                               unit="uatm"                            /> 
     232       <field id="pCO2sea"     long_name="surface ocean pCO2"                      unit="uatm"                            /> 
    232233       <field id="Dpo2"        long_name="Delta O2"                                unit="uatm"                            /> 
    233234       <field id="Heup"        long_name="Euphotic layer depth"                    unit="m"                               /> 
     235       <field id="AtmCo2"      long_name="Atmospheric CO2 concentration"           unit="ppm"                               /> 
    234236       <field id="Irondep"     long_name="Iron deposition from dust"               unit="mol/m2/s"                        /> 
    235237       <field id="Ironsed"     long_name="Iron deposition from sediment"           unit="mol/m2/s"   grid_ref="grid_T_3D" /> 
    236       <field id="FESCAV"      long_name="Scavenging of Iron"                      unit="mmol-Fe/m3/s"   grid_ref="grid_T_3D" /> 
     238       <field id="FESCAV"      long_name="Scavenging of Iron"                      unit="mmol-Fe/m3/s"   grid_ref="grid_T_3D" /> 
    237239       <field id="FECOLL"      long_name="Colloidal Pumping of FeL"                unit="mmol-FeL/m3/s"  grid_ref="grid_T_3D" /> 
    238240       <field id="LGWCOLL"     long_name="Coagulation loss of ligands"             unit="mmol-L/m3/s"  grid_ref="grid_T_3D" /> 
     
    250252       <field id="FEZOO2"      long_name="mesozooplankton iron recycling rate"     unit="nmol-FeL/m3/s"  grid_ref="grid_T_3D" /> 
    251253 
     254       <!-- PISCES tracers trends --> 
     255       <field id="INTdtAlk"    long_name="Vertically int. of change of alkalinity"             unit="mol/m2/s"                       /> 
     256       <field id="INTdtDIC"    long_name="Vertically int. of change of dissic    "             unit="mol/m2/s"                       /> 
     257       <field id="INTdtFer"    long_name="Vertically int. of change of iron      "             unit="mol/m2/s"                       /> 
     258       <field id="INTdtDIN"    long_name="Vertically int. of change of nitrogen  "             unit="mol/m2/s"                       /> 
     259       <field id="INTdtDIP"    long_name="Vertically int. of change of phophate  "             unit="mol/m2/s"                       /> 
     260       <field id="INTdtSil"    long_name="Vertically int. of change of silicon   "             unit="mol/m2/s"                       /> 
     261 
     262 
    252263       <!-- dbio_T on T grid : variables available with diaar5 --> 
    253264       <field id="TPP"         long_name="Total Primary production of phyto"                   unit="mol/m3/s"  grid_ref="grid_T_3D" /> 
     
    257268       <field id="O2MIN"       long_name="Oxygen minimum concentration"                        unit="mol/m3"                          /> 
    258269       <field id="ZO2MIN"      long_name="Depth of oxygen minimum concentration"               unit="m"                              /> 
    259        <field id="INTNFIX"     long_name="Nitrogen fixation rate : vert. integrated"           unit="mol/m2/s"                       /> 
    260        <field id="INTPPPHYN"   long_name="Vertically integrated primary production by nanophy" unit="mol/m2/s"                       /> 
    261        <field id="INTPPPHYD"   long_name="Vertically integrated primary production by diatom"  unit="mol/m2/s"                       /> 
    262        <field id="INTPP"       long_name="Vertically integrated primary production by phyto"   unit="mol/m2/s"                       /> 
    263        <field id="INTPNEW"     long_name="Vertically integrated new primary production"        unit="mol/m2/s"                       /> 
    264        <field id="INTPBFE"     long_name="Vertically integrated of biogenic iron production"   unit="mol/m2/s"                       /> 
    265        <field id="INTPBSI"     long_name="Vertically integrated of biogenic Si production"     unit="mol/m2/s"                       /> 
    266        <field id="INTPCAL"     long_name="Vertically integrated of calcite production"         unit="mol/m2/s"                       /> 
     270       <field id="INTNFIX"     long_name="Nitrogen fixation rate : vert. integrated"           unit="mol/m2/s"  grid_ref="grid_T_vsum"  detect_missing_value="true" > Nfix * e3t </field > 
     271       <field id="INTPPPHYN"   long_name="Vertically integrated primary production by nanophy" unit="mol/m2/s"  grid_ref="grid_T_vsum"  detect_missing_value="true" > PPPHYN * e3t </field > 
     272       <field id="INTPPPHYD"   long_name="Vertically integrated primary production by diatom"  unit="mol/m2/s"  grid_ref="grid_T_vsum"  detect_missing_value="true" > PPPHYD * e3t </field >                     
     273       <field id="INTPPPHYP"   long_name="Vertically integrated primary production by picophy" unit="mol/m2/s"  grid_ref="grid_T_vsum"  detect_missing_value="true" > PPPHYP * e3t </field >                     
     274       <field id="INTPP"       long_name="Vertically integrated primary production by phyto"   unit="mol/m2/s"  grid_ref="grid_T_vsum"  detect_missing_value="true" > TPP * e3t </field >                             
     275       <field id="INTPNEW"     long_name="Vertically integrated new primary production"        unit="mol/m2/s"  grid_ref="grid_T_vsum"  detect_missing_value="true" > TPNEW * e3t </field >                 
     276       <field id="INTPBFE"     long_name="Vertically integrated of biogenic iron production"   unit="mol/m2/s"  grid_ref="grid_T_vsum"  detect_missing_value="true" > TPBFE * e3t </field >                            
     277       <field id="INTPBSI"     long_name="Vertically integrated of biogenic Si production"     unit="mol/m2/s"  grid_ref="grid_T_vsum"  detect_missing_value="true" > PBSi * e3t </field >                            
    267278 
    268279       <!-- PISCES light : variables available with key_pisces_reduced --> 
  • NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/cfgs/SHARED/grid_def_nemo.xml

    r12315 r12928  
     1<?xml version="1.0"?> 
     2<!-- 
     3============================================================================================================ 
     4= grid definition = = DO NOT CHANGE = 
     5============================================================================================================ 
     6    --> 
    17 
    28    <grid_definition>     
     
    410        <!--  --> 
    511       <grid id="grid_T_2D" > 
    6          <domain id="grid_T" /> 
     12         <domain domain_ref="grid_T" /> 
    713       </grid> 
    814        <!--  --> 
    915       <grid id="grid_T_ncatice" > 
    10          <domain id="grid_T" /> 
    11          <axis id="ncatice" /> 
     16         <domain domain_ref="grid_T" /> 
     17         <axis axis_ref="ncatice" /> 
    1218       </grid> 
    1319        <!--  --> 
    1420       <grid id="grid_T_3D" > 
    15          <domain id="grid_T" /> 
    16          <axis id="deptht" /> 
     21         <domain domain_ref="grid_T" /> 
     22         <axis axis_ref="deptht" /> 
    1723       </grid> 
    1824        <!--  --> 
    1925       <grid id="grid_T_3DS" > 
    20          <domain id="grid_T" /> 
    21          <axis id="profsed" /> 
     26         <domain domain_ref="grid_T" /> 
     27         <axis axis_ref="profsed" /> 
    2228       </grid> 
    2329        <!--  --> 
    2430       <grid id="grid_U_2D" > 
    25          <domain id="grid_U" /> 
     31         <domain domain_ref="grid_U" /> 
    2632       </grid> 
    2733        <!--  --> 
    2834       <grid id="grid_U_3D" > 
    29          <domain id="grid_U" /> 
    30          <axis id="depthu" /> 
     35         <domain domain_ref="grid_U" /> 
     36         <axis axis_ref="depthu" /> 
    3137       </grid> 
    3238        <!--  --> 
    3339       <grid id="grid_V_2D" > 
    34          <domain id="grid_V" /> 
     40         <domain domain_ref="grid_V" /> 
    3541       </grid> 
    3642        <!--  --> 
    3743       <grid id="grid_V_3D" > 
    38          <domain id="grid_V" /> 
    39          <axis id="depthv" /> 
     44         <domain domain_ref="grid_V" /> 
     45         <axis axis_ref="depthv" /> 
    4046       </grid> 
    4147        <!--  --> 
    4248       <grid id="grid_W_2D" > 
    43          <domain id="grid_W" /> 
     49         <domain domain_ref="grid_W" /> 
    4450       </grid> 
    4551        <!--  --> 
    4652       <grid id="grid_W_3D" > 
    47          <domain id="grid_W" /> 
    48          <axis id="depthw" /> 
     53         <domain domain_ref="grid_W" /> 
     54         <axis axis_ref="depthw" /> 
    4955       </grid> 
    5056       <!--  --> 
     57       <grid id="grid_1point" > 
     58         <domain domain_ref="1point"/> 
     59       </grid> 
     60        <!--  --> 
    5161       <grid id="grid_T_nfloat" > 
    52          <domain id="grid_T" /> 
    53          <axis id="nfloat" /> 
    54        </grid> 
     62         <domain domain_ref="grid_T" /> 
     63         <axis axis_ref="nfloat" /> 
     64       </grid> 
     65        <!--  --> 
     66       <grid id="grid_EqT" > 
     67         <domain domain_ref="EqT" /> 
     68       </grid> 
     69        <!--  --> 
     70 
     71 
     72       <grid id="grid_znl_T_2D"> 
     73         <domain domain_ref="gznl" /> 
     74         <axis axis_ref="basin" /> 
     75       </grid> 
     76 
     77       <grid id="grid_znl_T_3D"> 
     78         <domain domain_ref="gznl" /> 
     79         <axis axis_ref="deptht"  /> 
     80         <axis axis_ref="basin" /> 
     81       </grid> 
     82 
     83       <grid id="grid_znl_W_3D"> 
     84         <domain domain_ref="gznl" /> 
     85         <axis axis_ref="depthw"  /> 
     86         <axis axis_ref="basin" /> 
     87       </grid> 
     88 
     89      <grid id="grid_ptr_T_2D"> 
     90         <domain domain_ref="ptr" /> 
     91         <axis axis_ref="basin" /> 
     92       </grid> 
     93 
     94       <grid id="grid_ptr_T_3D"> 
     95         <domain  domain_ref="ptr" /> 
     96         <axis axis_ref="deptht"  /> 
     97         <axis axis_ref="basin" /> 
     98       </grid> 
     99 
     100       <grid id="grid_ptr_W_3D"> 
     101         <domain  domain_ref="ptr" /> 
     102         <axis axis_ref="depthw"  /> 
     103         <axis axis_ref="basin" /> 
     104       </grid> 
     105 
     106       <grid id="grid_ptr_W_GLO"> 
     107         <domain  domain_ref="ptr" /> 
     108         <axis axis_ref="depthw"  /> 
     109         <scalar> 
     110           <extract_axis position="0" /> 
     111         </scalar> 
     112       </grid> 
     113 
     114       <grid id="grid_ptr_W_ATL"> 
     115         <domain  domain_ref="ptr" /> 
     116         <axis axis_ref="depthw"  /> 
     117         <scalar> 
     118           <extract_axis position="1" /> 
     119         </scalar> 
     120       </grid> 
     121 
     122       <grid id="grid_ptr_W_IND"> 
     123         <domain  domain_ref="ptr" /> 
     124         <axis axis_ref="depthw"  /> 
     125         <scalar> 
     126           <extract_axis position="2" /> 
     127         </scalar> 
     128       </grid> 
     129 
     130       <grid id="grid_T_SFC"> 
     131         <domain domain_ref="grid_T" /> 
     132         <scalar> 
     133           <extract_axis position="0" /> 
     134         </scalar> 
     135       </grid> 
     136 
     137       <grid id="grid_T_vsum"> 
     138         <domain domain_ref="grid_T"/> 
     139         <scalar> 
     140            <reduce_axis operation="sum" /> 
     141         </scalar> 
     142       </grid> 
     143 
     144       <grid id="grid_U_vsum"> 
     145         <domain domain_ref="grid_U"/> 
     146         <scalar> 
     147            <reduce_axis operation="sum" /> 
     148         </scalar> 
     149       </grid> 
     150 
     151       <grid id="grid_V_vsum"> 
     152         <domain domain_ref="grid_V"/> 
     153         <scalar> 
     154            <reduce_axis operation="sum" /> 
     155         </scalar> 
     156       </grid> 
     157 
     158       <!-- for ORCA2 grid  --> 
     159       <grid id="cumul_U"> 
     160         <axis axis_ref="cumul_U" n_glo="182" > 
     161           <reduce_domain local="true" operation="sum" direction="jDir" /> 
     162           <reduce_axis operation="sum" /> 
     163         </axis> 
     164         <axis axis_ref="depthu" /> 
     165       </grid> 
     166 
     167       <!-- for eORCA1 grid 
     168 
     169       <grid id="cumul_U"> 
     170         <axis axis_ref="cumul_U" n_glo="362" > 
     171           <reduce_domain local="true" operation="sum" direction="jDir" /> 
     172           <reduce_axis operation="sum" /> 
     173         </axis> 
     174         <axis axis_ref="depthu" /> 
     175       </grid> 
     176 
     177      --> 
     178 
     179 
     180       <grid id="grid_T_zoom_300"> 
     181         <domain domain_ref="grid_T" /> 
     182         <axis axis_ref="deptht300" /> 
     183       </grid> 
     184 
     185       <grid id="grid_U_scalar" > 
     186         <domain domain_ref="grid_U" /> 
     187         <scalar/> 
     188       </grid> 
     189 
     190       <grid id="grid_V_scalar" > 
     191         <domain domain_ref="grid_V" /> 
     192         <scalar/> 
     193       </grid> 
     194  
     195       <grid id="grid_U_4strait"> 
     196         <domain domain_ref="grid_U" /> 
     197         <axis axis_ref="section"> 
     198           <duplicate_scalar/> 
     199         </axis> 
     200       </grid> 
     201  
     202       <grid id="grid_V_4strait"> 
     203         <domain domain_ref="grid_V" /> 
     204         <axis axis_ref="section"> 
     205           <duplicate_scalar/> 
     206         </axis> 
     207       </grid> 
     208  
     209       <grid id="grid_U_4strait_hsum"> 
     210         <scalar > 
     211           <reduce_domain operation="sum" local="true"/> 
     212           <reduce_scalar operation="sum" /> 
     213         </scalar> 
     214         <axis axis_ref="section"/> 
     215       </grid> 
     216  
     217       <grid id="grid_V_4strait_hsum"> 
     218         <scalar > 
     219           <reduce_domain operation="sum" local="true"/> 
     220           <reduce_scalar operation="sum" /> 
     221         </scalar> 
     222         <axis axis_ref="section"/> 
     223       </grid> 
     224  
     225       <grid id="grid_4strait"> 
     226         <axis axis_ref="section"/> 
     227       </grid> 
     228 
     229      <grid id="grid_U_4strait_ice"> 
     230        <domain domain_ref="grid_U" /> 
     231        <axis axis_ref="section_ice"> 
     232          <duplicate_scalar/> 
     233        </axis> 
     234      </grid> 
     235 
     236      <grid id="grid_V_4strait_ice"> 
     237        <domain domain_ref="grid_V" /> 
     238        <axis axis_ref="section_ice"> 
     239         <duplicate_scalar/> 
     240        </axis> 
     241      </grid> 
     242 
     243      <grid id="grid_U_4strait_ice_hsum"> 
     244        <scalar > 
     245         <reduce_domain operation="sum" local="true"/> 
     246         <reduce_scalar operation="sum" /> 
     247        </scalar> 
     248        <axis axis_ref="section_ice"/> 
     249      </grid> 
     250 
     251      <grid id="grid_V_4strait_ice_hsum"> 
     252         <scalar > 
     253          <reduce_domain operation="sum" local="true"/> 
     254          <reduce_scalar operation="sum" /> 
     255         </scalar> 
     256         <axis axis_ref="section_ice"/> 
     257      </grid> 
     258 
     259      <grid id="grid_4strait_ice"> 
     260        <axis axis_ref="section_ice"/> 
     261      </grid> 
     262  
    55263      <!-- scalars --> 
    56264      <grid id="grid_scalar" > 
     
    137345     </grid_definition> 
    138346 
     347      <!-- ABL grid definition --> 
     348      <grid id="grid_TA_2D"> 
     349        <domain domain_ref="grid_T" /> 
     350      </grid> 
     351      <grid id="grid_TA_3D"> 
     352        <domain domain_ref="grid_T" /> 
     353        <axis  id="ght_abl" /> 
     354      </grid> 
     355      <grid id="grid_WA_3D"> 
     356        <domain domain_ref="grid_T" /> 
     357        <axis  id="ghw_abl" /> 
     358      </grid> 
     359      <!--  --> 
     360 
     361      <!-- grid definitions for multiple-linear-regression analysis (diamlr) --> 
     362      <grid id="diamlr_grid_scalar" > 
     363   <scalar /> 
     364   <scalar /> 
     365      </grid> 
     366      <grid id="diamlr_grid_T_2D" > 
     367   <domain domain_ref="grid_T" /> 
     368   <scalar /> 
     369      </grid> 
     370      <grid id="diamlr_grid_U_2D" > 
     371   <domain domain_ref="grid_U" /> 
     372   <scalar /> 
     373      </grid> 
     374      <grid id="diamlr_grid_V_2D" > 
     375   <domain domain_ref="grid_V" /> 
     376   <scalar /> 
     377      </grid> 
     378      <grid id="diamlr_grid_W_2D" > 
     379   <domain domain_ref="grid_W" /> 
     380   <scalar /> 
     381      </grid> 
     382      <grid id="diamlr_grid_2D_to_grid_T_3D" > 
     383   <domain domain_ref="grid_T" /> 
     384        <axis axis_ref="deptht"> 
     385          <duplicate_scalar /> 
     386        </axis> 
     387      </grid> 
     388      <grid id="diamlr_grid_2D_to_grid_U_3D" > 
     389   <domain domain_ref="grid_U" /> 
     390        <axis axis_ref="depthu"> 
     391          <duplicate_scalar /> 
     392        </axis> 
     393      </grid> 
     394      <grid id="diamlr_grid_2D_to_grid_V_3D" > 
     395   <domain domain_ref="grid_V" /> 
     396        <axis axis_ref="depthv"> 
     397          <duplicate_scalar /> 
     398        </axis> 
     399      </grid> 
     400      <grid id="diamlr_grid_2D_to_grid_W_3D" > 
     401   <domain domain_ref="grid_W" /> 
     402        <axis axis_ref="depthw"> 
     403          <duplicate_scalar /> 
     404        </axis> 
     405      </grid> 
     406      <grid id="diamlr_grid_2D_to_scalar" > 
     407        <scalar> 
     408          <reduce_domain operation="average" /> 
     409        </scalar> 
     410   <scalar /> 
     411      </grid> 
     412      <!-- grid definitions for the computation of daily detided model diagnostics (diadetide) --> 
     413      <grid id="diadetide_grid_T_2D" > 
     414   <domain domain_ref="grid_T" /> 
     415   <scalar /> 
     416      </grid> 
     417      <grid id="diadetide_grid_U_2D" > 
     418   <domain domain_ref="grid_U" /> 
     419   <scalar /> 
     420      </grid> 
     421      <grid id="diadetide_grid_V_2D" > 
     422   <domain domain_ref="grid_V" /> 
     423   <scalar /> 
     424      </grid> 
     425      <grid id="diadetide_grid_2D_to_grid_T_3D" > 
     426   <domain domain_ref="grid_T" /> 
     427        <axis axis_ref="deptht"> 
     428          <duplicate_scalar /> 
     429        </axis> 
     430      </grid> 
     431      <grid id="diadetide_grid_2D_to_grid_U_3D" > 
     432   <domain domain_ref="grid_U" /> 
     433        <axis axis_ref="depthu"> 
     434          <duplicate_scalar /> 
     435        </axis> 
     436      </grid> 
     437      <grid id="diadetide_grid_2D_to_grid_V_3D" > 
     438   <domain domain_ref="grid_V" /> 
     439        <axis axis_ref="depthv"> 
     440          <duplicate_scalar /> 
     441        </axis> 
     442      </grid> 
     443      <grid id="diadetide_grid_2D_to_grid_W_3D" > 
     444   <domain domain_ref="grid_W" /> 
     445        <axis axis_ref="depthw"> 
     446          <duplicate_scalar /> 
     447        </axis> 
     448      </grid> 
     449 
     450    </grid_definition> 
  • NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/cfgs/SHARED/namelist_ice_ref

    r12178 r12928  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! SI3 namelist:   
     2!! SI3 Reference namelist: 
    33!!              1 - Generic parameters                 (nampar) 
    44!!              2 - Ice thickness discretization       (namitd) 
     
    6868   ln_str_H79       = .true.          !  ice strength param.: Hibler_79   => P = pstar*<h>*exp(-c_rhg*A)                       
    6969      rn_pstar      =   2.0e+04       !     ice strength thickness parameter [N/m2] 
    70       rn_crhg       =   20.0          !     ice strength conc. parameter (-) 
     70      rn_crhg       =  20.0           !     ice strength conc. parameter (-) 
    7171                   ! -- ice_rdgrft -- ! 
    7272   rn_csrdg         =   0.5           !  fraction of shearing energy contributing to ridging 
  • NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/cfgs/SHARED/namelist_pisces_ref

    r12178 r12928  
    3434!              !  file name   ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    3535!              !              !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    36    sn_patm     = 'presatm'    ,     -1.           , 'patm'    ,  .true.      , .true. , 'yearly'  , ''       , ''       , '' 
    37    sn_atmco2   = 'presatmco2' ,     -1.           , 'xco2'    ,  .true.      , .true. , 'yearly'  , ''       , ''       , '' 
     36   sn_patm     = 'presatm'    ,     -1            , 'patm'    ,  .true.      , .true. , 'yearly'  , ''       , ''       , '' 
     37   sn_atmco2   = 'presatmco2' ,     -1            , 'xco2'    ,  .true.      , .true. , 'yearly'  , ''       , ''       , '' 
    3838   cn_dir      = './'     !  root directory for the location of the dynamical files 
    3939! 
     
    141141!              !  file name       ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    142142!              !                  !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    143    sn_par      = 'par.orca'       ,     24.           , 'fr_par'  ,  .true.      , .true. , 'yearly'  , ''       , ''       , '' 
     143   sn_par      = 'par.orca'       ,     24            , 'fr_par'  ,  .true.      , .true. , 'yearly'  , ''       , ''       , '' 
    144144   cn_dir      = './'      !  root directory for the location of the dynamical files 
    145145   ln_varpar   =  .true.   ! boolean for PAR variable 
     
    343343/ 
    344344!----------------------------------------------------------------------- 
    345 &nampissbc     !   parameters for inputs deposition 
     345&nampisbc     !   parameters for inputs deposition 
    346346!----------------------------------------------------------------------- 
    347347!              !  file name        ! frequency (hours) ! variable      ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    348348!              !                   !  (if <0  months)  !   name        !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    349    sn_dust     = 'dust.orca'       ,     -1.           , 'dust'        ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    350    sn_solub    = 'solubility.orca' ,    -12.           , 'solubility1' ,  .false.     , .true. , 'yearly'  , ''       , ''    , '' 
    351    sn_riverdic = 'river.orca'      ,    120.           , 'riverdic'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    352    sn_riverdoc = 'river.orca'      ,    120.           , 'riverdoc'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    353    sn_riverdin = 'river.orca'      ,    120.           , 'riverdin'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    354    sn_riverdon = 'river.orca'      ,    120.           , 'riverdon'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    355    sn_riverdip = 'river.orca'      ,    120.           , 'riverdip'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    356    sn_riverdop = 'river.orca'      ,    120.           , 'riverdop'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    357    sn_riverdsi = 'river.orca'      ,    120.           , 'riverdsi'    ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
    358    sn_ndepo    = 'ndeposition.orca',    -12.           , 'ndep'        ,  .false.     , .true. , 'yearly'  , ''       , ''    , '' 
    359    sn_ironsed  = 'bathy.orca'      ,    -12.           , 'bathy'       ,  .false.     , .true. , 'yearly'  , ''       , ''    , '' 
    360    sn_hydrofe  = 'hydrofe.orca'    ,    -12.           , 'epsdb'       ,  .false.     , .true. , 'yearly'  , ''       , ''    , '' 
     349   sn_dust     = 'dust.orca'       ,     -1            , 'dust'        ,  .true.      , .true. , 'yearly'  , ''       , ''    , '' 
     350   sn_ironsed  = 'bathy.orca'      ,    -12            , 'bathy'       ,  .false.     , .true. , 'yearly'  , ''       , ''    , '' 
     351   sn_hydrofe  = 'hydrofe.orca'    ,    -12            , 'epsdb'       ,  .false.     , .true. , 'yearly'  , ''       , ''    , '' 
    361352! 
    362353   cn_dir      = './'      !  root directory for the location of the dynamical files 
    363    ln_dust     =  .true.   ! boolean for dust input from the atmosphere 
    364    ln_solub    =  .true.   ! boolean for variable solubility of atm. Iron 
    365    ln_river    =  .true.   ! boolean for river input of nutrients 
    366    ln_ndepo    =  .true.   ! boolean for atmospheric deposition of N 
    367    ln_ironsed  =  .true.   ! boolean for Fe input from sediments 
    368    ln_ironice  =  .true.   ! boolean for Fe input from sea ice 
    369    ln_hydrofe  =  .true.   ! boolean for from hydrothermal vents 
     354   ln_ironsed  =  .false.   ! boolean for Fe input from sediments 
     355   ln_ironice  =  .false.   ! boolean for Fe input from sea ice 
     356   ln_hydrofe  =  .false.   ! boolean for from hydrothermal vents 
    370357   sedfeinput  =  2.e-9    ! Coastal release of Iron 
    371358   distcoast   =  5.e3     ! Distance off the coast for Iron from sediments 
    372    dustsolub   =  0.02     ! Solubility of the dusta 
    373359   mfrac       =  0.035    ! Fe mineral fraction of dust 
    374360   wdust       =  2.0      ! Dust sinking speed  
    375361   icefeinput  =  15.e-9   ! Iron concentration in sea ice 
     362   hratio      =  1.e+7    ! Fe to 3He ratio assumed for vent iron supply  
     363!                          ! ln_ligand 
     364   lgw_rath    =  0.5      ! Weak ligand ratio from sed hydro sources  
     365/ 
     366!----------------------------------------------------------------------- 
     367&nampissed     !   parameters for sediments mobilization 
     368!----------------------------------------------------------------------- 
    376369   nitrfix     =  1.e-7    ! Nitrogen fixation rate 
    377370   diazolight  =  50.      ! Diazotrophs sensitivity to light (W/m2) 
    378371   concfediaz  =  1.e-10   ! Diazotrophs half-saturation Cste for Iron 
    379    hratio      =  1.e+7    ! Fe to 3He ratio assumed for vent iron supply  
    380 !                          ! ln_ligand 
    381    lgw_rath    =  0.5      ! Weak ligand ratio from sed hydro sources  
    382372/ 
    383373!----------------------------------------------------------------------- 
  • NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/cfgs/SHARED/namelist_ref

    r12314 r12928  
    55!! namelists    2 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_cpl, 
    66!!                                    namsbc_sas, namtra_qsr, namsbc_rnf, 
    7 !!                                    namsbc_isf, namsbc_iscpl, namsbc_apr,  
     7!!                                    namisf, namsbc_apr,  
    88!!                                    namsbc_ssr, namsbc_wave, namberg) 
    99!!              3 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 
     
    4242   nn_leapy    =       0   !  Leap year calendar (1) or not (0) 
    4343   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T) 
    44       nn_euler    =    1      !  = 0 : start with forward time step if ln_rstart=T 
    45       nn_rstctl   =    0      !  restart control ==> activated only if ln_rstart=T 
     44      ln_1st_euler = .false.  !  =T force a start with forward time step (ln_rstart=T) 
     45      nn_rstctl    =    0     !  restart control ==> activated only if ln_rstart=T 
    4646      !                          !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
    4747      !                          !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 
     
    5151      cn_ocerst_out   = "restart"   !  suffix of ocean restart name (output) 
    5252      cn_ocerst_outdir = "."        !  directory in which to write output ocean restarts 
    53    ln_iscpl    = .false.   !  cavity evolution forcing or coupling to ice sheet model 
    5453   nn_istate   =       0   !  output the initial state (1) or not (0) 
    5554   ln_rst_list = .false.   !  output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F) 
     
    6564   ln_clobber  = .true.    !  clobber (overwrite) an existing file 
    6665   nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 
    67    ln_xios_read = .FALSE.  !  use XIOS to read restart file (only for a single file restart) 
     66   ln_xios_read = .false.  !  use XIOS to read restart file (only for a single file restart) 
    6867   nn_wxios = 0      !  use XIOS to write restart file 0 - no, 1 - single file output, 2 - multiple file output 
    6968/ 
     
    7271!----------------------------------------------------------------------- 
    7372   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    74    rn_isfhmin  =    1.00   !  treshold [m] to discriminate grounding ice from floating ice 
    75    ! 
    76    rn_rdt      = 5400.     !  time step for the dynamics and tracer 
     73   ! 
     74   rn_Dt       = 5400.     !  time step for the dynamics and tracer 
    7775   rn_atfp     =    0.1    !  asselin time filter parameter 
    7876   ! 
    7977   ln_crs      = .false.   !  Logical switch for coarsening module      (T => fill namcrs) 
    8078   ! 
    81    ln_meshmask = .false.   !  =T create a mesh file 
     79   ln_meshmask = .true.   !  =T create a mesh file 
    8280/ 
    8381!----------------------------------------------------------------------- 
    8482&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg) 
    8583!----------------------------------------------------------------------- 
    86    ln_read_cfg = .false.   !  (=T) read the domain configuration file 
    87       !                    !  (=F) user defined configuration           (F => create/check namusr_def) 
     84   ln_read_cfg = .false.     !  (=T) read the domain configuration file 
     85      !                      !  (=F) user defined configuration           (F => create/check namusr_def) 
    8886      cn_domcfg = "domain_cfg"  ! domain configuration filename 
    8987      ! 
    90       ln_closea    = .false.    !  T => keep closed seas (defined by closea_mask field) in the   
    91       !                         !       domain and apply special treatment of freshwater fluxes. 
    92       !                         !  F => suppress closed seas (defined by closea_mask field)  
    93       !                         !       from the bathymetry at runtime. 
    94       !                         !  If closea_mask field doesn't exist in the domain_cfg file 
    95       !                         !       then this logical does nothing. 
    96    ln_write_cfg = .false.   !  (=T) create the domain configuration file 
     88      ln_closea    = .false. !  (=T => fill namclo)  
     89      !                      !  (=F) no control of net precip/evap over closed sea 
     90      ! 
     91   ln_write_cfg = .false.    !  (=T) create the domain configuration file 
    9792      cn_domcfg_out = "domain_cfg_out" ! newly created domain configuration filename 
    9893      ! 
    99    ln_use_jattr = .false.  !  use (T) the file attribute: open_ocean_jstart, if present 
    100    !                       !  in netcdf input files, as the start j-row for reading 
     94   ln_use_jattr = .false.    !  use (T) the file attribute: open_ocean_jstart, if present 
     95   !                         !  in netcdf input files, as the start j-row for reading 
     96/ 
     97!----------------------------------------------------------------------- 
     98&namclo        !   parameters of the closed sea (cs) behavior                (default: OFF) 
     99!----------------------------------------------------------------------- 
     100   ln_maskcs = .false.        ! (=T) cs are masked ; So, in this case ln_mask_csundef and ln_clo_rnf have no effect. 
     101      !                       ! (=F => set ln_mask_csundef and ln_clo_rnf)  
     102      !                       ! cs masks are read and net evap/precip over closed sea spread out depending on domain_cfg.nc masks. 
     103      !                       ! See ln_mask_csundef and ln_clo_rnf for specific option related to this case 
     104      ! 
     105      ln_mask_csundef = .true.   ! (=T) undefined closed seas are masked ;  
     106      !                          ! (=F) undefined closed seas are kept and no specific treatment is done for these closed seas 
     107      ! 
     108      ln_clo_rnf = .true.        ! (=T) river mouth specified in domain_cfg.nc masks (rnf and emp case) are added to the runoff mask. 
     109      !                          !      allow the treatment of closed sea outflow grid-points to be the same as river mouth grid-points 
    101110/ 
    102111!----------------------------------------------------------------------- 
     
    106115   ln_tsd_init = .false.         !  ocean initialisation 
    107116   ln_tsd_dmp  = .false.         !  T-S restoring   (see namtra_dmp) 
    108     
     117 
    109118   cn_dir      = './'      !  root directory for the T-S data location 
    110119   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     
    184193!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T) 
    185194!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T) 
    186 !!   namsbc_isf      ice shelf melting/freezing                         (ln_isfcav  =T : read (ln_read_cfg=T) or set or usr_def_zgr ) 
    187 !!   namsbc_iscpl    coupling option between land ice model and ocean   (ln_isfcav  =T) 
    188195!!   namsbc_wave     external fields from wave model                    (ln_wave    =T) 
    189196!!   namberg         iceberg floats                                     (ln_icebergs=T) 
     
    195202   nn_fsbc     = 2         !  frequency of SBC module call 
    196203      !                    !  (control sea-ice & iceberg model call) 
    197                      ! Type of air-sea fluxes  
     204                     ! Type of air-sea fluxes 
    198205   ln_usr      = .false.   !  user defined formulation                  (T => check usrdef_sbc) 
    199206   ln_flx      = .false.   !  flux formulation                          (T => fill namsbc_flx ) 
    200207   ln_blk      = .false.   !  Bulk formulation                          (T => fill namsbc_blk ) 
     208   ln_abl      = .false.   !  ABL  formulation                          (T => fill namsbc_abl ) 
    201209      !              ! Type of coupling (Ocean/Ice/Atmosphere) : 
    202210   ln_cpl      = .false.   !  atmosphere coupled   formulation          ( requires key_oasis3 ) 
     
    205213      !                    !  =0 no opa-sas OASIS coupling: default single executable config. 
    206214      !                    !  =1 opa-sas OASIS coupling: multi executable config., OPA component 
    207       !                    !  =2 opa-sas OASIS coupling: multi executable config., SAS component  
     215      !                    !  =2 opa-sas OASIS coupling: multi executable config., SAS component 
    208216                     ! Sea-ice : 
    209    nn_ice      = 0         !  =0 no ice boundary condition     
     217   nn_ice      = 0         !  =0 no ice boundary condition 
    210218      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif ) 
    211       !                    !  =2 or 3 automatically for SI3 or CICE    ("key_si3" or "key_cice") 
    212       !                    !          except in AGRIF zoom where it has to be specified 
     219      !                    !  =2 or 3 for SI3 and CICE, respectively 
    213220   ln_ice_embd = .false.   !  =T embedded sea-ice (pressure + mass and salt exchanges) 
    214221      !                    !  =F levitating ice (no pressure, mass and salt exchanges) 
    215                      ! Misc. options of sbc :  
     222                     ! Misc. options of sbc : 
    216223   ln_traqsr   = .false.   !  Light penetration in the ocean            (T => fill namtra_qsr) 
    217224   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave 
     
    222229   ln_rnf      = .false.   !  runoffs                                   (T => fill namsbc_rnf) 
    223230   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr ) 
    224    ln_isf      = .false.   !  ice shelf                                 (T   => fill namsbc_isf & namsbc_iscpl) 
    225231   ln_wave     = .false.   !  Activate coupling with wave  (T => fill namsbc_wave) 
    226232   ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave) 
    227    ln_sdw      = .false.   !  Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave)  
     233   ln_sdw      = .false.   !  Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave) 
    228234   nn_sdrift   =  0        !  Parameterization for the calculation of 3D-Stokes drift from the surface Stokes drift 
    229235      !                    !   = 0 Breivik 2015 parameterization: v_z=v_0*[exp(2*k*z)/(1-8*k*z)] 
     
    250256/ 
    251257!----------------------------------------------------------------------- 
    252 &namsbc_blk    !   namsbc_blk  generic Bulk formula                     (ln_blk =T) 
     258&namsbc_blk    !   namsbc_blk  generic Bulk formula          (ln_blk =T) 
    253259!----------------------------------------------------------------------- 
    254260   !                    !  bulk algorithm : 
    255    ln_NCAR     = .false.   ! "NCAR"      algorithm   (Large and Yeager 2008) 
     261   ln_NCAR      = .true.    ! "NCAR"      algorithm   (Large and Yeager 2008) 
    256262   ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003) 
    257    ln_COARE_3p5 = .false.   ! "COARE 3.5" algorithm   (Edson et al. 2013) 
    258    ln_ECMWF    = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
     263   ln_COARE_3p6 = .false.   ! "COARE 3.6" algorithm   (Edson et al. 2013) 
     264   ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 45r1) 
    259265 
    260266   ln_humi_dpt    = .false.   !  Supply dewpoint tempearture instead of specific humidity (is true for ERA5) 
    261267 
    262268      ! 
    263       rn_zqt      = 10.       !  Air temperature & humidity reference height (m) 
    264       rn_zu       = 10.       !  Wind vector reference height (m) 
    265       ln_Cd_L12   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2012) 
    266       ln_Cd_L15   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2015) 
    267       ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data 
    268       rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    269       rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
    270       rn_vfac     = 0.        !  multiplicative factor for ocean & ice velocity used to 
    271       !                       !  calculate the wind stress (0.=absolute or 1.=relative winds) 
    272  
     269      rn_zqt     = 10.      !  Air temperature & humidity reference height (m) 
     270      rn_zu      = 10.      !  Wind vector reference height (m) 
     271      ln_Cd_L12  = .false.  !  air-ice drags = F(ice conc.) (Lupkes et al. 2012) 
     272      ln_Cd_L15  = .false.  !  air-ice drags = F(ice conc.) (Lupkes et al. 2015) 
     273      !                     !  - module of the mean stress" data 
     274      rn_pfac    = 1.       !  multipl. factor for precipitation (total & snow) 
     275      rn_efac    = 1.       !  multipl. factor for evaporation (0. or 1.) 
     276      rn_vfac    = 0.       !  multipl. factor for ocean & ice velocity 
     277      !                     !  used to calculate the wind stress 
     278      !                     ! (0. => absolute or 1. => relative winds) 
     279      ln_skin_cs = .false.  !  use the cool-skin parameterization 
     280      ln_skin_wl = .false.  !  use the warm-layer parameterization 
     281      !                     !   ==> only available in ECMWF and COARE algorithms 
     282      ln_humi_sph = .true.  !  humidity "sn_humi" is specific humidity  [kg/kg] 
     283      ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K] 
     284      ln_humi_rlh = .false. !  humidity "sn_humi" is relative humidity     [%] 
     285   ! 
    273286   cn_dir      = './'      !  root directory for the bulk data location 
    274287   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________! 
     
    281294   sn_tair     = 't_10.15JUNE2009_fill'       ,    6.        , 'T_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    282295   sn_humi     = 'q_10.15JUNE2009_fill'       ,    6.        , 'Q_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     296   sn_hpgi     = 'NOT USED'                   ,   24.        , 'uhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'UG'     , '' 
     297   sn_hpgj     = 'NOT USED'                   ,   24.        , 'vhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'VG'     , '' 
    283298   sn_prec     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'PRC_MOD1',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    284299   sn_snow     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'SNOW'    ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    285300   sn_slp      = 'slp.15JUNE2009_fill'        ,    6.        , 'SLP'     ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    286    sn_tdif     = 'taudif_core'                ,   24         , 'taudif'  ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     301/ 
     302!----------------------------------------------------------------------- 
     303&namsbc_abl    !   Atmospheric Boundary Layer formulation           (ln_abl = T) 
     304!----------------------------------------------------------------------- 
     305   cn_dir           = './'      !  root directory for the location of the ABL grid file 
     306   cn_dom           = 'dom_cfg_abl.nc' 
     307 
     308   cn_ablrst_in     = "restart_abl"   !  suffix of abl restart name (input) 
     309   cn_ablrst_out    = "restart_abl"   !  suffix of abl restart name (output) 
     310   cn_ablrst_indir  = "."             !  directory to read   input abl restarts 
     311   cn_ablrst_outdir = "."             !  directory to write output abl restarts 
     312 
     313   ln_hpgls_frc   = .false. 
     314   ln_geos_winds  = .false. 
     315   nn_dyn_restore = 2         ! restoring option for dynamical ABL variables: = 0 no restoring 
     316                              !                                               = 1 equatorial restoring 
     317                              !                                               = 2 global restoring 
     318   rn_ldyn_min   =  4.5       !  magnitude of the nudging on ABL dynamics at the bottom of the ABL   [hour] 
     319   rn_ldyn_max   =  1.5       !  magnitude of the nudging on ABL dynamics at the top of the ABL   [hour] 
     320   rn_ltra_min   =  4.5       !  magnitude of the nudging on ABL tracers  at the bottom of the ABL   [hour] 
     321   rn_ltra_max   =  1.5       !  magnitude of the nudging on ABL tracers  at the top of the ABL   [hour] 
     322   nn_amxl       =  0         ! mixing length: = 0 Deardorff 80 length-scale 
     323                              !                = 1 length-scale based on the distance to the PBL height 
     324                              !                = 2 Bougeault & Lacarrere 89 length-scale 
     325   rn_Cm         = 0.0667     ! 0.126 in MesoNH 
     326   rn_Ct         = 0.1667     ! 0.143 in MesoNH 
     327   rn_Ce         = 0.4        ! 0.4   in MesoNH 
     328   rn_Ceps       = 0.7        ! 0.85  in MesoNH 
     329   rn_Rod        = 0.15       ! c0 in RMCA17 mixing length formulation (not yet implemented) 
     330   rn_Ric        = 0.139      !  Critical Richardson number (to compute PBL height and diffusivities) 
    287331/ 
    288332!----------------------------------------------------------------------- 
     
    293337   !                       !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 
    294338   nn_cats_cpl   =     5   !  Number of sea ice categories over which coupling is to be carried out (if not 1) 
    295  
    296339   !_____________!__________________________!____________!_____________!______________________!________! 
    297340   !             !        description       !  multiple  !    vector   !       vector         ! vector ! 
     
    330373   sn_rcv_wper   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
    331374   sn_rcv_wnum   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
    332    sn_rcv_wstrf  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     375   sn_rcv_wfreq  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
    333376   sn_rcv_wdrag  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
    334377   sn_rcv_ts_ice =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     
    379422   nn_chldta   =      0       !  RGB : Chl data (=1) or cst value (=0) 
    380423   rn_si1      =   23.0       !  2BD : longest depth of extinction 
    381     
     424 
    382425   cn_dir      = './'      !  root directory for the chlorophyl data location 
    383426   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     
    396439      ln_sssr_bnd =  .true.   !  flag to bound erp term (associated with nn_sssr=2) 
    397440      rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day] 
     441      nn_sssr_ice =   1       ! control of sea surface restoring under sea-ice 
     442                              ! 0 = no restoration under ice : * (1-icefrac) 
     443                              ! 1 = restoration everywhere  
     444                              ! >1 = enhanced restoration under ice : 1+(nn_icedmp-1)*icefrac 
    398445 
    399446   cn_dir      = './'      !  root directory for the SST/SSS data location 
     
    443490/ 
    444491!----------------------------------------------------------------------- 
    445 &namsbc_isf    !  Top boundary layer (ISF)                              (ln_isfcav =T : read (ln_read_cfg=T)  
    446 !-----------------------------------------------------------------------             or set or usr_def_zgr ) 
    447    !                 ! type of top boundary layer  
    448    nn_isf      = 1         !  ice shelf melting/freezing 
    449                            !  1 = presence of ISF   ;  2 = bg03 parametrisation  
    450                            !  3 = rnf file for ISF  ;  4 = ISF specified freshwater flux 
    451                            !  options 1 and 4 need ln_isfcav = .true. (domzgr) 
    452       !              !  nn_isf = 1 or 2 cases: 
    453       rn_gammat0  = 1.e-4     ! gammat coefficient used in blk formula 
    454       rn_gammas0  = 1.e-4     ! gammas coefficient used in blk formula 
    455       !              !  nn_isf = 1 or 4 cases: 
    456       rn_hisf_tbl =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
    457       !                       ! 0 => thickness of the tbl = thickness of the first wet cell 
    458       !              ! nn_isf = 1 case 
    459       nn_isfblk   = 1         ! 1 ISOMIP  like: 2 equations formulation (Hunter et al., 2006) 
    460       !                       ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 
    461       nn_gammablk = 1         ! 0 = cst Gammat (= gammat/s) 
    462       !                       ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
    463       !                       ! 2 = velocity and stability dependent Gamma    (Holland et al. 1999) 
    464  
    465    !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
    466    !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    467    !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
    468 !* nn_isf = 4 case 
    469    sn_fwfisf   = 'rnfisf'    ,         -12.      ,'sowflisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    470 !* nn_isf = 3 case 
    471    sn_rnfisf   = 'rnfisf'    ,         -12.      ,'sofwfisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    472 !* nn_isf = 2 and 3 cases  
    473    sn_depmax_isf ='rnfisf'   ,         -12.      ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    474    sn_depmin_isf ='rnfisf'   ,         -12.      ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    475 !* nn_isf = 2 case 
    476    sn_Leff_isf = 'rnfisf'    ,         -12.      ,'Leff'     ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    477 / 
    478 !----------------------------------------------------------------------- 
    479 &namsbc_iscpl  !   land ice / ocean coupling option                     (ln_isfcav =T : read (ln_read_cfg=T)  
    480 !-----------------------------------------------------------------------             or set or usr_def_zgr ) 
    481    nn_drown    = 10        ! number of iteration of the extrapolation loop (fill the new wet cells) 
    482    ln_hsb      = .false.   ! activate conservation module (conservation exact after a time of rn_fiscpl) 
    483    nn_fiscpl   = 43800     ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency) 
     492&namisf       !  Top boundary layer (ISF)                               (default: OFF) 
     493!----------------------------------------------------------------------- 
     494   ! 
     495   ! ---------------- ice shelf load ------------------------------- 
     496   ! 
     497   cn_isfload = 'uniform'      ! scheme to compute ice shelf load (ln_isfcav = .true. in domain_cfg.nc) 
     498      rn_isfload_T = -1.9 
     499      rn_isfload_S = 34.4 
     500   ! 
     501   ! ---------------- ice shelf melt formulation ------------------------------- 
     502   ! 
     503   ln_isf = .false.           ! activate ice shelf module 
     504      ln_isfdebug = .false.      ! add debug print in ISF code (global min/max/sum of specific variable) 
     505      cn_isfdir   = './'         ! directory for all ice shelf input file 
     506      ! 
     507      ! ---------------- cavities opened ------------------------------- 
     508      ! 
     509      ln_isfcav_mlt = .false.    ! ice shelf melting into the cavity (need ln_isfcav = .true. in domain_cfg.nc) 
     510         cn_isfcav_mlt = '3eq'   ! ice shelf melting formulation (spe/2eq/3eq/oasis) 
     511         !                       ! spe = fwfisf is read from a forcing field 
     512         !                       ! 2eq = ISOMIP  like: 2 equations formulation (Hunter et al., 2006 for a short description) 
     513         !                       ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2016 for a short description) 
     514         !                       ! oasis = fwfisf is given by oasis and pattern by file sn_isfcav_fwf 
     515         !              !  cn_isfcav_mlt = 2eq or 3eq cases: 
     516         cn_gammablk = 'vel'     ! scheme to compute gammat/s (spe,ad15,hj99) 
     517         !                       ! spe      = constant transfert velocity (rn_gammat0, rn_gammas0) 
     518         !                       ! vel      = velocity dependent transfert velocity (u* * gammat/s) (Asay-Davis et al. 2016 for a short description) 
     519         !                       ! vel_stab = velocity and stability dependent transfert coeficient (Holland et al. 1999 for a complete description) 
     520         rn_gammat0  = 1.4e-2    ! gammat coefficient used in spe, vel and vel_stab gamma computation method 
     521         rn_gammas0  = 4.0e-4    ! gammas coefficient used in spe, vel and vel_stab gamma computation method 
     522         ! 
     523         rn_htbl     =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
     524         !                       ! 0 => thickness of the tbl = thickness of the first wet cell 
     525         ! 
     526         !* 'spe' and 'oasis' case 
     527         !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
     528         !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     529         !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     530         sn_isfcav_fwf = 'isfmlt_cav',      -12.      , 'fwflisf'  ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     531      ! 
     532      ! ---------------- cavities parametrised ------------------------------- 
     533      ! 
     534      ln_isfpar_mlt = .false.   ! ice shelf melting parametrised 
     535         cn_isfpar_mlt = 'spe'  ! ice shelf melting parametrisation (spe/bg03/oasis) 
     536         !                      ! spe   = fwfisf is read from a forcing field 
     537         !                      ! bg03  = melt computed using Beckmann and Goosse parametrisation 
     538         !                      ! oasis = fwfisf is given by oasis and pattern by file sn_isfpar_fwf 
     539         ! 
     540         !* all cases 
     541         !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
     542         !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     543         !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     544         sn_isfpar_zmax = 'isfmlt_par',       0        ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     545         sn_isfpar_zmin = 'isfmlt_par',       0        ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     546         !* 'spe' and 'oasis' case 
     547         sn_isfpar_fwf = 'isfmlt_par' ,      -12.      ,'sofwfisf' ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
     548         !* 'bg03' case 
     549         sn_isfpar_Leff = 'isfmlt_par',       0.       ,'Leff'     ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
     550      ! 
     551      ! ---------------- ice sheet coupling ------------------------------- 
     552      ! 
     553      ln_isfcpl = .false. 
     554         nn_drown       = 10        ! number of iteration of the extrapolation loop (fill the new wet cells) 
     555         ln_isfcpl_cons = .false. 
    484556/ 
    485557!----------------------------------------------------------------------- 
     
    506578   ln_icebergs = .false.      ! activate iceberg floats (force =F with "key_agrif") 
    507579   ! 
     580   !                          ! restart 
     581   cn_icbrst_in     = "restart_icb" !  suffix of iceberg restart name (input) 
     582   cn_icbrst_indir  = "./"          !  directory from which to read input ocean restarts 
     583   cn_icbrst_out    = "restart_icb" !  suffix of ocean restart name (output) 
     584   cn_icbrst_outdir = "./"          !  directory from which to read output ocean restarts           
     585   ! 
    508586   !                          ! diagnostics: 
    509587   ln_bergdia        = .true.        ! Calculate budgets 
     
    563641&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    564642!----------------------------------------------------------------------- 
     643   ln_agrif_2way = .true.  !  activate two way nesting 
    565644   ln_spc_dyn    = .true.  !  use 0 as special value for dynamics 
    566645   rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [m2/s] 
    567646   rn_sponge_dyn = 2880.   !  coefficient for dynamics sponge layer [m2/s] 
     647   rn_trelax_tra = 0.01    !  inverse of relaxation time (in steps) for tracers [] 
     648   rn_trelax_dyn = 0.01    !  inverse of relaxation time (in steps) for dynamics [] 
    568649   ln_chk_bathy  = .false. !  =T  check the parent bathymetry 
    569650/ 
     
    572653!----------------------------------------------------------------------- 
    573654   ln_tide     = .false.      ! Activate tides 
    574       ln_tide_pot   = .true.                !  use tidal potential forcing 
     655      nn_tide_var   = 1          !  Variant of tidal parameter set and tide-potential computation 
     656      !                          !     (1: default; 0: compatibility with previous versions) 
     657      ln_tide_dia   = .false.    !  Enable tidal diagnostic output 
     658      ln_tide_pot   = .false.               !  use tidal potential forcing 
     659         rn_tide_gamma = 0.7                   ! Tidal tilt factor 
    575660         ln_scal_load  = .false.               ! Use scalar approximation for 
    576661            rn_scal_load = 0.094               !     load potential 
    577662         ln_read_load  = .false.               ! Or read load potential from file 
    578663            cn_tide_load = 'tide_LOAD_grid_T.nc'  ! filename for load potential 
    579             !       
     664            ! 
    580665      ln_tide_ramp  = .false.               !  Use linear ramp for tides at startup 
    581          rdttideramp   =    0.                 !  ramp duration in days 
    582       clname(1)    = 'DUMMY'               !  name of constituent - all tidal components must be set in namelist_cfg 
     666         rn_tide_ramp_dt = 0.               !  ramp duration in days 
     667      sn_tide_cnames(1) = 'DUMMY'               !  name of constituent - all tidal components must be set in namelist_cfg 
    583668/ 
    584669!----------------------------------------------------------------------- 
     
    656741   filtide          = 'bdydta/amm12_bdytide_'   !  file name root of tidal forcing files 
    657742   ln_bdytide_2ddta = .false.                   ! 
    658    ln_bdytide_conj  = .false.                   !  
    659743/ 
    660744 
     
    683767!----------------------------------------------------------------------- 
    684768   rn_Cd0      =  1.e-3    !  drag coefficient [-] 
    685    rn_Uc0      =  0.4      !  ref. velocity [m/s] (linear drag=Cd0*Uc0)  
     769   rn_Uc0      =  0.4      !  ref. velocity [m/s] (linear drag=Cd0*Uc0) 
    686770   rn_Cdmax    =  0.1      !  drag value maximum [-] (logarithmic drag) 
    687771   rn_ke0      =  2.5e-3   !  background kinetic energy  [m2/s2] (non-linear cases) 
     
    694778!----------------------------------------------------------------------- 
    695779   rn_Cd0      =  1.e-3    !  drag coefficient [-] 
    696    rn_Uc0      =  0.4      !  ref. velocity [m/s] (linear drag=Cd0*Uc0)  
     780   rn_Uc0      =  0.4      !  ref. velocity [m/s] (linear drag=Cd0*Uc0) 
    697781   rn_Cdmax    =  0.1      !  drag value maximum [-] (logarithmic drag) 
    698782   rn_ke0      =  2.5e-3   !  background kinetic energy  [m2/s2] (non-linear cases) 
     
    761845      nn_cen_v   =  4            !  =2/4, vertical   2nd order CEN / 4th order COMPACT 
    762846   ln_traadv_fct = .false. !  FCT scheme 
    763       nn_fct_h   =  2            !  =2/4, horizontal 2nd / 4th order  
    764       nn_fct_v   =  2            !  =2/4, vertical   2nd / COMPACT 4th order  
     847      nn_fct_h   =  2            !  =2/4, horizontal 2nd / 4th order 
     848      nn_fct_v   =  2            !  =2/4, vertical   2nd / COMPACT 4th order 
    765849   ln_traadv_mus = .false. !  MUSCL scheme 
    766850      ln_mus_ups = .false.       !  use upstream scheme near river mouths 
     
    783867   ln_traldf_triad = .false.   !  iso-neutral (triad    operator) 
    784868   ! 
    785    !                       !  iso-neutral options:         
     869   !                       !  iso-neutral options: 
    786870   ln_traldf_msc   = .false.   !  Method of Stabilizing Correction      (both operators) 
    787871   rn_slpmax       =  0.01     !  slope limit                           (both operators) 
     
    793877   nn_aht_ijk_t    = 0         !  space/time variation of eddy coefficient: 
    794878      !                             !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
    795       !                             !   =  0           constant  
    796       !                             !   = 10 F(k)      =ldf_c1d  
    797       !                             !   = 20 F(i,j)    =ldf_c2d  
     879      !                             !   =  0           constant 
     880      !                             !   = 10 F(k)      =ldf_c1d 
     881      !                             !   = 20 F(i,j)    =ldf_c2d 
    798882      !                             !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    799883      !                             !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
    800884      !                             !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 
    801       !                        !  time invariant coefficients:  aht0 = 1/2  Ud*Ld   (lap case)  
     885      !                        !  time invariant coefficients:  aht0 = 1/2  Ud*Ld   (lap case) 
    802886      !                             !                           or   = 1/12 Ud*Ld^3 (blp case) 
    803887      rn_Ud        = 0.01           !  lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) 
     
    825909      nn_aei_ijk_t    = 0           !  space/time variation of eddy coefficient: 
    826910      !                             !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
    827       !                             !   =  0           constant  
    828       !                             !   = 10 F(k)      =ldf_c1d  
    829       !                             !   = 20 F(i,j)    =ldf_c2d  
     911      !                             !   =  0           constant 
     912      !                             !   = 10 F(k)      =ldf_c1d 
     913      !                             !   = 20 F(i,j)    =ldf_c2d 
    830914      !                             !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    831915      !                             !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
    832       !                        !  time invariant coefficients:  aei0 = 1/2  Ue*Le  
     916      !                        !  time invariant coefficients:  aei0 = 1/2  Ue*Le 
    833917      rn_Ue        = 0.02           !  lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) 
    834918      rn_Le        = 200.e+3        !  lateral diffusive length   [m]   (nn_aht_ijk_t= 0, 10) 
     
    870954   rn_lf_cutoff  =  5.0             !  cutoff frequency for low-pass filter  [days] 
    871955   rn_zdef_max   =  0.9             !  maximum fractional e3t deformation 
    872    ln_vvl_dbg    = .true.           !  debug prints    (T/F) 
     956   ln_vvl_dbg    = .false.          !  debug prints    (T/F) 
    873957/ 
    874958!----------------------------------------------------------------------- 
     
    890974   ln_dynvor_eeT = .false. !  energy conserving scheme (een using e3t) 
    891975   ln_dynvor_een = .false. !  energy & enstrophy scheme 
    892       nn_een_e3f = 0          ! =0  e3f = mi(mj(e3t))/4  
     976      nn_een_e3f = 0          ! =0  e3f = mi(mj(e3t))/4 
    893977      !                       ! =1  e3f = mi(mj(e3t))/mi(mj( tmask)) 
    894978   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T)        ==>>> PLEASE DO NOT ACTIVATE 
     
    913997      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    914998         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    915          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    916          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     999         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     1000         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    9171001      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    9181002         rn_bt_cmax   =  0.8        ! =T : the Maximum Courant Number allowed 
    919          nn_baro      = 30          ! =F : the number of sub-step in rn_rdt seconds 
     1003         nn_e         = 30          ! =F : the number of sub-step in rn_Dt seconds 
    9201004      rn_bt_alpha   = 0.         ! Temporal diffusion parameter (if ln_bt_av=F) 
    9211005/ 
     
    9351019      !                             !  =-30  read in eddy_viscosity_3D.nc file 
    9361020      !                             !  =-20  read in eddy_viscosity_2D.nc file 
    937       !                             !  =  0  constant  
     1021      !                             !  =  0  constant 
    9381022      !                             !  = 10  F(k)=c1d 
    9391023      !                             !  = 20  F(i,j)=F(grid spacing)=c2d 
     
    9411025      !                             !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
    9421026      !                             !  = 32  F(i,j,k)=F(local gridscale and deformation rate) 
    943       !                        !  time invariant coefficients :  ahm = 1/2  Uv*Lv   (lap case)  
     1027      !                        !  time invariant coefficients :  ahm = 1/2  Uv*Lv   (lap case) 
    9441028      !                             !                            or  = 1/12 Uv*Lv^3 (blp case) 
    9451029      rn_Uv      = 0.1              !  lateral viscous velocity [m/s] (nn_ahm_ijk_t= 0, 10, 20, 30) 
     
    10651149                              !        = 0  constant 10 m length scale 
    10661150                              !        = 1  0.5m at the equator to 30m poleward of 40 degrees 
    1067       rn_eice     =   4       !  below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4    
     1151      rn_eice     =   4       !  below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4 
    10681152/ 
    10691153!----------------------------------------------------------------------- 
     
    11241208   ln_mevar    = .true.    !  variable (T) or constant (F) mixing efficiency 
    11251209   ln_tsdiff   = .true.    !  account for differential T/S mixing (T) or not (F) 
    1126 / 
    1127  
     1210 
     1211   cn_dir      = './'      !  root directory for the iwm data location                                                                            
     1212   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     1213   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     1214   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     1215   sn_mpb      = 'NOT USED'              , -12               , 'mixing_power_bot' , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1216   sn_mpp      = 'NOT USED'              , -12               , 'mixing_power_pyc' , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1217   sn_mpc      = 'NOT USED'              , -12               , 'mixing_power_cri' , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1218   sn_dsb      = 'NOT USED'              , -12               , 'decay_scale_bot'  , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1219   sn_dsc      = 'NOT USED'              , -12               , 'decay_scale_cri'  , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1220/ 
    11281221!!====================================================================== 
    11291222!!                  ***  Diagnostics namelists  ***                   !! 
     
    11351228!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    11361229!!   namflo       float parameters                                      (default: OFF) 
    1137 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    11381230!!   nam_diadct   transports through some sections                      (default: OFF) 
    1139 !!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
    11401231!!   nam_dia25h   25h Mean Output                                       (default: OFF) 
    11411232!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4") 
     
    11621253!!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S 
    11631254!!gm 
    1164 !----------------------------------------------------------------------- 
    1165 &namptr        !   Poleward Transport Diagnostic                        (default: OFF) 
    1166 !----------------------------------------------------------------------- 
    1167    ln_diaptr   = .false.   !  Poleward heat and salt transport (T) or not (F) 
    1168    ln_subbas   = .false.   !  Atlantic/Pacific/Indian basins computation (T) or not 
    11691255/ 
    11701256!----------------------------------------------------------------------- 
     
    11951281/ 
    11961282!----------------------------------------------------------------------- 
    1197 &nam_diaharm   !   Harmonic analysis of tidal constituents              (default: OFF) 
    1198 !----------------------------------------------------------------------- 
    1199     ln_diaharm = .false.   ! Choose tidal harmonic output or not 
    1200        nit000_han = 1      !    First time step used for harmonic analysis 
    1201        nitend_han = 75     !    Last time step used for harmonic analysis 
    1202        nstep_han  = 15     !    Time step frequency for harmonic analysis 
    1203        tname(1)   = 'M2'   !    Name of tidal constituents 
    1204        tname(2)   = 'K1'   !              --- 
    1205 / 
    1206 !----------------------------------------------------------------------- 
    12071283&nam_diadct    !   transports through some sections                     (default: OFF) 
    12081284!----------------------------------------------------------------------- 
     
    12131289       !                   !     -1 : debug all section 
    12141290       !                   !  0 < n : debug section number n 
    1215 / 
    1216 !----------------------------------------------------------------------- 
    1217 &nam_diatmb    !  Top Middle Bottom Output                              (default: OFF) 
    1218 !----------------------------------------------------------------------- 
    1219    ln_diatmb   = .false.   !  Choose Top Middle and Bottom output or not 
    12201291/ 
    12211292!----------------------------------------------------------------------- 
     
    13391410&namctl        !   Control prints                                       (default: OFF) 
    13401411!----------------------------------------------------------------------- 
    1341    ln_ctl = .FALSE.                 ! Toggle all report printing on/off (T/F); Ignored if sn_cfctl%l_config is T 
     1412   sn_cfctl%l_glochk = .FALSE.    ! Range sanity checks are local (F) or global (T). Set T for debugging only 
     1413   sn_cfctl%l_allon  = .FALSE.    ! IF T activate all options. If F deactivate all unless l_config is T 
    13421414     sn_cfctl%l_config = .TRUE.     ! IF .true. then control which reports are written with the following 
    1343        sn_cfctl%l_runstat = .FALSE. ! switches and which areas produce reports with the proc integer settings. 
     1415       sn_cfctl%l_runstat = .TRUE. ! switches and which areas produce reports with the proc integer settings. 
    13441416       sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure 
    13451417       sn_cfctl%l_oceout  = .FALSE. ! that  all areas report. 
    13461418       sn_cfctl%l_layout  = .FALSE. ! 
    1347        sn_cfctl%l_mppout  = .FALSE. ! 
    1348        sn_cfctl%l_mpptop  = .FALSE. ! 
     1419       sn_cfctl%l_prtctl  = .FALSE. ! 
     1420       sn_cfctl%l_prttrc  = .FALSE. ! 
     1421       sn_cfctl%l_oasout  = .FALSE. ! 
    13491422       sn_cfctl%procmin   = 0       ! Minimum area number for reporting [default:0] 
    13501423       sn_cfctl%procmax   = 1000000 ! Maximum area number for reporting [default:1000000] 
  • NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/cfgs/SHARED/namelist_top_ref

    r12178 r12928  
    1414&namtrc_run      !   run information 
    1515!----------------------------------------------------------------------- 
    16    nn_dttrc      =  1        !  time step frequency for passive sn_tracers 
    1716   ln_top_euler  = .false.   !  use Euler time-stepping for TOP 
    1817   ln_rsttr      = .false.   !  start from a restart file (T) or not (F) 
     
    4140   ln_trcdmp     =  .false.  !  add a damping termn (T) or not (F) 
    4241   ln_trcdmp_clo =  .false.  !  damping term (T) or not (F) on closed seas 
     42   ln_trcbc      =  .false.  !  Surface, Lateral or Open Boundaries conditions 
    4343   ! 
    4444   jp_dia3d      = 0         ! Number of 3D diagnostic variables 
     
    135135   cn_dir_obc    =  './'     !  root directory for the location of OPEN data files 
    136136   ln_rnf_ctl    = .false.   !  Remove runoff dilution on tracers with absent river load 
    137    rn_bc_time    =  86400.   !  Time scaling factor for SBC and CBC data (seconds in a day) 
     137   rn_sbc_time   =  86400.   !  Time scaling factor for SBC data (seconds in a day) 
     138   rn_cbc_time   =  86400.   !  Time scaling factor for CBC data (seconds in a day) 
    138139/ 
    139140!---------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.