Skip to content
Snippets Groups Projects
field_def_nemo-oce.xml 156 KiB
Newer Older
Guillaume Samson's avatar
Guillaume Samson committed
<?xml version="1.0"?>
<!-- $id$ -->

<field_definition level="1" prec="4" operation="average" enabled=".TRUE." default_value="1.e20" > <!-- time step automaticaly defined -->

  <!--
      =====================================================================================================
      =                                   Configurable diagnostics                                        =
      =====================================================================================================
  -->

  <field_group id="diamlr_fields">

    <!--
        =====================================================================================================
        Configuration of multiple-linear-regression analysis (diamlr)
        =====================================================================================================

This field group configures diamlr for tidal harmonic analysis of field
ssh: in addition to a regressor for fitting the mean value (diamlr_r101),
it includes the regressors for the analysis of the tidal constituents
that are available in the tidal-forcing implementation (see
./src/OCE/SBC/tide.h90).

    -->

    <!-- Time -->
    <field id="diamlr_time" grid_ref="diamlr_grid_T_2D" prec="8" />

    <!-- Regressors for tidal harmonic analysis -->
    <field id="diamlr_r001"  field_ref="diamlr_time" expr="sin( __TDE_M2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:M2"   />
    <field id="diamlr_r002"  field_ref="diamlr_time" expr="cos( __TDE_M2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:M2"   />
    <field id="diamlr_r003"  field_ref="diamlr_time" expr="sin( __TDE_N2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:N2"   />
    <field id="diamlr_r004"  field_ref="diamlr_time" expr="cos( __TDE_N2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:N2"   />
    <field id="diamlr_r005"  field_ref="diamlr_time" expr="sin( __TDE_2N2_omega__  * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:2N2"  />
    <field id="diamlr_r006"  field_ref="diamlr_time" expr="cos( __TDE_2N2_omega__  * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:2N2"  />
    <field id="diamlr_r007"  field_ref="diamlr_time" expr="sin( __TDE_S2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:S2"   />
    <field id="diamlr_r008"  field_ref="diamlr_time" expr="cos( __TDE_S2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:S2"   />
    <field id="diamlr_r009"  field_ref="diamlr_time" expr="sin( __TDE_K2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:K2"   />
    <field id="diamlr_r010"  field_ref="diamlr_time" expr="cos( __TDE_K2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:K2"   />
    <field id="diamlr_r011"  field_ref="diamlr_time" expr="sin( __TDE_K1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:K1"   />
    <field id="diamlr_r012"  field_ref="diamlr_time" expr="cos( __TDE_K1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:K1"   />
    <field id="diamlr_r013"  field_ref="diamlr_time" expr="sin( __TDE_O1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:O1"   />
    <field id="diamlr_r014"  field_ref="diamlr_time" expr="cos( __TDE_O1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:O1"   />
    <field id="diamlr_r015"  field_ref="diamlr_time" expr="sin( __TDE_Q1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:Q1"   />
    <field id="diamlr_r016"  field_ref="diamlr_time" expr="cos( __TDE_Q1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:Q1"   />
    <field id="diamlr_r017"  field_ref="diamlr_time" expr="sin( __TDE_P1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:P1"   />
    <field id="diamlr_r018"  field_ref="diamlr_time" expr="cos( __TDE_P1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:P1"   />
    <field id="diamlr_r019"  field_ref="diamlr_time" expr="sin( __TDE_M4_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:M4"   />
    <field id="diamlr_r020"  field_ref="diamlr_time" expr="cos( __TDE_M4_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:M4"   />
    <field id="diamlr_r021"  field_ref="diamlr_time" expr="sin( __TDE_Mf_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:Mf"   />
    <field id="diamlr_r022"  field_ref="diamlr_time" expr="cos( __TDE_Mf_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:Mf"   />
    <field id="diamlr_r023"  field_ref="diamlr_time" expr="sin( __TDE_Mm_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:Mm"   />
    <field id="diamlr_r024"  field_ref="diamlr_time" expr="cos( __TDE_Mm_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:Mm"   />
    <field id="diamlr_r025"  field_ref="diamlr_time" expr="sin( __TDE_Msqm_omega__ * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:Msqm" />
    <field id="diamlr_r026"  field_ref="diamlr_time" expr="cos( __TDE_Msqm_omega__ * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:Msqm" />
    <field id="diamlr_r027"  field_ref="diamlr_time" expr="sin( __TDE_Mtm_omega__  * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:Mtm"  />
    <field id="diamlr_r028"  field_ref="diamlr_time" expr="cos( __TDE_Mtm_omega__  * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:Mtm"  />
    <field id="diamlr_r029"  field_ref="diamlr_time" expr="sin( __TDE_S1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:S1"   />
    <field id="diamlr_r030"  field_ref="diamlr_time" expr="cos( __TDE_S1_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:S1"   />
    <field id="diamlr_r031"  field_ref="diamlr_time" expr="sin( __TDE_MU2_omega__  * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:MU2"  />
    <field id="diamlr_r032"  field_ref="diamlr_time" expr="cos( __TDE_MU2_omega__  * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:MU2"  />
    <field id="diamlr_r033"  field_ref="diamlr_time" expr="sin( __TDE_NU2_omega__  * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:NU2"  />
    <field id="diamlr_r034"  field_ref="diamlr_time" expr="cos( __TDE_NU2_omega__  * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:NU2"  />
    <field id="diamlr_r035"  field_ref="diamlr_time" expr="sin( __TDE_L2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:L2"   />
    <field id="diamlr_r036"  field_ref="diamlr_time" expr="cos( __TDE_L2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:L2"   />
    <field id="diamlr_r037"  field_ref="diamlr_time" expr="sin( __TDE_T2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:sin:T2"   />
    <field id="diamlr_r038"  field_ref="diamlr_time" expr="cos( __TDE_T2_omega__   * diamlr_time )" enabled=".TRUE."  comment="harmonic:cos:T2"   />
    <field id="diamlr_r101"  field_ref="diamlr_time" expr="diamlr_time^0.0"                         enabled=".TRUE."  comment="mean"              />

    <!-- Fields selected for regression analysis -->
    <field id="diamlr_f001" field_ref="ssh" enabled=".TRUE." />

  </field_group>

  <!--
      ============================================================================================================
      =                                  definition of all existing variables                                    =
      =                                            DO NOT CHANGE                                                 =
      ============================================================================================================
  -->

  <field_group id="diadetide_fields">

    <!--
        =====================================================================================================
        Weight fields for the computation of daily detided model diagnostics (diadetide)
        =====================================================================================================

    -->

    <field id="diadetide_weight" grid_ref="diadetide_grid_T_2D" enabled=".TRUE." />
    <field id="diadetide_weight_grid_T_2D" field_ref="diadetide_weight" grid_ref="diadetide_grid_T_2D" enabled=".TRUE." > this </field>
    <field id="diadetide_weight_grid_U_2D" field_ref="diadetide_weight" grid_ref="diadetide_grid_U_2D" enabled=".TRUE." > this </field>
    <field id="diadetide_weight_grid_V_2D" field_ref="diadetide_weight" grid_ref="diadetide_grid_V_2D" enabled=".TRUE." > this </field>
    <field id="diadetide_weight_grid_T_3D" field_ref="diadetide_weight" grid_ref="diadetide_grid_2D_to_grid_T_3D" enabled=".TRUE." > this </field>
    <field id="diadetide_weight_grid_U_3D" field_ref="diadetide_weight" grid_ref="diadetide_grid_2D_to_grid_U_3D" enabled=".TRUE." > this </field>
    <field id="diadetide_weight_grid_V_3D" field_ref="diadetide_weight" grid_ref="diadetide_grid_2D_to_grid_V_3D" enabled=".TRUE." > this </field>
    <field id="diadetide_weight_grid_W_3D" field_ref="diadetide_weight" grid_ref="diadetide_grid_2D_to_grid_W_3D" enabled=".TRUE." > this </field>

  </field_group>

  <!--
      ============================================================================================================
      Physical ocean model variables
      ============================================================================================================
  -->

  <!-- T grid -->

  <field_group id="grid_T" grid_ref="grid_T_2D" >
    <field id="e3t"          long_name="T-cell thickness"                    standard_name="cell_thickness"        unit="m"   grid_ref="grid_T_3D"  />
Guillaume Samson's avatar
Guillaume Samson committed
    <field id="e3ts"         long_name="T-cell thickness"   field_ref="e3t"  standard_name="cell_thickness"        unit="m"   grid_ref="grid_T_SFC"/>
    <field id="e3t_0"        long_name="Initial T-cell thickness"            standard_name="ref_cell_thickness"    unit="m"   grid_ref="grid_T_3D" />
    <field id="e3tb"         long_name="bottom T-cell thickness"             standard_name="bottom_cell_thickness" unit="m"   grid_ref="grid_T_2D"/>
    <field id="e3t_300"      field_ref="e3t"                grid_ref="grid_T_zoom_300"       detect_missing_value="true" />
    <field id="e3t_vsum300"  field_ref="e3t_300"            grid_ref="grid_T_vsum"   detect_missing_value="true" />
    <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"/>
    <field id="volcello"     long_name="Ocean Volume"                   standard_name="ocean_volume"   unit="m3"       grid_ref="grid_T_3D" />
  
    <!-- EOS80 -->
    <field id="toce_pot"         long_name="potential temperature"           standard_name="sea_water_potential_temperature"      unit="degC"     grid_ref="grid_T_3D"/>
    <field id="toce_pot_e3t"     long_name="potential temperature (thickness weighted)"                                           unit="degC"     grid_ref="grid_T_3D" > toce_pot * e3t </field >
    <field id="soce_pra"         long_name="practical salinity"              standard_name="sea_water_practical_salinity"      unit="1e-3"     grid_ref="grid_T_3D"/>
    <field id="soce_pra_e3t"     long_name="practical salinity    (thickness weighted)"                                        unit="1e-3"     grid_ref="grid_T_3D" > soce_pra * e3t </field >
    <!-- TEOS10 -->
    <field id="toce_con"         long_name="conservative temperature"        standard_name="sea_water_conservative_temperature"   unit="degC"     grid_ref="grid_T_3D"/>
    <field id="toce_con_e3t"     long_name="conservative temperature (thickness weighted)"                                        unit="degC"     grid_ref="grid_T_3D" > toce_con * e3t </field >  
    <field id="soce_abs"         long_name="absolute salinity"               standard_name="sea_water_absolute_salinity"       unit="g/kg"     grid_ref="grid_T_3D"/>
    <field id="soce_abs_e3t"     long_name="absolute salinity    (thickness weighted)"                                         unit="g/kg"     grid_ref="grid_T_3D" > soce_abs * e3t </field >
    <!-- SEOS -->
    <field id="toce_seos"         long_name="temperature"        standard_name="sea_water_temperature"   unit="degC"     grid_ref="grid_T_3D"/>
    <field id="toce_seos_e3t"     long_name="temperature (thickness weighted)"                                        unit="degC"     grid_ref="grid_T_3D" > toce_seos * e3t </field >
    <field id="soce_seos"         long_name="salinity"               standard_name="sea_water_salinity"       unit="1e-3"     grid_ref="grid_T_3D"/>
    <field id="soce_seos_e3t"     long_name="salinity    (thickness weighted)"                                         unit="1e-3"     grid_ref="grid_T_3D" > soce_seos * e3t </field >
Guillaume Samson's avatar
Guillaume Samson committed

    <field id="toce_e3t_300"      field_ref="toce_e3t"          unit="degree_C"     grid_ref="grid_T_zoom_300"      detect_missing_value="true" />
    <field id="toce_e3t_vsum300"  field_ref="toce_e3t_300"      unit="degress_C*m"  grid_ref="grid_T_vsum"  detect_missing_value="true" />
    <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>

    <!-- AGRIF sponge -->
    <field id="agrif_spt"         long_name=" AGRIF t-sponge coefficient"   unit=" " />

    <!-- additions to diawri.F90 -->
    <field id="sssgrad"     long_name="module of surface salinity gradient"              unit="1e-3/m"   grid_ref="grid_T_2D_inner"/>
    <field id="sssgrad2"    long_name="square of module of surface salinity gradient"    unit="1e-6/m2"  grid_ref="grid_T_2D_inner"/>
    <field id="ke"          long_name="kinetic energy"          standard_name="specific_kinetic_energy_of_sea_water"   unit="m2/s2"  grid_ref="grid_T_3D" />
    <field id="ke_int"      long_name="vertical integration of kinetic energy"   unit="m3/s2"  grid_ref="grid_T_2D_inner" />
    <field id="salt2c"      long_name="Salt content vertically integrated"       unit="1e-3*kg/m2" />
Guillaume Samson's avatar
Guillaume Samson committed

    <!-- t-eddy viscosity coefficients (ldfdyn) -->
    <field id="ahmt_2d"      long_name=" surface t-eddy viscosity coefficient"   unit="m2/s or m4/s"                      />
    <field id="ahmt_3d"      long_name=" 3D      t-eddy viscosity coefficient"   unit="m2/s or m4/s"  grid_ref="grid_T_3D"/>

    
       <!-- EOS80 -->
        <field id="sst_pot"          long_name="sea surface potential temperature"                            standard_name="sea_surface_temperature"             unit="degC"     />
        <field id="sst2_pot"         long_name="square of sea surface potential temperature"                  standard_name="square_of_sea_surface_temperature"   unit="degC2"     > sst_pot * sst_pot </field >
        <field id="sstmax_pot"       long_name="max of sea surface potential temperature"   field_ref="sst_pot"   operation="maximum"                                                 />
        <field id="sstmin_pot"       long_name="min of sea surface potential temperature"   field_ref="sst_pot"   operation="minimum"                                                 />
        <field id="sstgrad_pot"      long_name="module of potential sst gradient"                                                                                 unit="degC/m"  grid_ref="grid_T_2D_inner" />
        <field id="sstgrad2_pot"     long_name="square of module of potential sst gradient"                                                                       unit="degC2/m2" grid_ref="grid_T_2D_inner" />
        <field id="sbt_pot"          long_name="sea bottom potential temperature"                                                                                 unit="degC"  grid_ref="grid_T_2D_inner" />
        <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" grid_ref="grid_T_2D_inner" />
        <field id="sst_wl_pot"       long_name="Delta potential SST of warm layer"                                                                                unit="degC"     />
        <field id="sst_cs_pot"       long_name="Delta potential SST of cool skin"                                                                                 unit="degC"     />
	<field id="temp_3m_pot"      long_name="potential temperature at 3m"                                                                                      unit="degC"     />

        <field id="sss_pra"          long_name="sea surface practical salinity"                               standard_name="sea_surface_practical_salinity"                unit="1e-3"     />
        <field id="sss2_pra"         long_name="square of sea surface practical salinity"                                                                         unit="1e-6"      > sss_pra * sss_pra </field >
        <field id="sssmax_pra"       long_name="max of sea surface practical salinity"      field_ref="sss_pra"   operation="maximum"                                                 />
        <field id="sssmin_pra"       long_name="min of sea surface practical salinity"      field_ref="sss_pra"   operation="minimum"                                                 />
        <field id="sbs_pra"          long_name="sea bottom practical salinity"                                                                                    unit="0.001"    grid_ref="grid_T_2D_inner" />
        <field id="somint_pra"       long_name="vertical integral of practical salinity times density"        standard_name="integral_wrt_depth_of_product_of_density_and_practical_salinity"  unit="(kg m2) x (1e-3)" grid_ref="grid_T_2D_inner" /> 
        
         <!-- TEOS10 -->
        <field id="sst_con"          long_name="sea surface conservative temperature"                            standard_name="sea_surface_conservative_temperature"             unit="degC"     />
        <field id="sst2_con"         long_name="square of sea surface conservative temperature"                  standard_name="square_of_sea_surface_temperature"   unit="degC2"     > sst_con * sst_con </field >
        <field id="sstmax_con"       long_name="max of sea surface conservative temperature"   field_ref="sst_con"   operation="maximum"                                                 />
        <field id="sstmin_con"       long_name="min of sea surface conservative temperature"   field_ref="sst_con"   operation="minimum"                                                 />
        <field id="sstgrad_con"      long_name="module of conservative sst gradient"                                                                                 unit="degC/m"  grid_ref="grid_T_2D_inner" />
        <field id="sstgrad2_con"     long_name="square of module of conservative sst gradient"                                                                       unit="degC2/m2" grid_ref="grid_T_2D_inner" />
        <field id="sbt_con"          long_name="sea bottom conservative temperature"                                                                                 unit="degC"    grid_ref="grid_T_2D_inner" />
        <field id="tosmint_con"      long_name="vertical integral of conservative temperature times density"     standard_name="integral_wrt_depth_of_product_of_density_and_conservative_temperature"  unit="(kg m2) degree_C" grid_ref="grid_T_2D_inner" />
        <field id="sst_wl_con"       long_name="Delta conservative SST of warm layer"                                                                                unit="degC"     />
        <field id="sst_cs_con"       long_name="Delta conservative SST of cool skin"                                                                                 unit="degC"     />
	<field id="temp_3m_con"      long_name="conservative temperature at 3m"                                                                                      unit="degC"     />

        <field id="sss_abs"          long_name="sea surface absolute salinity"                               standard_name="sea_surface_absolute_salinity"                unit="g/kg"     />
        <field id="sss2_abs"         long_name="square of sea surface absolute salinity"                                                                         unit="1e-6"      > sss_abs * sss_abs </field >
        <field id="sssmax_abs"       long_name="max of sea surface absolute salinity"      field_ref="sss_abs"   operation="maximum"                                                 />
        <field id="sssmin_abs"       long_name="min of sea surface absolute salinity"      field_ref="sss_abs"   operation="minimum"                                                 />
        <field id="sbs_abs"          long_name="sea bottom absolute salinity"                                                                                    unit="g/kg"   grid_ref="grid_T_2D_inner" />
        <field id="somint_abs"       long_name="vertical integral of absolute salinity times density"        standard_name="integral_wrt_depth_of_product_of_density_and_absolute_salinity"  unit="(kg m2) x (1e-3)" grid_ref="grid_T_2D_inner" />

        <!-- SEOS -->
        <field id="sst_seos"          long_name="sea surface temperature"                            standard_name="sea_surface_temperature"             unit="degC"     />
        <field id="sst2_seos"         long_name="square of sea surface temperature"                  standard_name="square_of_sea_surface_temperature"   unit="degC2"     > sst_seos * sst_seos </field >
        <field id="sstmax_seos"       long_name="max of sea surface temperature"   field_ref="sst_seos"   operation="maximum"                                                 />
        <field id="sstmin_seos"       long_name="min of sea surface temperature"   field_ref="sst_seos"   operation="minimum"                                                 />
        <field id="sstgrad_seos"      long_name="module of sst gradient"                                                                                 unit="degC/m"  grid_ref="grid_T_2D_inner" />
        <field id="sstgrad2_seos"     long_name="square of module of sst gradient"                                                                       unit="degC2/m2" grid_ref="grid_T_2D_inner" />
        <field id="sbt_seos"          long_name="sea bottom temperature"                                                                                 unit="degC"     grid_ref="grid_T_2D_inner" />
        <field id="tosmint_seos"      long_name="vertical integral of temperature times density"     standard_name="integral_wrt_depth_of_product_of_density_and_temperature"  unit="(kg m2) degree_C" grid_ref="grid_T_2D_inner" />
        <field id="sst_wl_seos"       long_name="Delta SST of warm layer"                                                                                unit="degC"     />
        <field id="sst_cs_seos"       long_name="Delta SST of cool skin"                                                                                 unit="degC"     />
	<field id="temp_3m_seos"      long_name="temperature at 3m"                                                                                      unit="degC"     />

        <field id="sss_seos"          long_name="sea surface salinity"                               standard_name="sea_surface_salinity"                unit="1e-3"     />
        <field id="sss2_seos"         long_name="square of sea surface salinity"                                                                         unit="1e-6"      > sss_seos * sss_seos </field >
        <field id="sssmax_seos"       long_name="max of sea surface salinity"      field_ref="sss_seos"   operation="maximum"                                                 />
        <field id="sssmin_seos"       long_name="min of sea surface salinity"      field_ref="sss_seos"   operation="minimum"                                                 />
        <field id="sbs_seos"          long_name="sea bottom salinity"                                                                                    unit="0.001"   grid_ref="grid_T_2D_inner" />
        <field id="somint_seos"       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)"  grid_ref="grid_T_2D_inner" />


    <field id="taubot"       long_name="bottom stress module"   unit="N/m2"   grid_ref="grid_T_2D_inner" />
    <field id="htau"         long_name="htau length scale" unit="m"/>
