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 13463 for NEMO/branches/2019/dev_r11351_fldread_with_XIOS/cfgs/SHARED – NEMO

Ignore:
Timestamp:
2020-09-14T17:40:34+02:00 (4 years ago)
Author:
andmirek
Message:

Ticket #2195:update to trunk 13461

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

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2019/dev_r11351_fldread_with_XIOS

    • Property svn:externals

      •  

        old new  
        33^/utils/build/mk@HEAD         mk 
        44^/utils/tools@HEAD            tools 
        5 ^/vendors/AGRIF/dev@HEAD      ext/AGRIF 
         5^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS      ext/AGRIF 
        66^/vendors/FCM@HEAD            ext/FCM 
        77^/vendors/IOIPSL@HEAD         ext/IOIPSL 
         8 
         9# SETTE 
         10^/utils/CI/sette@13382        sette 

  • NEMO/branches/2019/dev_r11351_fldread_with_XIOS/cfgs/SHARED/README.rst

    r10598 r13463  
    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_r11351_fldread_with_XIOS/cfgs/SHARED/domain_def_nemo.xml

    r9930 r13463  
    1010     </domain> 
    1111 
    12      <domain id="1point" domain_ref="grid_T" > 
    13        <zoom_domain ibegin="1" jbegin="1" ni="1" nj="1"/> 
    14      </domain> 
     12      <domain id="1point" domain_ref="grid_T" > 
     13         <zoom_domain ibegin="139" jbegin="119" ni="1" nj="1"/> 
     14      </domain> 
     15 
     16 
    1517     <!--   Eq section --> 
    1618     <domain id="EqT" domain_ref="grid_T" > <zoom_domain id="EqT"/> </domain> 
    … …  
    179181     <domain id="EqW" domain_ref="grid_W" > <zoom_domain id="EqW"/> </domain> 
    180182 
    181       
    182      <!--   zonal mean grid   --> 
    183      <domain id="gznl" long_name="gznl"/> 
    184      <domain id="znl_T" domain_ref="gznl" > <zoom_domain id="znl_T"/> </domain> 
    185      <domain id="znl_W" domain_ref="gznl" > <zoom_domain id="znl_W"/> </domain> 
     183              <!--   zonal mean grid   --> 
     184     <domain_group id="gznl"> 
     185        <domain id="gznl" long_name="gznl"/> 
     186        <domain id="ptr" domain_ref="gznl" >  
     187            <zoom_domain id="ptr" ibegin="0000" jbegin="0" ni="1" nj="0000" />  
     188        </domain> 
     189      </domain_group> 
    186190 
    187191      

  • NEMO/branches/2019/dev_r11351_fldread_with_XIOS/cfgs/SHARED/field_def_nemo-ice.xml

    r10911 r13463  
    3232  
    3333     <!-- general fields --> 
    34           <field id="icemass"      long_name="Sea-ice mass per area"                                   standard_name="sea_ice_amount"                            unit="kg/m2"/> 
     34          <field id="icemass"      long_name="Sea-ice mass per gridcell area"                          standard_name="sea_ice_amount"                            unit="kg/m2"/> 
    3535          <field id="snwmass"      long_name="Snow mass per area"                                      standard_name="liquid_water_content_of_surface_snow"      unit="kg/m2"/> 
    36           <field id="iceconc"      long_name="Sea-ice area fraction"                                   standard_name="sea_ice_area_fraction"                     unit=""    /> 
    37           <field id="icevolu"      long_name="Sea-ice volume per area"                                 standard_name="sea_ice_thickness"                         unit="m"   /> 
    38           <field id="icethic"      long_name="Sea-ice thickness per area"                              standard_name="sea_ice_thickness"                         unit="m"   /> 
    39           <field id="snwthic"      long_name="Snow thickness per area"                                 standard_name="snow_thickness"                            unit="m"   /> 
    40           <field id="icebrv"       long_name="brine volume"                                                                                                      unit="%"   /> 
    41           <field id="iceage"       long_name="ice age"                                                                                                           unit="days"/> 
    42           <field id="icehnew"      long_name="frazil ice collection thickness"                                                                                   unit="m"   /> 
    43           <field id="snwvolu"      long_name="snow volume"                                                                                                       unit="m"   /> 
    44           <field id="icefrb"       long_name="Sea-ice freeboard"                                       standard_name="sea_ice_freeboard"                         unit="m"   /> 
    45           <field id="icealb"       long_name="Sea-ice or snow albedo"                                  standard_name="sea_ice_albedo"                            unit=""    /> 
    46           <field id="tau_icebfr"   long_name="ice friction on ocean bottom for landfast ice"                                                                     unit="N/2" /> 
     36          <field id="iceconc"      long_name="Sea-ice area fraction"                                   standard_name="sea_ice_area_fraction"                     unit=""     /> 
     37          <field id="icevolu"      long_name="Sea-ice volume per area"                                 standard_name="sea_ice_thickness"                         unit="m"    /> 
     38          <field id="icethic"      long_name="Sea-ice thickness"                                       standard_name="sea_ice_thickness"                         unit="m"    /> 
     39          <field id="snwthic"      long_name="Snow thickness"                                          standard_name="snow_thickness"                            unit="m"    /> 
     40          <field id="icebrv"       long_name="brine volume"                                                                                                      unit="%"    /> 
     41          <field id="iceage"       long_name="ice age"                                                                                                           unit="days" detect_missing_value="true"/> 
     42          <field id="icehnew"      long_name="frazil ice collection thickness"                                                                                   unit="m"    /> 
     43          <field id="snwvolu"      long_name="snow volume"                                                                                                       unit="m"    /> 
     44          <field id="icefrb"       long_name="Sea-ice freeboard"                                       standard_name="sea_ice_freeboard"                         unit="m"    /> 
     45          <field id="icealb"       long_name="Sea-ice or snow albedo"                                  standard_name="sea_ice_albedo"                            unit=""    detect_missing_value="true" /> 
    4746      
    4847     <!-- melt ponds --> 
    49      <field id="iceapnd"      long_name="melt pond fraction"                                      standard_name="sea_ice_meltpond_fraction"                 unit="%" />  
     48     <field id="iceapnd"      long_name="melt pond concentration"                                 standard_name="sea_ice_meltpond_concentration"            unit=""  />  
     49          <field id="icehpnd"      long_name="melt pond depth"                                         standard_name="sea_ice_meltpond_depth"                    unit="m" />  
    5050          <field id="icevpnd"      long_name="melt pond volume"                                        standard_name="sea_ice_meltpond_volume"                   unit="m" />  
    5151      
    5252     <!-- heat --> 
    53           <field id="icetemp"      long_name="Mean ice temperature"                                                                                              unit="degC" /> 
    54           <field id="snwtemp"      long_name="Mean snow temperature"                                                                                             unit="degC" /> 
    55           <field id="icettop"      long_name="temperature at the ice surface"                                                                                    unit="degC" /> 
    56           <field id="icetbot"      long_name="temperature at the ice bottom"                                                                                     unit="degC" /> 
    57           <field id="icetsni"      long_name="temperature at the snow-ice interface"                                                                             unit="degC" /> 
     53          <field id="icetemp"      long_name="Mean ice temperature"                                                                                              unit="degC" detect_missing_value="true" /> 
     54          <field id="snwtemp"      long_name="Mean snow temperature"                                                                                             unit="degC" detect_missing_value="true" /> 
     55          <field id="icettop"      long_name="temperature at the ice surface"                                                                                    unit="degC" detect_missing_value="true" /> 
     56          <field id="icetbot"      long_name="temperature at the ice bottom"                                                                                     unit="degC" detect_missing_value="true" /> 
     57          <field id="icetsni"      long_name="temperature at the snow-ice interface"                                                                             unit="degC" detect_missing_value="true" /> 
    5858          <field id="icehc"        long_name="ice heat content"                                                                                                  unit="J/m2" />  
    5959          <field id="snwhc"        long_name="snow heat content"                                                                                                 unit="J/m2" /> 
    6060      
    6161     <!-- salt --> 
    62           <field id="icesalt"      long_name="Sea ice salinity"                                                                                                  unit="g/kg"  /> 
     62          <field id="icesalt"      long_name="Sea ice salinity"                                                                                                  unit="g/kg"  detect_missing_value="true" /> 
    6363          <field id="icesalm"      long_name="Mass of salt in sea ice per area"                        standard_name="sea_ice_salt_mass"                         unit="kg/m2" /> 
    6464      
    … …  
    7171          <field id="utau_oi"      long_name="X-component of ocean stress on sea ice"                  standard_name="sea_ice_base_upward_x_stress"              unit="N/m2" /> 
    7272          <field id="vtau_oi"      long_name="Y-component of ocean stress on sea ice"                  standard_name="sea_ice_base_upward_y_stress"              unit="N/m2" /> 
     73          <field id="utau_bi"      long_name="X-component of ocean bottom stress on sea ice -landfast" standard_name="ocean_bottom_upward_x_stress"              unit="N/m2" /> 
     74          <field id="vtau_bi"      long_name="Y-component of ocean bottom stress on sea ice -landfast" standard_name="ocean_bottom_upward_y_stress"              unit="N/m2" /> 
    7375          <field id="isig1"        long_name="1st principal stress component for EVP rhg"                                                                        unit=""     /> 
    7476          <field id="isig2"        long_name="2nd principal stress component for EVP rhg"                                                                        unit=""     /> 
    … …  
    9092          <field id="qemp_ice"     long_name="Downward Heat Flux from E-P over ice"                                                                              unit="W/m2" /> 
    9193          <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=""     /> 
    9296      
    9397     <!-- trends --> 
    94           <field id="afxthd"       long_name="sea-ice area fraction change from thermodynamics"  standard_name="tendency_of_sea_ice_area_fraction_due_to_dynamics" unit="s-1" /> 
    95           <field id="afxdyn"       long_name="sea-ice area fraction change from dynamics"        standard_name="tendency_of_sea_ice_area_fraction_due_to_dynamics" unit="s-1" /> 
    96      <field id="afxtot"       long_name="area tendency (total)"                                                                                               unit="s-1" /> 
     98          <field id="afxthd"       long_name="sea-ice area fraction change from thermodynamics" standard_name="tendency_of_sea_ice_area_fraction_due_to_dynamics" unit="s-1" /> 
     99          <field id="afxdyn"       long_name="sea-ice area fraction change from dynamics"       standard_name="tendency_of_sea_ice_area_fraction_due_to_dynamics" unit="s-1" /> 
     100     <field id="afxtot"       long_name="area tendency (total)"                                                                                              unit="s-1" /> 
    97101      
    98102     <!-- momentum (advection) --> 
    … …  
    168172          <field id="e3t_m"    unit="m"    /> 
    169173          <field id="frq_m"    unit="-"    /> 
    170  
    171      <!-- categories --> 
    172      <field id="iceconc_cat"  long_name="Sea-ice concentration per category"                unit=""       grid_ref="grid_T_3D_ncatice" /> 
    173           <field id="icethic_cat"  long_name="Sea-ice thickness per category"                    unit="m"      grid_ref="grid_T_3D_ncatice" /> 
    174           <field id="snwthic_cat"  long_name="Snow thickness per category"                       unit="m"      grid_ref="grid_T_3D_ncatice" /> 
    175           <field id="icesalt_cat"  long_name="Sea-Ice Bulk salinity per category"                unit="g/kg"   grid_ref="grid_T_3D_ncatice" /> 
    176           <field id="icetemp_cat"  long_name="Ice temperature per category"                      unit="degC"   grid_ref="grid_T_3D_ncatice" /> 
    177           <field id="snwtemp_cat"  long_name="Snow temperature per category"                     unit="degC"   grid_ref="grid_T_3D_ncatice" /> 
    178           <field id="icettop_cat"  long_name="Ice/snow surface temperature per category"         unit="degC"   grid_ref="grid_T_3D_ncatice" /> 
    179           <field id="iceapnd_cat"  long_name="Ice melt pond concentration per category"          unit="%"      grid_ref="grid_T_3D_ncatice" />  
    180           <field id="icehpnd_cat"  long_name="Ice melt pond thickness per category"              unit="m"      grid_ref="grid_T_3D_ncatice" />  
    181           <field id="iceafpnd_cat" long_name="Ice melt pond fraction per category"               unit="m"      grid_ref="grid_T_3D_ncatice" />  
    182           <field id="icemask_cat"  long_name="Fraction of time step with sea ice (per category)" unit=""       grid_ref="grid_T_3D_ncatice" /> 
    183           <field id="iceage_cat"   long_name="Ice age per category"                              unit="days"   grid_ref="grid_T_3D_ncatice" /> 
    184           <field id="icebrv_cat"   long_name="Brine volume per category"                         unit="%"      grid_ref="grid_T_3D_ncatice" /> 
    185174 
    186175     <!-- ================= --> 
    … …  
    241230 
    242231     <!-- momentum (advection) --> 
    243      <field id="xmtrpice"     long_name="X-component of ice mass transport"                       standard_name="ice_x_transport"                       unit="kg/s" /> 
    244           <field id="ymtrpice"     long_name="Y-component of ice mass transport"                       standard_name="ice_y_transport"                       unit="kg/s" /> 
    245           <field id="xmtrpsnw"     long_name="X-component of snw mass transport"                       standard_name="snw_x_transport"                       unit="kg/s" /> 
    246           <field id="ymtrpsnw"     long_name="Y-component of snw mass transport"                       standard_name="snw_y_transport"                       unit="kg/s" /> 
    247           <field id="xatrp"        long_name="X-component of ice area transport"                       standard_name="area_x_transport"                      unit="m2/s" /> 
    248           <field id="yatrp"        long_name="Y-component of ice area transport"                       standard_name="area_y_transport"                      unit="m2/s" /> 
    249           <field id="xmtrptot"     long_name="X-component of sea-ice mass transport"                   standard_name="sea_ice_x_transport"                   unit="kg/s" > xmtrpice + xmtrpsnw </field> 
    250           <field id="ymtrptot"     long_name="Y-component of sea-ice mass transport"                   standard_name="sea_ice_y_transport"                   unit="kg/s" > ymtrpice + ymtrpsnw </field> 
    251  
    252      <!-- categories --> 
    253           <field id="iceconc_cat_cmip"     long_name="Sea-ice area fractions in thickness categories"  standard_name="sea_ice_area_fraction_over_categories" detect_missing_value="true" unit=""  grid_ref="grid_T_3D_ncatice" > iceconc_cat      * icemask_cat + $missval * (1.-icemask_cat) </field> 
    254           <field id="icethic_cat_cmip"     long_name="Sea-ice thickness in thickness categories"       standard_name="sea_ice_thickness_over_categories"     detect_missing_value="true" unit="m" grid_ref="grid_T_3D_ncatice" > icethic_cat      * icemask_cat + $missval * (1.-icemask_cat) </field> 
    255           <field id="snwthic_cat_cmip"     long_name="Snow thickness in thickness categories"          standard_name="snow_thickness_over_categories"        detect_missing_value="true" unit="m" grid_ref="grid_T_3D_ncatice" > snwthic_cat      * icemask_cat + $missval * (1.-icemask_cat) </field> 
    256           <field id="iceconc_cat_pct_cmip" long_name="Sea-ice area fractions in thickness categories"  standard_name="sea_ice_area_fraction_over_categories" detect_missing_value="true" unit="%" grid_ref="grid_T_3D_ncatice" > iceconc_cat*100. * icemask_cat + $missval * (1.-icemask_cat) </field> 
    257      
     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 
    258276   </field_group> <!-- SBC_2D --> 
    259277 
     278   <!-- categories --> 
     279        <field_group id="SBC_3D" grid_ref="grid_T_ncatice" > 
     280 
     281          <!-- standard ice fields -->   
     282     <field id="iceconc_cat"  long_name="Sea-ice concentration per category"                unit=""        /> 
     283          <field id="icethic_cat"  long_name="Sea-ice thickness per category"                    unit="m"       detect_missing_value="true" /> 
     284          <field id="snwthic_cat"  long_name="Snow thickness per category"                       unit="m"       detect_missing_value="true" /> 
     285          <field id="icesalt_cat"  long_name="Sea-Ice Bulk salinity per category"                unit="g/kg"    detect_missing_value="true" /> 
     286          <field id="icetemp_cat"  long_name="Ice temperature per category"                      unit="degC"    detect_missing_value="true" /> 
     287          <field id="snwtemp_cat"  long_name="Snow temperature per category"                     unit="degC"    detect_missing_value="true" /> 
     288          <field id="icettop_cat"  long_name="Ice/snow surface temperature per category"         unit="degC"    detect_missing_value="true" /> 
     289          <field id="iceapnd_cat"  long_name="Ice melt pond concentration per category"          unit=""        />  
     290          <field id="icehpnd_cat"  long_name="Ice melt pond thickness per category"              unit="m"       detect_missing_value="true" />  
     291          <field id="iceafpnd_cat" long_name="Ice melt pond fraction per category"               unit=""        />  
     292          <field id="icemask_cat"  long_name="Fraction of time step with sea ice (per category)" unit=""        /> 
     293          <field id="iceage_cat"   long_name="Ice age per category"                              unit="days"    detect_missing_value="true" /> 
     294          <field id="icebrv_cat"   long_name="Brine volume per category"                         unit="%"       detect_missing_value="true" /> 
     295          <field id="icealb_cat"   long_name="Sea-ice or snow albedo"                            unit=""        detect_missing_value="true" /> 
     296 
     297          <!-- Add-ons for SIMIP --> 
     298          <field id="iceconc_cat_cmip"     long_name="Sea-ice area fractions in thickness categories"  standard_name="sea_ice_area_fraction_over_categories" detect_missing_value="true" unit=""   > iceconc_cat      * icemask_cat + $missval * (1.-icemask_cat) </field> 
     299          <field id="icethic_cat_cmip"     long_name="Sea-ice thickness in thickness categories"       standard_name="sea_ice_thickness_over_categories"     detect_missing_value="true" unit="m"  > icethic_cat      * icemask_cat + $missval * (1.-icemask_cat) </field> 
     300          <field id="snwthic_cat_cmip"     long_name="Snow thickness in thickness categories"          standard_name="snow_thickness_over_categories"        detect_missing_value="true" unit="m"  > snwthic_cat      * icemask_cat + $missval * (1.-icemask_cat) </field> 
     301          <field id="iceconc_cat_pct_cmip" long_name="Sea-ice area fractions in thickness categories"  standard_name="sea_ice_area_fraction_over_categories" detect_missing_value="true" unit="%"  > iceconc_cat*100. * icemask_cat + $missval * (1.-icemask_cat) </field> 
     302 
     303   </field_group> <!-- SBC_3D --> 
     304    
    260305   <!-- scalar variables --> 
    261    <field_group id="SBC_0D"  grid_ref="grid_1point" > 
     306   <field_group id="SBC_scalar"  grid_ref="grid_scalar" > 
     307          <field id="NH_iceextt"      long_name="Sea ice extent North"                   standard_name="sea_ice_extent_n"                   unit="1e6_km2"  /> 
     308          <field id="SH_iceextt"      long_name="Sea ice extent South"                   standard_name="sea_ice_extent_s"                   unit="1e6_km2"  /> 
     309          <field id="NH_icevolu"      long_name="Sea ice volume North"                   standard_name="sea_ice_volume_n"                   unit="1e3_km3"  /> 
     310          <field id="SH_icevolu"      long_name="Sea ice volume South"                   standard_name="sea_ice_volume_s"                   unit="1e3_km3"  /> 
     311          <field id="NH_icearea"      long_name="Sea ice area North"                     standard_name="sea_ice_area_n"                     unit="1e6_km2"  /> 
     312          <field id="SH_icearea"      long_name="Sea ice area South"                     standard_name="sea_ice_area_s"                     unit="1e6_km2"  /> 
     313 
    262314          <!-- available with ln_icediaout --> 
    263315          <field id="ibgfrcvoltop"    long_name="global mean ice/snow forcing at interface ice/snow-atm (volume equivalent ocean volume)"   unit="km3"      /> 
    … …  
    281333          <field id="sbgheat_tot"     long_name="global mean snow heat content"                                                             unit="1e20J"    /> 
    282334 
    283           <field id="NH_iceextt"      long_name="Sea ice extent North"                   standard_name="sea_ice_extent_n"                   unit="1e6_km2"  /> 
    284           <field id="SH_iceextt"      long_name="Sea ice extent South"                   standard_name="sea_ice_extent_s"                   unit="1e6_km2"  /> 
    285           <field id="NH_icevolu"      long_name="Sea ice volume North"                   standard_name="sea_ice_volume_n"                   unit="1e3_km3"  /> 
    286           <field id="SH_icevolu"      long_name="Sea ice volume South"                   standard_name="sea_ice_volume_s"                   unit="1e3_km3"  /> 
    287           <field id="NH_icearea"      long_name="Sea ice area North"                     standard_name="sea_ice_area_n"                     unit="1e6_km2"  /> 
    288           <field id="SH_icearea"      long_name="Sea ice area South"                     standard_name="sea_ice_area_s"                     unit="1e6_km2"  /> 
    289  
    290           <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" /> 
    291           <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" />   
    292           <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" /> 
    293    </field_group> <!-- SBC_0D --> 
     335   </field_group> 
    294336   
    295337   <!--  
    … …  
    299341       
    300342   <field_group id="myvarICE"        grid_ref="grid_T_2D" > 
    301           <field field_ref="icethic"          name="sithic"  /> 
    302           <field field_ref="icevolu"          name="sivolu"  /> 
    303           <field field_ref="iceconc"          name="siconc"  /> 
     343     <!-- ice mask --> 
     344     <field field_ref="icemask"          name="simsk"   /> 
     345     <field field_ref="icemask05"        name="simsk05" /> 
     346     <field field_ref="icemask15"        name="simsk15" /> 
     347      
     348     <!-- general --> 
     349     <field field_ref="snwvolu"          name="snvolu" /> 
     350     <field field_ref="snwthic"          name="snthic" /> 
     351     <field field_ref="icethic"          name="sithic" /> 
     352<!-- 
     353     <field field_ref="icethic"          name="sithic_max" operation="maximum" /> 
     354     <field field_ref="icethic"          name="sithic_min" operation="minimum" /> 
     355--> 
     356     <field field_ref="fasticepres"      name="fasticepres" /> 
     357     <field field_ref="icevolu"          name="sivolu" /> 
     358     <field field_ref="iceconc"          name="siconc" /> 
     359     <field field_ref="icesalt"          name="sisali" /> 
     360     <field field_ref="iceapnd"          name="siapnd" /> 
     361     <field field_ref="icehpnd"          name="sihpnd" /> 
     362     <field field_ref="icevpnd"          name="sivpnd" /> 
     363          <field field_ref="iceage"           name="siage"  /> 
     364     <field field_ref="sst_m"            name="sst_m"  /> 
     365     <field field_ref="sss_m"            name="sss_m"  /> 
     366      
     367     <!-- heat --> 
     368     <field field_ref="icetemp"          name="sitemp" /> 
     369     <field field_ref="snwtemp"          name="sntemp" /> 
     370     <field field_ref="icettop"          name="sittop" /> 
     371     <field field_ref="icetbot"          name="sitbot" /> 
     372     <field field_ref="icetsni"          name="sitsni" /> 
     373      
     374     <!-- momentum --> 
     375     <field field_ref="uice"             name="sivelu"  /> 
     376     <field field_ref="vice"             name="sivelv"  /> 
     377     <field field_ref="icevel"           name="sivelo"  /> 
     378     <field field_ref="utau_ai"          name="utau_ai" /> 
     379     <field field_ref="vtau_ai"          name="vtau_ai" /> 
     380     <field field_ref="utau_oi"          name="utau_oi" /> 
     381     <field field_ref="vtau_oi"          name="vtau_oi" /> 
     382      
     383     <!-- rheology --> 
     384     <field field_ref="icediv"           name="sidive"  /> 
     385     <field field_ref="iceshe"           name="sishea"  /> 
     386     <field field_ref="icestr"           name="sistre"  /> 
     387     <field field_ref="normstr"          name="normstr" /> 
     388     <field field_ref="sheastr"          name="sheastr" /> 
     389     <field field_ref="isig1"            name="isig1"   /> 
     390     <field field_ref="isig2"            name="isig2"   /> 
     391     <field field_ref="isig3"            name="isig3"   /> 
     392      
     393     <!-- heat fluxes --> 
     394     <field field_ref="qt_oce_ai"        name="qt_oce_ai"  /> 
     395     <field field_ref="qt_atm_oi"        name="qt_atm_oi"  /> 
     396     <field field_ref="qtr_ice_top"      name="qtr_ice_top"/> 
     397     <field field_ref="qtr_ice_bot"      name="qtr_ice_bot"/> 
     398     <field field_ref="qt_ice"           name="qt_ice"     /> 
     399     <field field_ref="qsr_ice"          name="qsr_ice"    /> 
     400     <field field_ref="qns_ice"          name="qns_ice"    /> 
     401     <field field_ref="qemp_ice"         name="qemp_ice"   /> 
     402     <field field_ref="albedo"           name="albedo"     /> 
     403          <field field_ref="icealb"           name="icealb"     /> 
     404      
     405     <field field_ref="hfxcndtop"        name="hfxcndtop"  /> 
     406     <field field_ref="hfxcndbot"        name="hfxcndbot"  /> 
     407     <field field_ref="hfxsensib"        name="hfxsensib"  /> 
     408 
     409     <field field_ref="icehc"         name="icehc"  /> 
     410     <field field_ref="snwhc"         name="snwhc"  /> 
     411     <field field_ref="hfxbog"        name="hfxbog"  /> 
     412     <field field_ref="hfxbom"        name="hfxbom"  /> 
     413     <field field_ref="hfxsum"        name="hfxsum"  /> 
     414     <field field_ref="hfxopw"        name="hfxopw"  /> 
     415     <field field_ref="hfxdif"        name="hfxdif"  /> 
     416     <field field_ref="hfxsnw"        name="hfxsnw"  /> 
     417     <field field_ref="hfxerr"        name="hfxerr"  /> 
     418     <field field_ref="hfxthd"        name="hfxthd"  /> 
     419     <field field_ref="hfxdyn"        name="hfxdyn"  /> 
     420     <field field_ref="hfxres"        name="hfxres"  /> 
     421     <field field_ref="hfxsub"        name="hfxsub"  /> 
     422     <field field_ref="hfxspr"        name="hfxspr"  /> 
     423 
     424<!-- 
     425     <field field_ref="hfxdhc"        name="hfxbudget"  operation="average" freq_op="1d" > @qt_atm_oi - @qt_oce_ai -@hfxdhc </field> 
     426--> 
     427      
     428     <!-- salt fluxes --> 
     429     <field field_ref="sfxice"           name="sfxice" /> 
     430      
     431     <!-- mass fluxes --> 
     432     <field field_ref="vfxice"           name="vfxice" /> 
     433     <field field_ref="vfxsnw"           name="vfxsnw" /> 
     434 
     435     <field field_ref="vfxbom"           name="vfxbom" /> 
     436     <field field_ref="vfxbog"           name="vfxbog" /> 
     437     <field field_ref="vfxsum"           name="vfxsum" /> 
     438     <field field_ref="vfxopw"           name="vfxopw" /> 
     439     <field field_ref="vfxlam"           name="vfxlam" /> 
     440     <field field_ref="vfxsni"           name="vfxsni" /> 
     441     <field field_ref="vfxdyn"           name="vfxdyn" /> 
     442     <field field_ref="vfxres"           name="vfxres" /> 
     443     <field field_ref="vfxpnd"           name="vfxpnd" /> 
     444     <field field_ref="vfxsub"           name="vfxsub" /> 
     445     <field field_ref="vfxsub_err"       name="vfxsub_err" /> 
     446          <field field_ref="vfxthin"          name="vfxthin" /> 
     447    
     448     <!-- diag error for negative ice volume after advection --> 
     449     <field field_ref="iceneg_pres"      name="sineg_pres" /> 
     450     <field field_ref="iceneg_volu"      name="sineg_volu" /> 
     451     <field field_ref="iceneg_hfx"       name="sineg_hfx"  /> 
     452        </field_group> 
     453 
     454 
     455        <field_group id="myvarICE_cat"        grid_ref="grid_T_ncatice" > 
     456 
     457     <!-- categories --> 
     458     <field field_ref="icemask_cat"      name="simskcat"/> 
     459     <field field_ref="snwthic_cat"      name="snthicat"/> 
     460     <field field_ref="iceconc_cat"      name="siconcat"/> 
     461     <field field_ref="icethic_cat"      name="sithicat"/> 
     462     <field field_ref="icesalt_cat"      name="sisalcat"/> 
     463     <field field_ref="icetemp_cat"      name="sitemcat"/> 
     464     <field field_ref="snwtemp_cat"      name="sntemcat"/> 
     465     <field field_ref="icettop_cat"      name="sitopcat"/> 
     466 
    304467   </field_group>     
    305468 
    306    <field_group id="ICE_globalbudget"  grid_ref="grid_1point" > 
     469   <field_group id="ICE_globalbudget"  grid_ref="grid_scalar" > 
     470     <!-- global contents --> 
    307471        <field field_ref="ibgvol_tot"       name="ibgvol_tot"   /> 
    308472     <field field_ref="sbgvol_tot"       name="sbgvol_tot"   /> 
    … …  
    312476        <field field_ref="sbgheat_tot"      name="sbgheat_tot"  /> 
    313477 
    314         <field field_ref="ibgvolume"        name="ibgvolume"    /> 
     478     <!-- global drifts (conservation checks) --> 
     479     <field field_ref="ibgvolume"        name="ibgvolume"    /> 
    315480        <field field_ref="ibgsaltco"        name="ibgsaltco"    /> 
    316481        <field field_ref="ibgheatco"        name="ibgheatco"    /> 
    317482          <field field_ref="ibgheatfx"        name="ibgheatfx"    /> 
    318483 
     484     <!-- global forcings  --> 
    319485        <field field_ref="ibgfrcvoltop"     name="ibgfrcvoltop" /> 
    320486        <field field_ref="ibgfrcvolbot"     name="ibgfrcvolbot" /> 
    … …  
    333499 
    334500   <!-- SIMIP monthly scalar variables --> 
    335    <field_group id="SImon_scalars"   grid_ref="grid_1point" > 
     501   <field_group id="SImon_scalars"   grid_ref="grid_scalar" > 
    336502          <!-- Integrated quantities --> 
    337503          <field field_ref="NH_iceextt"       name="siextentn"    /> 

  • NEMO/branches/2019/dev_r11351_fldread_with_XIOS/cfgs/SHARED/field_def_nemo-oce.xml

    r10824 r13463  
    1 <?xml version="1.0"?>  
     1<?xml version="1.0"?> 
    22    <!-- $id$ --> 
    3      
    4     <!--  
     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> 
     75 
     76    <!-- 
    577============================================================================================================ 
    678=                                  definition of all existing variables                                    = 
    … …  
    880============================================================================================================ 
    981    --> 
    10     <field_definition level="1" prec="4" operation="average" enabled=".TRUE." default_value="1.e20" > <!-- time step automaticaly defined --> 
    11  
    12     <!--  
     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> 
     102 
     103    <!-- 
    13104============================================================================================================ 
    14105                                  Physical ocean model variables 
    … …  
    17108 
    18109      <!-- T grid --> 
    19        
     110 
    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"                                                 /> 
    … …  
    43146        <field id="sst_cs"       long_name="Delta SST of cool skin"                                                                                 unit="degC"     /> 
    44147   <field id="temp_3m"      long_name="temperature at 3m"                                                                                      unit="degC"     /> 
    45          
     148 
    46149        <field id="sss"          long_name="sea surface salinity"                               standard_name="sea_surface_salinity"                unit="1e-3"     /> 
    47150        <field id="sss2"         long_name="square of sea surface salinity"                                                                         unit="1e-6"      > sss * sss </field > 
    … …  
    49152        <field id="sssmin"       long_name="min of sea surface salinity"      field_ref="sss"   operation="minimum"                                                 /> 
    50153        <field id="sbs"          long_name="sea bottom salinity"                                                                                    unit="0.001"    /> 
    51         <field id="somint"       long_name="vertical integral of salinity times density"        standard_name="integral_wrt_depth_of_product_of_density_and_salinity"  unit="(kg m2) x (1e-3)" />  
    52  
    53         <field id="taubot"       long_name="bottom stress module"                                                                                   unit="N/m2"     />  
     154        <field id="somint"       long_name="vertical integral of salinity times density"        standard_name="integral_wrt_depth_of_product_of_density_and_salinity"  unit="(kg m2) x (1e-3)" /> 
     155 
     156        <field id="taubot"       long_name="bottom stress module"                                                                                   unit="N/m2"     /> 
     157 
     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 
    54163 
    55164        <field id="ssh"          long_name="sea surface height"                                 standard_name="sea_surface_height_above_geoid"             unit="m" /> 
    … …  
    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 
    … …  
    235295        <field id="us_y"        long_name="j component of Stokes drift"                      unit="m/s"     /> 
    236296      </field_group> 
    237        
     297 
    238298      <!-- SBC --> 
    239299      <field_group id="SBC" > <!-- time step automaticaly defined based on nn_fsbc --> 
    … …  
    250310          <field id="runoffs"      long_name="River Runoffs"                        standard_name="water_flux_into_sea_water_from_rivers"                                unit="kg/m2/s"   /> 
    251311          <field id="precip"       long_name="Total precipitation"                  standard_name="precipitation_flux"                                                   unit="kg/m2/s"   /> 
    252       
     312          <field id="wclosea"      long_name="closed sea empmr correction"          standard_name="closea_empmr"                                                         unit="kg/m2/s"   /> 
     313 
    253314          <field id="qt"           long_name="Net Downward Heat Flux"                standard_name="surface_downward_heat_flux_in_sea_water"                              unit="W/m2"                           /> 
    254315          <field id="qns"          long_name="non solar Downward Heat Flux"                                                                                               unit="W/m2"                           /> 
    … …  
    256317          <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" /> 
    257318          <field id="qrp"          long_name="Surface Heat Flux: Damping"            standard_name="heat_flux_into_sea_water_due_to_newtonian_relaxation"                 unit="W/m2"                           /> 
     319          <field id="qclosea"      long_name="closed sea heat content flux"          standard_name="closea_heat_content_downward_flux"                                    unit="W/m2"     /> 
    258320          <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"                        /> 
    259321          <field id="taum"         long_name="wind stress module"                    standard_name="magnitude_of_surface_downward_stress"                                 unit="N/m2"                           /> 
    260322          <field id="wspd"         long_name="wind speed module"                     standard_name="wind_speed"                                                           unit="m/s"                            /> 
    261            
     323 
    262324          <!-- * variable relative to atmospheric pressure forcing : available with ln_apr_dyn --> 
    263325          <field id="ssh_ib"       long_name="Inverse barometer sea surface height"  standard_name="sea_surface_height_correction_due_to_air_pressure_at_low_frequency"   unit="m"        /> 
    264326 
    265327          <!-- * variable related to ice shelf forcing * --> 
    266           <field id="fwfisf"       long_name="Ice shelf melting"                             unit="kg/m2/s"  /> 
    267           <field id="fwfisf3d"     long_name="Ice shelf melting"                             unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
    268           <field id="qlatisf"      long_name="Ice shelf latent heat flux"                    unit="W/m2"     /> 
    269           <field id="qlatisf3d"    long_name="Ice shelf latent heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
    270           <field id="qhcisf"       long_name="Ice shelf heat content flux"                   unit="W/m2"     /> 
    271           <field id="qhcisf3d"     long_name="Ice shelf heat content flux"                   unit="W/m2"     grid_ref="grid_T_3D" /> 
    272           <field id="isfgammat"    long_name="transfert coefficient for isf (temperature) "  unit="m/s"      /> 
    273           <field id="isfgammas"    long_name="transfert coefficient for isf (salinity)    "  unit="m/s"      /> 
    274           <field id="stbl"         long_name="salinity in the Losh tbl                    "  unit="PSU"      /> 
    275           <field id="ttbl"         long_name="temperature in the Losh tbl                 "  unit="C"        /> 
    276           <field id="utbl"         long_name="zonal current in the Losh tbl at T point    "  unit="m/s"      /> 
    277           <field id="vtbl"         long_name="merid current in the Losh tbl at T point    "  unit="m/s"      /> 
    278           <field id="thermald"     long_name="thermal driving of ice shelf melting        "  unit="C"        /> 
    279           <field id="tfrz"         long_name="top freezing point (used to compute melt)   "  unit="C"        /> 
    280           <field id="tinsitu"      long_name="top insitu temperature (used to cmpt melt)  "  unit="C"        /> 
    281           <field id="ustar"        long_name="ustar at T point used in ice shelf melting  "  unit="m/s"      /> 
     328          <field id="isftfrz_cav"     long_name="freezing point temperature at ocean/isf interface"                unit="degC"     /> 
     329          <field id="isftfrz_par"     long_name="freezing point temperature in the parametrization boundary layer" unit="degC"     /> 
     330          <field id="fwfisf_cav"      long_name="Ice shelf melt rate"                           unit="kg/m2/s"  /> 
     331          <field id="fwfisf_par"      long_name="Ice shelf melt rate"                           unit="kg/m2/s"  /> 
     332          <field id="qoceisf_cav"     long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     /> 
     333          <field id="qoceisf_par"     long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     /> 
     334          <field id="qlatisf_cav"     long_name="Ice shelf latent heat flux"                    unit="W/m2"     /> 
     335          <field id="qlatisf_par"     long_name="Ice shelf latent heat flux"                    unit="W/m2"     /> 
     336          <field id="qhcisf_cav"      long_name="Ice shelf heat content flux of injected water" unit="W/m2"     /> 
     337          <field id="qhcisf_par"      long_name="Ice shelf heat content flux of injected water" unit="W/m2"     /> 
     338          <field id="fwfisf3d_cav"    long_name="Ice shelf melt rate"                           unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
     339          <field id="fwfisf3d_par"    long_name="Ice shelf melt rate"                           unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
     340          <field id="qoceisf3d_cav"   long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     341          <field id="qoceisf3d_par"   long_name="Ice shelf ocean  heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     342          <field id="qlatisf3d_cav"   long_name="Ice shelf latent heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     343          <field id="qlatisf3d_par"   long_name="Ice shelf latent heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     344          <field id="qhcisf3d_cav"    long_name="Ice shelf heat content flux of injected water" unit="W/m2"     grid_ref="grid_T_3D" /> 
     345          <field id="qhcisf3d_par"    long_name="Ice shelf heat content flux of injected water" unit="W/m2"     grid_ref="grid_T_3D" /> 
     346          <field id="ttbl_cav"        long_name="temperature in Losch tbl"                      unit="degC"     /> 
     347          <field id="ttbl_par"        long_name="temperature in the parametrisation boundary layer" unit="degC" /> 
     348          <field id="isfthermald_cav" long_name="thermal driving of ice shelf melting"          unit="degC"     /> 
     349          <field id="isfthermald_par" long_name="thermal driving of ice shelf melting"          unit="degC"     /> 
     350          <field id="isfgammat"       long_name="Ice shelf heat-transfert velocity"             unit="m/s"      /> 
     351          <field id="isfgammas"       long_name="Ice shelf salt-transfert velocity"             unit="m/s"      /> 
     352          <field id="stbl"            long_name="salinity in the Losh tbl"                      unit="1e-3"     /> 
     353          <field id="utbl"            long_name="zonal current in the Losh tbl at T point"      unit="m/s"      /> 
     354          <field id="vtbl"            long_name="merid current in the Losh tbl at T point"      unit="m/s"      /> 
     355          <field id="isfustar"        long_name="ustar at T point used in ice shelf melting"    unit="m/s"      /> 
     356          <field id="qconisf"         long_name="Conductive heat flux through the ice shelf"    unit="W/m2"     /> 
    282357 
    283358          <!-- *_oce variables available with ln_blk_clio or ln_blk_core --> 
     359          <field id="rho_air"      long_name="Air density at 10m above sea surface"         standard_name="rho_air_10m"                                        unit="kg/m3" /> 
     360          <field id="dt_skin"      long_name="SSST-SST temperature difference"              standard_name="SSST-SST"                                             unit="K"   /> 
    284361          <field id="qlw_oce"      long_name="Longwave Downward Heat Flux over open ocean"  standard_name="surface_net_downward_longwave_flux"                 unit="W/m2"  /> 
    285362          <field id="qsb_oce"      long_name="Sensible Downward Heat Flux over open ocean"  standard_name="surface_downward_sensible_heat_flux"                unit="W/m2"  /> 
    286363          <field id="qla_oce"      long_name="Latent Downward Heat Flux over open ocean"    standard_name="surface_downward_latent_heat_flux"                  unit="W/m2"  /> 
     364          <field id="evap_oce"     long_name="Evaporation over open ocean"                  standard_name="evaporation"                                        unit="kg/m2/s" /> 
    287365          <field id="qt_oce"       long_name="total flux at ocean surface"                  standard_name="surface_downward_heat_flux_in_sea_water"            unit="W/m2"  /> 
    288366          <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"  /> 
    … …  
    290368          <field id="qemp_oce"     long_name="Downward Heat Flux from E-P over open ocean"                                                                     unit="W/m2"  /> 
    291369          <field id="taum_oce"     long_name="wind stress module over open ocean"           standard_name="magnitude_of_surface_downward_stress"               unit="N/m2"  /> 
     370 
     371          <!-- variables computed by the bulk parameterization algorithms (ln_blk) --> 
     372          <field id="Cd_oce"      long_name="Drag coefficient over open ocean"              standard_name="drag_coefficient_water"                unit=""  /> 
     373          <field id="Ce_oce"      long_name="Evaporaion coefficient over open ocean"        standard_name="evap_coefficient_water"                unit=""  /> 
     374          <field id="Ch_oce"      long_name="Sensible heat coefficient over open ocean"     standard_name="sensible_heat_coefficient_water"       unit=""  /> 
     375          <field id="theta_zt"    long_name="Potential air temperature at z=zt"             standard_name="potential_air_temperature_at_zt"       unit="degC" /> 
     376          <field id="q_zt"        long_name="Specific air humidity at z=zt"                 standard_name="specific_air_humidity_at_zt"           unit="kg/kg" /> 
     377          <field id="theta_zu"    long_name="Potential air temperature at z=zu"             standard_name="potential_air_temperature_at_zu"       unit="degC" /> 
     378          <field id="q_zu"        long_name="Specific air humidity at z=zu"                 standard_name="specific_air_humidity_at_zu"           unit="kg/kg" /> 
     379          <field id="ssq"         long_name="Saturation specific humidity of air at z=0"    standard_name="surface_air_saturation_spec_humidity"  unit="kg/kg" /> 
     380          <field id="wspd_blk"    long_name="Bulk wind speed at z=zu"                       standard_name="bulk_wind_speed_at_zu"                 unit="m/s"   /> 
     381          <!-- ln_blk + key_si3 --> 
     382          <field id="Cd_ice"      long_name="Drag coefficient over ice"                     standard_name="drag_coefficient_ice"                 unit=""  /> 
     383          <field id="Ce_ice"      long_name="Evaporaion coefficient over ice"               standard_name="evap_coefficient_ice"                 unit=""  /> 
     384          <field id="Ch_ice"      long_name="Sensible heat coefficient over ice"            standard_name="sensible_heat_coefficient_ice"        unit=""  /> 
    292385 
    293386          <!-- available key_oasis3 --> 
    … …  
    302395 
    303396          <!-- available if key_oasis3 + conservative method --> 
    304           <field id="rain"          long_name="Liquid precipitation"                                     standard_name="rainfall_flux"                                                                 unit="kg/m2/s"  /> 
     397     <field id="rain"          long_name="Liquid precipitation"                                     standard_name="rainfall_flux"                                                                 unit="kg/m2/s"  /> 
     398          <field id="rain_ao_cea"   long_name="Liquid precipitation over ice-free ocean (cell average)"  standard_name="rainfall_flux"                                                                 unit="kg/m2/s"  /> 
    305399          <field id="evap_ao_cea"   long_name="Evaporation over ice-free ocean (cell average)"           standard_name="water_evaporation_flux"                                                        unit="kg/m2/s"  /> 
    306400          <field id="isnwmlt_cea"   long_name="Snow over Ice melting (cell average)"                     standard_name="surface_snow_melt_flux"                                                        unit="kg/m2/s"  /> 
    … …  
    309403          <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"     /> 
    310404          <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"     /> 
     405          <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"     /> 
    311406          <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"     /> 
    312407          <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"     /> 
    … …  
    325420          <!-- ice field (nn_ice=1)  --> 
    326421          <field id="ice_cover"    long_name="Ice fraction"                                                 standard_name="sea_ice_area_fraction"                              unit="1"            /> 
    327            
     422 
    328423          <!-- dilution --> 
    329424          <field id="emp_x_sst"    long_name="Concentration/Dilution term on SST"                                                                                              unit="kg*degC/m2/s" /> 
    330           <field id="emp_x_sss"    long_name="Concentration/Dilution term on SSS"                                                                                              unit="kg*1e-3/m2/s" />         
     425          <field id="emp_x_sss"    long_name="Concentration/Dilution term on SSS"                                                                                              unit="kg*1e-3/m2/s" /> 
    331426          <field id="rnf_x_sst"    long_name="Runoff term on SST"                                                                                                              unit="kg*degC/m2/s" /> 
    332427          <field id="rnf_x_sss"    long_name="Runoff term on SSS"                                                                                                              unit="kg*1e-3/m2/s" /> 
    333       
     428 
    334429     <!-- sbcssm variables --> 
    335430          <field id="sst_m"    unit="degC" /> 
    … …  
    342437 
    343438   </field_group> 
    344     
    345    <!-- scalar variables --> 
    346    <field_group id="SBC_0D"  grid_ref="grid_1point" > 
     439 
     440 
     441      </field_group> <!-- SBC --> 
     442 
     443      <!-- ABL --> 
     444      <field_group id="ABL" > <!-- time step automaticaly defined based on nn_fsbc --> 
     445 
     446   <!-- variables available with ABL on atmospheric T grid--> 
     447   <field_group id="grid_ABL3D" grid_ref="grid_TA_3D" > 
     448          <field id="u_abl"      long_name="ABL i-horizontal velocity"     standard_name="abl_x_velocity" unit="m/s"      /> 
     449          <field id="v_abl"      long_name="ABL j-horizontal velocity"     standard_name="abl_y_velocity" unit="m/s"      /> 
     450          <field id="t_abl"      long_name="ABL potential temperature"     standard_name="abl_theta"      unit="K"        /> 
     451          <field id="q_abl"      long_name="ABL specific humidity"         standard_name="abl_qspe"       unit="kg/kg"    /> 
     452          <!-- debug (to be removed) --> 
     453          <field id="u_dta"      long_name="DTA i-horizontal velocity"     standard_name="dta_x_velocity" unit="m/s"      /> 
     454          <field id="v_dta"      long_name="DTA j-horizontal velocity"     standard_name="dta_y_velocity" unit="m/s"      /> 
     455          <field id="t_dta"      long_name="DTA potential temperature"     standard_name="dta_theta"      unit="K"        /> 
     456          <field id="q_dta"      long_name="DTA specific humidity"         standard_name="dta_qspe"       unit="kg/kg"    /> 
     457          <field id="u_geo"      long_name="GEO i-horizontal velocity"     standard_name="geo_x_velocity" unit="m/s"      /> 
     458          <field id="v_geo"      long_name="GEO j-horizontal velocity"     standard_name="geo_y_velocity" unit="m/s"      /> 
     459          <field id="tke_abl"    long_name="ABL turbulent kinetic energy"  standard_name="abl_tke"        unit="m2/s2"    /> 
     460          <field id="avm_abl"    long_name="ABL turbulent viscosity"       standard_name="abl_avm"        unit="m2/s"     /> 
     461          <field id="avt_abl"    long_name="ABL turbulent diffusivity"     standard_name="abl_avt"        unit="m2/s"     /> 
     462          <field id="mxlm_abl"   long_name="ABL master mixing length"      standard_name="abl_mxlm"       unit="m"        /> 
     463          <field id="mxld_abl"   long_name="ABL dissipative mixing length" standard_name="abl_mxld"       unit="m"        /> 
    347464   </field_group> 
    348465 
    349       </field_group> <!-- SBC --> 
     466   <field_group id="grid_ABL2D" grid_ref="grid_TA_2D" > 
     467          <field id="pblh"       long_name="ABL height"                    standard_name="abl_height"     unit="m"        /> 
     468          <field id="uz1_abl"    long_name="ABL i-horizontal velocity"     standard_name="abl_x_velocity" unit="m/s"      /> 
     469          <field id="vz1_abl"    long_name="ABL j-horizontal velocity"     standard_name="abl_y_velocity" unit="m/s"      /> 
     470          <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> 
     471          <field id="tz1_abl"    long_name="ABL potential temperature"     standard_name="abl_theta"      unit="K"        /> 
     472          <field id="qz1_abl"    long_name="ABL specific humidity"         standard_name="abl_qspe"       unit="kg/kg"    /> 
     473          <field id="uz1_dta"    long_name="DTA i-horizontal velocity"     standard_name="dta_x_velocity" unit="m/s"      /> 
     474          <field id="vz1_dta"    long_name="DTA j-horizontal velocity"     standard_name="dta_y_velocity" unit="m/s"      /> 
     475          <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> 
     476          <field id="tz1_dta"    long_name="DTA potential temperature"     standard_name="dta_theta"      unit="K"        /> 
     477          <field id="qz1_dta"    long_name="DTA specific humidity"         standard_name="dta_qspe"       unit="kg/kg"    /> 
     478          <!-- debug (to be removed) --> 
     479          <field id="uz1_geo"    long_name="GEO i-horizontal velocity"     standard_name="geo_x_velocity" unit="m/s"      /> 
     480          <field id="vz1_geo"    long_name="GEO j-horizontal velocity"     standard_name="geo_y_velocity" unit="m/s"      /> 
     481          <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> 
     482   </field_group> 
     483 
     484      </field_group> <!-- ABL --> 
     485 
    350486 
    351487      <!-- U grid --> 
    352        
     488 
    353489      <field_group id="grid_U"   grid_ref="grid_U_2D"> 
    354         <field id="e3u"          long_name="U-cell thickness"                                       standard_name="cell_thickness"              unit="m"          grid_ref="grid_U_3D" /> 
    355         <field id="e3u_0"        long_name="Initial U-cell thickness"                               standard_name="ref_cell_thickness"          unit="m"          grid_ref="grid_U_3D"/> 
    356         <field id="utau"         long_name="Wind Stress along i-axis"                               standard_name="surface_downward_x_stress"   unit="N/m2"                            /> 
    357         <field id="uoce"         long_name="ocean current along i-axis"                             standard_name="sea_water_x_velocity"        unit="m/s"        grid_ref="grid_U_3D" /> 
    358         <field id="uoce_e3u"     long_name="ocean current along i-axis  (thickness weighted)"                                                   unit="m/s"        grid_ref="grid_U_3D"  > uoce * e3u </field> 
     490        <field id="e2u"           long_name="U-cell width in meridional direction"                   standard_name="cell_width"                  unit="m"                               /> 
     491        <field id="e3u"           long_name="U-cell thickness"                                       standard_name="cell_thickness"              unit="m"          grid_ref="grid_U_3D" /> 
     492        <field id="e3u_0"         long_name="Initial U-cell thickness"                               standard_name="ref_cell_thickness"          unit="m"          grid_ref="grid_U_3D"/> 
     493        <field id="utau"          long_name="Wind Stress along i-axis"                               standard_name="surface_downward_x_stress"   unit="N/m2"                            /> 
     494        <field id="uoce"          long_name="ocean current along i-axis"                             standard_name="sea_water_x_velocity"        unit="m/s"        grid_ref="grid_U_3D" /> 
     495      <field id="uoce_e3u"      long_name="ocean current along i-axis  (thickness weighted)"                                                   unit="m/s"        grid_ref="grid_U_3D"  > uoce * e3u </field> 
     496      <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"/> 
     497        <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> 
     498 
     499        <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> 
     500        <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" /> 
     501        <field id="msftbarot"         long_name="ocean_barotropic_mass_streamfunction"   unit="kg s-1" > uocetr_vsum_cumul * $rau0 </field> 
     502 
     503 
    359504        <field id="ssu"          long_name="ocean surface current along i-axis"                                                                 unit="m/s"                             /> 
    360505        <field id="sbu"          long_name="ocean bottom current along i-axis"                                                                  unit="m/s"                             /> 
    … …  
    363508        <field id="uocet"        long_name="ocean transport along i-axis times temperature (CRS)"                                               unit="degC*m/s"   grid_ref="grid_U_3D" /> 
    364509        <field id="uoces"        long_name="ocean transport along i-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_U_3D" /> 
    365  
     510        <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> 
     511        <!-- AGRIF sponge --> 
     512        <field id="agrif_spu"    long_name=" AGRIF u-sponge coefficient"   unit=" " /> 
    366513        <!-- u-eddy diffusivity coefficients (available if ln_traldf_OFF=F) --> 
    367514        <field id="ahtu_2d"      long_name=" surface u-eddy diffusivity coefficient"   unit="m2/s or m4/s" /> 
    … …  
    375522 
    376523        <!-- uoce_eiv: available EIV (ln_ldfeiv=T and ln_ldfeiv_dia=T) --> 
    377         <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" /> 
     524      <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" /> 
     525        <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" /> 
     526        <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"                         /> 
     527        <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"                        /> 
     528        <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" /> 
     529        <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" /> 
    378530 
    379531        <!-- uoce_bbl: available with ln_trabbl=T and nn_bbl_adv=1 --> 
    … …  
    389541        <field id="utbl"         long_name="zonal current in the Losh tbl"     unit="m/s" /> 
    390542 
    391         <field id="u_masstr"      long_name="Ocean Mass X Transport"                                          standard_name="ocean_mass_x_transport"                         unit="kg/s" grid_ref="grid_U_3D" /> 
     543        <!-- variables available with diaar5 --> 
     544        <field id="u_masstr"      long_name="Ocean Mass X Transport"                                          standard_name="ocean_mass_x_transport"                         unit="kg/s"  grid_ref="grid_U_3D" /> 
    392545        <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"                      /> 
    393546        <field id="u_heattr"      long_name="ocean eulerian heat transport along i-axis"                      standard_name="ocean_heat_x_transport"                         unit="W"                         /> 
    … …  
    395548        <field id="uadv_heattr"   long_name="ocean advective heat transport along i-axis"                     standard_name="advectice_ocean_heat_x_transport"               unit="W"                         /> 
    396549        <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"                 /> 
    397         <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"                         /> 
    398         <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"                        /> 
    399         <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" /> 
    400         <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" /> 
    401550        <field id="udiff_heattr"  long_name="ocean diffusion heat transport along i-axis"                     standard_name="ocean_heat_x_transport_due_to_diffusion"        unit="W"                         /> 
    402551        <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"                 /> 
    403552      </field_group> 
    404        
     553 
    405554      <!-- V grid --> 
    406        
     555 
    407556      <field_group id="grid_V"   grid_ref="grid_V_2D"> 
     557        <field id="e1v"          long_name="V-cell width in longitudinal direction"                 standard_name="cell_width"                  unit="m"                              /> 
    408558        <field id="e3v"          long_name="V-cell thickness"                                       standard_name="cell_thickness"              unit="m"          grid_ref="grid_V_3D" /> 
    409559        <field id="e3v_0"        long_name="Initial V-cell thickness"                               standard_name="ref_cell_thickness"          unit="m"          grid_ref="grid_V_3D" /> 
    … …  
    417567        <field id="vocet"        long_name="ocean transport along j-axis times temperature (CRS)"                                               unit="degC*m/s"   grid_ref="grid_V_3D" /> 
    418568        <field id="voces"        long_name="ocean transport along j-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_V_3D" /> 
    419  
     569        <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> 
     570        <!-- AGRIF sponge --> 
     571        <field id="agrif_spv"    long_name=" AGRIF v-sponge coefficient"   unit=" " /> 
    420572        <!-- v-eddy diffusivity coefficients (available if ln_traldf_OFF=F) --> 
    421573        <field id="ahtv_2d"      long_name=" surface v-eddy diffusivity coefficient"     unit="m2/s or (m4/s)^1/2" /> 
    … …  
    429581 
    430582        <!-- voce_eiv: available EIV (ln_ldfeiv=T and ln_ldfeiv_dia=T)  --> 
    431         <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" /> 
     583   <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" /> 
     584        <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" /> 
     585        <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"                         /> 
     586        <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"                        /> 
     587        <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" /> 
     588        <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" /> 
     589 
    432590 
    433591        <!-- voce_bbl: available with ln_trabbl=T and nn_bbl_adv=1 --> 
    … …  
    449607        <field id="vadv_heattr"   long_name="ocean advective heat transport along j-axis"   standard_name="advectice_ocean_heat_y_transport"                unit="W"                         /> 
    450608        <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"                 /> 
    451         <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"                         /> 
    452         <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"                        /> 
    453         <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" /> 
    454         <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" /> 
    455609        <field id="vdiff_heattr"  long_name="ocean diffusion heat transport along j-axis"   standard_name="ocean_heat_y_transport_due_to_diffusion"         unit="W"                         /> 
    456610        <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"                 /> 
    457611      </field_group> 
    458        
     612 
    459613      <!-- W grid --> 
    460        
     614 
    461615      <field_group id="grid_W" grid_ref="grid_W_3D"> 
    462616        <field id="e3w"          long_name="W-cell thickness"                              standard_name="cell_thickness"                         unit="m"    /> 
    463         <field id="woce"         long_name="ocean vertical velocity"                       standard_name="upward_sea_water_velocity"              unit="m/s"  /> 
     617   <field id="woce"         long_name="ocean vertical velocity"                       standard_name="upward_sea_water_velocity"              unit="m/s"  /> 
     618   <field id="woce_e3w"     long_name="ocean vertical velocity * e3w"                                                                        unit="m2/s"  > woce * e3w </field> 
    464619        <field id="wocetr_eff"   long_name="effective ocean vertical transport"                                                                   unit="m3/s" /> 
    465620 
    466621        <!-- woce_eiv: available with EIV  (ln_ldfeiv=T and ln_ldfeiv_dia=T)  --> 
    467622        <field id="woce_eiv"     long_name="EIV ocean vertical velocity"                    standard_name="bolus_upward_sea_water_velocity"       unit="m/s"  /> 
    468  
    469         <field id="avt"          long_name="vertical eddy diffusivity"                      standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
     623        <field id="weiv_masstr"  long_name="EIV Upward Ocean Mass Transport"  standard_name="bolus_upward_ocean_mass_transport"             unit="kg/s"   /> 
     624        <field id="weiv_heattr3d" long_name="ocean bolus heat transport"    standard_name="ocean_heat_z_transport_due_to_bolus_advection"   unit="W"    /> 
     625        <field id="weiv_salttr3d" long_name="ocean bolus salt transport"    standard_name="ocean_salt_z_transport_due_to_bolus_advection"   unit="kg"   /> 
     626 
     627   <field id="avt"          long_name="vertical eddy diffusivity"                      standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
     628        <field id="avt_e3w"      long_name="vertical heat diffusivity * e3w"                unit="m3/s" > avt * e3w </field> 
    470629        <field id="logavt"       long_name="logarithm of vertical eddy diffusivity"         standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
    471630        <field id="avm"          long_name="vertical eddy viscosity"                        standard_name="ocean_vertical_momentum_diffusivity"   unit="m2/s" /> 
     631        <field id="avm_e3w"      long_name="vertical eddy viscosity * e3w"   unit="m3/s" > avm * e3w </field> 
    472632 
    473633        <!-- avs: /= avt with ln_zdfddm=T --> 
    474634        <field id="avs"          long_name="salt vertical eddy diffusivity"                 standard_name="ocean_vertical_salt_diffusivity"       unit="m2/s" /> 
    475         <field id="logavs"       long_name="logarithm of salt vertical eddy diffusivity"    standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
     635        <field id="avs_e3w"      long_name="vertical salt diffusivity * e3w"   unit="m3/s" > avs * e3w </field> 
     636   <field id="logavs"       long_name="logarithm of salt vertical eddy diffusivity"    standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
    476637 
    477638        <!-- avt_evd and avm_evd: available with ln_zdfevd --> 
    478639        <field id="avt_evd"      long_name="convective enhancement of vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_convection"     unit="m2/s" /> 
    479         <field id="avm_evd"      long_name="convective enhancement of vertical viscosity"   standard_name="ocean_vertical_momentum_diffusivity_due_to_convection"   unit="m2/s" /> 
     640        <field id="avt_evd_e3w"  long_name="convective enhancement to vertical diffusivity * e3w "    unit="m3/s" > avt_evd * e3w </field> 
     641   <field id="avm_evd"      long_name="convective enhancement of vertical viscosity"   standard_name="ocean_vertical_momentum_diffusivity_due_to_convection"   unit="m2/s" /> 
    480642 
    481643        <!-- avt_tide: available with ln_zdfiwm=T --> 
    … …  
    489651        <field id="wstokes"      long_name="Stokes Drift vertical velocity"                 standard_name="upward_StokesDrift_velocity"   unit="m/s" /> 
    490652 
    491         <!-- variables available with diaar5 -->    
     653        <!-- variables available with diaar5 --> 
    492654        <field id="w_masstr"     long_name="vertical mass transport"                        standard_name="upward_ocean_mass_transport"             unit="kg/s"   /> 
    493655        <field id="w_masstr2"    long_name="square of vertical mass transport"              standard_name="square_of_upward_ocean_mass_transport"   unit="kg2/s2" /> 
    494656 
    495657      </field_group> 
    496        
     658 
    497659      <!-- F grid --> 
     660      <!-- AGRIF sponge --> 
     661      <field id="agrif_spf"    long_name=" AGRIF f-sponge coefficient"   unit=" " /> 
    498662      <!-- f-eddy viscosity coefficients (ldfdyn) --> 
    499663      <field id="ahmf_2d"      long_name=" surface f-eddy viscosity coefficient"   unit="m2/s or m4/s" /> 
    500664      <field id="ahmf_3d"      long_name=" 3D      f-eddy viscosity coefficient"   unit="m2/s or m4/s"  grid_ref="grid_T_3D"/> 
    501665 
    502       <field_group id="scalar"  grid_ref="grid_T_2D"  > 
     666      <field_group id="scalar"  grid_ref="grid_scalar"  > 
    503667         <!-- Need to have a salinity reference climatological file : sali_ref_clim_monthly --> 
    504668        <field id="voltot"     long_name="global total volume"                          standard_name="sea_water_volume"                               unit="m3"   /> 
    … …  
    508672        <field id="masstot"    long_name="global total mass"                            standard_name="sea_water_mass"                                 unit="kg"   /> 
    509673        <field id="temptot"    long_name="global mean temperature"                      standard_name="sea_water_potential_temperature"                unit="degC" /> 
    510         <field id="saltot"     long_name="global mean salinity"                         standard_name="sea_water_salinity"                             unit="1e-3" /> 
    511         <field id="fram_trans" long_name="Sea Ice Mass Transport Through Fram Strait"   standard_name="sea_ice_transport_across_line"                  unit="kg/s" /> 
     674   <field id="saltot"     long_name="global mean salinity"                         standard_name="sea_water_salinity"                             unit="1e-3" /> 
     675        <field id="ssttot"     long_name="global mean sea surface temperature"          standard_name="sea_water_conservative_temperature"             unit="degC" /> 
     676         <!-- EOS = TEOS-10 --> 
     677        <field id="temptot_pot" long_name="global mean potential temperature"            standard_name="sea_water_potential_temperature"               unit="degC" /> 
    512678 
    513679         <!-- available with ln_diahsb --> 
    … …  
    526692        <field id="bgmissal"     long_name="global mean salinity error due to free surface (linssh true)"                                                                          unit="1e-3"     /> 
    527693      </field_group> 
    528        
    529       <!-- variables available with key_float --> 
     694 
     695 
     696      <!-- transects --> 
     697      <field_group id="oce_straits"> 
     698         <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> 
     699         <field id="uoce_e3u_ave_vsum"    long_name="Vertical sum of u*e3u"                           field_ref="uoce_e3u_ave"         grid_ref="grid_U_vsum"     /> 
     700         <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> 
     701         <field id="uocetr_strait"        long_name="Total transport across lines in i-direction"     field_ref="uocetr_vsum_section"  grid_ref="grid_U_4strait" /> 
     702         <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> 
     703 
     704         <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> 
     705         <field id="voce_e3v_ave_vsum"    long_name="Vertical sum of v*e3v"                           field_ref="voce_e3v_ave"         grid_ref="grid_V_vsum"      /> 
     706         <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> 
     707         <field id="vocetr_strait"        long_name="Total transport across lines in j-direction"     field_ref="vocetr_vsum_section"  grid_ref="grid_V_4strait"  /> 
     708         <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> 
     709 
     710         <field id="masstr_strait"        long_name="Sea water transport across line"                                                  grid_ref="grid_4strait"  > u_masstr_strait + v_masstr_strait </field> 
     711      </field_group> 
     712 
     713 
     714      <!-- variables available with ln_floats --> 
    530715 
    531716      <field_group id="floatvar" grid_ref="grid_T_nfloat"  operation="instant" > 
    … …  
    541726      <!-- variables available with iceberg trajectories --> 
    542727 
    543       <field_group id="icbvar" domain_ref="grid_T"  >  
     728      <field_group id="icbvar" domain_ref="grid_T"  > 
    544729        <field id="berg_melt"          long_name="icb melt rate of icebergs"                       unit="kg/m2/s"                    /> 
    545730        <field id="berg_melt_hcflx"    long_name="icb heat flux to ocean due to melting heat content"   unit="J/m2/s"                /> 
    … …  
    559744      </field_group> 
    560745 
    561       <!-- Poleward transport : ptr -->      
    562       <field_group id="diaptr" >   
    563         <field id="zomsfglo"          long_name="Meridional Stream-Function: Global"           unit="Sv"        grid_ref="gznl_W_3D" /> 
    564         <field id="zomsfatl"          long_name="Meridional Stream-Function: Atlantic"         unit="Sv"        grid_ref="gznl_W_3D" /> 
    565         <field id="zomsfpac"          long_name="Meridional Stream-Function: Pacific"          unit="Sv"        grid_ref="gznl_W_3D" /> 
    566         <field id="zomsfind"          long_name="Meridional Stream-Function: Indian"           unit="Sv"        grid_ref="gznl_W_3D" /> 
    567         <field id="zomsfipc"          long_name="Meridional Stream-Function: Pacific+Indian"   unit="Sv"        grid_ref="gznl_W_3D" /> 
    568         <field id="zotemglo"          long_name="Zonal Mean Temperature : Global"              unit="degree_C"  grid_ref="gznl_T_3D" /> 
    569         <field id="zotematl"          long_name="Zonal Mean Temperature : Atlantic"            unit="degree_C"  grid_ref="gznl_T_3D" /> 
    570         <field id="zotempac"          long_name="Zonal Mean Temperature : Pacific"             unit="degree_C"  grid_ref="gznl_T_3D" /> 
    571         <field id="zotemind"          long_name="Zonal Mean Temperature : Indian"              unit="degree_C"  grid_ref="gznl_T_3D" /> 
    572         <field id="zotemipc"          long_name="Zonal Mean Temperature : Pacific+Indian"      unit="degree_C"  grid_ref="gznl_T_3D" /> 
    573         <field id="zosalglo"          long_name="Zonal Mean Salinity : Global"                 unit="0.001"     grid_ref="gznl_T_3D" /> 
    574         <field id="zosalatl"          long_name="Zonal Mean Salinity : Atlantic"               unit="0.001"     grid_ref="gznl_T_3D" /> 
    575         <field id="zosalpac"          long_name="Zonal Mean Salinity : Pacific"                unit="0.001"     grid_ref="gznl_T_3D" /> 
    576         <field id="zosalind"          long_name="Zonal Mean Salinity : Indian"                 unit="0.001"     grid_ref="gznl_T_3D" /> 
    577         <field id="zosalipc"          long_name="Zonal Mean Salinity : Pacific+Indian"         unit="0.001"     grid_ref="gznl_T_3D" /> 
    578         <field id="zosrfglo"          long_name="Zonal Mean Surface"                           unit="m2"        grid_ref="gznl_T_3D" /> 
    579         <field id="zosrfatl"          long_name="Zonal Mean Surface : Atlantic"                unit="m2"        grid_ref="gznl_T_3D" /> 
    580         <field id="zosrfpac"          long_name="Zonal Mean Surface : Pacific"                 unit="m2"        grid_ref="gznl_T_3D" /> 
    581         <field id="zosrfind"          long_name="Zonal Mean Surface : Indian"                  unit="m2"        grid_ref="gznl_T_3D" /> 
    582         <field id="zosrfipc"          long_name="Zonal Mean Surface : Pacific+Indian"          unit="m2"        grid_ref="gznl_T_3D" /> 
    583         <field id="sophtadv"          long_name="Advective Heat Transport"                     unit="PW"        grid_ref="gznl_T_2D" /> 
    584         <field id="sophtadv_atl"      long_name="Advective Heat Transport: Atlantic"           unit="PW"        grid_ref="gznl_T_2D" /> 
    585         <field id="sophtadv_pac"      long_name="Advective Heat Transport: Pacific"            unit="PW"        grid_ref="gznl_T_2D" /> 
    586         <field id="sophtadv_ind"      long_name="Advective Heat Transport: Indian"             unit="PW"        grid_ref="gznl_T_2D" /> 
    587         <field id="sophtadv_ipc"      long_name="Advective Heat Transport: Pacific+Indian"     unit="PW"        grid_ref="gznl_T_2D" /> 
    588         <field id="sophtldf"          long_name="Diffusive Heat Transport"                     unit="PW"        grid_ref="gznl_T_2D" /> 
    589         <field id="sophtldf_atl"      long_name="Diffusive Heat Transport: Atlantic"           unit="PW"        grid_ref="gznl_T_2D" /> 
    590         <field id="sophtldf_pac"      long_name="Diffusive Heat Transport: Pacific"            unit="PW"        grid_ref="gznl_T_2D" /> 
    591         <field id="sophtldf_ind"      long_name="Diffusive Heat Transport: Indian"             unit="PW"        grid_ref="gznl_T_2D" /> 
    592         <field id="sophtldf_ipc"      long_name="Diffusive Heat Transport: Pacific+Indian"     unit="PW"        grid_ref="gznl_T_2D" /> 
    593         <field id="sophtove"          long_name="Overturning Heat Transport"                   unit="PW"        grid_ref="gznl_T_2D" /> 
    594         <field id="sophtove_atl"      long_name="Overturning Heat Transport: Atlantic"         unit="PW"        grid_ref="gznl_T_2D" /> 
    595         <field id="sophtove_pac"      long_name="Overturning Heat Transport: Pacific"          unit="PW"        grid_ref="gznl_T_2D" /> 
    596         <field id="sophtove_ind"      long_name="Overturning Heat Transport: Indian"           unit="PW"        grid_ref="gznl_T_2D" /> 
    597         <field id="sophtove_ipc"      long_name="Overturning Heat Transport: Pacific+Indian"   unit="PW"        grid_ref="gznl_T_2D" /> 
    598         <field id="sophtbtr"          long_name="Barotropic Heat Transport"                    unit="PW"        grid_ref="gznl_T_2D" /> 
    599         <field id="sophtbtr_atl"      long_name="Barotropic Heat Transport: Atlantic"          unit="PW"        grid_ref="gznl_T_2D" /> 
    600         <field id="sophtbtr_pac"      long_name="Barotropic Heat Transport: Pacific"           unit="PW"        grid_ref="gznl_T_2D" /> 
    601         <field id="sophtbtr_ind"      long_name="Barotropic Heat Transport: Indian"            unit="PW"        grid_ref="gznl_T_2D" /> 
    602         <field id="sophtbtr_ipc"      long_name="Barotropic Heat Transport: Pacific+Indian"    unit="PW"        grid_ref="gznl_T_2D" /> 
    603         <field id="sophteiv"          long_name="Heat Transport from mesoscale eddy advection"                     unit="PW"       grid_ref="gznl_T_2D" /> 
    604         <field id="sophteiv_atl"      long_name="Heat Transport from mesoscale eddy advection: Atlantic"           unit="PW"       grid_ref="gznl_T_2D" /> 
    605         <field id="sophteiv_pac"      long_name="Heat Transport from mesoscale eddy advection: Pacific"            unit="PW"       grid_ref="gznl_T_2D" /> 
    606         <field id="sophteiv_ind"      long_name="Heat Transport from mesoscale eddy advection: Indian"             unit="PW"       grid_ref="gznl_T_2D" /> 
    607         <field id="sophteiv_ipc"      long_name="Heat Transport from mesoscale eddy advection: Pacific+Indian"     unit="PW"       grid_ref="gznl_T_2D" /> 
    608         <field id="sopstadv"          long_name="Advective Salt Transport"                     unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    609         <field id="sopstadv_atl"      long_name="Advective Salt Transport: Atlantic"           unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    610         <field id="sopstadv_pac"      long_name="Advective Salt Transport: Pacific"            unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    611         <field id="sopstadv_ind"      long_name="Advective Salt Transport: Indian"             unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    612         <field id="sopstadv_ipc"      long_name="Advective Salt Transport: Pacific+Indian"     unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    613         <field id="sopstove"          long_name="Overturning Salt Transport"                   unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    614         <field id="sopstove_atl"      long_name="Overturning Salt Transport: Atlantic"         unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    615         <field id="sopstove_pac"      long_name="Overturning Salt Transport: Pacific"          unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    616         <field id="sopstove_ind"      long_name="Overturning Salt Transport: Indian"           unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    617         <field id="sopstove_ipc"      long_name="Overturning Salt Transport: Pacific+Indian"   unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    618         <field id="sopstbtr"          long_name="Barotropic Salt Transport"                    unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    619         <field id="sopstbtr_atl"      long_name="Barotropic Salt Transport: Atlantic"          unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    620         <field id="sopstbtr_pac"      long_name="Barotropic Salt Transport: Pacific"           unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    621         <field id="sopstbtr_ind"      long_name="Barotropic Salt Transport: Indian"            unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    622         <field id="sopstbtr_ipc"      long_name="Barotropic Salt Transport: Pacific+Indian"    unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    623         <field id="sopstldf"          long_name="Diffusive Salt Transport"                     unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    624         <field id="sopstldf_atl"      long_name="Diffusive Salt Transport: Atlantic"           unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    625         <field id="sopstldf_pac"      long_name="Diffusive Salt Transport: Pacific"            unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    626         <field id="sopstldf_ind"      long_name="Diffusive Salt Transport: Indian"             unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    627         <field id="sopstldf_ipc"      long_name="Diffusive Salt Transport: Pacific+Indian"     unit="Giga g/s"  grid_ref="gznl_T_2D" /> 
    628         <field id="sopsteiv"          long_name="Salt Transport from mesoscale eddy advection"                     unit="Giga g/s"       grid_ref="gznl_T_2D" /> 
    629         <field id="sopsteiv_atl"      long_name="Salt Transport from mesoscale eddy advection: Atlantic"           unit="Giga g/s"       grid_ref="gznl_T_2D" /> 
    630         <field id="sopsteiv_pac"      long_name="Salt Transport from mesoscale eddy advection: Pacific"            unit="Giga g/s"       grid_ref="gznl_T_2D" /> 
    631         <field id="sopsteiv_ind"      long_name="Salt Transport from mesoscale eddy advection: Indian"             unit="Giga g/s"       grid_ref="gznl_T_2D" /> 
    632         <field id="sopsteiv_ipc"      long_name="Salt Transport from mesoscale eddy advection: Pacific+Indian"     unit="Giga g/s"       grid_ref="gznl_T_2D" />        
    633       </field_group> 
    634  
    635     <!--  
     746      <!-- Poleward transport : ptr --> 
     747      <field_group id="diaptr" > 
     748        <field id="zomsf"         long_name="Overturning Stream-Function : All basins"                     unit="Sv"         grid_ref="grid_znl_W_3D" /> 
     749        <field id="zotem"         long_name="Zonal Mean Temperature : All basins"                          unit="degree_C"   grid_ref="grid_znl_T_3D" /> 
     750        <field id="zosal"         long_name="Zonal Mean Salinity : All basins"                             unit="0.001"      grid_ref="grid_znl_T_3D" /> 
     751        <field id="zosrf"         long_name="Zonal Mean Surface : All basins"                              unit="m2"         grid_ref="grid_znl_T_3D" /> 
     752        <field id="sophtove"      long_name="Overturning Heat Transport: All basins"                       unit="PW"         grid_ref="grid_znl_T_2D" /> 
     753        <field id="sopstove"      long_name="Overturning Salt Transport: All basins"                       unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
     754        <field id="sophtbtr"      long_name="Barotropic Heat Transport: All basins"                        unit="PW"         grid_ref="grid_znl_T_2D" /> 
     755        <field id="sopstbtr"      long_name="Barotropic Salt Transport: All basins"                        unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
     756        <field id="sophtadv"      long_name="Advective Heat Transport: All basins"                         unit="PW"         grid_ref="grid_znl_T_2D" /> 
     757        <field id="sopstadv"      long_name="Advective Salt Transport: All basins"                         unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
     758        <field id="sophtldf"      long_name="Diffusive Heat Transport: All basins"                         unit="PW"         grid_ref="grid_znl_T_2D" /> 
     759        <field id="sopstldf"      long_name="Diffusive Salt Transport: All basins"                         unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
     760        <field id="sophtvtr"      long_name="Heat Transport : All basins"                                  unit="PW"         grid_ref="grid_znl_T_2D" /> 
     761        <field id="sopstvtr"      long_name="Salt Transport : All basins"                                  unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
     762        <field id="sophteiv"      long_name="Heat Transport from mesoscale eddy advection: All basins"     unit="PW"         grid_ref="grid_znl_T_2D" /> 
     763        <field id="sopsteiv"      long_name="Salt Transport from mesoscale eddy advection : All basins"    unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
     764   <field id="sopstadv"      long_name="Advective Salt Transport"                                     unit="Giga g/s"   grid_ref="grid_znl_T_2D" /> 
     765        <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> 
     766   <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> 
     767       </field_group> 
     768 
     769       <field_group id="constant_fields"  grid_ref="grid_T_2D"  operation="once" > 
     770         <field id="bathy"       long_name="Sea floor depth below geoid"               standard_name="sea_floor_depth_below_geoid"               unit="m"/> 
     771         <field id="areacello"   long_name="Horizontal area of ocean grid cells"       standard_name="cell_area"                                 unit="m2" /> 
     772         <field id="hfgeou"      long_name="Upward geothermal heat flux at sea floor"  standard_name="upward_geothermal_heat_flux_at_sea_floor"  unit="W/m2"/> 
     773         <field id="basins"      long_name="ocean tracer region masks"                 standard_name="ocean_tracer_region_masks"   unit="none" grid_ref="grid_basin"  /> 
     774       </field_group> 
     775 
     776 
     777    <!-- 
    636778============================================================================================================ 
    637779                  Physical ocean model trend diagnostics : temperature, KE, PE, momentum 
    … …  
    664806     <field id="ttrd_totad"    long_name="temperature-trend: total advection"         standard_name="tendency_of_sea_water_salinity_due_to_advection"              unit="degC/s"                        /> 
    665807     <field id="strd_totad"    long_name="salinity   -trend: total advection"         standard_name="tendency_of_sea_water_salinity_due_to_advection"              unit="1e-3/s"                        /> 
    666      <field id="ttrd_sad"      long_name="temperature-trend: surface adv. (linssh true)"                                                                                unit="degC/s"   grid_ref="grid_T_2D" /> 
    667      <field id="strd_sad"      long_name="salinity   -trend: surface adv. (linssh true)"                                                                                unit="1e-3/s"   grid_ref="grid_T_2D" /> 
     808     <field id="ttrd_sad"      long_name="temperature-trend: surface adv. (linssh true)"                                                                           unit="degC/s"   grid_ref="grid_T_2D" /> 
     809     <field id="strd_sad"      long_name="salinity   -trend: surface adv. (linssh true)"                                                                           unit="1e-3/s"   grid_ref="grid_T_2D" /> 
    668810     <field id="ttrd_ldf"      long_name="temperature-trend: lateral  diffusion"      standard_name="tendency_of_sea_water_temperature_due_to_horizontal_mixing"   unit="degC/s"                        /> 
    669811     <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"                        /> 
    … …  
    722864     <field id="ttrd_npc_e3t"      unit="degC/s * m"  >  ttrd_npc * e3t </field> 
    723865     <field id="strd_npc_e3t"      unit="1e-3/s * m"  >  strd_npc * e3t </field> 
    724      <field id="ttrd_qns_e3t"      unit="degC/s * m"  >  ttrd_qns * e3t_surf </field> 
    725      <field id="strd_cdt_e3t"      unit="degC/s * m"  >  strd_cdt * e3t_surf </field> 
     866     <field id="ttrd_qns_e3t"      unit="degC/s * m"  >  ttrd_qns * e3ts </field> 
     867     <field id="strd_cdt_e3t"      unit="degC/s * m"  >  strd_cdt * e3ts </field> 
    726868     <field id="ttrd_qsr_e3t"      unit="degC/s * m"  >  ttrd_qsr * e3t </field> 
    727869     <field id="ttrd_bbc_e3t"      unit="degC/s * m"  >  ttrd_bbc * e3t </field> 
    … …  
    774916     <field id="ketrd_zdf"     long_name="ke-trend: vertical  diffusion"                    unit="W/s^3"                        /> 
    775917     <field id="ketrd_tau"     long_name="ke-trend: wind stress "                           unit="W/s^3"   grid_ref="grid_T_2D" /> 
    776      <field id="ketrd_bfr"     long_name="ke-trend: bottom friction (explicit)"             unit="W/s^3"                        />    
    777      <field id="ketrd_bfri"    long_name="ke-trend: bottom friction (implicit)"             unit="W/s^3"                        />    
    778      <field id="ketrd_atf"     long_name="ke-trend: asselin time filter trend"              unit="W/s^3"                        />   
     918     <field id="ketrd_bfr"     long_name="ke-trend: bottom friction (explicit)"             unit="W/s^3"                        /> 
     919     <field id="ketrd_bfri"    long_name="ke-trend: bottom friction (implicit)"             unit="W/s^3"                        /> 
     920     <field id="ketrd_atf"     long_name="ke-trend: asselin time filter trend"              unit="W/s^3"                        /> 
    779921     <field id="ketrd_convP2K" long_name="ke-trend: conversion (potential to kinetic)"      unit="W/s^3"                        /> 
    780      <field id="KE"            long_name="kinetic energy: u(n)*u(n+1)/2"                    unit="W/s^2"                        />    
     922     <field id="KE"            long_name="kinetic energy: u(n)*u(n+1)/2"                    unit="W/s^2"                        /> 
    781923 
    782924     <!-- variables available when explicit lateral mixing is used (ln_dynldf_OFF=F) --> 
    783      <field id="dispkexyfo"    long_name="KE-trend: lateral  mixing induced dissipation"   standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_xy_friction"                   unit="W/m^2" grid_ref="grid_T_2D" />    
    784      <field id="dispkevfo"     long_name="KE-trend: vertical mixing induced dissipation"   standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_vertical_friction"             unit="W/m^2" grid_ref="grid_T_2D" />    
     925     <field id="dispkexyfo"    long_name="KE-trend: lateral  mixing induced dissipation"   standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_xy_friction"                   unit="W/m^2" grid_ref="grid_T_2D" /> 
     926     <field id="dispkevfo"     long_name="KE-trend: vertical mixing induced dissipation"   standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_vertical_friction"             unit="W/m^2" grid_ref="grid_T_2D" /> 
    785927     <!-- variables available with ln_traadv_eiv=T and ln_diaeiv=T --> 
    786      <field id="eketrd_eiv"    long_name="EKE-trend due to parameterized eddy advection"   standard_name="tendency_of_ocean_eddy_kinetic_energy_content_due_to_parameterized_eddy_advection"   unit="W/m^2" grid_ref="grid_T_2D" />    
     928     <field id="eketrd_eiv"    long_name="EKE-trend due to parameterized eddy advection"   standard_name="tendency_of_ocean_eddy_kinetic_energy_content_due_to_parameterized_eddy_advection"   unit="W/m^2" grid_ref="grid_T_2D" /> 
    787929 
    788930     <!-- variables available with ln_PE_trd --> 
    … …  
    801943     <field id="petrd_bbc"     long_name="pe-trend: geothermal heating"         unit="W/m^3"                        /> 
    802944     <field id="petrd_atf"     long_name="pe-trend: asselin time filter"        unit="W/m^3"                        /> 
    803      <field id="PEanom"        long_name="potential energy anomaly"             unit="1"                            />    
    804      <field id="alphaPE"       long_name="partial deriv. of PEanom wrt T"       unit="degC-1"                       />    
    805      <field id="betaPE"        long_name="partial deriv. of PEanom wrt S"       unit="1e3"                          />    
     945     <field id="PEanom"        long_name="potential energy anomaly"             unit="1"                            /> 
     946     <field id="alphaPE"       long_name="partial deriv. of PEanom wrt T"       unit="degC-1"                       /> 
     947     <field id="betaPE"        long_name="partial deriv. of PEanom wrt S"       unit="1e3"                          /> 
    806948   </field_group> 
    807949 
    … …  
    820962     <field id="utrd_zdf"       long_name="i-trend: vertical  diffusion"                    unit="m/s^2"                        /> 
    821963     <field id="utrd_tau"       long_name="i-trend: wind stress "                           unit="m/s^2"   grid_ref="grid_U_2D" /> 
    822      <field id="utrd_bfr"       long_name="i-trend: bottom friction (explicit)"             unit="m/s^2"                        />    
    823      <field id="utrd_bfri"      long_name="i-trend: bottom friction (implicit)"             unit="m/s^2"                        />    
    824      <field id="utrd_tot"       long_name="i-trend: total momentum trend before atf"        unit="m/s^2"                        />    
    825      <field id="utrd_atf"       long_name="i-trend: asselin time filter trend"              unit="m/s^2"                        />    
     964     <field id="utrd_bfr"       long_name="i-trend: bottom friction (explicit)"             unit="m/s^2"                        /> 
     965     <field id="utrd_bfri"      long_name="i-trend: bottom friction (implicit)"             unit="m/s^2"                        /> 
     966     <field id="utrd_tot"       long_name="i-trend: total momentum trend before atf"        unit="m/s^2"                        /> 
     967     <field id="utrd_atf"       long_name="i-trend: asselin time filter trend"              unit="m/s^2"                        /> 
    826968   </field_group> 
    827969 
    … …  
    840982     <field id="vtrd_zdf"       long_name="j-trend: vertical  diffusion"                    unit="m/s^2"                        /> 
    841983     <field id="vtrd_tau"       long_name="j-trend: wind stress "                           unit="m/s^2"   grid_ref="grid_V_2D" /> 
    842      <field id="vtrd_bfr"       long_name="j-trend: bottom friction (explicit)"             unit="m/s^2"                        />    
    843      <field id="vtrd_bfri"      long_name="j-trend: bottom friction (implicit)"             unit="m/s^2"                        />    
    844      <field id="vtrd_tot"       long_name="j-trend: total momentum trend before atf"        unit="m/s^2"                        />    
    845      <field id="vtrd_atf"       long_name="j-trend: asselin time filter trend"              unit="m/s^2"                        />    
     984     <field id="vtrd_bfr"       long_name="j-trend: bottom friction (explicit)"             unit="m/s^2"                        /> 
     985     <field id="vtrd_bfri"      long_name="j-trend: bottom friction (implicit)"             unit="m/s^2"                        /> 
     986     <field id="vtrd_tot"       long_name="j-trend: total momentum trend before atf"        unit="m/s^2"                        /> 
     987     <field id="vtrd_atf"       long_name="j-trend: asselin time filter trend"              unit="m/s^2"                        /> 
    846988   </field_group> 
    847989 
    848990 
    849     <!--  
     991    <!-- 
    850992============================================================================================================ 
    851993                                        Definitions for iodef_demo.xml 
    … …  
    8651007      <field field_ref="strd_zdfp_li"    name="osaltdiff" /> 
    8661008    </field_group> 
    867      
     1009 
    8681010    <field_group id="mooring" > 
    8691011      <field field_ref="toce"         name="thetao"   long_name="sea_water_potential_temperature"      /> 
    … …  
    8741016      <field field_ref="avt"          name="difvho"   long_name="ocean_vertical_heat_diffusivity"      /> 
    8751017      <field field_ref="avm"          name="difvmo"   long_name="ocean_vertical_momentum_diffusivity"  /> 
    876        
     1018 
    8771019      <field field_ref="sst"          name="tos"      long_name="sea_surface_temperature"                       /> 
    8781020      <field field_ref="sst2"         name="tossq"    long_name="square_of_sea_surface_temperature"             /> 
    … …  
    9221064      <field field_ref="BLT"          name="blt"      long_name="Barrier Layer Thickness"                       /> 
    9231065    </field_group> 
    924      
     1066 
    9251067    <field_group id="groupU" > 
    9261068      <field field_ref="uoce"         name="uo"      long_name="sea_water_x_velocity"      /> 
    927       <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"    /> 
    9281069      <field field_ref="utau"         name="tauuo"   long_name="surface_downward_x_stress" /> 
    9291070    </field_group> 
    930      
     1071 
    9311072    <field_group id="groupV" > 
    9321073      <field field_ref="voce"         name="vo"      long_name="sea_water_y_velocity"      /> 
    933       <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"    /> 
    9341074      <field field_ref="vtau"         name="tauvo"   long_name="surface_downward_y_stress" /> 
    9351075    </field_group> 
    936      
     1076 
    9371077    <field_group id="groupW" > 
    9381078      <field field_ref="woce"         name="wo"       long_name="ocean vertical velocity"  /> 
    939     </field_group> 
    940  
    941     <!-- TMB diagnostic output --> 
    942     <field_group  id="1h_grid_T_tmb" grid_ref="grid_T_2D" operation="instant"> 
    943       <field id="top_temp"           name="votemper_top"  unit="degC"  /> 
    944       <field id="mid_temp"           name="votemper_mid"  unit="degC"  /> 
    945       <field id="bot_temp"           name="votemper_bot"  unit="degC"  /> 
    946       <field id="top_sal"            name="vosaline_top"  unit="psu"   /> 
    947       <field id="mid_sal"            name="vosaline_mid"  unit="psu"   /> 
    948       <field id="bot_sal"            name="vosaline_bot"  unit="psu"   /> 
    949       <field id="sshnmasked"         name="sossheig"      unit="m"     />  
    9501079    </field_group> 
    9511080 
    … …  
    9881117    </field_group> 
    9891118 
    990     <!--  
     1119    <!-- 
    9911120============================================================================================================ 
    9921121    --> 
    993     <!-- output variables for my configuration (example) -->  
    994      
     1122    <!-- output variables for my configuration (example) --> 
     1123 
    9951124    <field_group id="myvarOCE" > 
    996       <!-- grid T -->  
     1125      <!-- grid T --> 
    9971126      <field field_ref="e3t"          name="e3t"      long_name="vertical scale factor"           /> 
    9981127      <field field_ref="sst"          name="tos"      long_name="sea_surface_temperature"         /> 
    9991128      <field field_ref="sss"          name="sos"      long_name="sea_surface_salinity"            /> 
    10001129      <field field_ref="ssh"          name="zos"      long_name="sea_surface_height_above_geoid"  /> 
    1001        
    1002       <!-- grid U -->  
     1130 
     1131      <!-- grid U --> 
    10031132      <field field_ref="e3u"          name="e3u"     long_name="vertical scale factor"            /> 
    10041133      <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"           /> 
    1005        
    1006       <!-- grid V -->  
     1134 
     1135      <!-- grid V --> 
    10071136      <field field_ref="e3v"          name="e3v"     long_name="vertical scale factor"            /> 
    1008       <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"           />      
    1009     </field_group>     
     1137      <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"           /> 
     1138    </field_group> 
    10101139 
    10111140   </field_definition> 

  • NEMO/branches/2019/dev_r11351_fldread_with_XIOS/cfgs/SHARED/field_def_nemo-pisces.xml

    r10416 r13463  
    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 --> 
    … …  
    290301     </field_group> 
    291302 
    292      <field_group id="tracer_scalar"  grid_ref="grid_T_2D"  > 
     303     <field_group id="tracer_scalar"  grid_ref="grid_scalar"  > 
    293304     <!-- PISCES scalar  --> 
    294305       <field id="pno3tot"         long_name="Global mean nitrate concentration"                  unit="mol/m3"   /> 

  • NEMO/branches/2019/dev_r11351_fldread_with_XIOS/cfgs/SHARED/grid_def_nemo.xml

    r10226 r13463  
     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" /> 
    7        </grid> 
    8         <!--  --> 
    9        <grid id="grid_T_3D_ncatice" > 
    10          <domain id="grid_T" /> 
    11          <axis id="ncatice" /> 
     12         <domain domain_ref="grid_T" /> 
     13       </grid> 
     14        <!--  --> 
     15       <grid id="grid_T_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        <!--  --> 
    … …  
    5460        <!--  --> 
    5561       <grid id="grid_T_nfloat" > 
    56          <domain id="grid_T" /> 
    57          <axis id="nfloat" /> 
    58        </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  
     263      <!-- scalars --> 
     264      <grid id="grid_scalar" > 
     265        <scalar/> 
     266      </grid> 
     267 
     268      <!-- ABL grid definition --> 
     269      <grid id="grid_TA_2D"> 
     270        <domain domain_ref="grid_T" /> 
     271      </grid> 
     272      <grid id="grid_TA_3D"> 
     273        <domain domain_ref="grid_T" /> 
     274        <axis  id="ght_abl" /> 
     275      </grid> 
     276      <grid id="grid_WA_3D"> 
     277        <domain domain_ref="grid_T" /> 
     278        <axis  id="ghw_abl" /> 
     279      </grid> 
     280      <!--  --> 
     281 
     282      <!-- grid definitions for multiple-linear-regression analysis (diamlr) --> 
     283      <grid id="diamlr_grid_scalar" > 
     284   <scalar /> 
     285   <scalar /> 
     286      </grid> 
     287      <grid id="diamlr_grid_T_2D" > 
     288   <domain domain_ref="grid_T" /> 
     289   <scalar /> 
     290      </grid> 
     291      <grid id="diamlr_grid_U_2D" > 
     292   <domain domain_ref="grid_U" /> 
     293   <scalar /> 
     294      </grid> 
     295      <grid id="diamlr_grid_V_2D" > 
     296   <domain domain_ref="grid_V" /> 
     297   <scalar /> 
     298      </grid> 
     299      <grid id="diamlr_grid_W_2D" > 
     300   <domain domain_ref="grid_W" /> 
     301   <scalar /> 
     302      </grid> 
     303      <grid id="diamlr_grid_2D_to_grid_T_3D" > 
     304   <domain domain_ref="grid_T" /> 
     305        <axis axis_ref="deptht"> 
     306          <duplicate_scalar /> 
     307        </axis> 
     308      </grid> 
     309      <grid id="diamlr_grid_2D_to_grid_U_3D" > 
     310   <domain domain_ref="grid_U" /> 
     311        <axis axis_ref="depthu"> 
     312          <duplicate_scalar /> 
     313        </axis> 
     314      </grid> 
     315      <grid id="diamlr_grid_2D_to_grid_V_3D" > 
     316   <domain domain_ref="grid_V" /> 
     317        <axis axis_ref="depthv"> 
     318          <duplicate_scalar /> 
     319        </axis> 
     320      </grid> 
     321      <grid id="diamlr_grid_2D_to_grid_W_3D" > 
     322   <domain domain_ref="grid_W" /> 
     323        <axis axis_ref="depthw"> 
     324          <duplicate_scalar /> 
     325        </axis> 
     326      </grid> 
     327      <grid id="diamlr_grid_2D_to_scalar" > 
     328        <scalar> 
     329          <reduce_domain operation="average" /> 
     330        </scalar> 
     331   <scalar /> 
     332      </grid> 
     333      <!-- grid definitions for the computation of daily detided model diagnostics (diadetide) --> 
     334      <grid id="diadetide_grid_T_2D" > 
     335   <domain domain_ref="grid_T" /> 
     336   <scalar /> 
     337      </grid> 
     338      <grid id="diadetide_grid_U_2D" > 
     339   <domain domain_ref="grid_U" /> 
     340   <scalar /> 
     341      </grid> 
     342      <grid id="diadetide_grid_V_2D" > 
     343   <domain domain_ref="grid_V" /> 
     344   <scalar /> 
     345      </grid> 
     346      <grid id="diadetide_grid_2D_to_grid_T_3D" > 
     347   <domain domain_ref="grid_T" /> 
     348        <axis axis_ref="deptht"> 
     349          <duplicate_scalar /> 
     350        </axis> 
     351      </grid> 
     352      <grid id="diadetide_grid_2D_to_grid_U_3D" > 
     353   <domain domain_ref="grid_U" /> 
     354        <axis axis_ref="depthu"> 
     355          <duplicate_scalar /> 
     356        </axis> 
     357      </grid> 
     358      <grid id="diadetide_grid_2D_to_grid_V_3D" > 
     359   <domain domain_ref="grid_V" /> 
     360        <axis axis_ref="depthv"> 
     361          <duplicate_scalar /> 
     362        </axis> 
     363      </grid> 
     364      <grid id="diadetide_grid_2D_to_grid_W_3D" > 
     365   <domain domain_ref="grid_W" /> 
     366        <axis axis_ref="depthw"> 
     367          <duplicate_scalar /> 
     368        </axis> 
     369      </grid> 
    59370 
    60371    </grid_definition> 
    61      

  • NEMO/branches/2019/dev_r11351_fldread_with_XIOS/cfgs/SHARED/namelist_ice_ref

    r10911 r13463  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! SI3 namelist:   
     2!! SI3 Reference namelist: 
    33!!              1 - Generic parameters                 (nampar) 
    44!!              2 - Ice thickness discretization       (namitd) 
    … …  
    5656   rn_ishlat        =   2.            !  lbc : free slip (0) ; partial slip (0-2) ; no slip (2) ; strong slip (>2) 
    5757   ln_landfast_L16  = .false.         !  landfast: parameterization from Lemieux 2016 
    58    ln_landfast_home = .false.         !  landfast: parameterization from "home made" 
    5958      rn_depfra     =   0.125         !        fraction of ocean depth that ice must reach to initiate landfast 
    60                                       !          recommended range: [0.1 ; 0.25] - L16=0.125 - home=0.15 
    61       rn_icebfr     =  15.            !        ln_landfast_L16:  maximum bottom stress per unit volume [N/m3] 
    62                                       !        ln_landfast_home: maximum bottom stress per unit area of contact [N/m2] 
    63                                       !          recommended range: ?? L16=15 - home=10 
     59                                      !          recommended range: [0.1 ; 0.25] 
     60      rn_icebfr     =  15.            !        maximum bottom stress per unit volume [N/m3] 
    6461      rn_lfrelax    =   1.e-5         !        relaxation time scale to reach static friction [s-1] 
    65       rn_tensile    =   0.2           !        ln_landfast_L16: isotropic tensile strength 
     62      rn_tensile    =   0.05          !        isotropic tensile strength [0-0.5??] 
    6663/ 
    6764!------------------------------------------------------------------------------ 
    … …  
    7168   ln_str_H79       = .true.          !  ice strength param.: Hibler_79   => P = pstar*<h>*exp(-c_rhg*A)                       
    7269      rn_pstar      =   2.0e+04       !     ice strength thickness parameter [N/m2] 
    73       rn_crhg       =   20.0          !     ice strength conc. parameter (-) 
     70      rn_crhg       =  20.0           !     ice strength conc. parameter (-) 
    7471                   ! -- ice_rdgrft -- ! 
    7572   rn_csrdg         =   0.5           !  fraction of shearing energy contributing to ridging 
    … …  
    104101&namdyn_adv     !   Ice advection 
    105102!------------------------------------------------------------------------------ 
    106    ln_adv_Pra       = .false.         !  Advection scheme (Prather) 
    107    ln_adv_UMx       = .true.          !  Advection scheme (Ultimate-Macho) 
     103   ln_adv_Pra       = .true.         !  Advection scheme (Prather) 
     104   ln_adv_UMx       = .false.          !  Advection scheme (Ultimate-Macho) 
    108105      nn_UMx        =   5             !     order of the scheme for UMx (1-5 ; 20=centered 2nd order) 
    109106/ 
    … …  
    178175&namthd_pnd     !   Melt ponds 
    179176!------------------------------------------------------------------------------ 
    180    ln_pnd_H12       = .false.         !  activate evolutive melt ponds (from Holland et al 2012) 
    181    ln_pnd_CST       = .false.         !  activate constant  melt ponds 
    182       rn_apnd       =   0.2           !     prescribed pond fraction, at Tsu=0 degC 
    183       rn_hpnd       =   0.05          !     prescribed pond depth,    at Tsu=0 degC 
    184    ln_pnd_alb       = .false.         !  melt ponds affect albedo or not 
     177   ln_pnd           = .false.         !  activate melt ponds or not 
     178     ln_pnd_H12     = .false.         !  activate evolutive melt ponds (from Holland et al 2012) 
     179     ln_pnd_CST     = .false.         !  activate constant  melt ponds 
     180       rn_apnd      =   0.2           !     prescribed pond fraction, at Tsu=0 degC 
     181       rn_hpnd      =   0.05          !     prescribed pond depth,    at Tsu=0 degC 
     182     ln_pnd_alb     = .false.         !  melt ponds affect albedo or not 
    185183/ 
    186184!------------------------------------------------------------------------------ 
    … …  
    189187   ln_iceini        = .true.          !  activate ice initialization (T) or not (F) 
    190188   ln_iceini_file   = .false.         !  netcdf file provided for initialization (T) or not (F) 
    191    rn_thres_sst     =   2.0           !  max delta temp. above Tfreeze with initial ice = (sst - tfreeze) 
    192    rn_hts_ini_n     =   0.3           !  initial real snow thickness (m), North 
     189   rn_thres_sst     =   2.0           !  max temp. above Tfreeze with initial ice = (sst - tfreeze) 
     190   rn_hti_ini_n     =   3.0           !  initial ice thickness       (m), North 
     191   rn_hti_ini_s     =   1.0           !        "            "             South 
     192   rn_hts_ini_n     =   0.3           !  initial snow thickness      (m), North 
    193193   rn_hts_ini_s     =   0.3           !        "            "             South 
    194    rn_hti_ini_n     =   3.0           !  initial real ice thickness  (m), North 
    195    rn_hti_ini_s     =   1.0           !        "            "             South 
    196194   rn_ati_ini_n     =   0.9           !  initial ice concentration   (-), North 
    197195   rn_ati_ini_s     =   0.9           !        "            "             South 
    198196   rn_smi_ini_n     =   6.3           !  initial ice salinity     (g/kg), North 
    199197   rn_smi_ini_s     =   6.3           !        "            "             South 
    200    rn_tmi_ini_n     = 270.            !  initial ice/snw temperature (K), North 
     198   rn_tmi_ini_n     = 270.            !  initial ice temperature    (K), North 
    201199   rn_tmi_ini_s     = 270.            !        "            "             South 
    202  
     200   rn_tsu_ini_n     = 270.            !  initial surface temperature (K), North 
     201   rn_tsu_ini_s     = 270.            !        "            "             South 
     202   rn_tms_ini_n     = 270.            !  initial snw temperature     (K), North 
     203   rn_tms_ini_s     = 270.            !        "            "             South 
     204   rn_apd_ini_n     =   0.2           !  initial pond fraction       (-), North 
     205   rn_apd_ini_s     =   0.2           !        "            "             South 
     206   rn_hpd_ini_n     =   0.05          !  initial pond depth          (m), North 
     207   rn_hpd_ini_s     =   0.05          !        "            "             South 
     208   ! -- for ln_iceini_file = T 
    203209   sn_hti = 'Ice_initialization'    , -12 ,'hti'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    204210   sn_hts = 'Ice_initialization'    , -12 ,'hts'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    205211   sn_ati = 'Ice_initialization'    , -12 ,'ati'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     212   sn_smi = 'Ice_initialization'    , -12 ,'smi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     213   sn_tmi = 'Ice_initialization'    , -12 ,'tmi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    206214   sn_tsu = 'Ice_initialization'    , -12 ,'tsu'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    207    sn_tmi = 'Ice_initialization'    , -12 ,'tmi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    208    sn_smi = 'Ice_initialization'    , -12 ,'smi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     215   sn_tms = 'NOT USED'              , -12 ,'tms'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     216   !      melt ponds (be careful, sn_apd is the pond concentration (not fraction), so it differs from rn_apd) 
     217   sn_apd = 'NOT USED'              , -12 ,'apd'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     218   sn_hpd = 'NOT USED'              , -12 ,'hpd'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
    209219   cn_dir='./' 
    210220/ 
    … …  
    222232&namdia         !   Diagnostics 
    223233!------------------------------------------------------------------------------ 
    224    ln_icediachk     = .false.         !  check online the heat, mass & salt budgets (T) or not (F) 
     234   ln_icediachk     = .false.         !  check online heat, mass & salt budgets 
     235      !                               !   rate of ice spuriously gained/lost at each time step => rn_icechk=1 <=> 1.e-6 m/hour 
     236      rn_icechk_cel =  100.           !     check at each gridcell          (1.e-4m/h)=> stops the code if violated (and writes a file) 
     237      rn_icechk_glo =  1.             !     check over the entire ice cover (1.e-6m/h)=> only prints warnings 
    225238   ln_icediahsb     = .false.         !  output the heat, mass & salt budgets (T) or not (F) 
    226239   ln_icectl        = .false.         !  ice points output for debug (T or F) 

  • NEMO/branches/2019/dev_r11351_fldread_with_XIOS/cfgs/SHARED/namelist_pisces_ref

    r10721 r13463  
    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      ! 
    349349   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'  , ''       , ''    , '' 
    359350   sn_ironsed  = 'bathy.orca'      ,    -12            , 'bathy'       ,  .false.     , .true. , 'yearly'  , ''       , ''    , '' 
    360351   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_r11351_fldread_with_XIOS/cfgs/SHARED/namelist_ref

    r11485 r13463  
    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 
    … …  
    5050      cn_ocerst_indir = "."         !  directory from which to read input ocean restarts 
    5151      cn_ocerst_out   = "restart"   !  suffix of ocean restart name (output) 
    52       cn_ocerst_outdir = "."         !  directory in which to write output ocean restarts 
    53    ln_iscpl    = .false.   !  cavity evolution forcing or coupling to ice sheet model 
     52      cn_ocerst_outdir = "."        !  directory in which to write output ocean restarts 
    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) 
    56    nn_stock    =    5840   !  frequency of creation of a restart file (modulo referenced to 1) 
     55   nn_stock    =       0   !  used only if ln_rst_list = F: output restart freqeuncy (modulo referenced to 1) 
     56      !                          !    =  0 force to write restart files only at the end of the run 
     57      !                          !    = -1 do not do any restart 
    5758   nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written 
    58    nn_write    =    5840   !  frequency of write in the output file   (modulo referenced to nn_it000) 
    59    ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
     59   nn_write    =       0   !  used only if key_iomput is not defined: output frequency (modulo referenced to nn_it000) 
     60      !                          !    =  0 force to write output files only at the end of the run 
     61      !                          !    = -1 do not do any output file 
     62   ln_mskland  = .false.   !  mask land points in NetCDF outputs 
    6063   ln_cfmeta   = .false.   !  output additional data to netCDF files required for compliance with the CF metadata standard 
    6164   ln_clobber  = .true.    !  clobber (overwrite) an existing file 
    6265   nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 
    63    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) 
    6467   nn_wxios = 0      !  use XIOS to write restart file 0 - no, 1 - single file output, 2 - multiple file output 
    6568/ 
    … …  
    6871!----------------------------------------------------------------------- 
    6972   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    70    rn_isfhmin  =    1.00   !  treshold [m] to discriminate grounding ice from floating ice 
    71    ! 
    72    rn_rdt      = 5400.     !  time step for the dynamics and tracer 
     73   ! 
     74   rn_Dt       = 5400.     !  time step for the dynamics and tracer 
    7375   rn_atfp     =    0.1    !  asselin time filter parameter 
    7476   ! 
    7577   ln_crs      = .false.   !  Logical switch for coarsening module      (T => fill namcrs) 
    7678   ! 
    77    ln_meshmask = .false.   !  =T create a mesh file 
     79   ln_meshmask = .true.   !  =T create a mesh file 
    7880/ 
    7981!----------------------------------------------------------------------- 
    8082&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg) 
    8183!----------------------------------------------------------------------- 
    82    ln_read_cfg = .false.   !  (=T) read the domain configuration file 
    83       !                    !  (=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) 
    8486      cn_domcfg = "domain_cfg"  ! domain configuration filename 
    8587      ! 
    86       ln_closea    = .false.    !  T => keep closed seas (defined by closea_mask field) in the   
    87       !                         !       domain and apply special treatment of freshwater fluxes. 
    88       !                         !  F => suppress closed seas (defined by closea_mask field)  
    89       !                         !       from the bathymetry at runtime. 
    90       !                         !  If closea_mask field doesn't exist in the domain_cfg file 
    91       !                         !       then this logical does nothing. 
    92    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 
    9392      cn_domcfg_out = "domain_cfg_out" ! newly created domain configuration filename 
    9493      ! 
    95    ln_use_jattr = .false.  !  use (T) the file attribute: open_ocean_jstart, if present 
    96    !                       !  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 
    97110/ 
    98111!----------------------------------------------------------------------- 
    … …  
    103116   ln_tsd_dmp  = .false.         !  T-S restoring   (see namtra_dmp) 
    104117   ln_tsd_xios = .false.         !  read tsd using XIOS 
    105     
    106118   cn_dir      = './'      !  root directory for the T-S data location 
    107119   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
    108120   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
    109121   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
    110    sn_tem = 'data_1m_potential_temperature_nomask',  -1      , 'votemper',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
    111    sn_sal = 'data_1m_salinity_nomask'             ,  -1      , 'vosaline',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     122   sn_tem = 'data_1m_potential_temperature_nomask',  -1.     , 'votemper',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     123   sn_sal = 'data_1m_salinity_nomask'             ,  -1.     , 'vosaline',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
    112124/ 
    113125!----------------------------------------------------------------------- 
    … …  
    164176   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
    165177   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
    166    sn_ucur     = 'ucurrent'              ,         -1        ,'u_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Ume'   , '' 
    167    sn_vcur     = 'vcurrent'              ,         -1        ,'v_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Vme'   , '' 
     178   sn_ucur     = 'ucurrent'              ,         -1.       ,'u_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Ume'   , '' 
     179   sn_vcur     = 'vcurrent'              ,         -1.       ,'v_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Vme'   , '' 
    168180/ 
    169181 
    … …  
    181193!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T) 
    182194!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T) 
    183 !!   namsbc_isf      ice shelf melting/freezing                         (ln_isfcav  =T : read (ln_read_cfg=T) or set or usr_def_zgr ) 
    184 !!   namsbc_iscpl    coupling option between land ice model and ocean   (ln_isfcav  =T) 
    185195!!   namsbc_wave     external fields from wave model                    (ln_wave    =T) 
    186196!!   namberg         iceberg floats                                     (ln_icebergs=T) 
    … …  
    192202   nn_fsbc     = 2         !  frequency of SBC module call 
    193203      !                    !  (control sea-ice & iceberg model call) 
    194                      ! Type of air-sea fluxes  
     204                     ! Type of air-sea fluxes 
    195205   ln_usr      = .false.   !  user defined formulation                  (T => check usrdef_sbc) 
    196206   ln_flx      = .false.   !  flux formulation                          (T => fill namsbc_flx ) 
    197207   ln_blk      = .false.   !  Bulk formulation                          (T => fill namsbc_blk ) 
     208   ln_abl      = .false.   !  ABL  formulation                          (T => fill namsbc_abl ) 
    198209      !              ! Type of coupling (Ocean/Ice/Atmosphere) : 
    199210   ln_cpl      = .false.   !  atmosphere coupled   formulation          ( requires key_oasis3 ) 
    … …  
    202213      !                    !  =0 no opa-sas OASIS coupling: default single executable config. 
    203214      !                    !  =1 opa-sas OASIS coupling: multi executable config., OPA component 
    204       !                    !  =2 opa-sas OASIS coupling: multi executable config., SAS component  
     215      !                    !  =2 opa-sas OASIS coupling: multi executable config., SAS component 
    205216                     ! Sea-ice : 
    206    nn_ice      = 0         !  =0 no ice boundary condition     
     217   nn_ice      = 0         !  =0 no ice boundary condition 
    207218      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif ) 
    208       !                    !  =2 or 3 automatically for SI3 or CICE    ("key_si3" or "key_cice") 
    209       !                    !          except in AGRIF zoom where it has to be specified 
     219      !                    !  =2 or 3 for SI3 and CICE, respectively 
    210220   ln_ice_embd = .false.   !  =T embedded sea-ice (pressure + mass and salt exchanges) 
    211221      !                    !  =F levitating ice (no pressure, mass and salt exchanges) 
    212                      ! Misc. options of sbc :  
     222                     ! Misc. options of sbc : 
    213223   ln_traqsr   = .false.   !  Light penetration in the ocean            (T => fill namtra_qsr) 
    214224   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave 
    … …  
    219229   ln_rnf      = .false.   !  runoffs                                   (T => fill namsbc_rnf) 
    220230   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr ) 
    221    ln_isf      = .false.   !  ice shelf                                 (T   => fill namsbc_isf & namsbc_iscpl) 
    222231   ln_wave     = .false.   !  Activate coupling with wave  (T => fill namsbc_wave) 
    223232   ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave) 
    224    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) 
    225234   nn_sdrift   =  0        !  Parameterization for the calculation of 3D-Stokes drift from the surface Stokes drift 
    226235      !                    !   = 0 Breivik 2015 parameterization: v_z=v_0*[exp(2*k*z)/(1-8*k*z)] 
    … …  
    240249   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
    241250   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
    242    sn_utau     = 'utau'                  ,        24         , 'utau'    , .false.     , .false., 'yearly'  , ''               , ''       , '' 
    243    sn_vtau     = 'vtau'                  ,        24         , 'vtau'    , .false.     , .false., 'yearly'  , ''               , ''       , '' 
    244    sn_qtot     = 'qtot'                  ,        24         , 'qtot'    , .false.     , .false., 'yearly'  , ''               , ''       , '' 
    245    sn_qsr      = 'qsr'                   ,        24         , 'qsr'     , .false.     , .false., 'yearly'  , ''               , ''       , '' 
    246    sn_emp      = 'emp'                   ,        24         , 'emp'     , .false.     , .false., 'yearly'  , ''               , ''       , '' 
    247 / 
    248 !----------------------------------------------------------------------- 
    249 &namsbc_blk    !   namsbc_blk  generic Bulk formula                     (ln_blk =T) 
     251   sn_utau     = 'utau'                  ,        24.        , 'utau'    , .false.     , .false., 'yearly'  , ''               , ''       , '' 
     252   sn_vtau     = 'vtau'                  ,        24.        , 'vtau'    , .false.     , .false., 'yearly'  , ''               , ''       , '' 
     253   sn_qtot     = 'qtot'                  ,        24.        , 'qtot'    , .false.     , .false., 'yearly'  , ''               , ''       , '' 
     254   sn_qsr      = 'qsr'                   ,        24.        , 'qsr'     , .false.     , .false., 'yearly'  , ''               , ''       , '' 
     255   sn_emp      = 'emp'                   ,        24.        , 'emp'     , .false.     , .false., 'yearly'  , ''               , ''       , '' 
     256/ 
     257!----------------------------------------------------------------------- 
     258&namsbc_blk    !   namsbc_blk  generic Bulk formula          (ln_blk =T) 
    250259!----------------------------------------------------------------------- 
    251260   !                    !  bulk algorithm : 
    252    ln_NCAR     = .false.   ! "NCAR"      algorithm   (Large and Yeager 2008) 
     261   ln_NCAR      = .true.    ! "NCAR"      algorithm   (Large and Yeager 2008) 
    253262   ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003) 
    254    ln_COARE_3p5 = .false.   ! "COARE 3.5" algorithm   (Edson et al. 2013) 
    255    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) 
    256265      ! 
    257       rn_zqt      = 10.       !  Air temperature & humidity reference height (m) 
    258       rn_zu       = 10.       !  Wind vector reference height (m) 
    259       ln_Cd_L12   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2012) 
    260       ln_Cd_L15   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2015) 
    261       ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data 
    262       rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    263       rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
    264       rn_vfac     = 0.        !  multiplicative factor for ocean & ice velocity used to 
    265       !                       !  calculate the wind stress (0.=absolute or 1.=relative winds) 
    266  
     266      rn_zqt     = 10.      !  Air temperature & humidity reference height (m) 
     267      rn_zu      = 10.      !  Wind vector reference height (m) 
     268      ln_Cd_L12  = .false.  !  air-ice drags = F(ice conc.) (Lupkes et al. 2012) 
     269      ln_Cd_L15  = .false.  !  air-ice drags = F(ice conc.) (Lupkes et al. 2015) 
     270      ln_crt_fbk = .false.     !  Add surface current feedback to the wind stress (Renault et al. 2020, doi: 10.1029/2019MS001715) 
     271         rn_stau_a = -2.9e-3   !     Alpha from eq. 10: Stau = Alpha * Wnd + Beta 
     272         rn_stau_b =  8.0e-3   !     Beta  
     273      rn_pfac    = 1.       !  multipl. factor for precipitation (total & snow) 
     274      rn_efac    = 1.       !  multipl. factor for evaporation (0. or 1.) 
     275      ln_skin_cs = .false.  !  use the cool-skin parameterization 
     276      ln_skin_wl = .false.  !  use the warm-layer parameterization 
     277      !                     !   ==> only available in ECMWF and COARE algorithms 
     278      ln_humi_sph = .true.  !  humidity "sn_humi" is specific humidity  [kg/kg] 
     279      ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K] 
     280      ln_humi_rlh = .false. !  humidity "sn_humi" is relative humidity     [%] 
     281      ln_tpot     = .true.  !!GS: compute potential temperature or not 
     282   ! 
    267283   cn_dir      = './'      !  root directory for the bulk data location 
    268284   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________! 
    269285   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename               ! rotation ! land/sea mask ! 
    270286   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                                      ! pairing  !    filename   ! 
    271    sn_wndi     = 'u_10.15JUNE2009_fill'       ,    6         , 'U_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'  , 'Uwnd'   , '' 
    272    sn_wndj     = 'v_10.15JUNE2009_fill'       ,    6         , 'V_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'  , 'Vwnd'   , '' 
    273    sn_qsr      = 'ncar_rad.15JUNE2009_fill'   ,   24         , 'SWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    274    sn_qlw      = 'ncar_rad.15JUNE2009_fill'   ,   24         , 'LWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    275    sn_tair     = 't_10.15JUNE2009_fill'       ,    6         , 'T_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    276    sn_humi     = 'q_10.15JUNE2009_fill'       ,    6         , 'Q_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    277    sn_prec     = 'ncar_precip.15JUNE2009_fill',   -1         , 'PRC_MOD1',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    278    sn_snow     = 'ncar_precip.15JUNE2009_fill',   -1         , 'SNOW'    ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    279    sn_slp      = 'slp.15JUNE2009_fill'        ,    6         , 'SLP'     ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    280    sn_tdif     = 'taudif_core'                ,   24         , 'taudif'  ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     287   sn_wndi     = 'u_10.15JUNE2009_fill'       ,    6.        , 'U_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'  , 'Uwnd'   , '' 
     288   sn_wndj     = 'v_10.15JUNE2009_fill'       ,    6.        , 'V_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'  , 'Vwnd'   , '' 
     289   sn_qsr      = 'ncar_rad.15JUNE2009_fill'   ,   24.        , 'SWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     290   sn_qlw      = 'ncar_rad.15JUNE2009_fill'   ,   24.        , 'LWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     291   sn_tair     = 't_10.15JUNE2009_fill'       ,    6.        , 'T_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     292   sn_humi     = 'q_10.15JUNE2009_fill'       ,    6.        , 'Q_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     293   sn_prec     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'PRC_MOD1',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     294   sn_snow     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'SNOW'    ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     295   sn_slp      = 'slp.15JUNE2009_fill'        ,    6.        , 'SLP'     ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     296   sn_uoatm    = 'NOT USED'                   ,    6.        , 'UOATM'   ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , 'Uoceatm', '' 
     297   sn_voatm    = 'NOT USED'                   ,    6.        , 'VOATM'   ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , 'Voceatm', '' 
     298   sn_hpgi     = 'NOT USED'                   ,   24.        , 'uhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'UG'     , '' 
     299   sn_hpgj     = 'NOT USED'                   ,   24.        , 'vhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'VG'     , '' 
     300/ 
     301!----------------------------------------------------------------------- 
     302&namsbc_abl    !   Atmospheric Boundary Layer formulation           (ln_abl = T) 
     303!----------------------------------------------------------------------- 
     304   cn_dir           = './'      !  root directory for the location of the ABL grid file 
     305   cn_dom           = 'dom_cfg_abl.nc' 
     306 
     307   cn_ablrst_in     = "restart_abl"   !  suffix of abl restart name (input) 
     308   cn_ablrst_out    = "restart_abl"   !  suffix of abl restart name (output) 
     309   cn_ablrst_indir  = "."             !  directory to read   input abl restarts 
     310   cn_ablrst_outdir = "."             !  directory to write output abl restarts 
     311 
     312   ln_rstart_abl  = .false. 
     313   ln_hpgls_frc   = .false. 
     314   ln_geos_winds  = .false. 
     315   ln_smth_pblh   = .false. 
     316   nn_dyn_restore = 0         ! restoring option for dynamical ABL variables: = 0 no restoring 
     317                              !                                               = 1 equatorial restoring 
     318                              !                                               = 2 global restoring 
     319   rn_ldyn_min   =  4.5       ! dynamics nudging magnitude inside the ABL [hour] (~3 rn_Dt) 
     320   rn_ldyn_max   =  1.5       ! dynamics nudging magnitude above  the ABL [hour] (~1 rn_Dt) 
     321   rn_ltra_min   =  4.5       ! tracers  nudging magnitude inside the ABL [hour] (~3 rn_Dt) 
     322   rn_ltra_max   =  1.5       ! tracers  nudging magnitude above  the ABL [hour] (~1 rn_Dt) 
     323   nn_amxl       =  0         ! mixing length: = 0 Deardorff 80 length-scale 
     324                              !                = 1 length-scale based on the distance to the PBL height 
     325                              !                = 2 Bougeault & Lacarrere 89 length-scale 
     326                              ! CBR00  ! CCH02  ! MesoNH ! 
     327   rn_Cm          = 0.0667    ! 0.0667 ! 0.1260 ! 0.1260 ! 
     328   rn_Ct          = 0.1667    ! 0.1667 ! 0.1430 ! 0.1430 ! 
     329   rn_Ce          = 0.40      ! 0.40   ! 0.34   ! 0.40   ! 
     330   rn_Ceps        = 0.700     ! 0.700  ! 0.845  ! 0.850  ! 
     331   rn_Ric         = 0.139     ! 0.139  ! 0.143  !   ?    ! Critical Richardson number (to compute PBL height and diffusivities) 
     332   rn_Rod         = 0.15      ! c0 in RMCA17 mixing length formulation (not yet implemented) 
    281333/ 
    282334!----------------------------------------------------------------------- 
    … …  
    287339   !                       !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 
    288340   nn_cats_cpl   =     5   !  Number of sea ice categories over which coupling is to be carried out (if not 1) 
    289  
    290341   !_____________!__________________________!____________!_____________!______________________!________! 
    291342   !             !        description       !  multiple  !    vector   !       vector         ! vector ! 
    … …  
    324375   sn_rcv_wper   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
    325376   sn_rcv_wnum   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
    326    sn_rcv_wstrf  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     377   sn_rcv_wfreq  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
    327378   sn_rcv_wdrag  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
    328379   sn_rcv_ts_ice =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
    … …  
    344395   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
    345396   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
    346    sn_usp      = 'sas_grid_U'            ,       120         , 'uos'     ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
    347    sn_vsp      = 'sas_grid_V'            ,       120         , 'vos'     ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
    348    sn_tem      = 'sas_grid_T'            ,       120         , 'sosstsst',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
    349    sn_sal      = 'sas_grid_T'            ,       120         , 'sosaline',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
    350    sn_ssh      = 'sas_grid_T'            ,       120         , 'sossheig',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
    351    sn_e3t      = 'sas_grid_T'            ,       120         , 'e3t_m'   ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
    352    sn_frq      = 'sas_grid_T'            ,       120         , 'frq_m'   ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     397   sn_usp      = 'sas_grid_U'            ,       120.        , 'uos'     ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     398   sn_vsp      = 'sas_grid_V'            ,       120.        , 'vos'     ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     399   sn_tem      = 'sas_grid_T'            ,       120.        , 'sosstsst',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     400   sn_sal      = 'sas_grid_T'            ,       120.        , 'sosaline',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     401   sn_ssh      = 'sas_grid_T'            ,       120.        , 'sossheig',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     402   sn_e3t      = 'sas_grid_T'            ,       120.        , 'e3t_m'   ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     403   sn_frq      = 'sas_grid_T'            ,       120.        , 'frq_m'   ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
    353404/ 
    354405!----------------------------------------------------------------------- 
    … …  
    373424   nn_chldta   =      0       !  RGB : Chl data (=1) or cst value (=0) 
    374425   rn_si1      =   23.0       !  2BD : longest depth of extinction 
    375     
     426 
    376427   cn_dir      = './'      !  root directory for the chlorophyl data location 
    377428   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
    378429   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
    379430   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
    380    sn_chl      ='chlorophyll'            ,        -1         , 'CHLA'    ,   .true.    , .true. , 'yearly'  , ''               , ''       , '' 
     431   sn_chl      ='chlorophyll'            ,        -1.        , 'CHLA'    ,   .true.    , .true. , 'yearly'  , ''               , ''       , '' 
    381432/ 
    382433!----------------------------------------------------------------------- 
    … …  
    390441      ln_sssr_bnd =  .true.   !  flag to bound erp term (associated with nn_sssr=2) 
    391442      rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day] 
     443      nn_sssr_ice =   1       ! control of sea surface restoring under sea-ice 
     444                              ! 0 = no restoration under ice : * (1-icefrac) 
     445                              ! 1 = restoration everywhere  
     446                              ! >1 = enhanced restoration under ice : 1+(nn_icedmp-1)*icefrac 
    392447 
    393448   cn_dir      = './'      !  root directory for the SST/SSS data location 
    … …  
    395450   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
    396451   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
    397    sn_sst      = 'sst_data'              ,        24         ,  'sst'    ,    .false.  , .false., 'yearly'  ,    ''            ,    ''    ,     '' 
    398    sn_sss      = 'sss_data'              ,        -1         ,  'sss'    ,    .true.   , .true. , 'yearly'  ,    ''            ,    ''    ,     '' 
     452   sn_sst      = 'sst_data'              ,        24.        ,  'sst'    ,    .false.  , .false., 'yearly'  ,    ''            ,    ''    ,     '' 
     453   sn_sss      = 'sss_data'              ,        -1.        ,  'sss'    ,    .true.   , .true. , 'yearly'  ,    ''            ,    ''    ,     '' 
    399454/ 
    400455!----------------------------------------------------------------------- 
    … …  
    417472   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
    418473   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
    419    sn_rnf      = 'runoff_core_monthly'   ,        -1         , 'sorunoff',   .true.    , .true. , 'yearly'  , ''               , ''       , '' 
    420    sn_cnf      = 'runoff_core_monthly'   ,         0         , 'socoefr0',   .false.   , .true. , 'yearly'  , ''               , ''       , '' 
    421    sn_s_rnf    = 'runoffs'               ,        24         , 'rosaline',   .true.    , .true. , 'yearly'  , ''               , ''       , '' 
    422    sn_t_rnf    = 'runoffs'               ,        24         , 'rotemper',   .true.    , .true. , 'yearly'  , ''               , ''       , '' 
    423    sn_dep_rnf  = 'runoffs'               ,         0         , 'rodepth' ,   .false.   , .true. , 'yearly'  , ''               , ''       , '' 
     474   sn_rnf      = 'runoff_core_monthly'   ,        -1.        , 'sorunoff',   .true.    , .true. , 'yearly'  , ''               , ''       , '' 
     475   sn_cnf      = 'runoff_core_monthly'   ,         0.        , 'socoefr0',   .false.   , .true. , 'yearly'  , ''               , ''       , '' 
     476   sn_s_rnf    = 'runoffs'               ,        24.        , 'rosaline',   .true.    , .true. , 'yearly'  , ''               , ''       , '' 
     477   sn_t_rnf    = 'runoffs'               ,        24.        , 'rotemper',   .true.    , .true. , 'yearly'  , ''               , ''       , '' 
     478   sn_dep_rnf  = 'runoffs'               ,         0.        , 'rodepth' ,   .false.   , .true. , 'yearly'  , ''               , ''       , '' 
    424479/ 
    425480!----------------------------------------------------------------------- 
    … …  
    434489   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
    435490   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
    436    sn_apr      = 'patm'                  ,         -1        ,'somslpre' ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,      '' 
    437 / 
    438 !----------------------------------------------------------------------- 
    439 &namsbc_isf    !  Top boundary layer (ISF)                              (ln_isfcav =T : read (ln_read_cfg=T)  
    440 !-----------------------------------------------------------------------             or set or usr_def_zgr ) 
    441    !                 ! type of top boundary layer  
    442    nn_isf      = 1         !  ice shelf melting/freezing 
    443                            !  1 = presence of ISF   ;  2 = bg03 parametrisation  
    444                            !  3 = rnf file for ISF  ;  4 = ISF specified freshwater flux 
    445                            !  options 1 and 4 need ln_isfcav = .true. (domzgr) 
    446       !              !  nn_isf = 1 or 2 cases: 
    447       rn_gammat0  = 1.e-4     ! gammat coefficient used in blk formula 
    448       rn_gammas0  = 1.e-4     ! gammas coefficient used in blk formula 
    449       !              !  nn_isf = 1 or 4 cases: 
    450       rn_hisf_tbl =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
    451       !                       ! 0 => thickness of the tbl = thickness of the first wet cell 
    452       !              ! nn_isf = 1 case 
    453       nn_isfblk   = 1         ! 1 ISOMIP  like: 2 equations formulation (Hunter et al., 2006) 
    454       !                       ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 
    455       nn_gammablk = 1         ! 0 = cst Gammat (= gammat/s) 
    456       !                       ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
    457       !                       ! 2 = velocity and stability dependent Gamma    (Holland et al. 1999) 
    458  
    459    !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
    460    !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    461    !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
    462 !* nn_isf = 4 case 
    463    sn_fwfisf   = 'rnfisf'    ,         -12       ,'sowflisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    464 !* nn_isf = 3 case 
    465    sn_rnfisf   = 'rnfisf'    ,         -12       ,'sofwfisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    466 !* nn_isf = 2 and 3 cases  
    467    sn_depmax_isf ='rnfisf'    ,         -12       ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    468    sn_depmin_isf ='rnfisf'    ,         -12       ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    469 !* nn_isf = 2 case 
    470    sn_Leff_isf = 'rnfisf'    ,         -12       ,'Leff'     ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    471 / 
    472 !----------------------------------------------------------------------- 
    473 &namsbc_iscpl  !   land ice / ocean coupling option                     (ln_isfcav =T : read (ln_read_cfg=T)  
    474 !-----------------------------------------------------------------------             or set or usr_def_zgr ) 
    475    nn_drown    = 10        ! number of iteration of the extrapolation loop (fill the new wet cells) 
    476    ln_hsb      = .false.   ! activate conservation module (conservation exact after a time of rn_fiscpl) 
    477    nn_fiscpl   = 43800     ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency) 
     491   sn_apr      = 'patm'                  ,         -1.       ,'somslpre' ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,      '' 
     492/ 
     493!----------------------------------------------------------------------- 
     494&namisf       !  Top boundary layer (ISF)                               (default: OFF) 
     495!----------------------------------------------------------------------- 
     496   ! 
     497   ! ---------------- ice shelf load ------------------------------- 
     498   ! 
     499   cn_isfload = 'uniform'      ! scheme to compute ice shelf load (ln_isfcav = .true. in domain_cfg.nc) 
     500      rn_isfload_T = -1.9 
     501      rn_isfload_S = 34.4 
     502   ! 
     503   ! ---------------- ice shelf melt formulation ------------------------------- 
     504   ! 
     505   ln_isf = .false.           ! activate ice shelf module 
     506      ln_isfdebug = .false.      ! add debug print in ISF code (global min/max/sum of specific variable) 
     507      cn_isfdir   = './'         ! directory for all ice shelf input file 
     508      ! 
     509      ! ---------------- cavities opened ------------------------------- 
     510      ! 
     511      ln_isfcav_mlt = .false.    ! ice shelf melting into the cavity (need ln_isfcav = .true. in domain_cfg.nc) 
     512         cn_isfcav_mlt = '3eq'   ! ice shelf melting formulation (spe/2eq/3eq/oasis) 
     513         !                       ! spe = fwfisf is read from a forcing field 
     514         !                       ! 2eq = ISOMIP  like: 2 equations formulation (Hunter et al., 2006 for a short description) 
     515         !                       ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2016 for a short description) 
     516         !                       ! oasis = fwfisf is given by oasis and pattern by file sn_isfcav_fwf 
     517         !              !  cn_isfcav_mlt = 2eq or 3eq cases: 
     518         cn_gammablk = 'vel'     ! scheme to compute gammat/s (spe,ad15,hj99) 
     519         !                       ! spe      = constant transfert velocity (rn_gammat0, rn_gammas0) 
     520         !                       ! vel      = velocity dependent transfert velocity (u* * gammat/s) (Asay-Davis et al. 2016 for a short description) 
     521         !                       ! vel_stab = velocity and stability dependent transfert coeficient (Holland et al. 1999 for a complete description) 
     522         rn_gammat0  = 1.4e-2    ! gammat coefficient used in spe, vel and vel_stab gamma computation method 
     523         rn_gammas0  = 4.0e-4    ! gammas coefficient used in spe, vel and vel_stab gamma computation method 
     524         ! 
     525         rn_htbl     =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
     526         !                       ! 0 => thickness of the tbl = thickness of the first wet cell 
     527         ! 
     528         !* 'spe' and 'oasis' case 
     529         !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
     530         !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     531         !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     532         sn_isfcav_fwf = 'isfmlt_cav',      -12.      , 'fwflisf'  ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     533      ! 
     534      ! ---------------- cavities parametrised ------------------------------- 
     535      ! 
     536      ln_isfpar_mlt = .false.   ! ice shelf melting parametrised 
     537         cn_isfpar_mlt = 'spe'  ! ice shelf melting parametrisation (spe/bg03/oasis) 
     538         !                      ! spe   = fwfisf is read from a forcing field 
     539         !                      ! bg03  = melt computed using Beckmann and Goosse parametrisation 
     540         !                      ! oasis = fwfisf is given by oasis and pattern by file sn_isfpar_fwf 
     541         ! 
     542         !* all cases 
     543         !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
     544         !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     545         !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     546         sn_isfpar_zmax = 'isfmlt_par',       0.       ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     547         sn_isfpar_zmin = 'isfmlt_par',       0.       ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     548         !* 'spe' and 'oasis' case 
     549         sn_isfpar_fwf = 'isfmlt_par' ,      -12.      ,'sofwfisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     550         !* 'bg03' case 
     551         sn_isfpar_Leff = 'isfmlt_par',       0.       ,'Leff'     ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     552      ! 
     553      ! ---------------- ice sheet coupling ------------------------------- 
     554      ! 
     555      ln_isfcpl = .false. 
     556         nn_drown       = 10        ! number of iteration of the extrapolation loop (fill the new wet cells) 
     557         ln_isfcpl_cons = .false. 
    478558/ 
    479559!----------------------------------------------------------------------- 
    … …  
    484564   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
    485565   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
    486    sn_cdg      =  'sdw_ecwaves_orca2'    ,        6          , 'drag_coeff' ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
    487    sn_usd      =  'sdw_ecwaves_orca2'    ,        6          , 'u_sd2d'     ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
    488    sn_vsd      =  'sdw_ecwaves_orca2'    ,        6          , 'v_sd2d'     ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
    489    sn_hsw      =  'sdw_ecwaves_orca2'    ,        6          , 'hs'         ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
    490    sn_wmp      =  'sdw_ecwaves_orca2'    ,        6          , 'wmp'        ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
    491    sn_wfr      =  'sdw_ecwaves_orca2'    ,        6          , 'wfr'        ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
    492    sn_wnum     =  'sdw_ecwaves_orca2'    ,        6          , 'wave_num'   ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
    493    sn_tauwoc   =  'sdw_ecwaves_orca2'    ,        6          , 'wave_stress',  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
    494    sn_tauwx    =  'sdw_ecwaves_orca2'    ,        6          , 'wave_stress',  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
    495    sn_tauwy    =  'sdw_ecwaves_orca2'    ,        6          , 'wave_stress',  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
     566   sn_cdg      =  'sdw_ecwaves_orca2'    ,        6.         , 'drag_coeff' ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
     567   sn_usd      =  'sdw_ecwaves_orca2'    ,        6.         , 'u_sd2d'     ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
     568   sn_vsd      =  'sdw_ecwaves_orca2'    ,        6.         , 'v_sd2d'     ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
     569   sn_hsw      =  'sdw_ecwaves_orca2'    ,        6.         , 'hs'         ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
     570   sn_wmp      =  'sdw_ecwaves_orca2'    ,        6.         , 'wmp'        ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
     571   sn_wfr      =  'sdw_ecwaves_orca2'    ,        6.         , 'wfr'        ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
     572   sn_wnum     =  'sdw_ecwaves_orca2'    ,        6.         , 'wave_num'   ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
     573   sn_tauwoc   =  'sdw_ecwaves_orca2'    ,        6.         , 'wave_stress',  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
     574   sn_tauwx    =  'sdw_ecwaves_orca2'    ,        6.         , 'wave_stress',  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
     575   sn_tauwy    =  'sdw_ecwaves_orca2'    ,        6.         , 'wave_stress',  .true.  , .true. , 'yearly'  ,  ''              , ''       , '' 
    496576/ 
    497577!----------------------------------------------------------------------- 
    … …  
    499579!----------------------------------------------------------------------- 
    500580   ln_icebergs = .false.      ! activate iceberg floats (force =F with "key_agrif") 
     581   ! 
     582   !                          ! restart 
     583   cn_icbrst_in     = "restart_icb" !  suffix of iceberg restart name (input) 
     584   cn_icbrst_indir  = "./"          !  directory from which to read input ocean restarts 
     585   cn_icbrst_out    = "restart_icb" !  suffix of ocean restart name (output) 
     586   cn_icbrst_outdir = "./"          !  directory from which to read output ocean restarts           
    501587   ! 
    502588   !                          ! diagnostics: 
    … …  
    533619   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
    534620   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
    535    sn_icb     =  'calving'              ,         -1         ,'calvingmask',  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     621   sn_icb     =  'calving'              ,         -1.        ,'calvingmask',  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
    536622/ 
    537623 
    … …  
    557643&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    558644!----------------------------------------------------------------------- 
    559    ln_spc_dyn    = .true.  !  use 0 as special value for dynamics 
    560    rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [m2/s] 
    561    rn_sponge_dyn = 2880.   !  coefficient for dynamics sponge layer [m2/s] 
    562    ln_chk_bathy  = .false. !  =T  check the parent bathymetry 
     645   ln_agrif_2way   = .true.  !  activate two way nesting 
     646   ln_init_chfrpar = .false. !  initialize child grids from parent 
     647   ln_spc_dyn      = .true.  !  use 0 as special value for dynamics 
     648   rn_sponge_tra   = 0.002   !  coefficient for tracer   sponge layer [] 
     649   rn_sponge_dyn   = 0.002   !  coefficient for dynamics sponge layer [] 
     650   rn_trelax_tra   = 0.01    !  inverse of relaxation time (in steps) for tracers [] 
     651   rn_trelax_dyn   = 0.01    !  inverse of relaxation time (in steps) for dynamics [] 
     652   ln_chk_bathy    = .false. !  =T  check the parent bathymetry 
    563653/ 
    564654!----------------------------------------------------------------------- 
    … …  
    566656!----------------------------------------------------------------------- 
    567657   ln_tide     = .false.      ! Activate tides 
    568       ln_tide_pot   = .true.                !  use tidal potential forcing 
     658      nn_tide_var   = 1          !  Variant of tidal parameter set and tide-potential computation 
     659      !                          !     (1: default; 0: compatibility with previous versions) 
     660      ln_tide_dia   = .false.    !  Enable tidal diagnostic output 
     661      ln_tide_pot   = .false.               !  use tidal potential forcing 
     662         rn_tide_gamma = 0.7                   ! Tidal tilt factor 
    569663         ln_scal_load  = .false.               ! Use scalar approximation for 
    570664            rn_scal_load = 0.094               !     load potential 
    571665         ln_read_load  = .false.               ! Or read load potential from file 
    572666            cn_tide_load = 'tide_LOAD_grid_T.nc'  ! filename for load potential 
    573             !       
     667            ! 
    574668      ln_tide_ramp  = .false.               !  Use linear ramp for tides at startup 
    575          rdttideramp   =    0.                 !  ramp duration in days 
    576       clname(1)    = 'DUMMY'               !  name of constituent - all tidal components must be set in namelist_cfg 
     669         rn_tide_ramp_dt = 0.               !  ramp duration in days 
     670      sn_tide_cnames(1) = 'DUMMY'               !  name of constituent - all tidal components must be set in namelist_cfg 
    577671/ 
    578672!----------------------------------------------------------------------- 
    … …  
    599693   nn_ice_dta    =  0         !  = 0, bdy data are equal to the initial state 
    600694   !                          !  = 1, bdy data are read in 'bdydata   .nc' files 
    601    rn_ice_tem    = 270.       !  si3 only: arbitrary temperature of incoming sea ice 
    602    rn_ice_sal    = 10.        !  si3 only:      --   salinity           -- 
    603    rn_ice_age    = 30.        !  si3 only:      --   age                -- 
    604695   ! 
    605696   ln_tra_dmp    =.false.     !  open boudaries conditions for tracers 
    606697   ln_dyn3d_dmp  =.false.     !  open boundary condition for baroclinic velocities 
    607698   rn_time_dmp   =  1.        !  Damping time scale in days 
    608    rn_time_dmp_out = 1.        !  Outflow damping time scale 
     699   rn_time_dmp_out = 1.       !  Outflow damping time scale 
    609700   nn_rimwidth   = 10         !  width of the relaxation zone 
    610701   ln_vol        = .false.    !  total volume correction (see nn_volctl parameter) 
    611702   nn_volctl     =  1         !  = 0, the total water flux across open boundaries is zero 
    612    nb_jpk_bdy    = -1         ! number of levels in the bdy data (set < 0 if consistent with planned run) 
    613703/ 
    614704!----------------------------------------------------------------------- 
    615705&nambdy_dta    !  open boundaries - external data                       (see nam_bdy) 
    616706!----------------------------------------------------------------------- 
    617    ln_full_vel = .false.      !  ??? 
    618  
     707   ln_zinterp  = .false.      !  T if a vertical interpolation is required. Variables gdep[tuv] and e3[tuv] must exist in the file 
     708   !                          !  automatically defined to T if the number of vertical levels in bdy dta /= jpk 
     709   ln_full_vel = .false.      !  T if [uv]3d are "full" velocities and not only its baroclinic components 
     710   !                          !  in this case, baroclinic and barotropic velocities will be recomputed -> [uv]2d not needed 
     711   ! 
    619712   cn_dir      = 'bdydta/'    !  root directory for the BDY data location 
    620713   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
    621714   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
    622715   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
    623    bn_ssh      = 'amm12_bdyT_u2d'        ,         24        , 'sossheig',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
    624    bn_u2d      = 'amm12_bdyU_u2d'        ,         24        , 'vobtcrtx',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
    625    bn_v2d      = 'amm12_bdyV_u2d'        ,         24        , 'vobtcrty',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
    626    bn_u3d      = 'amm12_bdyU_u3d'        ,         24        , 'vozocrtx',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
    627    bn_v3d      = 'amm12_bdyV_u3d'        ,         24        , 'vomecrty',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
    628    bn_tem      = 'amm12_bdyT_tra'        ,         24        , 'votemper',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
    629    bn_sal      = 'amm12_bdyT_tra'        ,         24        , 'vosaline',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     716   bn_ssh      = 'amm12_bdyT_u2d'        ,         24.       , 'sossheig',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     717   bn_u2d      = 'amm12_bdyU_u2d'        ,         24.       , 'vobtcrtx',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     718   bn_v2d      = 'amm12_bdyV_u2d'        ,         24.       , 'vobtcrty',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     719   bn_u3d      = 'amm12_bdyU_u3d'        ,         24.       , 'vozocrtx',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     720   bn_v3d      = 'amm12_bdyV_u3d'        ,         24.       , 'vomecrty',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     721   bn_tem      = 'amm12_bdyT_tra'        ,         24.       , 'votemper',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     722   bn_sal      = 'amm12_bdyT_tra'        ,         24.       , 'vosaline',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
    630723!* for si3 
    631 !   bn_a_i     = 'amm12_bdyT_ice'        ,         24        , 'ileadfra',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
    632 !   bn_h_i     = 'amm12_bdyT_ice'        ,         24        , 'iicethic',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
    633 !   bn_h_s     = 'amm12_bdyT_ice'        ,         24        , 'isnowthi',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     724   bn_a_i      = 'amm12_bdyT_ice'        ,         24.       , 'siconc'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     725   bn_h_i      = 'amm12_bdyT_ice'        ,         24.       , 'sithic'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     726   bn_h_s      = 'amm12_bdyT_ice'        ,         24.       , 'snthic'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     727   bn_t_i      = 'NOT USED'              ,         24.       , 'sitemp'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     728   bn_t_s      = 'NOT USED'              ,         24.       , 'sntemp'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     729   bn_tsu      = 'NOT USED'              ,         24.       , 'sittop'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     730   bn_s_i      = 'NOT USED'              ,         24.       , 'sisalt'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     731   ! melt ponds (be careful, bn_aip is the pond concentration (not fraction), so it differs from rn_iceapnd) 
     732   bn_aip      = 'NOT USED'              ,         24.       , 'siapnd'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     733   bn_hip      = 'NOT USED'              ,         24.       , 'sihpnd'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     734   ! if bn_t_i etc are "not used", then define arbitrary temperatures and salinity and ponds 
     735   rn_ice_tem  = 270.         !  arbitrary temperature               of incoming sea ice 
     736   rn_ice_sal  = 10.          !       --   salinity                            -- 
     737   rn_ice_age  = 30.          !       --   age                                 -- 
     738   rn_ice_apnd = 0.2          !       --   pond fraction = a_ip/a_i            -- 
     739   rn_ice_hpnd = 0.05         !       --   pond depth                          -- 
    634740/ 
    635741!----------------------------------------------------------------------- 
    … …  
    638744   filtide          = 'bdydta/amm12_bdytide_'   !  file name root of tidal forcing files 
    639745   ln_bdytide_2ddta = .false.                   ! 
    640    ln_bdytide_conj  = .false.                   !  
    641746/ 
    642747 
    … …  
    645750!!                                                                    !! 
    646751!!   namdrg        top/bottom drag coefficient                          (default: NO selection) 
    647 !!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T) 
    648 !!   namdrg_bot    bottom friction                                      (ln_OFF=F) 
     752!!   namdrg_top    top    friction                                      (ln_drg_OFF=F & ln_isfcav=T) 
     753!!   namdrg_bot    bottom friction                                      (ln_drg_OFF=F) 
    649754!!   nambbc        bottom temperature boundary condition                (default: OFF) 
    650755!!   nambbl        bottom boundary layer scheme                         (default: OFF) 
    … …  
    654759&namdrg        !   top/bottom drag coefficient                          (default: NO selection) 
    655760!----------------------------------------------------------------------- 
    656    ln_OFF      = .false.   !  free-slip       : Cd = 0                  (F => fill namdrg_bot 
     761   ln_drg_OFF  = .false.   !  free-slip       : Cd = 0                  (F => fill namdrg_bot 
    657762   ln_lin      = .false.   !      linear  drag: Cd = Cd0 Uc0                   &   namdrg_top) 
    658763   ln_non_lin  = .false.   !  non-linear  drag: Cd = Cd0 |U| 
    … …  
    662767/ 
    663768!----------------------------------------------------------------------- 
    664 &namdrg_top    !   TOP friction                                         (ln_OFF =F & ln_isfcav=T) 
     769&namdrg_top    !   TOP friction                                         (ln_drg_OFF =F & ln_isfcav=T) 
    665770!----------------------------------------------------------------------- 
    666771   rn_Cd0      =  1.e-3    !  drag coefficient [-] 
    667    rn_Uc0      =  0.4      !  ref. velocity [m/s] (linear drag=Cd0*Uc0)  
     772   rn_Uc0      =  0.4      !  ref. velocity [m/s] (linear drag=Cd0*Uc0) 
    668773   rn_Cdmax    =  0.1      !  drag value maximum [-] (logarithmic drag) 
    669774   rn_ke0      =  2.5e-3   !  background kinetic energy  [m2/s2] (non-linear cases) 
    … …  
    673778/ 
    674779!----------------------------------------------------------------------- 
    675 &namdrg_bot    !   BOTTOM friction                                      (ln_OFF =F) 
     780&namdrg_bot    !   BOTTOM friction                                      (ln_drg_OFF =F) 
    676781!----------------------------------------------------------------------- 
    677782   rn_Cd0      =  1.e-3    !  drag coefficient [-] 
    678    rn_Uc0      =  0.4      !  ref. velocity [m/s] (linear drag=Cd0*Uc0)  
     783   rn_Uc0      =  0.4      !  ref. velocity [m/s] (linear drag=Cd0*Uc0) 
    679784   rn_Cdmax    =  0.1      !  drag value maximum [-] (logarithmic drag) 
    680785   rn_ke0      =  2.5e-3   !  background kinetic energy  [m2/s2] (non-linear cases) 
    … …  
    743848      nn_cen_v   =  4            !  =2/4, vertical   2nd order CEN / 4th order COMPACT 
    744849   ln_traadv_fct = .false. !  FCT scheme 
    745       nn_fct_h   =  2            !  =2/4, horizontal 2nd / 4th order  
    746       nn_fct_v   =  2            !  =2/4, vertical   2nd / COMPACT 4th order  
     850      nn_fct_h   =  2            !  =2/4, horizontal 2nd / 4th order 
     851      nn_fct_v   =  2            !  =2/4, vertical   2nd / COMPACT 4th order 
    747852   ln_traadv_mus = .false. !  MUSCL scheme 
    748853      ln_mus_ups = .false.       !  use upstream scheme near river mouths 
    … …  
    765870   ln_traldf_triad = .false.   !  iso-neutral (triad    operator) 
    766871   ! 
    767    !                       !  iso-neutral options:         
     872   !                       !  iso-neutral options: 
    768873   ln_traldf_msc   = .false.   !  Method of Stabilizing Correction      (both operators) 
    769874   rn_slpmax       =  0.01     !  slope limit                           (both operators) 
    … …  
    775880   nn_aht_ijk_t    = 0         !  space/time variation of eddy coefficient: 
    776881      !                             !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
    777       !                             !   =  0           constant  
    778       !                             !   = 10 F(k)      =ldf_c1d  
    779       !                             !   = 20 F(i,j)    =ldf_c2d  
     882      !                             !   =  0           constant 
     883      !                             !   = 10 F(k)      =ldf_c1d 
     884      !                             !   = 20 F(i,j)    =ldf_c2d 
    780885      !                             !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    781886      !                             !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
    782887      !                             !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 
    783       !                        !  time invariant coefficients:  aht0 = 1/2  Ud*Ld   (lap case)  
     888      !                        !  time invariant coefficients:  aht0 = 1/2  Ud*Ld   (lap case) 
    784889      !                             !                           or   = 1/12 Ud*Ld^3 (blp case) 
    785890      rn_Ud        = 0.01           !  lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) 
    … …  
    807912      nn_aei_ijk_t    = 0           !  space/time variation of eddy coefficient: 
    808913      !                             !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
    809       !                             !   =  0           constant  
    810       !                             !   = 10 F(k)      =ldf_c1d  
    811       !                             !   = 20 F(i,j)    =ldf_c2d  
     914      !                             !   =  0           constant 
     915      !                             !   = 10 F(k)      =ldf_c1d 
     916      !                             !   = 20 F(i,j)    =ldf_c2d 
    812917      !                             !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    813918      !                             !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
    814       !                        !  time invariant coefficients:  aei0 = 1/2  Ue*Le  
     919      !                        !  time invariant coefficients:  aei0 = 1/2  Ue*Le 
    815920      rn_Ue        = 0.02           !  lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) 
    816921      rn_Le        = 200.e+3        !  lateral diffusive length   [m]   (nn_aht_ijk_t= 0, 10) 
    … …  
    852957   rn_lf_cutoff  =  5.0             !  cutoff frequency for low-pass filter  [days] 
    853958   rn_zdef_max   =  0.9             !  maximum fractional e3t deformation 
    854    ln_vvl_dbg    = .true.           !  debug prints    (T/F) 
     959   ln_vvl_dbg    = .false.          !  debug prints    (T/F) 
    855960/ 
    856961!----------------------------------------------------------------------- 
    … …  
    872977   ln_dynvor_eeT = .false. !  energy conserving scheme (een using e3t) 
    873978   ln_dynvor_een = .false. !  energy & enstrophy scheme 
    874       nn_een_e3f = 0          ! =0  e3f = mi(mj(e3t))/4  
     979      nn_een_e3f = 0          ! =0  e3f = mi(mj(e3t))/4 
    875980      !                       ! =1  e3f = mi(mj(e3t))/mi(mj( tmask)) 
    876981   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T)        ==>>> PLEASE DO NOT ACTIVATE 
    … …  
    8951000      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    8961001         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    897          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    898          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     1002         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     1003         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    8991004      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    9001005         rn_bt_cmax   =  0.8        ! =T : the Maximum Courant Number allowed 
    901          nn_baro      = 30          ! =F : the number of sub-step in rn_rdt seconds 
     1006         nn_e         = 30          ! =F : the number of sub-step in rn_Dt seconds 
    9021007      rn_bt_alpha   = 0.         ! Temporal diffusion parameter (if ln_bt_av=F) 
    9031008/ 
    … …  
    9171022      !                             !  =-30  read in eddy_viscosity_3D.nc file 
    9181023      !                             !  =-20  read in eddy_viscosity_2D.nc file 
    919       !                             !  =  0  constant  
     1024      !                             !  =  0  constant 
    9201025      !                             !  = 10  F(k)=c1d 
    9211026      !                             !  = 20  F(i,j)=F(grid spacing)=c2d 
    … …  
    9231028      !                             !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
    9241029      !                             !  = 32  F(i,j,k)=F(local gridscale and deformation rate) 
    925       !                        !  time invariant coefficients :  ahm = 1/2  Uv*Lv   (lap case)  
     1030      !                        !  time invariant coefficients :  ahm = 1/2  Uv*Lv   (lap case) 
    9261031      !                             !                            or  = 1/12 Uv*Lv^3 (blp case) 
    9271032      rn_Uv      = 0.1              !  lateral viscous velocity [m/s] (nn_ahm_ijk_t= 0, 10, 20, 30) 
    … …  
    9451050   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
    9461051   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
    947    sn_tem      = 'dyna_grid_T'           ,       120         , 'votemper'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
    948    sn_sal      = 'dyna_grid_T'           ,       120         , 'vosaline'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
    949    sn_mld      = 'dyna_grid_T'           ,       120         , 'somixhgt'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
    950    sn_emp      = 'dyna_grid_T'           ,       120         , 'sowaflup'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
    951    sn_fmf      = 'dyna_grid_T'           ,       120         , 'iowaflup'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
    952    sn_ice      = 'dyna_grid_T'           ,       120         , 'soicecov'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
    953    sn_qsr      = 'dyna_grid_T'           ,       120         , 'soshfldo'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
    954    sn_wnd      = 'dyna_grid_T'           ,       120         , 'sowindsp'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
    955    sn_uwd      = 'dyna_grid_U'           ,       120         , 'uocetr_eff',  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
    956    sn_vwd      = 'dyna_grid_V'           ,       120         , 'vocetr_eff',  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
    957    sn_wwd      = 'dyna_grid_W'           ,       120         , 'wocetr_eff',  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
    958    sn_avt      = 'dyna_grid_W'           ,       120         , 'voddmavs'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
    959    sn_ubl      = 'dyna_grid_U'           ,       120         , 'sobblcox'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
    960    sn_vbl      = 'dyna_grid_V'           ,       120         , 'sobblcoy'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     1052   sn_tem      = 'dyna_grid_T'           ,       120.        , 'votemper'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     1053   sn_sal      = 'dyna_grid_T'           ,       120.        , 'vosaline'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     1054   sn_mld      = 'dyna_grid_T'           ,       120.        , 'somixhgt'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     1055   sn_emp      = 'dyna_grid_T'           ,       120.        , 'sowaflup'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     1056   sn_fmf      = 'dyna_grid_T'           ,       120.        , 'iowaflup'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     1057   sn_ice      = 'dyna_grid_T'           ,       120.        , 'soicecov'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     1058   sn_qsr      = 'dyna_grid_T'           ,       120.        , 'soshfldo'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     1059   sn_wnd      = 'dyna_grid_T'           ,       120.        , 'sowindsp'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     1060   sn_uwd      = 'dyna_grid_U'           ,       120.        , 'uocetr_eff',  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     1061   sn_vwd      = 'dyna_grid_V'           ,       120.        , 'vocetr_eff',  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     1062   sn_wwd      = 'dyna_grid_W'           ,       120.        , 'wocetr_eff',  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     1063   sn_avt      = 'dyna_grid_W'           ,       120.        , 'voddmavs'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     1064   sn_ubl      = 'dyna_grid_U'           ,       120.        , 'sobblcox'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     1065   sn_vbl      = 'dyna_grid_V'           ,       120.        , 'sobblcoy'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
    9611066/ 
    9621067 
    … …  
    10351140   !                       !                 = 3 as =2 with distinct dissipative an mixing length scale 
    10361141   ln_mxl0     = .true.    !  surface mixing length scale = F(wind stress) (T) or not (F) 
     1142      nn_mxlice    = 0        ! type of scaling under sea-ice 
     1143                              !    = 0 no scaling under sea-ice 
     1144                              !    = 1 scaling with constant sea-ice thickness 
     1145                              !    = 2  scaling with mean sea-ice thickness ( only with SI3 sea-ice model ) 
     1146                              !    = 3  scaling with maximum sea-ice thickness 
     1147      rn_mxlice   = 10.       ! max constant ice thickness value when scaling under sea-ice ( nn_mxlice=1) 
    10371148   rn_mxl0     =   0.04    !  surface  buoyancy lenght scale minimum value 
    1038    ln_drg      = .false.   !  top/bottom friction added as boundary condition of TKE 
    10391149   ln_lc       = .true.    !  Langmuir cell parameterisation (Axell 2002) 
    10401150      rn_lc       =   0.15    !  coef. associated to Langmuir cells 
    … …  
    10471157                              !        = 0  constant 10 m length scale 
    10481158                              !        = 1  0.5m at the equator to 30m poleward of 40 degrees 
    1049       rn_eice     =   4       !  below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4    
     1159      rn_eice     =   4       !  below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4 
    10501160/ 
    10511161!----------------------------------------------------------------------- 
    … …  
    10931203   ln_mevar    = .true.    !  variable (T) or constant (F) mixing efficiency 
    10941204   ln_tsdiff   = .true.    !  account for differential T/S mixing (T) or not (F) 
    1095 / 
    1096  
     1205 
     1206   cn_dir      = './'      !  root directory for the iwm data location                                                                            
     1207   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     1208   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     1209   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     1210   sn_mpb      = 'NOT USED'              , -12.              , 'mixing_power_bot' , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1211   sn_mpp      = 'NOT USED'              , -12.              , 'mixing_power_pyc' , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1212   sn_mpc      = 'NOT USED'              , -12.              , 'mixing_power_cri' , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1213   sn_dsb      = 'NOT USED'              , -12.              , 'decay_scale_bot'  , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1214   sn_dsc      = 'NOT USED'              , -12.              , 'decay_scale_cri'  , .false.  , .true. , 'yearly' , '' , ''  , '' 
     1215/ 
    10971216!!====================================================================== 
    10981217!!                  ***  Diagnostics namelists  ***                   !! 
    … …  
    11031222!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    11041223!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    1105 !!   namflo       float parameters                                      ("key_float") 
    1106 !!   nam_diaharm  Harmonic analysis of tidal constituents               ("key_diaharm") 
    1107 !!   namdct       transports through some sections                      ("key_diadct") 
    1108 !!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
     1224!!   namflo       float parameters                                      (default: OFF) 
     1225!!   nam_diadct   transports through some sections                      (default: OFF) 
    11091226!!   nam_dia25h   25h Mean Output                                       (default: OFF) 
    11101227!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4") 
    … …  
    11311248!!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S 
    11321249!!gm 
    1133 !----------------------------------------------------------------------- 
    1134 &namptr        !   Poleward Transport Diagnostic                        (default: OFF) 
    1135 !----------------------------------------------------------------------- 
    1136    ln_diaptr   = .false.   !  Poleward heat and salt transport (T) or not (F) 
    1137    ln_subbas   = .false.   !  Atlantic/Pacific/Indian basins computation (T) or not 
    11381250/ 
    11391251!----------------------------------------------------------------------- 
    11401252&namhsb        !  Heat and salt budgets                                 (default: OFF) 
    11411253!----------------------------------------------------------------------- 
    1142    ln_diahsb   = .false.   !  check the heat and salt budgets (T) or not (F) 
     1254   ln_diahsb   = .false.   !  output the heat and salt budgets (T) or not (F) 
    11431255/ 
    11441256!----------------------------------------------------------------------- 
    … …  
    11491261/ 
    11501262!----------------------------------------------------------------------- 
    1151 &namflo        !   float parameters                                     ("key_float") 
    1152 !----------------------------------------------------------------------- 
    1153    jpnfl       = 1         !  total number of floats during the run 
    1154    jpnnewflo   = 0         !  number of floats for the restart 
    1155    ln_rstflo   = .false.   !  float restart (T) or not (F) 
    1156    nn_writefl  =      75   !  frequency of writing in float output file 
    1157    nn_stockfl  =    5475   !  frequency of creation of the float restart file 
    1158    ln_argo     = .false.   !  Argo type floats (stay at the surface each 10 days) 
    1159    ln_flork4   = .false.   !  trajectories computed with a 4th order Runge-Kutta (T) 
    1160    !                       !  or computed with Blanke' scheme (F) 
    1161    ln_ariane   = .true.    !  Input with Ariane tool convention(T) 
    1162    ln_flo_ascii = .true.    !  Output with Ariane tool netcdf convention(F) or ascii file (T) 
    1163 / 
    1164 !----------------------------------------------------------------------- 
    1165 &nam_diaharm   !   Harmonic analysis of tidal constituents              ("key_diaharm") 
    1166 !----------------------------------------------------------------------- 
    1167     nit000_han = 1         ! First time step used for harmonic analysis 
    1168     nitend_han = 75        ! Last time step used for harmonic analysis 
    1169     nstep_han  = 15        ! Time step frequency for harmonic analysis 
    1170     tname(1)   = 'M2'      ! Name of tidal constituents 
    1171     tname(2)   = 'K1' 
    1172 / 
    1173 !----------------------------------------------------------------------- 
    1174 &namdct        ! transports through some sections                       ("key_diadct") 
    1175 !----------------------------------------------------------------------- 
    1176     nn_dct     = 15        !  time step frequency for transports computing 
    1177     nn_dctwri  = 15        !  time step frequency for transports writing 
    1178     nn_secdebug = 112       !      0 : no section to debug 
    1179     !                      !     -1 : debug all section 
    1180     !                      !  0 < n : debug section number n 
    1181 / 
    1182 !----------------------------------------------------------------------- 
    1183 &nam_diatmb    !  Top Middle Bottom Output                              (default: OFF) 
    1184 !----------------------------------------------------------------------- 
    1185    ln_diatmb   = .false.   !  Choose Top Middle and Bottom output or not 
     1263&namflo        !   float parameters                                     (default: OFF) 
     1264!----------------------------------------------------------------------- 
     1265   ln_floats   = .false.      ! activate floats or not 
     1266      jpnfl       = 1         !    total number of floats during the run 
     1267      jpnnewflo   = 0         !    number of floats for the restart 
     1268      ln_rstflo   = .false.   !    float restart (T) or not (F) 
     1269      nn_writefl  =      75   !    frequency of writing in float output file 
     1270      nn_stockfl  =    5475   !    frequency of creation of the float restart file 
     1271      ln_argo     = .false.   !    Argo type floats (stay at the surface each 10 days) 
     1272      ln_flork4   = .false.   !    trajectories computed with a 4th order Runge-Kutta (T) 
     1273      !                       !    or computed with Blanke' scheme (F) 
     1274      ln_ariane   = .true.    !    Input with Ariane tool convention(T) 
     1275      ln_flo_ascii= .true.    !    Output with Ariane tool netcdf convention(F) or ascii file (T) 
     1276/ 
     1277!----------------------------------------------------------------------- 
     1278&nam_diadct    !   transports through some sections                     (default: OFF) 
     1279!----------------------------------------------------------------------- 
     1280    ln_diadct  = .false.   ! Calculate transport thru sections or not 
     1281       nn_dct     = 15     !  time step frequency for transports computing 
     1282       nn_dctwri  = 15     !  time step frequency for transports writing 
     1283       nn_secdebug = 112   !      0 : no section to debug 
     1284       !                   !     -1 : debug all section 
     1285       !                   !  0 < n : debug section number n 
    11861286/ 
    11871287!----------------------------------------------------------------------- 
    … …  
    12961396&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi") 
    12971397!----------------------------------------------------------------------- 
    1298    cn_mpi_send =  'I'      !  mpi send/recieve type   ='S', 'B', or 'I' for standard send, 
    1299    !                       !  buffer blocking send or immediate non-blocking sends, resp. 
    1300    nn_buffer   =   0       !  size in bytes of exported buffer ('B' case), 0 no exportation 
     1398   ln_listonly =  .false.  !  do nothing else than listing the best domain decompositions (with land domains suppression) 
     1399   !                       !  if T: the largest number of cores tested is defined by max(mppsize, jpni*jpnj) 
    13011400   ln_nnogather =  .true.  !  activate code to avoid mpi_allgather use at the northfold 
    1302    jpni        =   0       !  jpni   number of processors following i (set automatically if < 1) 
    1303    jpnj        =   0       !  jpnj   number of processors following j (set automatically if < 1) 
     1401   jpni        =   0       !  number of processors following i (set automatically if < 1), see also ln_listonly = T 
     1402   jpnj        =   0       !  number of processors following j (set automatically if < 1), see also ln_listonly = T 
     1403   nn_hls      =   1       !  halo width (applies to both rows and columns) 
    13041404/ 
    13051405!----------------------------------------------------------------------- 
    13061406&namctl        !   Control prints                                       (default: OFF) 
    13071407!----------------------------------------------------------------------- 
    1308    ln_ctl = .FALSE.                 ! Toggle all report printing on/off (T/F); Ignored if sn_cfctl%l_config is T 
    1309      sn_cfctl%l_config = .TRUE.     ! IF .true. then control which reports are written with the following 
    1310        sn_cfctl%l_runstat = .FALSE. ! switches and which areas produce reports with the proc integer settings. 
    1311        sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure 
    1312        sn_cfctl%l_oceout  = .FALSE. ! that  all areas report. 
    1313        sn_cfctl%l_layout  = .FALSE. ! 
    1314        sn_cfctl%l_mppout  = .FALSE. ! 
    1315        sn_cfctl%l_mpptop  = .FALSE. ! 
    1316        sn_cfctl%procmin   = 0       ! Minimum area number for reporting [default:0] 
    1317        sn_cfctl%procmax   = 1000000 ! Maximum area number for reporting [default:1000000] 
    1318        sn_cfctl%procincr  = 1       ! Increment for optional subsetting of areas [default:1] 
    1319        sn_cfctl%ptimincr  = 1       ! Timestep increment for writing time step progress info 
    1320    nn_print    =    0      !  level of print (0 no extra print) 
    1321    nn_ictls    =    0      !  start i indice of control sum (use to compare mono versus 
    1322    nn_ictle    =    0      !  end   i indice of control sum        multi processor runs 
    1323    nn_jctls    =    0      !  start j indice of control               over a subdomain) 
    1324    nn_jctle    =    0      !  end   j indice of control 
    1325    nn_isplt    =    1      !  number of processors in i-direction 
    1326    nn_jsplt    =    1      !  number of processors in j-direction 
    1327    ln_timing   = .false.   !  timing by routine write out in timing.output file 
    1328    ln_diacfl   = .false.   !  CFL diagnostics write out in cfl_diagnostics.ascii 
     1408   sn_cfctl%l_runstat = .TRUE.    ! switches and which areas produce reports with the proc integer settings. 
     1409   sn_cfctl%l_trcstat = .FALSE.   ! The default settings for the proc integers should ensure 
     1410   sn_cfctl%l_oceout  = .FALSE.   ! that  all areas report. 
     1411   sn_cfctl%l_layout  = .FALSE.   ! 
     1412   sn_cfctl%l_prtctl  = .FALSE.   ! 
     1413   sn_cfctl%l_prttrc  = .FALSE.   ! 
     1414   sn_cfctl%l_oasout  = .FALSE.   ! 
     1415   sn_cfctl%procmin   = 0         ! Minimum area number for reporting [default:0] 
     1416   sn_cfctl%procmax   = 1000000   ! Maximum area number for reporting [default:1000000] 
     1417   sn_cfctl%procincr  = 1         ! Increment for optional subsetting of areas [default:1] 
     1418   sn_cfctl%ptimincr  = 1         ! Timestep increment for writing time step progress info 
     1419   nn_ictls    =    0             !  start i indice of control sum (use to compare mono versus 
     1420   nn_ictle    =    0             !  end   i indice of control sum        multi processor runs 
     1421   nn_jctls    =    0             !  start j indice of control               over a subdomain) 
     1422   nn_jctle    =    0             !  end   j indice of control 
     1423   nn_isplt    =    1             !  number of processors in i-direction 
     1424   nn_jsplt    =    1             !  number of processors in j-direction 
     1425   ln_timing   = .false.          !  timing by routine write out in timing.output file 
     1426   ln_diacfl   = .false.          !  CFL diagnostics write out in cfl_diagnostics.ascii 
    13291427/ 
    13301428!----------------------------------------------------------------------- 

  • NEMO/branches/2019/dev_r11351_fldread_with_XIOS/cfgs/SHARED/namelist_top_ref

    r10375 r13463  
    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 
    … …  
    5959!                !  file name        ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    6060!                !                   !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    61    sn_trcdta(1)  = 'data_TRC_nomask' ,        -12        ,  'TRC'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , '' 
     61   sn_trcdta(1)  = 'data_TRC_nomask' ,        -12.       ,  'TRC'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , '' 
    6262   ! 
    6363   cn_dir        =  './'     !  root directory for the location of the data files 
    … …  
    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.