Guillaume Samson's avatar
Guillaume Samson committed
    <!-- Case EOS = TEOS-10 : output potential temperature -->
    <field id="toce_pot"     long_name="Sea Water Potential Temperature"              standard_name="sea_water_potential_temperature"   unit="degC"     grid_ref="grid_T_3D"/>
    <field id="sst_pot"      long_name="potential sea surface temperature"             standard_name="sea_surface_temperature"             unit="degC"     />
    <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" />

    <field id="ht"           long_name="water column height at T point"                     standard_name="water_column_height_T"                      unit="m" />
    <field id="ssh"          long_name="sea surface height"                                 standard_name="sea_surface_height_above_geoid"             unit="m" />
    <field id="ssh2"         long_name="square of sea surface height"                       standard_name="square_of_sea_surface_height_above_geoid"   unit="m2" > ssh * ssh </field >
    <field id="wetdep"       long_name="wet depth"                                          standard_name="wet_depth"                                  unit="m" />
    <field id="sshmax"       long_name="max of sea surface height"        field_ref="ssh"   operation="maximum"                                                 />

    <field id="mldkz5"       long_name="Turbocline depth (Kz = 5e-4)"                       standard_name="ocean_mixed_layer_thickness_defined_by_vertical_tracer_diffusivity"                unit="m"          />
    <field id="mldr10_1"     long_name="Mixed Layer Depth (dsigma = 0.01 wrt 10m)"          standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta"                                unit="m"          />
    <field id="mldr10_1max"  long_name="Max of Mixed Layer Depth (dsigma = 0.01 wrt 10m)"   field_ref="mldr10_1"   operation="maximum"                                                                          />
    <field id="mldr10_1min"  long_name="Min of Mixed Layer Depth (dsigma = 0.01 wrt 10m)"   field_ref="mldr10_1"   operation="minimum"                                                                          />
    <field id="mldzint_1"    long_name="Mixed Layer Depth interpolated"                     standard_name="ocean_mixed_layer_thickness"                                                       unit="m"      grid_ref="grid_T_2D"    />
    <field id="mldzint_2"    long_name="Mixed Layer Depth interpolated"                     standard_name="ocean_mixed_layer_thickness"                                                       unit="m"      grid_ref="grid_T_2D"    />
    <field id="mldzint_3"    long_name="Mixed Layer Depth interpolated"                     standard_name="ocean_mixed_layer_thickness"                                                       unit="m"      grid_ref="grid_T_2D"    />
    <field id="mldzint_4"    long_name="Mixed Layer Depth interpolated"                     standard_name="ocean_mixed_layer_thickness"                                                       unit="m"      grid_ref="grid_T_2D"    />
    <field id="mldzint_5"    long_name="Mixed Layer Depth interpolated"                     standard_name="ocean_mixed_layer_thickness"                                                       unit="m"      grid_ref="grid_T_2D"    />
    <field id="mldhtc_1"     long_name="Mixed Layer Depth integrated heat content"          standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content"   unit="J/m2"   grid_ref="grid_T_2D"    />
    <field id="mldhtc_2"     long_name="Mixed Layer Depth integrated heat content"          standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content"   unit="J/m2"   grid_ref="grid_T_2D"   />
    <field id="mldhtc_3"     long_name="Mixed Layer Depth integrated heat content"          standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content"   unit="J/m2"   grid_ref="grid_T_2D"    />
    <field id="mldhtc_4"     long_name="Mixed Layer Depth integrated heat content"          standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content"   unit="J/m2"   grid_ref="grid_T_2D"    />
    <field id="mldhtc_5"     long_name="Mixed Layer Depth integrated heat content"          standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content"   unit="J/m2"   grid_ref="grid_T_2D"    />
    <field id="heatc"        long_name="Heat content vertically integrated"                 standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content"   unit="J/m2"   grid_ref="grid_T_2D_inner" />
    <field id="saltc"        long_name="Salt content vertically integrated"                                                                                                                   unit="1e-3.m"    grid_ref="grid_T_2D_inner" />
    <field id="saltc2"       long_name="square of Salt content vertically integrated"                                                                                                         unit="1e-6.m2"   grid_ref="grid_T_2D_inner" />
Guillaume Samson's avatar
Guillaume Samson committed

    <!-- EOS -->
    <field id="alpha"        long_name="thermal expansion"                                                         unit="degC-1" grid_ref="grid_T_3D" />
    <field id="beta"         long_name="haline contraction"                                                        unit="1e3"    grid_ref="grid_T_3D" />
    <field id="rhop"         long_name="potential density (sigma0)"        standard_name="sea_water_sigma_theta"   unit="kg/m3"  grid_ref="grid_T_3D" />
    <field id="rhop_e3t"     long_name="potential density (sigma0) thickness weighted"        standard_name="sea_water_sigma_theta"   unit="kg/m3"  grid_ref="grid_T_3D"> rhop * e3t </field>
    <field id="N_2d"         long_name="Depth-mean N"                                                              unit="m/s" />
    <field id="modslp"       long_name="sqrt( slpi^2 + slpj^2 )"                                                   unit="1"      grid_ref="grid_T_3D" />
    <field id="RossRad"      long_name="Rossby radius"                                                             unit="m" />
    <field id="RossRadlim"   long_name="Rossby radius (limited)"                                                   unit="m" />
    <field id="Tclinic_recip" long_name="recip of baroclinic timescale"                                            unit="s-1" />
    <field id="RR_GS"        long_name="Rossby radius / min(dx,dy)"                                                unit="1" />
Guillaume Samson's avatar
Guillaume Samson committed

    <!-- Energy - horizontal divergence -->
    <field id="eken"         long_name="kinetic energy"          standard_name="specific_kinetic_energy_of_sea_water"   unit="m2/s2"  grid_ref="grid_T_3D" />
Guillaume Samson's avatar
Guillaume Samson committed
    <field id="sKE"          long_name="surface kinetic energy"  standard_name="specific_kinetic_energy_of_sea_water"   unit="m2/s2"  grid_ref="grid_T_2D_inner" />
    <field id="hdiv"         long_name="horizontal divergence"                                                          unit="s-1"    grid_ref="grid_T_3D" />

    <!-- variables available with MLE -->
    <field id="Lf_NHpf"      long_name="MLE: Lf = N H / f"   unit="m" />

    <!-- next variables available with ln_zad_Aimp=.true. -->
    <field id="Courant"      long_name="Courant number"                                                                 unit="#"   grid_ref="grid_T_3D" />
    <field id="wimp"         long_name="Implicit vertical velocity"                                                     unit="m/s" grid_ref="grid_T_3D" />
    <field id="wexp"         long_name="Explicit vertical velocity"                                                     unit="m/s" grid_ref="grid_T_3D" />
    <field id="wi_cff"       long_name="Fraction of implicit vertical velocity"                                         unit="#"   grid_ref="grid_T_3D" />

    <!-- next variables available with key_diahth -->
    <field id="mlddzt"       long_name="Thermocline Depth (depth of max dT/dz)"         standard_name="depth_at_maximum_upward_derivative_of_sea_water_potential_temperature"             unit="m"                         />
    <field id="mldr10_3"     long_name="Mixed Layer Depth (dsigma = 0.03 wrt 10m)"      standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta"                                unit="m"                         />
    <field id="mldr0_1"      long_name="Mixed Layer Depth (dsigma = 0.01 wrt sfc)"      standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta"                                unit="m"                         />
    <field id="mldr0_3"      long_name="Mixed Layer Depth (dsigma = 0.03 wrt sfc)"      standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta"                                unit="m"                         />
    <field id="mld_dt02"     long_name="Mixed Layer Depth (|dT| = 0.2 wrt 10m)"         standard_name="ocean_mixed_layer_thickness_defined_by_temperature"                                unit="m"                         />
    <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"                         />
    <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"                         />
    <field id="BLT"          long_name="Barrier Layer Thickness"                                                                                                                          unit="m"   > topthdep - pycndep </field>
    <field id="tinv"         long_name="Max of vertical invertion of temperature"                                                                                                         unit="degC"                      />
    <field id="depti"        long_name="Depth of max. vert. inv. of temperature"                                                                                                          unit="m"                         />
    <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" />
    <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"   />
    <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" />
    <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"                      />
    <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"                      />
    <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"                      />

    <!-- variables available with diaar5 -->
    <field id="botpres"      long_name="Sea Water Pressure at Sea Floor"          standard_name="sea_water_pressure_at_sea_floor"                    unit="dbar" />
    <field id="sshdyn"       long_name="dynamic sea surface height"               standard_name="dynamic_sea_surface_height_above_geoid"             unit="m"    />
    <field id="sshdyn2"      long_name="square of dynamic sea surface height"     standard_name="dynamic_sea_surface_height_above_geoid_squared"     unit="m2"    > sshdyn * sshdyn </field>
    <field id="tnpeo"      long_name="Tendency of ocean potential energy content"                                                                    unit="W/m2" />

    <!-- variables available ln_linssh=.FALSE. -->
    <field id="tpt_dep"      long_name="T-point depth"                  standard_name="depth_below_geoid"   unit="m"   grid_ref="grid_T_3D" />
    <field id="e3tdef"       long_name="T-cell thickness deformation"                                       unit="%"   grid_ref="grid_T_3D" />

    <!-- variables available with ln_diacfl=.true. -->
    <field id="cfl_cu"       long_name="u-courant number"   unit="#" />
    <field id="cfl_cv"       long_name="v-courant number"   unit="#" />
    <field id="cfl_cw"       long_name="w-courant number"   unit="#" />

    <!-- variables available with ln_zdfmfc=.true. -->
    <field id="mf_Tp"       long_name="plume_temperature"      standard_name="plume_temperature"     unit="degC"   grid_ref="grid_T_3D" />
    <field id="mf_Sp"       long_name="plume_salinity"         standard_name="plume_salinity"        unit="1e-3"   grid_ref="grid_T_3D" />
    <field id="mf_mf"       long_name="mass flux"              standard_name="mf_mass_flux"          unit="m"      grid_ref="grid_T_3D" />

    <!-- fluxes from damping -->
    <field id="sflx_dmp_cea"  long_name="salt flux due to damping"  standard_name="salt_flux_due_to_damping"          unit="g/m2/s"   />
    <field id="hflx_dmp_cea"  long_name="heat flux due to damping"  standard_name="heat_flux_due_to_damping"          unit="W/m2"     />

    <!-- * variable related to ice shelf forcing * -->

    <!-- * fwf * -->
    <field id="fwfisf_cav"      long_name="Ice shelf fresh water flux ( from isf to oce )"                    unit="kg/m2/s"  />
    <field id="fwfisf_par"      long_name="Ice shelf fresh water flux ( from isf to oce )"                    unit="kg/m2/s"  />
    <field id="fwfisf3d_cav"    long_name="3d Ice shelf fresh water flux ( from isf to oce )"                 unit="kg/m2/s"  grid_ref="grid_T_3D" />
    <field id="fwfisf3d_par"    long_name="3d Ice shelf fresh water flux ( from isf to oce )"                 unit="kg/m2/s"  grid_ref="grid_T_3D" />

    <!-- * heat fluxes * -->
    <field id="qoceisf_cav"     long_name="Ice shelf ocean  heat flux ( from isf to oce )"                    unit="W/m2"     />
    <field id="qoceisf_par"     long_name="Ice shelf ocean  heat flux ( from isf to oce )"                    unit="W/m2"     />
    <field id="qlatisf_cav"     long_name="Ice shelf latent heat flux ( from isf to oce )"                    unit="W/m2"     />
    <field id="qlatisf_par"     long_name="Ice shelf latent heat flux ( from isf to oce )"                    unit="W/m2"     />
    <field id="qhcisf_cav"      long_name="Ice shelf heat content flux of injected water ( from isf to oce )" unit="W/m2"     />
    <field id="qhcisf_par"      long_name="Ice shelf heat content flux of injected water ( from isf to oce )" unit="W/m2"     />
    <field id="qoceisf3d_cav"   long_name="Ice shelf ocean  heat flux ( from isf to oce )"                    unit="W/m2"     grid_ref="grid_T_3D" />
    <field id="qoceisf3d_par"   long_name="Ice shelf ocean  heat flux ( from isf to oce )"                    unit="W/m2"     grid_ref="grid_T_3D" />
    <field id="qlatisf3d_cav"   long_name="Ice shelf latent heat flux ( from isf to oce )"                    unit="W/m2"     grid_ref="grid_T_3D" />
    <field id="qlatisf3d_par"   long_name="Ice shelf latent heat flux ( from isf to oce )"                    unit="W/m2"     grid_ref="grid_T_3D" />
    <field id="qhcisf3d_cav"    long_name="Ice shelf heat content flux of injected water ( from isf to oce )" unit="W/m2"     grid_ref="grid_T_3D" />
    <field id="qhcisf3d_par"    long_name="Ice shelf heat content flux of injected water ( from isf to oce )" unit="W/m2"     grid_ref="grid_T_3D" />
    <field id="qconisf"         long_name="Conductive heat flux through the ice shelf ( from isf to oce )"    unit="W/m2"     />

    <field id="qlatisf"      long_name="Ice shelf latent heat flux"                    unit="W/m2"     />
    <field id="qhcisf"       long_name="Ice shelf heat content flux"                   unit="W/m2"     />
    <field id="fwfisf"       long_name="Ice shelf melting"                             unit="kg/m2/s"  />

Guillaume Samson's avatar
Guillaume Samson committed
    <!-- top boundary layer properties -->
    <field id="isftfrz_cav"     long_name="freezing point temperature at ocean/isf interface"                unit="degC"     />
    <field id="isftfrz_par"     long_name="freezing point temperature in the parametrization boundary layer" unit="degC"     />
    <field id="isfthermald_cav" long_name="thermal driving of ice shelf melting"          unit="degC"     />
    <field id="isfthermald_par" long_name="thermal driving of ice shelf melting"          unit="degC"     />
    <field id="isfgammat"       long_name="Ice shelf heat-transfert velocity"             unit="m/s"      />
    <field id="isfgammas"       long_name="Ice shelf salt-transfert velocity"             unit="m/s"      />
    <field id="ttbl_cav"        long_name="temperature in Losch tbl"                      unit="degC"     />
    <field id="ttbl_par"        long_name="temperature in the parametrisation boundary layer" unit="degC" />
    <field id="stbl"            long_name="salinity in the Losh tbl"                      unit="1e-3"     />
    <field id="utbl"            long_name="zonal current in the Losh tbl at T point"      unit="m/s"      />
    <field id="vtbl"            long_name="merid current in the Losh tbl at T point"      unit="m/s"      />
    <field id="isfustar"        long_name="ustar at T point used in ice shelf melting"    unit="m/s"      />

  </field_group> <!-- grid_T -->

  <!-- Tides -->

  <field_group id="Tides_T" grid_ref="grid_T_2D" operation="once" >
    <!-- Tidal potential -->
    <field id="tide_pot"      long_name="Total tidal potential" unit="m" />
    <field id="tide_pot_M2"   long_name="M2 tidal potential"    unit="m" />
    <field id="tide_pot_N2"   long_name="N2 tidal potential"    unit="m" />
    <field id="tide_pot_2N2"  long_name="2N2 tidal potential"   unit="m" />
    <field id="tide_pot_S2"   long_name="S2 tidal potential"    unit="m" />
    <field id="tide_pot_K2"   long_name="K2 tidal potential"    unit="m" />
    <field id="tide_pot_K1"   long_name="K1 tidal potential"    unit="m" />
    <field id="tide_pot_O1"   long_name="O1 tidal potential"    unit="m" />
    <field id="tide_pot_Q1"   long_name="Q1 tidal potential"    unit="m" />
    <field id="tide_pot_P1"   long_name="P1 tidal potential"    unit="m" />
    <field id="tide_pot_M4"   long_name="M4 tidal potential"    unit="m" />
    <field id="tide_pot_Mf"   long_name="Mf tidal potential"    unit="m" />
    <field id="tide_pot_Mm"   long_name="Mm tidal potential"    unit="m" />
    <field id="tide_pot_Msqm" long_name="Msqm tidal potential"  unit="m" />
    <field id="tide_pot_Mtm"  long_name="Mtm tidal potential"   unit="m" />
    <field id="tide_pot_S1"   long_name="S1 tidal potential"    unit="m" />
    <field id="tide_pot_MU2"  long_name="MU2 tidal potential"   unit="m" />
    <field id="tide_pot_NU2"  long_name="NU2 tidal potential"   unit="m" />
    <field id="tide_pot_L2"   long_name="L2 tidal potential"    unit="m" />
    <field id="tide_pot_T2"   long_name="T2 tidal potential"    unit="m" />
  </field_group>

  <!-- OSMOSIS variables (available with ln_zdfosm=.true.) -->

  <field_group id="OSMOSIS_T" grid_ref="grid_T_2D">
    <field id="hml"                 long_name="mixed layr depth"                         unit="m"       />
    <field id="hbl"                 long_name="boundary layer depth"                     unit="m"       />
    <field id="dh"                  long_name="Pycnocline thickness"                     unit=" m"      />
    <field id="ibld"                long_name="index of boundary layer depth"            unit="#"       />
    <field id="imld"                long_name="index of mixed layer depth"               unit="#"       />
    <field id="jp_ext"              long_name="flag =1 if pycnocline well resolved"      unit="#"       />
    <field id="j_ddh"               long_name="index of mixed layer depth"               unit="#"       />
    <field id="zshear"              long_name="shear production of TKE "                 unit="m^3/s^3" />
    <field id="zhbl"                long_name="boundary layer depth -grid"               unit="m"       />
    <field id="zhml"                long_name="mixed layer depth - grid"                 unit="m"       />
    <field id="zdh"                 long_name="Pycnocline  depth - grid"                 unit=" m"      />
    <field id="zustke"              long_name="magnitude of stokes drift  at T-points"   unit="m/s"     />
    <field id="us_x"                long_name="i component of active Stokes drift"       unit="m/s"     />
    <field id="us_y"                long_name="j component of active Stokes drift"       unit="m/s"     />
    <field id="dstokes"             long_name="stokes drift  depth scale"                unit="m"       />
    <field id="zwth0"               long_name="surface non-local temperature flux"       unit="deg m/s" />
    <field id="zws0"                long_name="surface non-local salinity flux"          unit="psu m/s" />
    <field id="zwb0"                long_name="surface non-local buoyancy flux"          unit="m^2/s^3" />
    <field id="zwstrc"              long_name="convective velocity scale"                unit="m/s"     />
    <field id="zustar"              long_name="friction velocity"                        unit="m/s"     />
    <field id="zwstrl"              long_name="langmuir velocity scale"                  unit="m/s"     />
    <field id="zvstr"               long_name="mixed velocity scale"                     unit="m/s"     />
    <field id="zla"                 long_name="langmuir number"                          unit="#"       />
    <field id="wmp"                 long_name="wave mean period"                         unit="s"       />
    <field id="hsw"                 long_name="significant wave height"                  unit="m"       />
    <field id="wmp_NP"              long_name="wave mean period from NP spectrum"        unit="s"       />
    <field id="hsw_NP"              long_name="significant wave height from NP spectrum" unit="m"       />
    <field id="wndm"                long_name="U_10"                                     unit="m/s"     />
Guillaume Samson's avatar
Guillaume Samson committed
    <field id="wind_wave_abs_power" long_name="\rho |U_s| x u*^2"                        unit="mW"      />
    <field id="wind_wave_power"     long_name="U_s \dot  tau"                            unit="mW"      />
    <field id="wind_power"          long_name="\rho  u*^3"                               unit="mW"      />

    <!-- interior BL OSMOSIS diagnostics -->
    <field id="zwbav"               long_name="av turb flux of buoyancy in ml"           unit="m^2/s^3" />
    <field id="zt_ml"               long_name="av T in ml"                               unit="deg"     />
    <field id="zhol"                long_name="Hoenekker number"                         unit="#"       />
    <field id="zws_ent"            long_name="entrainment turb flux of S"                unit="10^-3 m/s" />
    <field id="zwth_ent"            long_name="entrainment turb flux of T"               unit="deg m/s" />
    <field id="zwb_ent"            long_name="entrainment turb flux of buoyancy"         unit="m^2/s^-3" />

    <field id="zdt_bl"             long_name="temperature jump at base of BL"            unit="deg"      />
    <field id="zds_bl"             long_name="salinity jump at base of BL"               unit="10^-3"    />
    <field id="zdb_bl"             long_name="buoyancy jump at base of BL"               unit="m/s^2"    />
    <field id="zdu_bl"             long_name="u jump at base of BL"                      unit="m/s"      />
    <field id="zdv_bl"             long_name="v jump at base of BL"                      unit="m/s"      />
    <field id="zdt_ml"             long_name="temperature jump at base of ML"            unit="deg"      />
    <field id="zds_ml"             long_name="salinity jump at base of ML"               unit="10^-3"    />
    <field id="zdb_ml"             long_name="buoyancy jump at base of ML"               unit="m/s^2"    />
    <field id="pb_coup"            long_name="bottom coupling velocity"                  unit="m/s"      />
    <!-- extra OSMOSIS diagnostics for debugging -->
    <field id="zsc_uw_1_0"       long_name="zsc u-momentum flux on T after Stokes"                       unit="m^2/s^2" />
    <field id="zsc_uw_1_f"       long_name="zsc u-momentum flux on T after Coriolis"                       unit="m^2/s^2" />
    <field id="zsc_vw_1_f"       long_name="zsc v-momentum flux on T after Coriolis"                       unit="m^2/s^2" />
    <field id="zsc_uw_2_f"       long_name="2nd zsc u-momentum flux on T after Coriolis"                       unit="m^2/s^2" />
    <field id="zsc_vw_2_f"       long_name="2nd zsc v-momentum flux on T after Coriolis"                       unit="m^2/s^2" />

    <!-- FK_OSM OSMOSIS diagnostics (require also ln_osm_mle=.true.-->
    <field id="hmle"          long_name="OBL FK-layer thickness"                                     unit="m"        />
    <field id="mld_prof"              long_name="FK-layer depth index"                  unit="#" />
    <field id="zmld"          long_name="target FK-layer thickness"                                     unit="m"        />
    <field id="zwb_fk"          long_name="FK b-flux"                                     unit="m^2 s^-3"        />
    <field id="zwb_fk_b"          long_name="layer averaged FK b-flux"                 unit="m^2 s^-3"       />
    <field id="zdiff_mle"          long_name="max FK diffusivity in MLE"       unit=" 10^-4 m^2 s^-1"       />
    <field id="zvel_mle"          long_name="FK velocity scale in MLE"       unit=" m s^-1"       />
  </field_group>

  <field_group id="OSMOSIS_W" grid_ref="grid_W_3D" >
    <field id="zviscos"       long_name="BL viscosity"   unit="m^2/s" />
    <field id="ghamt"       long_name="non-local temperature flux"                       unit="deg m/s" />
    <field id="ghams"       long_name="non-local salinity flux"                          unit="psu m/s" />
    <field id="zdtdz_pyc"   long_name="Pycnocline temperature gradient"                  unit=" deg/m"  />
    <field id="zdsdz_pyc"   long_name="Pycnocline salinity gradient"                  unit=" 10^-3/m"  />
    <field id="zdbdz_pyc"   long_name="Pycnocline buoyancy gradient"                  unit=" s^-2"  />
    <field id="zdudz_pyc"   long_name="Pycnocline u gradient"                  unit=" s^-2"  />
    <field id="zdvdz_pyc"   long_name="Pycnocline v gradient"                  unit=" s^-2"  />

    <!-- extra OSMOSIS diagnostics for debugging -->
    <field id="ghamu_00"       long_name="initial non-local u-momentum flux"   unit="m^2/s^2" />
    <field id="ghamv_00"       long_name="initial non-local v-momentum flux"   unit="m^2/s^2" />
    <field id="ghamu_0"       long_name="after dstokes non-local u-momentum flux"   unit="m^2/s^2" />
    <field id="ghamu_f"       long_name="after Coriolis non-local u-momentum flux"   unit="m^2/s^2" />
    <field id="ghamv_f"       long_name="after Coriolis  non-local v-momentum flux"   unit="m^2/s^2" />
    <field id="ghamu_b"       long_name="after buoyancy added non-local u-momentum flux"   unit="m^2/s^2" />
    <field id="ghamv_b"       long_name="after buoyancy added  non-local v-momentum flux"  unit="m^2/s^2" />
    <field id="ghamu_1"       long_name="after entrainment non-local u-momentum flux"   unit="m^2/s^2" />
    <field id="ghamv_1"       long_name="after entrainment  non-local v-momentum flux"  unit="m^2/s^2" />
  </field_group>

  <field_group id="OSMOSIS_U" grid_ref="grid_U_2D" >
    <field id="ghamu"       long_name="non-local u-momentum flux"   grid_ref="grid_U_3D" unit="m^2/s^2" />
    <!-- FK_OSM OSMOSIS diagnostics (require also ln_osm_mle=.true.-->
    <field id="zdtdx"          long_name="FK  T x-gradient"                                     unit=" deg C m^-1"        />
    <field id="zdsdx"          long_name="FK  S x-gradient"                                     unit=" 10^-3 m^-1"        />
    <field id="dbdx_mle"          long_name="FK  B x-gradient"                                     unit=" s^-2"        />
  </field_group>

  <field_group id="OSMOSIS_V" grid_ref="grid_V_2D" >
    <field id="ghamv"       long_name="non-local v-momentum flux"   grid_ref="grid_V_3D" unit="m^2/s^2" />
    <!-- FK_OSM OSMOSIS diagnostics (require also ln_osm_mle=.true.-->
    <field id="zdtdy"          long_name="FK T y-gradient"                                     unit=" deg C m^-1"        />
    <field id="zdsdy"          long_name="FK S y-gradient"                                     unit=" 10^-3 m^-1"        />
    <field id="dbdy_mle"          long_name="FK B y-gradient"                                     unit=" s^-2"        />
  </field_group>

  <!-- SBC -->
  <field_group id="SBC" > <!-- time step automaticaly defined based on nn_fsbc -->

    <field_group id="SBC_2D" grid_ref="grid_T_2D" >

      <field id="empmr"        long_name="Net Upward Water Flux"                standard_name="water_flux_out_of_sea_ice_and_sea_water"                              unit="kg/m2/s"   />
      <field id="empbmr"       long_name="Net Upward Water Flux at pre. tstep"  standard_name="water_flux_out_of_sea_ice_and_sea_water"                              unit="kg/m2/s"   />
      <field id="emp_oce"      long_name="Evap minus Precip over ocean"         standard_name="evap_minus_precip_over_sea_water"                                     unit="kg/m2/s"   />
      <field id="emp_ice"      long_name="Evap minus Precip over ice"           standard_name="evap_minus_precip_over_sea_ice"                                       unit="kg/m2/s"   />
      <field id="saltflx"      long_name="Downward salt flux"                                                                                                        unit="g/m2/s"    />
      <field id="fmmflx"       long_name="Water flux due to freezing/melting"                                                                                        unit="kg/m2/s"   />
      <field id="snowpre"      long_name="Snow precipitation"                   standard_name="snowfall_flux"                                                        unit="kg/m2/s"   />
      <field id="runoffs"      long_name="River Runoffs"                        standard_name="water_flux_into_sea_water_from_rivers"                                unit="kg/m2/s"   />
      <field id="precip"       long_name="Total precipitation"                  standard_name="precipitation_flux"                                                   unit="kg/m2/s"   />
      <field id="wclosea"      long_name="closed sea empmr correction"          standard_name="closea_empmr"                                                         unit="kg/m2/s"   />

      <field id="qt"           long_name="Net Downward Heat Flux"                standard_name="surface_downward_heat_flux_in_sea_water"                              unit="W/m2"                           />
      <field id="qns"          long_name="non solar Downward Heat Flux"                                                                                               unit="W/m2"                           />
      <field id="qsr"          long_name="Shortwave Radiation"                   standard_name="net_downward_shortwave_flux_at_sea_water_surface"                     unit="W/m2"                           />
      <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_inner" />
Guillaume Samson's avatar
Guillaume Samson committed
      <field id="qrp"          long_name="Surface Heat Flux: Damping"            standard_name="heat_flux_into_sea_water_due_to_newtonian_relaxation"                 unit="W/m2"                           />
      <field id="qclosea"      long_name="closed sea heat content flux"          standard_name="closea_heat_content_downward_flux"                                    unit="W/m2"     />
      <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"                        />
      <field id="taum"         long_name="wind stress module"                    standard_name="magnitude_of_surface_downward_stress"                                 unit="N/m2"                           />
      <field id="wspd"         long_name="wind speed module"                     standard_name="wind_speed"                                                           unit="m/s"                            />

      <!-- * variable relative to atmospheric pressure forcing : available with ln_apr_dyn -->
      <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"        />

      <!-- *_oce variables available with ln_blk_clio or ln_blk_core -->
      <field id="rho_air"      long_name="Air density at 10m above sea surface"         standard_name="rho_air_10m"                                        unit="kg/m3" />
      <field id="dt_skin"      long_name="SSST-SST temperature difference"              standard_name="SSST-SST"                                             unit="K"   />
      <field id="qlw_oce"      long_name="Longwave Downward Heat Flux over open ocean"  standard_name="surface_net_downward_longwave_flux"                 unit="W/m2"  />
      <field id="qsb_oce"      long_name="Sensible Downward Heat Flux over open ocean"  standard_name="surface_downward_sensible_heat_flux"                unit="W/m2"  />
      <field id="qla_oce"      long_name="Latent Downward Heat Flux over open ocean"    standard_name="surface_downward_latent_heat_flux"                  unit="W/m2"  />
      <field id="evap_oce"     long_name="Evaporation over open ocean"                  standard_name="evaporation"                                        unit="kg/m2/s" />
      <field id="qt_oce"       long_name="total flux at ocean surface"                  standard_name="surface_downward_heat_flux_in_sea_water"            unit="W/m2"  />
      <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"  />
      <field id="qns_oce"      long_name="non-solar heat flux at ocean surface (including E-P)"                                                            unit="W/m2"  />
      <field id="qemp_oce"     long_name="Downward Heat Flux from E-P over open ocean"                                                                     unit="W/m2"  />
      <field id="taum_oce"     long_name="wind stress module over open ocean"           standard_name="magnitude_of_surface_downward_stress"               unit="N/m2"  />
      <field id="utau_oce"     long_name="Wind Stress along i-axis over open ocean (T-points)"  standard_name="surf_down_x_stress_open_oce_Tpoints"        unit="N/m2"   />
      <field id="vtau_oce"     long_name="Wind Stress along j-axis over open ocean (T-points)"  standard_name="surf_down_y_stress_open_oce_Tpoints"        unit="N/m2"    />

      <!-- variables computed by the bulk parameterization algorithms (ln_blk) -->
      <field id="Cd_oce"      long_name="Drag coefficient over open ocean"              standard_name="drag_coefficient_water"                unit=""  />
      <field id="Ce_oce"      long_name="Evaporaion coefficient over open ocean"        standard_name="evap_coefficient_water"                unit=""  />
      <field id="Ch_oce"      long_name="Sensible heat coefficient over open ocean"     standard_name="sensible_heat_coefficient_water"       unit=""  />
      <field id="theta_zt"    long_name="Potential air temperature at z=zt"             standard_name="potential_air_temperature_at_zt"       unit="degC" />
      <field id="q_zt"        long_name="Specific air humidity at z=zt"                 standard_name="specific_air_humidity_at_zt"           unit="kg/kg" />
      <field id="theta_zu"    long_name="Potential air temperature at z=zu"             standard_name="potential_air_temperature_at_zu"       unit="degC" />
      <field id="q_zu"        long_name="Specific air humidity at z=zu"                 standard_name="specific_air_humidity_at_zu"           unit="kg/kg" />
      <field id="ssq"         long_name="Saturation specific humidity of air at z=0"    standard_name="surface_air_saturation_spec_humidity"  unit="kg/kg" />
      <field id="wspd_blk"    long_name="Bulk wind speed at z=zu"                       standard_name="bulk_wind_speed_at_zu"                 unit="m/s"   />
      <!-- ln_blk + key_si3 -->
      <field id="Cd_ice"      long_name="Drag coefficient over ice"                     standard_name="drag_coefficient_ice"                 unit=""  />
      <field id="Ce_ice"      long_name="Evaporaion coefficient over ice"               standard_name="evap_coefficient_ice"                 unit=""  />
      <field id="Ch_ice"      long_name="Sensible heat coefficient over ice"            standard_name="sensible_heat_coefficient_ice"        unit=""  />

      <!-- available key_oasis3 -->
      <field id="snow_ao_cea"  long_name="Snow over ice-free ocean (cell average)"   standard_name="snowfall_flux"                             unit="kg/m2/s"  />
      <field id="snow_ai_cea"  long_name="Snow over sea-ice (cell average)"          standard_name="snowfall_flux"                             unit="kg/m2/s"  />
      <field id="subl_ai_cea"  long_name="Sublimation over sea-ice (cell average)"   standard_name="surface_snow_and_ice_sublimation_flux"     unit="kg/m2/s"  />
      <field id="icealb_cea"   long_name="Ice albedo (cell average)"                 standard_name="sea_ice_albedo"                            unit="1"        />
      <field id="calving_cea"  long_name="Calving"                                   standard_name="water_flux_into_sea_water_from_icebergs"   unit="kg/m2/s"  />
      <field id="iceberg_cea"  long_name="Iceberg"                                   standard_name="water_flux_into_sea_water_from_icebergs"   unit="kg/m2/s"  />
      <field id="iceshelf_cea" long_name="Iceshelf"                                  standard_name="water_flux_into_sea_water_from_iceshelf"   unit="kg/m2/s"  />


      <!-- available if key_oasis3 + conservative method -->
      <field id="rain"          long_name="Liquid precipitation"                                     standard_name="rainfall_flux"                                                                 unit="kg/m2/s"  />
      <field id="rain_ao_cea"   long_name="Liquid precipitation over ice-free ocean (cell average)"  standard_name="rainfall_flux"                                                                 unit="kg/m2/s"  />
      <field id="evap_ao_cea"   long_name="Evaporation over ice-free ocean (cell average)"           standard_name="water_evaporation_flux"                                                        unit="kg/m2/s"  />
      <field id="isnwmlt_cea"   long_name="Snow over Ice melting (cell average)"                     standard_name="surface_snow_melt_flux"                                                        unit="kg/m2/s"  />
      <field id="fsal_virt_cea" long_name="Virtual salt flux due to ice formation (cell average)"    standard_name="virtual_salt_flux_into_sea_water_due_to_sea_ice_thermodynamics"                unit="kg/m2/s"  />
      <field id="fsal_real_cea" long_name="Real salt flux due to ice formation (cell average)"       standard_name="downward_sea_ice_basal_salt_flux"                                              unit="kg/m2/s"  />
      <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"     />
      <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"     />
      <field id="hflx_subl_cea" long_name="heat flux due to sublimation (from atm. forcings)"        standard_name="temperature_flux_due_to_sublimation_expressed_as_heat_flux_out_of_sea_ice"     unit="W/m2"     />
      <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"     />
      <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"     />
      <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"     />
      <field id="hflx_snow_ao_cea" long_name="heat flux due to snow falling over ice-free ocean"     standard_name="heat_flux_onto_sea_water_due_to_snow_thermodynamics"                           unit="W/m2"     />
      <field id="hflx_ice_cea"  long_name="heat flux due to ice thermodynamics"                      standard_name="heat_flux_into_sea_water_due_to_sea_ice_thermodynamics"                        unit="W/m2"     />
      <field id="hflx_rnf_cea"  long_name="heat flux due to runoffs"                                 standard_name="temperature_flux_due_to_runoff_expressed_as_heat_flux_into_sea_water"          unit="W/m2"     />
      <field id="sflx_rnf_cea"  long_name="salt flux due to runoffs"                                 standard_name="salt_flux_due_to_runoffs"                                                      unit="g/m2/s"   />
      <field id="hflx_cal_cea"  long_name="heat flux due to calving"                                 standard_name="heat_flux_into_sea_water_due_to_calving"                                       unit="W/m2"     />
      <field id="hflx_icb_cea"  long_name="heat flux due to iceberg"                                 standard_name="heat_flux_into_sea_water_due_to_icebergs"                                      unit="W/m2"     />
      <field id="hflx_isf_cea"  long_name="heat flux due to iceshelf"                                standard_name="heat_flux_into_sea_water_due_to_iceshelf"                                      unit="W/m2"     />
      <field id="bicemel_cea"   long_name="Rate of Melt at Sea Ice Base (cell average)"              standard_name="tendency_of_sea_ice_amount_due_to_basal_melting"                               unit="kg/m2/s"  />
      <field id="licepro_cea"   long_name="Lateral Sea Ice Growth Rate (cell average)"               standard_name="tendency_of_sea_ice_amount_due_to_lateral_growth_of_ice_floes"                 unit="kg/m2/s"  />
      <field id="snowmel_cea"   long_name="Snow Melt Rate (cell average)"                            standard_name="surface_snow_melt_flux"                                                        unit="kg/m2/s"  />
      <field id="sntoice_cea"   long_name="Snow-Ice Formation Rate (cell average)"                   standard_name="tendency_of_sea_ice_amount_due_to_snow_conversion"                             unit="kg/m2/s"  />
      <field id="ticemel_cea"   long_name="Rate of Melt at Upper Surface of Sea Ice (cell average)"  standard_name="tendency_of_sea_ice_amount_due_to_surface_melting"                             unit="kg/m2/s"  />

      <!-- fluxes from relaxation and freshwater budget -->
      <field id="sflx_ssr_cea"  long_name="salt flux due to restoring"    standard_name="salt_flux_due_to_restoring"    unit="g/m2/s"   />
      <field id="hflx_ssr_cea"  long_name="heat flux due to restoring"    standard_name="heat_flux_due_to_restoring"    unit="W/m2"     />
      <field id="vflx_ssr_cea"  long_name="volume flux due to restoring"  standard_name="volume_flux_due_to_restoring"  unit="kg/m2/s"  />
      <field id="hflx_fwb_cea"  long_name="heat flux due to fwb"          standard_name="heat_flux_due_to_fwb"          unit="W/m2"     />
      <field id="vflx_fwb_cea"  long_name="volume flux due to fwb"        standard_name="volume_flux_due_to_fwb"        unit="kg/m2/s"  />
      
      <!-- ice field (nn_ice=1)  -->
      <field id="ice_cover"    long_name="Ice fraction"                                                 standard_name="sea_ice_area_fraction"                              unit="1"            />

      <!-- dilution -->
      <field id="emp_x_sst"    long_name="Concentration/Dilution term on SST"                                                                                              unit="kg*degC/m2/s" />
      <field id="emp_x_sss"    long_name="Concentration/Dilution term on SSS"                                                                                              unit="kg*1e-3/m2/s" />
      <field id="rnf_x_sst"    long_name="Runoff term on SST"                                                                                                              unit="kg*degC/m2/s" />
      <field id="rnf_x_sss"    long_name="Runoff term on SSS"                                                                                                              unit="kg*1e-3/m2/s" />

      <!-- sbcssm variables -->
      <field id="sst_m"    unit="degC" />
      <field id="sss_m"    unit="psu"  />

      <field id="sst_m_pot"    unit="degC" /> 
      <!-- EOS-80 -->
      <field id="sss_m_pra"    unit="psu"  />
      <!-- TEOS-10 -->
      <field id="sss_m_abs"    unit="g/kg"  />
      <!-- SEOS -->
      <field id="sss_m_seos"   unit="psu"  />

Guillaume Samson's avatar
Guillaume Samson committed
      <field id="ssu_m"    unit="m/s"  />
      <field id="ssv_m"    unit="m/s"  />
      <field id="ssh_m"    unit="m"    />
      <field id="e3t_m"    unit="m"    />
      <field id="frq_m"    unit="-"    />

    </field_group>


  </field_group> <!-- SBC -->

  <!-- ABL -->
  <field_group id="ABL" > <!-- time step automaticaly defined based on nn_fsbc -->

    <!-- variables available with ABL on atmospheric T grid-->
    <field_group id="grid_ABL3D" grid_ref="grid_TA_3D" >
      <field id="u_abl"      long_name="ABL i-horizontal velocity"     standard_name="abl_x_velocity" unit="m/s"      />
      <field id="v_abl"      long_name="ABL j-horizontal velocity"     standard_name="abl_y_velocity" unit="m/s"      />
      <field id="t_abl"      long_name="ABL potential temperature"     standard_name="abl_theta"      unit="K"        />
      <field id="q_abl"      long_name="ABL specific humidity"         standard_name="abl_qspe"       unit="kg/kg"    />
      <!-- debug (to be removed) -->
      <field id="u_dta"      long_name="DTA i-horizontal velocity"     standard_name="dta_x_velocity" unit="m/s"      />
      <field id="v_dta"      long_name="DTA j-horizontal velocity"     standard_name="dta_y_velocity" unit="m/s"      />
      <field id="t_dta"      long_name="DTA potential temperature"     standard_name="dta_theta"      unit="K"        />
      <field id="q_dta"      long_name="DTA specific humidity"         standard_name="dta_qspe"       unit="kg/kg"    />
      <field id="u_geo"      long_name="GEO i-horizontal velocity"     standard_name="geo_x_velocity" unit="m/s"      />
      <field id="v_geo"      long_name="GEO j-horizontal velocity"     standard_name="geo_y_velocity" unit="m/s"      />
      <field id="tke_abl"    long_name="ABL turbulent kinetic energy"  standard_name="abl_tke"        unit="m2/s2"    />
      <field id="avm_abl"    long_name="ABL turbulent viscosity"       standard_name="abl_avm"        unit="m2/s"     />
      <field id="avt_abl"    long_name="ABL turbulent diffusivity"     standard_name="abl_avt"        unit="m2/s"     />
      <field id="mxlm_abl"   long_name="ABL master mixing length"      standard_name="abl_mxlm"       unit="m"        />
      <field id="mxld_abl"   long_name="ABL dissipative mixing length" standard_name="abl_mxld"       unit="m"        />
    </field_group>

    <field_group id="grid_ABL2D" grid_ref="grid_TA_2D" >
      <field id="pblh"       long_name="ABL height"                    standard_name="abl_height"     unit="m"        />
      <field id="uz1_abl"    long_name="ABL i-horizontal velocity"     standard_name="abl_x_velocity" unit="m/s"      />
      <field id="vz1_abl"    long_name="ABL j-horizontal velocity"     standard_name="abl_y_velocity" unit="m/s"      />
      <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>
      <field id="tz1_abl"    long_name="ABL potential temperature"     standard_name="abl_theta"      unit="K"        />
      <field id="qz1_abl"    long_name="ABL specific humidity"         standard_name="abl_qspe"       unit="kg/kg"    />
      <field id="uz1_dta"    long_name="DTA i-horizontal velocity"     standard_name="dta_x_velocity" unit="m/s"      />
      <field id="vz1_dta"    long_name="DTA j-horizontal velocity"     standard_name="dta_y_velocity" unit="m/s"      />
      <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>
      <field id="tz1_dta"    long_name="DTA potential temperature"     standard_name="dta_theta"      unit="K"        />
      <field id="qz1_dta"    long_name="DTA specific humidity"         standard_name="dta_qspe"       unit="kg/kg"    />
      <!-- debug (to be removed) -->
      <field id="uz1_geo"    long_name="GEO i-horizontal velocity"     standard_name="geo_x_velocity" unit="m/s"      />
      <field id="vz1_geo"    long_name="GEO j-horizontal velocity"     standard_name="geo_y_velocity" unit="m/s"      />
      <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>
    </field_group>

  </field_group> <!-- ABL -->


  <!-- U grid -->

  <field_group id="grid_U"   grid_ref="grid_U_2D">
    <field id="hu"            long_name="water column height at U point"                         standard_name="water_column_height_U"       unit="m" />
    <field id="e2u"           long_name="U-cell width in meridional direction"                   standard_name="cell_width"                  unit="m"                               />
    <field id="e3u"           long_name="U-cell thickness"                                       standard_name="cell_thickness"              unit="m"          grid_ref="grid_U_3D" />
    <field id="e3u_0"         long_name="Initial U-cell thickness"                               standard_name="ref_cell_thickness"          unit="m"          grid_ref="grid_U_3D"/>
    <field id="utau"          long_name="Wind Stress along i-axis"                               standard_name="surface_downward_x_stress"   unit="N/m2"                            />
    <field id="uoce"          long_name="ocean current along i-axis"                             standard_name="sea_water_x_velocity"        unit="m/s"        grid_ref="grid_U_3D" />
    <field id="uoce_e3u"      long_name="ocean current along i-axis  (thickness weighted)"                                                   unit="m/s"        grid_ref="grid_U_3D"  > uoce * e3u </field>
    <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"/>
    <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>

    <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>
    <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" />
    <field id="msftbarot"         long_name="ocean_barotropic_mass_streamfunction"   unit="kg s-1" > uocetr_vsum_cumul * $rho0 </field>
    
    <field id="uoce2_e3u"         long_name="ocean current along i-axis squared (thickness weighted)"                                       unit="m3/s2"      grid_ref="grid_U_3D"  > uoce * uoce * e3u </field>
Guillaume Samson's avatar
Guillaume Samson committed
    <field id="ssu"          long_name="ocean surface current along i-axis"                                                                 unit="m/s"                             />
    <field id="sbu"          long_name="ocean bottom current along i-axis"                                                                  unit="m/s"  grid_ref="grid_U_2D_inner" />
    <field id="ubar"         long_name="ocean barotropic current along i-axis"                                                              unit="m/s"                             />
    <field id="uocetr_eff"   long_name="Effective ocean transport along i-axis"                 standard_name="ocean_volume_x_transport"    unit="m3/s"       grid_ref="grid_U_3D" />
    <field id="uocet"        long_name="ocean transport along i-axis times temperature (CRS)"                                               unit="degC*m/s"   grid_ref="grid_U_3D" />
    <field id="uoces"        long_name="ocean transport along i-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_U_3D" />
    <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>
    <!-- AGRIF sponge -->
    <field id="agrif_spu"    long_name=" AGRIF u-sponge coefficient"   unit=" " />
    <!-- u-eddy diffusivity coefficients (available if ln_traldf_OFF=F) -->
    <field id="ahtu_2d"      long_name=" surface u-eddy diffusivity coefficient"   unit="m2/s or m4/s" />
    <field id="ahtu_3d"      long_name=" 3D u-EIV coefficient"                     unit="m2/s or m4/s"      grid_ref="grid_U_3D"/>
    <!-- u-eiv diffusivity coefficients (available if ln_ldfeiv=F) -->
    <field id="aeiu_2d"      long_name=" surface u-EIV coefficient"                unit="m2/s" />
    <field id="aeiu_3d"      long_name=" 3D u-EIV coefficient"                     unit="m2/s"              grid_ref="grid_U_3D"/>

    <!-- variables available with MLE (ln_mle=T) -->
    <field id="psiu_mle"     long_name="MLE streamfunction along i-axis"   unit="m3/s"   grid_ref="grid_U_3D" />

    <!-- uoce_eiv: available EIV (ln_ldfeiv=T and ln_ldfeiv_dia=T) -->
    <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" />
    <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" />
    <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"                         />
    <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"                        />
    <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" />
    <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" />

    <!-- uoce_bbl: available with ln_trabbl=T and nn_bbl_adv=1 -->
    <field id="uoce_bbl"     long_name="BBL ocean current along i-axis"    unit="m/s"  />
    <!-- ahu_bbl : available with ln_trabbl=T and nn_bbl_ldf=1 -->
    <field id="ahu_bbl"      long_name="BBL diffusive flux along i-axis"   unit="m3/s" />

    <!-- variables available with WAVE (ln_wave=T) -->
    <field id="ustokes"      long_name="Stokes Drift Velocity i-axis"      standard_name="StokesDrift_x_velocity"      unit="m/s"        grid_ref="grid_U_3D" />
    <field id="ustokes_e3u"  long_name="Stokes Drift Velocity i-axis  (thickness weighted)"                            unit="m/s"        grid_ref="grid_U_3D"  > ustokes * e3u </field>

    <!-- variable for ice shelves -->
    <field id="utbl"         long_name="zonal current in the Losh tbl"     unit="m/s" />

    <!-- variables available with diaar5 -->
    <field id="u_masstr"      long_name="Ocean Mass X Transport"                                          standard_name="ocean_mass_x_transport"                         unit="kg/s"   grid_ref="grid_U_3D" />
    <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"   grid_ref="grid_U_2D_inner" />
    <field id="u_heattr"      long_name="ocean eulerian heat transport along i-axis"                      standard_name="ocean_heat_x_transport"                         unit="W"      grid_ref="grid_U_2D_inner" />
    <field id="u_salttr"      long_name="ocean eulerian salt transport along i-axis"                      standard_name="ocean_salt_x_transport"                         unit="1e-3*kg/s" grid_ref="grid_U_2D_inner" />
    <field id="uadv_heattr"   long_name="ocean advective heat transport along i-axis"                     standard_name="advectice_ocean_heat_x_transport"               unit="W"                         />
    <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"                 />
    <field id="udiff_heattr"  long_name="ocean diffusion heat transport along i-axis"                     standard_name="ocean_heat_x_transport_due_to_diffusion"        unit="W"                         />
    <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"                 />
  </field_group>

  <!-- V grid -->

  <field_group id="grid_V"   grid_ref="grid_V_2D">
    <field id="e1v"          long_name="V-cell width in longitudinal direction"                 standard_name="cell_width"                  unit="m"                              />
    <field id="e3v"          long_name="V-cell thickness"                                       standard_name="cell_thickness"              unit="m"          grid_ref="grid_V_3D" />
    <field id="e3v_0"        long_name="Initial V-cell thickness"                               standard_name="ref_cell_thickness"          unit="m"          grid_ref="grid_V_3D" />
    <field id="hv"            long_name="water column height at V point"                        standard_name="water_column_height_V"       unit="m" />
    <field id="vtau"         long_name="Wind Stress along j-axis"                               standard_name="surface_downward_y_stress"   unit="N/m2"                            />
    <field id="voce"         long_name="ocean current along j-axis"                             standard_name="sea_water_y_velocity"        unit="m/s"        grid_ref="grid_V_3D" />
    <field id="voce_e3v"     long_name="ocean current along j-axis  (thickness weighted)"                                                   unit="m/s"        grid_ref="grid_V_3D"  > voce * e3v </field>

    <field id="voce2_e3v"    long_name="ocean current along j-axis squared (thickness weighted)"                                            unit="m3/s2"      grid_ref="grid_V_3D"  > voce * voce * e3v </field>
Guillaume Samson's avatar
Guillaume Samson committed
    <field id="ssv"          long_name="ocean surface current along j-axis"                                                                 unit="m/s"                             />
    <field id="sbv"          long_name="ocean bottom current along j-axis"                                                                  unit="m/s"  grid_ref="grid_V_2D_inner" />
    <field id="vbar"         long_name="ocean barotropic current along j-axis"                                                              unit="m/s"                             />
    <field id="vocetr_eff"   long_name="Effective ocean transport along j-axis"                 standard_name="ocean_volume_y_transport"    unit="m3/s"       grid_ref="grid_V_3D" />
    <field id="vocet"        long_name="ocean transport along j-axis times temperature (CRS)"                                               unit="degC*m/s"   grid_ref="grid_V_3D" />
    <field id="voces"        long_name="ocean transport along j-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_V_3D" />
    <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>
    <!-- AGRIF sponge -->
    <field id="agrif_spv"    long_name=" AGRIF v-sponge coefficient"   unit=" " />
    <!-- v-eddy diffusivity coefficients (available if ln_traldf_OFF=F) -->
    <field id="ahtv_2d"      long_name=" surface v-eddy diffusivity coefficient"     unit="m2/s or (m4/s)^1/2" />
    <field id="ahtv_3d"      long_name=" 3D v-eddy diffusivity coefficient"          unit="m2/s or (m4/s)^1/2"           grid_ref="grid_V_3D"/>
    <!-- v-eiv diffusivity coefficients (available if ln_ldfeiv=F) -->
    <field id="aeiv_2d"      long_name=" surface v-EIV coefficient"                  unit="m2/s" />
    <field id="aeiv_3d"      long_name=" 3D v-EIV coefficient"                       unit="m2/s"                         grid_ref="grid_V_3D" />

    <!-- variables available with MLE (ln_mle=T) -->
    <field id="psiv_mle"     long_name="MLE streamfunction along j-axis"   unit="m3/s"   grid_ref="grid_V_3D" />

    <!-- voce_eiv: available EIV (ln_ldfeiv=T and ln_ldfeiv_dia=T)  -->
    <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" />
    <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" />
    <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"                         />
    <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"                        />
    <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" />
    <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" />


    <!-- voce_bbl: available with ln_trabbl=T and nn_bbl_adv=1 -->
    <field id="voce_bbl"     long_name="BBL ocean current along j-axis"    unit="m/s"  />
    <!-- ahu_bbl : available with ln_trabbl=T and nn_bbl_ldf=1 -->
    <field id="ahv_bbl"      long_name="BBL diffusive flux along j-axis"   unit="m3/s" />

    <!-- variables available with WAVE (ln_wave=T) -->
    <field id="vstokes"      long_name="Stokes Drift Velocity j-axis"      standard_name="StokesDrift_y_velocity"      unit="m/s"        grid_ref="grid_V_3D" />
    <field id="vstokes_e3v"  long_name="Stokes Drift Velocity j-axis  (thickness weighted)"                            unit="m/s"        grid_ref="grid_V_3D"  > vstokes * e3v </field>

    <!-- variable for ice shelves -->
    <field id="vtbl"         long_name="meridional current in the Losh tbl"   unit="m/s" />

    <!-- variables available with diaar5 -->
    <field id="v_masstr"      long_name="ocean eulerian mass transport along j-axis"    standard_name="ocean_mass_y_transport"                          unit="kg/s" grid_ref="grid_V_3D" />
    <field id="v_heattr"      long_name="ocean eulerian heat transport along j-axis"    standard_name="ocean_heat_y_transport"                          unit="W"          grid_ref="grid_V_2D_inner" />
    <field id="v_salttr"      long_name="ocean eulerian salt transport along i-axis"    standard_name="ocean_salt_y_transport"                          unit="1e-3*kg/s"  grid_ref="grid_V_2D_inner" />
    <field id="vadv_heattr"   long_name="ocean advective heat transport along j-axis"   standard_name="advectice_ocean_heat_y_transport"                unit="W"                         />
    <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"                 />
    <field id="vdiff_heattr"  long_name="ocean diffusion heat transport along j-axis"   standard_name="ocean_heat_y_transport_due_to_diffusion"         unit="W"                         />
    <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"                 />
  </field_group>

  <!-- W grid -->

  <field_group id="grid_W" grid_ref="grid_W_3D">
    <field id="e3w"          long_name="W-cell thickness"                              standard_name="cell_thickness"                         unit="m"    />
    <field id="woce"         long_name="ocean vertical velocity"                       standard_name="upward_sea_water_velocity"              unit="m/s"  />
    <field id="woce_e3w"     long_name="ocean vertical velocity * e3w"                                                                        unit="m2/s"  > woce * e3w </field>
    <field id="wocetr_eff"   long_name="effective ocean vertical transport"                                                                   unit="m3/s" />

    <!-- woce_eiv: available with EIV  (ln_ldfeiv=T and ln_ldfeiv_dia=T)  -->
    <field id="woce_eiv"     long_name="EIV ocean vertical velocity"                    standard_name="bolus_upward_sea_water_velocity"       unit="m/s"  />
    <field id="weiv_masstr"  long_name="EIV Upward Ocean Mass Transport"  standard_name="bolus_upward_ocean_mass_transport"             unit="kg/s"   />
    <field id="weiv_heattr3d" long_name="ocean bolus heat transport"    standard_name="ocean_heat_z_transport_due_to_bolus_advection"   unit="W"    />
    <field id="weiv_salttr3d" long_name="ocean bolus salt transport"    standard_name="ocean_salt_z_transport_due_to_bolus_advection"   unit="kg"   />

    <field id="avt"          long_name="vertical eddy diffusivity"                      standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" />
    <field id="avt_e3w"      long_name="vertical heat diffusivity * e3w"                unit="m3/s" > avt * e3w </field>
    <field id="logavt"       long_name="logarithm of vertical eddy diffusivity"         standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" />
    <field id="avm"          long_name="vertical eddy viscosity"                        standard_name="ocean_vertical_momentum_diffusivity"   unit="m2/s" />
    <field id="avm_e3w"      long_name="vertical eddy viscosity * e3w"   unit="m3/s" > avm * e3w </field>

    <!-- avs: /= avt with ln_zdfddm=T -->
    <field id="avs"          long_name="salt vertical eddy diffusivity"                 standard_name="ocean_vertical_salt_diffusivity"       unit="m2/s" />
    <field id="avs_e3w"      long_name="vertical salt diffusivity * e3w"   unit="m3/s" > avs * e3w </field>
    <field id="logavs"       long_name="logarithm of salt vertical eddy diffusivity"    standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" />

    <!-- avt_evd and avm_evd: available with ln_zdfevd -->
    <field id="avt_evd"      long_name="convective enhancement of vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_convection"     unit="m2/s" />
    <field id="avt_evd_e3w"  long_name="convective enhancement to vertical diffusivity * e3w "    unit="m3/s" > avt_evd * e3w </field>
    <field id="avm_evd"      long_name="convective enhancement of vertical viscosity"   standard_name="ocean_vertical_momentum_diffusivity_due_to_convection"   unit="m2/s" />

    <!-- mf_app and mf_wp: available with ln_zdfmfc -->
    <field id="mf_app"      long_name="convective area"        standard_name="mf_convective_area"    unit="%"      grid_ref="grid_W_3D" />
    <field id="mf_wp"       long_name="convective velocity"    standard_name="mf_convective_velo"    unit="m/s"    grid_ref="grid_W_3D" />

    <!-- av_tmx: available with ln_zdftmx -->
    <field id="av_tmx"      long_name="vertical diffusivity due to tidal mixing" unit="m2/s" />
Guillaume Samson's avatar
Guillaume Samson committed

    <!-- avt_tide: available with ln_zdfiwm=T -->
    <field id="av_ratio"     long_name="S over T diffusivity ratio"                     standard_name="salinity_over_temperature_diffusivity_ratio"                     unit="1"    />
    <field id="av_wave"      long_name="internal wave-induced vertical diffusivity"     standard_name="ocean_vertical_tracer_diffusivity_due_to_internal_waves"         unit="m2/s" />
    <field id="bflx_iwm"     long_name="internal wave-induced buoyancy flux"            standard_name="buoyancy_flux_due_to_internal_waves"                             unit="W/kg" />
    <field id="pcmap_iwm"    long_name="power consumed by wave-driven mixing"           standard_name="vertically_integrated_power_consumption_by_wave_driven_mixing"   unit="W/m2"      grid_ref="grid_W_2D" />
    <field id="emix_iwm"     long_name="power density available for mixing"             standard_name="power_available_for_mixing_from_breaking_internal_waves"         unit="W/kg" />

    <!-- variables available with WAVE (ln_wave=T) -->
    <field id="wstokes"      long_name="Stokes Drift vertical velocity"                 standard_name="upward_StokesDrift_velocity"   unit="m/s" />

    <!-- variables available with diaar5 -->
    <field id="w_masstr"     long_name="vertical mass transport"                        standard_name="upward_ocean_mass_transport"             unit="kg/s"   />
    <field id="w_masstr2"    long_name="square of vertical mass transport"              standard_name="square_of_upward_ocean_mass_transport"   unit="kg2/s2" />

    <!-- EOS -->
    <field id="bn2"          long_name="squared Brunt-Vaisala frequency"                unit="s-2" />

    <!-- dissipation diagnostics (note: ediss_k is only available with tke scheme) -->   
    <field id="avt_k"        long_name="vertical eddy diffusivity from closure schemes" standard_name="ocean_vertical_eddy_diffusivity"       unit="m2/s" />
    <field id="avm_k"        long_name="vertical eddy viscosity from closure schemes"   standard_name="ocean_vertical_eddy_viscosity"         unit="m2/s" />
    <field id="ediss_k"      long_name="Kolmogorov energy dissipation (tke scheme)"     standard_name="Kolmogorov_energy_dissipation"         unit="W/kg" />
    <field id="eshear_k"     long_name="energy source from vertical shear"              standard_name="energy_source_from_shear"                                                      unit="W/kg"   grid_ref="grid_W_3D_inner" />
Guillaume Samson's avatar
Guillaume Samson committed
    <field id="estrat_k"     long_name="energy sink from stratification"                standard_name="energy_sink_from_stratification"       unit="W/kg" />
    
  </field_group>

  <!-- F grid -->
  <field_group id="grid_F" grid_ref="grid_F_2D">
    <field id="e3f"          long_name="F-cell thickness"                      standard_name="cell_thickness"         unit="m"   grid_ref="grid_F_3D" />
    <field id="e3f_0"        long_name="F-cell thickness"                      standard_name="cell_thickness"         unit="m"   grid_ref="grid_F_3D" />
    <field id="hf"           long_name="water column height at F point"        standard_name="water_column_height_F"  unit="m"                     />
    <field id="ssKEf"        long_name="surface kinetic energy at F point"     standard_name="specific_kinetic_energy_of_sea_water"   unit="m2/s2" grid_ref="grid_F_2D_inner" />
    <field id="ssrelvor"     long_name="surface relative vorticity"            standard_name="relative_vorticity"     unit="1/s"       grid_ref="grid_F_2D_inner" />
    <field id="ssplavor"     long_name="surface planetary vorticity"           standard_name="planetary_vorticity"    unit="1/s"       />
    <field id="ssrelpotvor"  long_name="surface relative potential vorticity"  standard_name="relpot_vorticity"       unit="1/m.s"     grid_ref="grid_F_2D_inner" />
    <field id="ssabspotvor"  long_name="surface absolute potential vorticity"  standard_name="abspot_vorticity"       unit="1/m.s"     grid_ref="grid_F_2D_inner" />
    <field id="ssEns"        long_name="surface enstrophy"                     standard_name="enstrophy"              unit="1/m2.s2"   grid_ref="grid_F_2D_inner" />
  </field_group>

  <!-- AGRIF sponge -->
  <field id="agrif_spf"    long_name=" AGRIF f-sponge coefficient"   unit=" " />
  <!-- f-eddy viscosity coefficients (ldfdyn) -->
  <field id="ahmf_2d"      long_name=" surface f-eddy viscosity coefficient"   unit="m2/s or m4/s" />
  <field id="ahmf_3d"      long_name=" 3D      f-eddy viscosity coefficient"   unit="m2/s or m4/s"  grid_ref="grid_T_3D"/>
  
  
  <!-- product fields -->
  <field_group id="diaprod">
  
    <field id="ut"           long_name="product_of_sea_water_x_velocity_and_potential_temperature"      unit="degree_C m/s"      grid_ref="grid_U_3D"   />
    <field id="ut_e3u"       long_name="product_of_sea_water_x_velocity_and_potential_temperature * e3u"  unit="degree_C m2/s"   grid_ref="grid_U_3D" > ut * e3u </field >
    <field id="us"           long_name="product_of_sea_water_x_velocity_and_salinity"                   unit="PSU m/s"       grid_ref="grid_U_3D"   />
    <field id="us_e3u"       long_name="product_of_sea_water_x_velocity_and_salinity * e3u"             unit="PSU m2/s"      grid_ref="grid_U_3D" > us * e3u </field >
    <field id="urhop"        long_name="product_of_sea_water_x_velocity_and_potential_density"          unit="(kg/m3).(m/s)" grid_ref="grid_U_3D"   />
    <field id="urhop_e3u"    long_name="product_of_sea_water_x_velocity_and_potential_density * e3u"    unit="(kg/m3).(m2/s)"   grid_ref="grid_U_3D" > urhop * e3u </field >
    <field id="vt"           long_name="product_of_sea_water_y_velocity_and_potential_temperature"      unit="degree_C m/s"      grid_ref="grid_V_3D"   />
    <field id="vt_e3v"       long_name="product_of_sea_water_y_velocity_and_potential_temperature * e3v"  unit="degree_C m2/s"   grid_ref="grid_V_3D" > vt * e3v </field >
    <field id="vs"           long_name="product_of_sea_water_y_velocity_and_salinity"                   unit="PSU m/s"       grid_ref="grid_V_3D"   />
    <field id="vs_e3v"       long_name="product_of_sea_water_y_velocity_and_salinity * e3t"             unit="PSU m2/s"      grid_ref="grid_V_3D" > vs * e3v </field >
    <field id="vrhop"        long_name="product_of_sea_water_y_velocity_and_potential_density"          unit="(kg/m3).(m/s)" grid_ref="grid_V_3D"   />
    <field id="vrhop_e3v"    long_name="product_of_sea_water_y_velocity_and_potential_density * e3t"    unit="(kg/m3).(m2/s)"  grid_ref="grid_V_3D" > vrhop * e3v </field >
    <field id="wt"           long_name="product_of_upward_sea_water_velocity_and_potential_temperature" unit="degree_C m/s"      grid_ref="grid_W_3D"   />
    <field id="ws"           long_name="product_of_upward_sea_water_velocity_and_salinity"              unit="PSU m/s"       grid_ref="grid_W_3D"   />
    <field id="wrhop"        long_name="product_of_upward_sea_water_velocity_and_potential_density"     unit="(kg/m3).(m/s)" grid_ref="grid_W_3D"   />
    <field id="uv"           long_name="product_of_sea_water_x_velocity_and_sea_water_y_velocity"       unit="m2/s2   "      grid_ref="grid_T_3D"   />
    <field id="uw"           long_name="product_of_upward_sea_water_velocity_and_sea_water_x_velocity"  unit="m2/s2   "      grid_ref="grid_W_3D"   />
    <field id="vw"           long_name="product_of_upward_sea_water_velocity_and_sea_water_y_velocity"  unit="m2/s2"         grid_ref="grid_W_3D"   />
  
  </field_group>
Guillaume Samson's avatar
Guillaume Samson committed
908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936 937 938 939 940 941 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 973 974 975 976 977 978 979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 1000 1001 1002 1003 1004 1005 1006 1007 1008 1009 1010 1011 1012 1013 1014 1015 1016 1017 1018 1019 1020 1021 1022 1023 1024 1025 1026 1027 1028 1029 1030 1031 1032 1033 1034 1035 1036 1037 1038 1039 1040 1041 1042 1043 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062 1063 1064 1065 1066 1067 1068 1069 1070 1071 1072 1073 1074 1075 1076 1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089 1090 1091 1092 1093 1094 1095 1096 1097 1098 1099 1100 1101 1102 1103 1104 1105 1106 1107 1108 1109 1110 1111 1112 1113 1114 1115 1116 1117 1118 1119 1120 1121 1122 1123 1124 1125 1126 1127 1128 1129 1130 1131 1132 1133 1134 1135 1136 1137 1138 1139 1140 1141 1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184 1185 1186 1187 1188 1189 1190 1191 1192 1193 1194 1195 1196 1197 1198 1199 1200 1201 1202 1203 1204 1205 1206 1207 1208 1209 1210 1211 1212 1213 1214 1215 1216 1217 1218 1219 1220 1221 1222 1223 1224 1225 1226 1227 1228 1229 1230 1231 1232 1233 1234 1235 1236 1237 1238 1239 1240 1241 1242 1243 1244 1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 1255 1256 1257 1258 1259 1260 1261 1262 1263 1264 1265 1266 1267 1268 1269

  <field_group id="scalar"  grid_ref="grid_scalar"  >
    <!-- Need to have a salinity reference climatological file : sali_ref_clim_monthly -->
    <field id="voltot"     long_name="global total volume"                          standard_name="sea_water_volume"                               unit="m3"   />
    <field id="sshtot"     long_name="global mean ssh"                              standard_name="global_average_sea_level_change"                unit="m"    />
    <field id="sshsteric"  long_name="global mean ssh steric"                       standard_name="global_average_steric_sea_level_change"         unit="m"    />
    <field id="sshthster"  long_name="global mean ssh thermosteric"                 standard_name="global_average_thermosteric_sea_level_change"   unit="m"    />
    <field id="masstot"    long_name="global total mass"                            standard_name="sea_water_mass"                                 unit="kg"   />
    <field id="temptot"    long_name="global mean temperature"                      standard_name="sea_water_potential_temperature"                unit="degC" />
    <field id="saltot"     long_name="global mean salinity"                         standard_name="sea_water_salinity"                             unit="1e-3" />
    <field id="ssttot"     long_name="global mean sea surface temperature"          standard_name="sea_water_conservative_temperature"             unit="degC" />
    <!-- EOS = TEOS-10 -->
    <field id="temptot_pot" long_name="global mean potential temperature"            standard_name="sea_water_potential_temperature"               unit="degC" />

    <!-- available with ln_diahsb -->
    <field id="bgtemper"     long_name="drift in global mean temperature wrt timestep 1"                 standard_name="change_over_time_in_sea_water_potential_temperature"   unit="degC"     />
    <field id="bgsaline"     long_name="drift in global mean salinity wrt timestep 1"                    standard_name="change_over_time_in_sea_water_practical_salinity"      unit="1e-3"     />
    <field id="bgheatco"     long_name="drift in global mean heat content wrt timestep 1"                                                                                      unit="1.e20J"   />
    <field id="bgheatfx"     long_name="drift in global mean heat flux    wrt timestep 1"                                                                                      unit="W/m2"     />
    <field id="bgsaltco"     long_name="drift in global mean salt content wrt timestep 1"                                                                                      unit="1e-3*km3" />
    <field id="bgvolssh"     long_name="drift in global mean ssh volume wrt timestep 1"                                                                                        unit="km3"      />
    <field id="bgvole3t"     long_name="drift in global mean volume variation (e3t) wrt timestep 1"                                                                            unit="km3"      />
    <field id="bgfrcvol"     long_name="global mean volume from forcing"                                                                                                       unit="km3"      />
    <field id="bgfrctem"     long_name="global mean heat content from forcing"                                                                                                 unit="1.e20J"   />
    <field id="bgfrchfx"     long_name="global mean heat flux from forcing"                                                                                                    unit="W/m2"     />
    <field id="bgfrcsal"     long_name="global mean salt content from forcing"                                                                                                 unit="1e-3*km3" />
    <field id="bgmistem"     long_name="global mean temperature error due to free surface (linssh true)"                                                                       unit="degC"     />
    <field id="bgmissal"     long_name="global mean salinity error due to free surface (linssh true)"                                                                          unit="1e-3"     />
  </field_group>


  <!-- transects -->
  <field_group id="oce_straits">
    <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>
    <field id="uoce_e3u_ave_vsum"    long_name="Vertical sum of u*e3u"                           field_ref="uoce_e3u_ave"         grid_ref="grid_U_vsum"     />
    <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>
    <field id="uocetr_strait"        long_name="Total transport across lines in i-direction"     field_ref="uocetr_vsum_section"  grid_ref="grid_U_4strait" />
    <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 * $rho0 </field>

    <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>
    <field id="voce_e3v_ave_vsum"    long_name="Vertical sum of v*e3v"                           field_ref="voce_e3v_ave"         grid_ref="grid_V_vsum"      />
    <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>
    <field id="vocetr_strait"        long_name="Total transport across lines in j-direction"     field_ref="vocetr_vsum_section"  grid_ref="grid_V_4strait"  />
    <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 * $rho0 </field>

    <field id="masstr_strait"        long_name="Sea water transport across line"                                                  grid_ref="grid_4strait"  > u_masstr_strait + v_masstr_strait </field>
  </field_group>

  <!-- variables available with ln_floats -->

  <field_group id="floatvar" grid_ref="grid_T_nfloat"  operation="instant" >
    <field id="traj_lon"      long_name="floats longitude"                                                           unit="degrees_east"  />
    <field id="traj_lat"      long_name="floats latitude"                                                            unit="degrees_north" />
    <field id="traj_dep"      long_name="floats depth"                                                               unit="m"             />
    <field id="traj_temp"     long_name="floats temperature"       standard_name="sea_water_potential_temperature"   unit="degC"          />
    <field id="traj_salt"     long_name="floats salinity"          standard_name="sea_water_practical_salinity"      unit="1e-3"          />
    <field id="traj_dens"     long_name="floats in-situ density"   standard_name="sea_water_density"                 unit="kg/m3"         />
    <field id="traj_group"    long_name="floats group"                                                               unit="1"             />
  </field_group>

  <!-- variables available with iceberg trajectories -->

  <field_group id="icbvar" domain_ref="grid_T"  >
    <field id="berg_melt"          long_name="icb melt rate of icebergs"                       unit="kg/m2/s"                    />
    <field id="berg_melt_hcflx"    long_name="icb heat flux to ocean due to melting heat content"   unit="J/m2/s"                />
    <field id="berg_melt_qlat"     long_name="icb heat flux to ocean due to melting latent heat"    unit="J/m2/s"                />
    <field id="berg_buoy_melt"     long_name="icb buoyancy component of iceberg melt rate"     unit="kg/m2/s"                    />
    <field id="berg_eros_melt"     long_name="icb erosion component of iceberg melt rate"      unit="kg/m2/s"                    />
    <field id="berg_conv_melt"     long_name="icb convective component of iceberg melt rate"   unit="kg/m2/s"                    />
    <field id="berg_virtual_area"  long_name="icb virtual coverage by icebergs"                unit="m2"                         />
    <field id="bits_src"           long_name="icb mass source of bergy bits"                   unit="kg/m2/s"                    />
    <field id="bits_melt"          long_name="icb melt rate of bergy bits"                     unit="kg/m2/s"                    />
    <field id="bits_mass"          long_name="icb bergy bit density field"                     unit="kg/m2"                      />
    <field id="berg_mass"          long_name="icb iceberg density field"                       unit="kg/m2"                      />
    <field id="calving"            long_name="icb calving mass input"                          unit="kg/s"                       />
    <field id="berg_floating_melt" long_name="icb melt rate of icebergs + bits"                unit="kg/m2/s"                    />
    <field id="berg_real_calving"  long_name="icb calving into iceberg class"                  unit="kg/s"     axis_ref="icbcla" />
    <field id="berg_stored_ice"    long_name="icb accumulated ice mass by class"               unit="kg"       axis_ref="icbcla" />
  </field_group>

  <!-- Poleward transport : ptr -->
  <field_group id="diaptr" >
    <field id="zomsf"         long_name="Overturning Stream-Function : All basins"                     unit="Sv"         grid_ref="grid_znl_W_3D" />
    <field id="zotem"         long_name="Zonal Mean Temperature : All basins"                          unit="degree_C"   grid_ref="grid_znl_T_3D" />
    <field id="zosal"         long_name="Zonal Mean Salinity : All basins"                             unit="0.001"      grid_ref="grid_znl_T_3D" />
    <field id="zosrf"         long_name="Zonal Mean Surface : All basins"                              unit="m2"         grid_ref="grid_znl_T_3D" />
    <field id="sophtove"      long_name="Overturning Heat Transport: All basins"                       unit="PW"         grid_ref="grid_znl_T_2D" />
    <field id="sopstove"      long_name="Overturning Salt Transport: All basins"                       unit="Giga g/s"   grid_ref="grid_znl_T_2D" />
    <field id="sophtbtr"      long_name="Barotropic Heat Transport: All basins"                        unit="PW"         grid_ref="grid_znl_T_2D" />
    <field id="sopstbtr"      long_name="Barotropic Salt Transport: All basins"                        unit="Giga g/s"   grid_ref="grid_znl_T_2D" />
    <field id="sophtadv"      long_name="Advective Heat Transport: All basins"                         unit="PW"         grid_ref="grid_znl_T_2D" />
    <field id="sopstadv"      long_name="Advective Salt Transport: All basins"                         unit="Giga g/s"   grid_ref="grid_znl_T_2D" />
    <field id="sophtldf"      long_name="Diffusive Heat Transport: All basins"                         unit="PW"         grid_ref="grid_znl_T_2D" />
    <field id="sopstldf"      long_name="Diffusive Salt Transport: All basins"                         unit="Giga g/s"   grid_ref="grid_znl_T_2D" />
    <field id="sophtvtr"      long_name="Heat Transport : All basins"                                  unit="PW"         grid_ref="grid_znl_T_2D" />
    <field id="sopstvtr"      long_name="Salt Transport : All basins"                                  unit="Giga g/s"   grid_ref="grid_znl_T_2D" />
    <field id="sophteiv"      long_name="Heat Transport from mesoscale eddy advection: All basins"     unit="PW"         grid_ref="grid_znl_T_2D" />
    <field id="sopsteiv"      long_name="Salt Transport from mesoscale eddy advection : All basins"    unit="Giga g/s"   grid_ref="grid_znl_T_2D" />
    <field id="sopstadv"      long_name="Advective Salt Transport"                                     unit="Giga g/s"   grid_ref="grid_znl_T_2D" />
    <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>
    <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>
  </field_group>

  <field_group id="constant_fields"  grid_ref="grid_T_2D"  operation="once" >
    <field id="bathy"       long_name="Sea floor depth below geoid"               standard_name="sea_floor_depth_below_geoid"               unit="m"/>
    <field id="areacello"   long_name="Horizontal area of ocean grid cells"       standard_name="cell_area"                                 unit="m2" />
    <field id="hfgeou"      long_name="Upward geothermal heat flux at sea floor"  standard_name="upward_geothermal_heat_flux_at_sea_floor"  unit="W/m2"/>
    <field id="basins"      long_name="ocean tracer region masks"                 standard_name="ocean_tracer_region_masks"   unit="none" grid_ref="grid_basin"  />
  </field_group>


  <!--
      ============================================================================================================
      Physical ocean model trend diagnostics : temperature, KE, PE, momentum
      ============================================================================================================
  -->

  <!-- variables available with ln_tra_trd -->
  <!-- Asselin trends  calculated on odd time steps-->
  <field_group id="trendT_odd"  grid_ref="grid_T_3D">
    <field id="ttrd_atf"      long_name="temperature-trend: asselin time filter"       unit="degree_C/s" />
    <field id="strd_atf"      long_name="salinity   -trend: asselin time filter"       unit="0.001/s" />
    <!-- Thickness weighted versions: -->
    <field id="ttrd_atf_e3t"      unit="degC/s * m"  >  ttrd_atf * e3t </field>
    <field id="strd_atf_e3t"      unit="1e-3/s * m"  >  strd_atf * e3t </field>
    <!-- OMIP  layer-integrated trends -->
    <field id="ttrd_atf_li"      long_name="layer integrated heat-trend: asselin time filter "       unit="W/m^2" > ttrd_atf_e3t * 1026.0 * 3991.86795711963  </field>
    <field id="strd_atf_li"      long_name="layer integrated salt   -trend: asselin time filter "       unit="kg/(m^2 s)" > strd_atf_e3t * 1026.0 * 0.001 </field>
  </field_group>

  <!-- Other trends  calculated on even time steps-->
  <field_group id="trendT_even" grid_ref="grid_T_3D">
    <field id="ttrd_xad"      long_name="temperature-trend: i-advection"                                                                                          unit="degC/s"                        />
    <field id="strd_xad"      long_name="salinity   -trend: i-advection"                                                                                          unit="1e-3/s"                        />
    <field id="ttrd_yad"      long_name="temperature-trend: j-advection"                                                                                          unit="degC/s"                        />
    <field id="strd_yad"      long_name="salinity   -trend: j-advection"                                                                                          unit="1e-3/s"                        />
    <field id="ttrd_zad"      long_name="temperature-trend: k-advection"                                                                                          unit="degC/s"                        />
    <field id="strd_zad"      long_name="salinity   -trend: k-advection"                                                                                          unit="1e-3/s"                        />
    <field id="ttrd_ad"       long_name="temperature-trend: advection"               standard_name="tendency_of_sea_water_temperature_due_to_advection"           unit="degC/s"                         > sqrt( ttrd_xad^2 + ttrd_yad^2 + ttrd_zad^2 ) </field>
    <field id="strd_ad"       long_name="salinity   -trend: advection"               standard_name="tendency_of_sea_water_salinity_due_to_advection"              unit="1e-3/s"                         > sqrt( strd_xad^2 + strd_yad^2 + strd_zad^2 ) </field>
    <field id="ttrd_totad"    long_name="temperature-trend: total advection"         standard_name="tendency_of_sea_water_salinity_due_to_advection"              unit="degC/s"                        />
    <field id="strd_totad"    long_name="salinity   -trend: total advection"         standard_name="tendency_of_sea_water_salinity_due_to_advection"              unit="1e-3/s"                        />
    <field id="ttrd_sad"      long_name="temperature-trend: surface adv. (linssh true)"                                                                           unit="degC/s"   grid_ref="grid_T_2D" />
    <field id="strd_sad"      long_name="salinity   -trend: surface adv. (linssh true)"                                                                           unit="1e-3/s"   grid_ref="grid_T_2D" />
    <field id="ttrd_ldf"      long_name="temperature-trend: lateral  diffusion"      standard_name="tendency_of_sea_water_temperature_due_to_horizontal_mixing"   unit="degC/s"                        />
    <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"                        />
    <field id="ttrd_zdf"      long_name="temperature-trend: vertical diffusion"      standard_name="tendency_of_sea_water_temperature_due_to_vertical_mixing"     unit="degC/s"                        />
    <field id="strd_zdf"      long_name="salinity   -trend: vertical diffusion"      standard_name="tendency_of_sea_water_salinity_due_to_vertical_mixing"        unit="1e-3/s"                        />
    <field id="ttrd_evd"      long_name="temperature-trend: EVD convection"                                                                                       unit="degC/s"                        />
    <field id="strd_evd"      long_name="salinity   -trend: EVD convection"                                                                                       unit="1e-3/s"                        />

    <!-- ln_traldf_iso=T only (iso-neutral diffusion) -->
    <field id="ttrd_iso"      long_name="temperature-trend: isopycnal diffusion"                             unit="degC/s" > ttrd_ldf + ttrd_zdf - ttrd_zdfp </field>
    <field id="strd_iso"      long_name="salinity   -trend: isopycnal diffusion"                             unit="1e-3/s" > strd_ldf + strd_zdf - strd_zdfp </field>
    <field id="ttrd_zdfp"     long_name="temperature-trend: pure vert. diffusion"   unit="degC/s" />
    <field id="strd_zdfp"     long_name="salinity   -trend: pure vert. diffusion"   unit="1e-3/s" />

    <!-- ln_zdfosm=T only (OSMOSIS-OBL) -->
    <field id="ttrd_osm"      long_name="temperature-trend: OSM-OSBL non-local forcing"                             unit="degC/s" />
    <field id="strd_osm"      long_name="salinity   -trend: OSM-OSBL non-local forcing"                             unit="1e-3/s" />


    <!-- -->
    <field id="ttrd_dmp"      long_name="temperature-trend: interior restoring"        unit="degC/s" />
    <field id="strd_dmp"      long_name="salinity   -trend: interior restoring"        unit="1e-3/s" />
    <field id="ttrd_bbl"      long_name="temperature-trend: bottom boundary layer"     unit="degC/s" />
    <field id="strd_bbl"      long_name="salinity   -trend: bottom boundary layer"     unit="1e-3/s" />
    <field id="ttrd_npc"      long_name="temperature-trend: non-penetrative conv."     unit="degC/s" />
    <field id="strd_npc"      long_name="salinity   -trend: non-penetrative conv."     unit="1e-3/s" />
    <field id="ttrd_qns"      long_name="temperature-trend: non-solar flux + runoff"   unit="degC/s" grid_ref="grid_T_2D" />
    <field id="strd_cdt"      long_name="salinity   -trend: C/D term       + runoff"   unit="degC/s" grid_ref="grid_T_2D" />
    <field id="ttrd_qsr"      long_name="temperature-trend: solar penetr. heating"     unit="degC/s" />
    <field id="ttrd_bbc"      long_name="temperature-trend: geothermal heating"        unit="degC/s" />

    <!-- Thickness weighted versions: -->
    <field id="ttrd_xad_e3t"      unit="degC/s * m" >  ttrd_xad * e3t </field>
    <field id="strd_xad_e3t"      unit="1e-3/s * m" >  strd_xad * e3t </field>
    <field id="ttrd_yad_e3t"      unit="degC/s * m" >  ttrd_yad * e3t </field>
    <field id="strd_yad_e3t"      unit="1e-3/s * m" >  strd_yad * e3t </field>
    <field id="ttrd_zad_e3t"      unit="degC/s * m" >  ttrd_zad * e3t </field>
    <field id="strd_zad_e3t"      unit="1e-3/s * m" >  strd_zad * e3t </field>
    <field id="ttrd_ad_e3t"       unit="degC/s * m" >  ttrd_ad  * e3t </field>
    <field id="strd_ad_e3t"       unit="1e-3/s * m" >  strd_ad  * e3t </field>
    <field id="ttrd_totad_e3t"    unit="degC/s * m" >  ttrd_totad  * e3t </field>
    <field id="strd_totad_e3t"    unit="1e-3/s * m" >  strd_totad  * e3t </field>
    <field id="ttrd_ldf_e3t"      unit="degC/s * m" >  ttrd_ldf * e3t </field>
    <field id="strd_ldf_e3t"      unit="1e-3/s * m" >  strd_ldf * e3t </field>
    <field id="ttrd_zdf_e3t"      unit="degC/s * m" >  ttrd_zdf * e3t </field>
    <field id="strd_zdf_e3t"      unit="1e-3/s * m" >  strd_zdf * e3t </field>
    <field id="ttrd_evd_e3t"      unit="degC/s * m" >  ttrd_evd * e3t </field>
    <field id="strd_evd_e3t"      unit="1e-3/s * m" >  strd_evd * e3t </field>

    <!-- ln_traldf_iso=T only (iso-neutral diffusion) -->
    <field id="ttrd_iso_e3t"      unit="degC/s * m"  >  ttrd_iso * e3t </field>
    <field id="strd_iso_e3t"      unit="1e-3/s * m"  >  strd_iso * e3t </field>
    <field id="ttrd_zdfp_e3t"     unit="degC/s * m"  >  ttrd_zdfp * e3t </field>
    <field id="strd_zdfp_e3t"     unit="1e-3/s * m"  >  strd_zdfp * e3t </field>

    <!-- ln_zdfosm=T only (OSMOSIS-OBL) -->
    <field id="ttrd_osm_e3t"      long_name="temperature-trend: OSM-OSBL non-local forcing"                             unit="degC/s * m" >  ttrd_osm * e3t </field>
    <field id="strd_osm_e3t"      long_name="salinity   -trend: OSM-OSBL non-local forcing"                             unit="1e-3/s * m" >  strd_osm * e3t </field>

    <!-- -->
    <field id="ttrd_dmp_e3t"      unit="degC/s * m"  >  ttrd_dmp * e3t </field>
    <field id="strd_dmp_e3t"      unit="1e-3/s * m"  >  strd_dmp * e3t </field>
    <field id="ttrd_bbl_e3t"      unit="degC/s * m"  >  ttrd_bbl * e3t </field>
    <field id="strd_bbl_e3t"      unit="1e-3/s * m"  >  strd_bbl * e3t </field>
    <field id="ttrd_npc_e3t"      unit="degC/s * m"  >  ttrd_npc * e3t </field>
    <field id="strd_npc_e3t"      unit="1e-3/s * m"  >  strd_npc * e3t </field>
    <field id="ttrd_qns_e3t"      unit="degC/s * m"  >  ttrd_qns * e3ts </field>
    <field id="strd_cdt_e3t"      unit="degC/s * m"  >  strd_cdt * e3ts </field>
    <field id="ttrd_qsr_e3t"      unit="degC/s * m"  >  ttrd_qsr * e3t </field>
    <field id="ttrd_bbc_e3t"      unit="degC/s * m"  >  ttrd_bbc * e3t </field>

    <!-- OMIP  layer-integrated trends -->
    <field id="ttrd_totad_li"    long_name="layer integrated heat-trend: total advection"         unit="W/m^2"     > ttrd_totad_e3t * 1026.0 * 3991.86795711963 </field>
    <field id="strd_totad_li"    long_name="layer integrated salt-trend: total advection"         unit="kg/(m^2 s)"    > strd_totad_e3t * 1026.0 * 0.001  </field>
    <field id="ttrd_osm_li"    long_name="layer integrated heat-trend: non-local OSM"         unit="W/m^2"     > ttrd_osm_e3t * 1026.0 * 3991.86795711963 </field>
    <field id="strd_osm_li"    long_name="layer integrated salt-trend: non-local OSM"         unit="kg/(m^2 s)"    > strd_osm_e3t * 1026.0 * 0.001  </field>
    <field id="ttrd_evd_li"      long_name="layer integrated heat-trend: EVD convection"          unit="W/m^2"    > ttrd_evd_e3t * 1026.0 * 3991.86795711963 </field>
    <field id="strd_evd_li"      long_name="layer integrated salt-trend: EVD convection"          unit="kg/(m^2 s)"  > strd_evd_e3t * 1026.0 * 0.001  </field>
    <field id="ttrd_iso_li"      long_name="layer integrated heat-trend: isopycnal diffusion"     unit="W/m^2" > ttrd_iso_e3t * 1026.0 * 3991.86795711963 </field>
    <field id="strd_iso_li"      long_name="layer integrated salt-trend: isopycnal diffusion"     unit="kg/(m^2 s)" > strd_iso_e3t * 1026.0 * 0.001  </field>
    <field id="ttrd_zdfp_li"     long_name="layer integrated heat-trend: pure vert. diffusion"    unit="W/m^2" > ttrd_zdfp_e3t * 1026.0 * 3991.86795711963 </field>
    <field id="strd_zdfp_li"     long_name="layer integrated salt-trend: pure vert. diffusion"    unit="kg/(m^2 s)" > strd_zdfp_e3t * 1026.0 * 0.001  </field>
    <field id="ttrd_qns_li"      long_name="layer integrated heat-trend: non-solar flux + runoff" unit="W/m^2" grid_ref="grid_T_2D"> ttrd_qns_e3t * 1026.0 * 3991.86795711963 </field>
    <field id="ttrd_qsr_li"      long_name="layer integrated heat-trend: solar flux"              unit="W/m^2"  grid_ref="grid_T_3D"> ttrd_qsr_e3t * 1026.0 * 3991.86795711963 </field>
    <field id="ttrd_bbl_li"      long_name="layer integrated heat-trend: bottom boundary layer "  unit="W/m^2" > ttrd_bbl_e3t * 1026.0 * 3991.86795711963 </field>
    <field id="strd_bbl_li"      long_name="layer integrated salt-trend: bottom boundary layer "  unit="kg/(m^2 s)" > strd_bbl_e3t * 1026.0 * 0.001  </field>
    <field id="ttrd_evd_li"      long_name="layer integrated heat-trend: evd convection "         unit="W/m^2" >ttrd_evd_e3t * 1026.0 * 3991.86795711963  </field>
    <field id="strd_evd_li"      long_name="layer integrated salt-trend: evd convection "         unit="kg/(m^2 s)" > strd_evd_e3t * 1026.0 * 0.001  </field>

  </field_group>

  <!--  Total trends calculated every time step-->
  <field_group id="trendT" grid_ref="grid_T_3D">
    <field id="ttrd_tot"      long_name="temperature-trend: total model trend"         unit="degC/s" />
    <field id="strd_tot"      long_name="salinity   -trend: total model trend"         unit="1e-3/s" />
    <!-- Thickness weighted versions: -->
    <field id="ttrd_tot_e3t"      unit="degC/s * m"  >  ttrd_tot * e3t </field>
    <field id="strd_tot_e3t"      unit="1e-3/s * m"  >  strd_tot * e3t </field>
    <!-- OMIP  layer-integrated total trends -->
    <field id="ttrd_tot_li"      long_name="layer integrated heat-trend: total model trend :"         unit="W/m^2" > ttrd_tot_e3t * 1026.0 * 3991.86795711963 </field>
    <field id="strd_tot_li"      long_name="layer integrated salt-trend: total model trend :"         unit="kg/(m^2 s)" > strd_tot_e3t * 1026.0 * 0.001  </field>

    <!-- **** these trends have not been apportioned to all/even/odd ts yet **** -->
    <!-- variables available with ln_KE_trd -->
    <field id="ketrd_hpg"     long_name="ke-trend: hydrostatic pressure gradient"          unit="W/s^3"                        />
    <field id="ketrd_spg"     long_name="ke-trend: surface     pressure gradient"          unit="W/s^3"                        />
    <field id="ketrd_spgexp"  long_name="ke-trend: surface pressure gradient (explicit)"   unit="W/s^3"                        />
    <field id="ketrd_spgflt"  long_name="ke-trend: surface pressure gradient (filter)"     unit="W/s^3"                        />
    <field id="ssh_flt"       long_name="filtered contribution to ssh (dynspg_flt)"        unit="m"       grid_ref="grid_T_2D" />
    <field id="w0"            long_name="surface vertical velocity"                        unit="m/s"     grid_ref="grid_T_2D" />
    <field id="pw0_exp"       long_name="surface pressure flux due to ssh"                 unit="W/s^2"   grid_ref="grid_T_2D" />
    <field id="pw0_flt"       long_name="surface pressure flux due to filtered ssh"        unit="W/s^2"   grid_ref="grid_T_2D" />
    <field id="ketrd_keg"     long_name="ke-trend: KE gradient         or hor. adv."       unit="W/s^3"                        />
    <field id="ketrd_rvo"     long_name="ke-trend: relative  vorticity or metric term"     unit="W/s^3"                        />
    <field id="ketrd_pvo"     long_name="ke-trend: planetary vorticity"                    unit="W/s^3"                        />
    <field id="ketrd_zad"     long_name="ke-trend: vertical  advection"                    unit="W/s^3"                        />
    <field id="ketrd_udx"     long_name="ke-trend: U.dx[U]"                                unit="W/s^3"                        />
    <field id="ketrd_ldf"     long_name="ke-trend: lateral   diffusion"                    unit="W/s^3"                        />
    <field id="ketrd_zdf"     long_name="ke-trend: vertical  diffusion"                    unit="W/s^3"                        />
    <field id="ketrd_tau"     long_name="ke-trend: wind stress "                           unit="W/s^3"   grid_ref="grid_T_2D" />
    <field id="ketrd_bfr"     long_name="ke-trend: bottom friction (explicit)"             unit="W/s^3"                        />
    <field id="ketrd_bfri"    long_name="ke-trend: bottom friction (implicit)"             unit="W/s^3"                        />
    <field id="ketrd_atf"     long_name="ke-trend: asselin time filter trend"              unit="W/s^3"                        />
    <field id="ketrd_convP2K" long_name="ke-trend: conversion (potential to kinetic)"      unit="W/s^3"                        />
    <field id="KE"            long_name="kinetic energy: u(n)*u(n+1)/2"                    unit="W/s^2"                        />

    <!-- variables available when explicit lateral mixing is used (ln_dynldf_OFF=F) -->
    <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" />
    <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" />
    <!-- variables available with ln_traadv_eiv=T and ln_diaeiv=T -->
    <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" />

    <!-- variables available with ln_PE_trd -->
    <field id="petrd_xad"     long_name="pe-trend: i-advection"                unit="W/m^3"                        />
    <field id="petrd_yad"     long_name="pe-trend: j-advection"                unit="W/m^3"                        />
    <field id="petrd_zad"     long_name="pe-trend: k-advection"                unit="W/m^3"                        />
    <field id="petrd_sad"     long_name="pe-trend: surface adv. (linssh true)" unit="W/m^3"   grid_ref="grid_T_2D" />
    <field id="petrd_ldf"     long_name="pe-trend: lateral  diffusion"         unit="W/m^3"                        />
    <field id="petrd_zdf"     long_name="pe-trend: vertical diffusion"         unit="W/m^3"                        />
    <field id="petrd_zdfp"    long_name="pe-trend: pure vert. diffusion"       unit="W/m^3"                        />
    <field id="petrd_dmp"     long_name="pe-trend: interior restoring"         unit="W/m^3"                        />
    <field id="petrd_bbl"     long_name="pe-trend: bottom boundary layer"      unit="W/m^3"                        />
    <field id="petrd_npc"     long_name="pe-trend: non-penetrative conv."      unit="W/m^3"                        />
    <field id="petrd_nsr"     long_name="pe-trend: surface forcing + runoff"   unit="W/m^3"                        />
    <field id="petrd_qsr"     long_name="pe-trend: solar penetr. heating"      unit="W/m^3"                        />
    <field id="petrd_bbc"     long_name="pe-trend: geothermal heating"         unit="W/m^3"                        />
    <field id="petrd_atf"     long_name="pe-trend: asselin time filter"        unit="W/m^3"                        />
    <field id="PEanom"        long_name="potential energy anomaly"             unit="1"                            />
    <field id="alphaPE"       long_name="partial deriv. of PEanom wrt T"       unit="degC-1"                       />
    <field id="betaPE"        long_name="partial deriv. of PEanom wrt S"       unit="1e3"                          />
  </field_group>

  <field_group id="trendU" grid_ref="grid_U_3D">
    <!-- variables available with ln_dyn_trd -->
    <field id="utrd_hpg"       long_name="i-trend: hydrostatic pressure gradient"          unit="m/s^2"                        />
    <field id="utrd_spg"       long_name="i-trend: surface     pressure gradient"          unit="m/s^2"                        />
    <field id="utrd_spgexp"    long_name="i-trend: surface pressure gradient (explicit)"   unit="m/s^2"                        />
    <field id="utrd_spgflt"    long_name="i-trend: surface pressure gradient (filtered)"   unit="m/s^2"                        />
    <field id="utrd_keg"       long_name="i-trend: KE gradient         or hor. adv."       unit="m/s^2"                        />
    <field id="utrd_rvo"       long_name="i-trend: relative  vorticity or metric term"     unit="m/s^2"                        />
    <field id="utrd_pvo"       long_name="i-trend: planetary vorticity"                    unit="m/s^2"                        />
    <field id="utrd_zad"       long_name="i-trend: vertical  advection"                    unit="m/s^2"                        />
    <field id="utrd_udx"       long_name="i-trend: U.dx[U]"                                unit="m/s^2"                        />
    <field id="utrd_ldf"       long_name="i-trend: lateral   diffusion"                    unit="m/s^2"                        />
    <field id="utrd_zdf"       long_name="i-trend: vertical  diffusion"                    unit="m/s^2"                        />
    <field id="utrd_tau"       long_name="i-trend: wind stress "                           unit="m/s^2"   grid_ref="grid_U_2D" />
    <field id="utrd_bfr"       long_name="i-trend: bottom friction (explicit)"             unit="m/s^2"                        />
    <field id="utrd_bfri"      long_name="i-trend: bottom friction (implicit)"             unit="m/s^2"                        />
    <field id="utrd_tot"       long_name="i-trend: total momentum trend before atf"        unit="m/s^2"                        />
    <field id="utrd_atf"       long_name="i-trend: asselin time filter trend"              unit="m/s^2"                        />
  </field_group>

  <field_group id="trendV" grid_ref="grid_V_3D">
    <!-- variables available with ln_dyn_trd -->
    <field id="vtrd_hpg"       long_name="j-trend: hydrostatic pressure gradient"          unit="m/s^2"                        />
    <field id="vtrd_spg"       long_name="j-trend: surface     pressure gradient"          unit="m/s^2"                        />
    <field id="vtrd_spgexp"    long_name="j-trend: surface pressure gradient (explicit)"   unit="m/s^2"                        />
    <field id="vtrd_spgflt"    long_name="j-trend: surface pressure gradient (filtered)"   unit="m/s^2"                        />
    <field id="vtrd_keg"       long_name="j-trend: KE gradient         or hor. adv."       unit="m/s^2"                        />
    <field id="vtrd_rvo"       long_name="j-trend: relative  vorticity or metric term"     unit="m/s^2"                        />
    <field id="vtrd_pvo"       long_name="j-trend: planetary vorticity"                    unit="m/s^2"                        />
    <field id="vtrd_zad"       long_name="j-trend: vertical  advection"                    unit="m/s^2"                        />
    <field id="vtrd_vdy"       long_name="i-trend: V.dx[V]"                                unit="m/s^2"                        />
    <field id="vtrd_ldf"       long_name="j-trend: lateral   diffusion"                    unit="m/s^2"                        />
    <field id="vtrd_zdf"       long_name="j-trend: vertical  diffusion"                    unit="m/s^2"                        />
    <field id="vtrd_tau"       long_name="j-trend: wind stress "                           unit="m/s^2"   grid_ref="grid_V_2D" />
    <field id="vtrd_bfr"       long_name="j-trend: bottom friction (explicit)"             unit="m/s^2"                        />
    <field id="vtrd_bfri"      long_name="j-trend: bottom friction (implicit)"             unit="m/s^2"                        />
    <field id="vtrd_tot"       long_name="j-trend: total momentum trend before atf"        unit="m/s^2"                        />
    <field id="vtrd_atf"       long_name="j-trend: asselin time filter trend"              unit="m/s^2"                        />
  </field_group>

  <!-- shared variables available with TOP interface -->
  <field_group id="top_shared" grid_ref="grid_T_3D">
    <field id="xeps"           long_name="Broadband light attenuation"                     unit="-"                            />
    <field id="Heup"           long_name="Euphotic layer depth"                            unit="m"     grid_ref="grid_T_2D"   />
  </field_group>


  <!--
      ============================================================================================================
      Definitions for iodef_demo.xml
      ============================================================================================================
  -->

  <field_group id="TRD" >
    <field field_ref="ttrd_totad_li"   name="opottempadvect"  />
    <field field_ref="ttrd_iso_li"     name="opottemppmdiff"  />
    <field field_ref="ttrd_zdfp_li"    name="opottempdiff"  />
    <field field_ref="ttrd_evd_li"     name="opottempevd" />
    <field field_ref="strd_evd_li"     name="osaltevd" />
    <field field_ref="ttrd_qns_li"     name="opottempqns"  />
    <field field_ref="ttrd_qsr_li"     name="rsdoabsorb" operation="accumulate" />
    <field field_ref="strd_totad_li"   name="osaltadvect" />
    <field field_ref="strd_iso_li"     name="osaltpmdiff"  />
    <field field_ref="strd_zdfp_li"    name="osaltdiff" />
  </field_group>

  <field_group id="mooring" >

    <!-- EOS80 -->
    <field field_ref="toce_pot"     name="thetao_pot"   long_name="sea_water_potential_temperature"  />
    <field field_ref="soce_pra"     name="so_pra"       long_name="sea_water_practical_salinity"     />
    <!-- TEOS10 -->
    <field field_ref="toce_con"     name="thetao_con"   long_name="sea_water_conservative_temperature" />
    <field field_ref="soce_abs"     name="so_con"       long_name="sea_water_absolute_salinity"      />
    <!-- SEOS -->
    <field field_ref="toce_seos"    name="thetao_seos"   long_name="sea_water_temperature"           />
    <field field_ref="soce_seos"    name="so_seos"       long_name="sea_water_salinity"              />

Guillaume Samson's avatar
Guillaume Samson committed
    <field field_ref="uoce"         name="uo"       long_name="sea_water_x_velocity"                 />
    <field field_ref="voce"         name="vo"       long_name="sea_water_y_velocity"                 />
    <field field_ref="woce"         name="wo"       long_name="sea_water_z_velocity"                 />
    <field field_ref="avt"          name="difvho"   long_name="ocean_vertical_heat_diffusivity"      />
    <field field_ref="avm"          name="difvmo"   long_name="ocean_vertical_momentum_diffusivity"  />

      <!-- EOS80 -->
    <field field_ref="sst_pot"          name="tos_pot"      long_name="sea_surface_potential_temperature"                       />
    <field field_ref="sst2_pot"         name="tossq_pot"    long_name="square_of_sea_surface_potential_temperature"             />
    <field field_ref="sstgrad_pot"      name="tosgrad_pot"  long_name="module_of_sea_surface_potential_temperature_gradient"    />
    <field field_ref="sss_pra"          name="sos_pra"      long_name="sea_surface_absolute_salinity"                          />
    <!-- TEOS10 -->
    <field field_ref="sst_con"          name="tos_con"      long_name="sea_surface_conservative_temperature"                       />
    <field field_ref="sst2_con"         name="tossq_con"    long_name="square_of_sea_surface_conservative_temperature"             /> <field field_ref="sstgrad_con"      name="tosgrad_con" long_name="module_of_sea_surface_conservative_temperature_gradient"    />
    <field field_ref="sss_abs"          name="sos_abs"      long_name="sea_surface_absolute_salinity"                          />
    <!-- SEOS -->
    <field field_ref="sst_seos"          name="tos_seos"      long_name="sea_surface_temperature"                       />    
    <field field_ref="sst2_seos"         name="tossq_seos"    long_name="square_of_sea_surface_temperature"             />
    <field field_ref="sstgrad_seos"      name="tosgrad_seos"  long_name="module_of_sea_surface_temperature_gradient"    />
    <field field_ref="sss_seos"          name="sos_seos"      long_name="sea_surface_salinity"                          />

Guillaume Samson's avatar
Guillaume Samson committed
    <field field_ref="ssh"          name="zos"      long_name="sea_surface_height_above_geoid"                />
    <field field_ref="empmr"        name="wfo"      long_name="water_flux_into_sea_water"                     />
    <field field_ref="qsr"          name="rsntds"   long_name="surface_net_downward_shortwave_flux"           />
    <field field_ref="qt"           name="tohfls"   long_name="surface_net_downward_total_heat_flux"          />
    <field field_ref="taum"                                                                                   />
    <field field_ref="20d"                                                                                    />
    <field field_ref="mldkz5"                                                                                 />
    <field field_ref="mldr10_1"                                                                               />
    <field field_ref="mldr10_3"                                                                               />
    <field field_ref="mldr0_1"                                                                                />
    <field field_ref="mldr0_3"                                                                                />
    <field field_ref="mld_dt02"                                                                               />
    <field field_ref="topthdep"                                                                               />
    <field field_ref="pycndep"                                                                                />
    <field field_ref="tinv"                                                                                   />
    <field field_ref="depti"                                                                                  />
    <field field_ref="BLT"          name="blt"      long_name="barrier_layer_thickness"                       />
    <field field_ref="utau"         name="tauuo"    long_name="surface_downward_x_stress"                     />
    <field field_ref="vtau"         name="tauvo"    long_name="surface_downward_y_stress"                     />
  </field_group>

  <field_group id="groupT" >
  <!-- EOS80 -->
    <field id="toce_pot"           long_name="sea_water_potential_temperature"                grid_ref="grid_T_3D"      />
    <field id="soce_pra"           long_name="sea_water_practical_salinity"                   grid_ref="grid_T_3D"      />
    <field id="sst_pot"            long_name="sea_surface_potential_temperature"              grid_ref="grid_T_2D"      />
    <field id="sst2_pot"           long_name="square_of_sea_surface_potential_temperature"    grid_ref="grid_T_2D"      />
    <field id="sss_pra"            long_name="sea_surface_practical_salinity"                 grid_ref="grid_T_2D"      />
  <!-- TEOS10 -->
    <field id="toce_con"           long_name="sea_water_conservative_temperature"             grid_ref="grid_T_3D"      />
    <field id="soce_abs"           long_name="sea_water_absolute_salinity"                    grid_ref="grid_T_3D"      />
    <field id="sst_con"            long_name="sea_surface_conservative_temperature"           grid_ref="grid_T_2D"   unit="degC"   />                
    <field id="sst2_con"           long_name="square_of_sea_surface_conservative_temperature" grid_ref="grid_T_2D"      /> 
    <field id="sss_abs"            long_name="sea_surface_absolute_salinity"                  grid_ref="grid_T_2D"      />
  <!-- SEOS -->
    <field id="toce_seos"          long_name="sea_water_temperature"                        grid_ref="grid_T_3D"      />
    <field id="soce_seos"          long_name="sea_water_salinity"                           grid_ref="grid_T_3D"      />
    <field id="sst_seos"           long_name="sea_surface_temperature"                      grid_ref="grid_T_2D"      />
    <field id="sst2_seos"          long_name="square_of_sea_surface_temperature"            grid_ref="grid_T_2D"      />
    <field id="sss_seos"           long_name="sea_surface_salinity"                         grid_ref="grid_T_2D"      />
    
    <field field_ref="rhop"         name="sigma0"   long_name="potential density"                             />
Guillaume Samson's avatar
Guillaume Samson committed
    <field field_ref="ssh"          name="zos"      long_name="sea_surface_height_above_geoid"                />
    <field field_ref="empmr"        name="wfo"      long_name="water_flux_into_sea_water"                     />
    <field field_ref="qsr"          name="rsntds"   long_name="surface_net_downward_shortwave_flux"           />
    <field field_ref="qt"           name="tohfls"   long_name="surface_net_downward_total_heat_flux"          />
    <field field_ref="taum"                                                                                   />
    <field field_ref="20d"                                                                                    />
    <field field_ref="mldkz5"                                                                                 />
    <field field_ref="mldr10_1"                                                                               />
    <field field_ref="mldr10_3"                                                                               />
    <field field_ref="mld_dt02"                                                                               />
    <field field_ref="topthdep"                                                                               />
    <field field_ref="pycndep"                                                                                />
    <field field_ref="tinv"                                                                                   />
    <field field_ref="depti"                                                                                  />
    <field field_ref="BLT"          name="blt"      long_name="Barrier Layer Thickness"                       />
  </field_group>

  <field_group id="groupU" >
    <field field_ref="uoce"         name="uo"      long_name="sea_water_x_velocity"      />
    <field field_ref="utau"         name="tauuo"   long_name="surface_downward_x_stress" />
  </field_group>

  <field_group id="groupV" >
    <field field_ref="voce"         name="vo"      long_name="sea_water_y_velocity"      />
    <field field_ref="vtau"         name="tauvo"   long_name="surface_downward_y_stress" />
  </field_group>

  <field_group id="groupW" >
    <field field_ref="woce"         name="wo"       long_name="ocean vertical velocity"  />
  </field_group>

  <!-- TMB diagnostic output -->
  <field_group  id="1h_grid_T_tmb" grid_ref="grid_T_2D" operation="instant">


    <!-- EOS80 -->
    <field id="top_temp_pot"           name="votemper_top_pot"  unit="degC"  />
    <field id="mid_temp_pot"           name="votemper_mid_pot"  unit="degC"  />
    <field id="bot_temp_pot"           name="votemper_bot_pot"  unit="degC"  />
    <field id="top_sal_pra"            name="vosaline_top_pra"  unit="psu"   />
    <field id="mid_sal_pra"            name="vosaline_mid_pra"  unit="psu"   />
    <field id="bot_sal_pra"            name="vosaline_bot_pra"  unit="psu"   />
    <!-- TEOS10 -->
    <field id="top_temp_con"           name="votemper_top_con"  unit="degC"  />
    <field id="mid_temp_con"           name="votemper_mid_con"  unit="degC"  />
    <field id="bot_temp_con"           name="votemper_bot_con"  unit="degC"  />
    <field id="top_sal_abs"            name="vosaline_top_abs"  unit="g/kg"   />
    <field id="mid_sal_abs"            name="vosaline_mid_abs"  unit="g/kg"   />
    <field id="bot_sal_abs"            name="vosaline_bot_abs"  unit="g/kg"   />
    <!-- SEOS -->
    <field id="top_temp_seos"           name="votemper_top_seos"  unit="degC"  />
    <field id="mid_temp_seos"           name="votemper_mid_seos"  unit="degC"  />
    <field id="bot_temp_seos"           name="votemper_bot_seos"  unit="degC"  />
    <field id="top_sal_seos"            name="vosaline_top_seos"  unit="psu"   />
    <field id="mid_sal_seos"            name="vosaline_mid_seos"  unit="psu"   />
    <field id="bot_sal_seos"            name="vosaline_bot_seos"  unit="psu"   />

Guillaume Samson's avatar
Guillaume Samson committed
    <field id="top_temp"           name="votemper_top"  unit="degC"  />
    <field id="mid_temp"           name="votemper_mid"  unit="degC"  />
    <field id="bot_temp"           name="votemper_bot"  unit="degC"  />
    <field id="top_sal"            name="vosaline_top"  unit="psu"   />
    <field id="mid_sal"            name="vosaline_mid"  unit="psu"   />
    <field id="bot_sal"            name="vosaline_bot"  unit="psu"   />
    <field id="sshnmasked"         name="sossheig"      unit="m"     />
  </field_group>

  <field_group  id="1h_grid_U_tmb" grid_ref="grid_U_2D" operation="instant">
    <field id="top_u"           name="vozocrtx_top"  unit="m/s"  />
    <field id="mid_u"           name="vozocrtx_mid"  unit="m/s"  />
    <field id="bot_u"           name="vozocrtx_bot"  unit="m/s"  />
    <field id="baro_u"          name="vobtcrtx"      unit="m/s"  />
  </field_group>

  <field_group  id="1h_grid_V_tmb" grid_ref="grid_V_2D" operation="instant">
    <field id="top_v"           name="vomecrty_top"  unit="m/s"  />
    <field id="mid_v"           name="vomecrty_mid"  unit="m/s"  />
    <field id="bot_v"           name="vomecrty_bot"  unit="m/s"  />
    <field id="baro_v"          name="vobtcrty"      unit="m/s"  />
  </field_group>

  <!-- 25h diagnostic output -->
  <field_group id="25h_grid_T" grid_ref="grid_T_3D_inner" operation="instant">
    <field id="tempis25h"         name="insitu temperature 25h mean"       unit="degC" />

    <!-- EOS80 -->
    <field id="temper25h_pot"     name="potential temperature 25h mean"    unit="degC" />
    <field id="salin25h_pra"      name="practical salinity 25h mean"       unit="psu"  />
    <!-- TEOS10 -->
    <field id="temper25h_con"     name="conservative temperature 25h mean" unit="degC" />
    <field id="salin25h_abs"      name="absolute salinity 25h mean"        unit="g/kg"  />
    <!-- SEOS -->
    <field id="temper25h_seos"    name="temperature 25h mean"              unit="degC" />
    <field id="salin25h_seos"     name="salinity 25h mean"                 unit="psu"  />
    
    <field id="ssh25h"            name="sea surface height 25h mean"       unit="m" grid_ref="grid_T_2D"   />
Guillaume Samson's avatar
Guillaume Samson committed
  </field_group>

  <field_group id="25h_grid_U" grid_ref="grid_U_3D_inner" operation="instant" >
    <field id="vozocrtx25h"         name="i current 25h mean"    unit="m/s"   />
  </field_group>

  <field_group id="25h_grid_V" grid_ref="grid_V_3D_inner" operation="instant">
    <field id="vomecrty25h"         name="j current 25h mean"    unit="m/s"    />
  </field_group>

  <field_group id="25h_grid_W" grid_ref="grid_W_3D_inner" operation="instant">
    <field id="vovecrtz25h"         name="k current 25h mean"                 unit="m/s"      />
    <field id="avt25h"              name="vertical diffusivity25h mean"       unit="m2/s" />
    <field id="avm25h"              name="vertical viscosity 25h mean"        unit="m2/s" />
    <field id="tke25h"              name="turbulent kinetic energy 25h mean" />
    <field id="mxln25h"             name="mixing length 25h mean"             unit="m" />
  </field_group>

  <!--
      ============================================================================================================
  -->
  <!-- output variables for my configuration (example) -->

  <field_group id="myvarOCE" >
    <!-- grid T -->
    <field field_ref="e3t"          name="e3t"      long_name="vertical scale factor"           />
    <field field_ref="sst"          name="tos"      long_name="sea_surface_temperature"         />
    <field field_ref="sss"          name="sos"      long_name="sea_surface_salinity"            />
    <field field_ref="ssh"          name="zos"      long_name="sea_surface_height_above_geoid"  />

    <!-- grid U -->
    <field field_ref="e3u"          name="e3u"     long_name="vertical scale factor"            />
    <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"           />

    <!-- grid V -->
    <field field_ref="e3v"          name="e3v"     long_name="vertical scale factor"            />
    <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"           />
  </field_group>

  <!--============================-->
  <!--  CONSERVATION diagnostics  -->
  <!--============================-->
  <!-- BE CAREFUL: this group (OCE_budget) cannot be called in file_def.xml as such (unless nn_fsbc=1)
                   If doing so, the last output (in time) of the netcdf file 
		   would be corrupted (NaN values). However calling each of these
		   variables directly in the file_def.xml works. It is probably 
		   because there is a mix up of sbc variables with other variables
    -->
  <field_group id="OCE_budget"        grid_ref="grid_T_2D" >
    <field field_ref="sst"                 name="tos"          />
    <field field_ref="sss"                 name="sos"          />
    <field field_ref="ssh"                 name="zos"          />
    <!-- mass flux -->
    <field field_ref="empmr"               name="empmr"        />
    <field field_ref="runoffs"             name="runoffs"      />
    <field field_ref="emp_ice"             name="emp_ice"      />
    <field field_ref="emp_oce"             name="emp_oce"      />
    <field field_ref="iceshelf_cea"        name="iceshelf"     />
    <field field_ref="iceberg_cea"         name="iceberg"      />
    <field field_ref="calving_cea"         name="calving"      />
    <!-- <field field_ref="berg_floating_melt"  name="calving" /> -->
    <field field_ref="precip"              name="precip"       />
    <field field_ref="snowpre"             name="snowpre"      />
    <field field_ref="rain"                name="rain"         />
    <field field_ref="evap_ao_cea"         name="evap_ao"      />
    <field field_ref="subl_ai_cea"         name="subl_ai"      />
    <field field_ref="snow_ai_cea"         name="snow_ai"      />
    <field field_ref="snow_ao_cea"         name="snow_ao"      />
    <!-- heat flux -->
    <field field_ref="qsr"                 name="qsr"          />
    <field field_ref="qns"                 name="qns"          />
    <field field_ref="qt_oce"              name="qt_oce"       />
    <field field_ref="qemp_oce"            name="qemp_oce"     />
    <field field_ref="hflx_rain_cea"       name="hflx_rain"    />
    <field field_ref="hflx_evap_cea"       name="hflx_evap"    />
    <field field_ref="hflx_snow_cea"       name="hflx_snow"    />
    <field field_ref="hflx_snow_ao_cea"    name="hflx_snow_ao" />
    <field field_ref="hflx_snow_ai_cea"    name="hflx_snow_ai" />
    <field field_ref="hflx_rnf_cea"        name="hflx_rnf"     />
    <field field_ref="hflx_icb_cea"        name="hflx_icb"     />
    <field field_ref="hflx_isf_cea"        name="hflx_isf"     />
    <!-- salt flux (includes ssr) -->
    <field field_ref="saltflx"             name="saltflx"      />
    <field field_ref="sflx_rnf_cea"        name="sflx_rnf"     />
    <!-- relaxation and damping -->
    <field field_ref="hflx_ssr_cea"        name="hflx_ssr"     />
    <field field_ref="vflx_ssr_cea"        name="vflx_ssr"     />
    <field field_ref="sflx_ssr_cea"        name="sflx_ssr"     />
    <field field_ref="hflx_dmp_cea"        name="hflx_dmp"     />
    <field field_ref="sflx_dmp_cea"        name="sflx_dmp"     />
    <field field_ref="hflx_fwb_cea"        name="hflx_fwb"     />
    <field field_ref="vflx_fwb_cea"        name="vflx_fwb"     />
  </field_group>

  <field_group id="OCE_globalbudget"  grid_ref="grid_scalar" >
    <field field_ref="voltot"              name="scvoltot"     />
    <field field_ref="saltot"              name="scsaltot"     />
    <field field_ref="temptot"             name="sctemtot"     />
  </field_group>

  
</field_definition>