<!--
<?xml version="1.0"?>
       -->
<!-- $id$ -->
<!-- 
============================================================================================================
=                                  definition of all existing variables                                    =
=                                            DO NOT CHANGE                                                 =
============================================================================================================
    -->
<field_definition level="1" prec="4" operation="average" enabled=".TRUE." default_value="1.e20"> <!-- time step automaticaly defined -->

    <!-- 
============================================================================================================
                                  Physical ocean and ice 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"/>
         <field id="e3t_surf" long_name="T-cell thickness" field_ref="e3t" standard_name="cell_thickness" unit="m" grid_ref="grid_T_surface_extract"/>
         <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"/> 
         <field id="toce" long_name="Sea Water Potential Temperature" standard_name="sea_water_potential_temperature" unit="degC" grid_ref="grid_T_3D" name="thetao"/>
         <field id="toce_e3t" long_name="temperature * e3t" unit="degree_C*m" grid_ref="grid_T_3D"> toce * e3t </field>
         <field id="soce" long_name="Sea Water Salinity" standard_name="sea_water_salinity" unit="0.001" grid_ref="grid_T_3D" name="so"/>
         <field id="soce_e3t" long_name="salinity * e3t" unit="0.001*m" grid_ref="grid_T_3D"> soce * e3t </field>

         <field id="sst" long_name="Sea Surface Temperature" standard_name="sea_surface_temperature" unit="degC" name="tos"/>
         <field id="sst2" long_name="Square of Sea Surface Temperature" standard_name="square_of_sea_surface_temperature" unit="degree_C2"> sst * sst </field>
         <field id="sstmax" long_name="max of sea surface temperature" field_ref="sst" operation="maximum"/>
         <field id="sstmin" long_name="min of sea surface temperature" field_ref="sst" operation="minimum"/>
         <field id="sstgrad" long_name="module of sst gradient" unit="degree_C/m"/>
         <field id="sstgrad2" long_name="square of module of sst gradient" unit="degree_C2/m2"/>
         <field id="sbt" long_name="sea bottom temperature" unit="degree_C"/>
         <field id="tosmint" long_name="vertical integral of temperature times density" standard_name="integral_wrt_depth_of_product_of_density_and_potential_temperature" unit="(kg m2) degree_C"/>
	 
         <field id="sss" long_name="Sea Surface Salinity" standard_name="sea_surface_salinity" unit="0.001" name="sos"/>
         <field id="sss2" long_name="square of sea surface salinity" unit="1e-6"> sss * sss </field>
         <field id="sssmax" long_name="max of sea surface salinity" field_ref="sss" operation="maximum"/>
         <field id="sssmin" long_name="min of sea surface salinity" field_ref="sss" operation="minimum"/>
         <field id="sbs" long_name="sea bottom salinity" unit="0.001"/>
         <field id="somint" long_name="vertical integral of salinity times density" standard_name="integral_wrt_depth_of_product_of_density_and_salinity" unit="(kg m2) x (1e-3)"/> 

         <field id="taubot" long_name="bottom stress module" unit="N/m2"/> 

         <field id="ssh" long_name="Sea Surface Height Above Geoid" standard_name="sea_surface_height_above_geoid" unit="m" name="zos"/>
         <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" name="somixhgt"/>
         <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" name="mlots010"/>
         <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="mldkr125" long_name="Mixed Layer Depth (rhoc=.125)" unit="m" name="mlots125" standard_name="Mixed Layer Depth 0.125"/>
         <field id="mldkr03" long_name="Mixed Layer Depth (rhoc=.03)" unit="m" name="mlots030" standard_name="Mixed Layer Depth 0.03"/>
         <field id="mldt05" long_name="Mixed Layer Depth (t_c=0.5)" unit="m" name="mlott05" standard_name="Mixed Layer Depth 0.5C"/>
         <field id="mldt02" long_name="Mixed Layer Depth (t_c=0.2)" unit="m" name="mlott02" standard_name="Mixed Layer Depth 0.2C"/>
         <field id="mldkara" long_name="Mixed Layer Depth (kara)" unit="m"/>
         <field id="mldzint_1" long_name="Mixed Layer Depth (mldzint_1)" unit="m"/>
         <field id="mldzint_2" long_name="Mixed Layer Depth (mldzint_2)" unit="m"/>
         <field id="mldzint_3" long_name="Mixed Layer Depth (mldzint_3)" unit="m"/>
         <field id="mldzint_4" long_name="Mixed Layer Depth (mldzint_4)" unit="m"/>
         <field id="mldzint_5" long_name="Mixed Layer Depth (mldzint_5)" unit="m"/>
         <field id="mldhtc_1" long_name="Mixed Layer Depth (mldhtc_1)" unit="m"/>
         <field id="mldhtc_2" long_name="Mixed Layer Depth (mldhtc_2)" unit="m"/>
         <field id="mldhtc_3" long_name="Mixed Layer Depth (mldhtc_3)" unit="m"/>
         <field id="mldhtc_4" long_name="Mixed Layer Depth (mldhtc_4)" unit="m"/>
         <field id="mldhtc_5" long_name="Mixed Layer Depth (mldhtc_5)" unit="m"/>
         <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" name="hcbtm"/>
         <field id="saltc" long_name="Salt content vertically integrated" unit="0.001*kg/m2" name="scbtm" standard_name="Integrated Salinity column"/>
         <field id="cdn10" long_name="Neutral drag coefficient at 10m" unit=""/>
         <field id="cfl" long_name="CFL diagnostics" unit="" grid_ref="grid_T_3D"/>
         <field id="div" long_name="Divergence" unit="1/x" grid_ref="grid_T_3D"/>
         <field id="rot" long_name="Rotation" unit="1/x" grid_ref="grid_T_3D"/>
         <field id="uocee" long_name="ocean current along east direction" standard_name="eastward_sea_water_velocity" unit="m/s" grid_ref="grid_T_3D" name="uoe"/>
         <field id="vocen" long_name="ocean current along north direction" standard_name="northward_sea_water_velocity" unit="m/s" grid_ref="grid_T_3D" name="von"/>
         <field id="uocees" long_name="ocean surface current along east direction" standard_name="surface_eastward_sea_water_velocity" unit="m/s" name="uose"/>
         <field id="vocens" long_name="ocean surface current along north direction" standard_name="surface_northward_sea_water_velocity" unit="m/s" name="uosn"/>
         <field id="utaue" long_name="Wind Stress along east direction" standard_name="surface_downward_eastward_stress" unit="N/m2" name="taueo"/>
         <field id="vtaun" long_name="Wind Stress along north direction" standard_name="surface_downward_northward_stress" unit="N/m2" name="tauno"/>
         <field id="hst" long_name="Steric height" unit="m" name="sostheig" standard_name="Steric Height"/>
         <field id="hbp" long_name="Bottom Pressure Equivalent Height" unit="m" name="sobpheig" standard_name="Bottom Pressure Equivalent Height"/>
         <field id="sigmat" long_name="Sigmat" unit="N/A" grid_ref="grid_T_3D" name="sigmat" standard_name="Sigmat"/>
         <field id="vottrdmp" long_name="3D-damping T" unit="C/s" grid_ref="grid_T_3D" name="vottrdmp" standard_name="3D-damping T"/>
         <field id="vostrdmp" long_name="3D-damping S" unit="psu/s" grid_ref="grid_T_3D" name="vostrdmp" standard_name="3D-damping S"/>
         <field id="votbiasd" long_name="3D-bias T direct" unit="C" grid_ref="grid_T_3D" name="votbiasd" standard_name="3D-bias T direct"/>
         <field id="vosbiasd" long_name="3D-bias S direct" unit="psu" grid_ref="grid_T_3D" name="vosbiasd" standard_name="3D-bias S direct"/>
         <field id="votbiasp" long_name="3D-bias T pressure" unit="C" grid_ref="grid_T_3D" name="votbiasp" standard_name="3D-bias T pressure"/>
         <field id="vosbiasp" long_name="3D-bias S pressure" unit="psu" grid_ref="grid_T_3D" name="vosbiasp" standard_name="3D-bias S pressure"/>
         <field id="vorbiasp" long_name="3D-bias rho pressure" unit="" grid_ref="grid_T_3D" name="vorbiasp" standard_name="3D-bias rho pressure"/>
         <field id="sodsstmo" long_name="SST mod - obs monthly" unit="K"/>
         <field id="sodsstyy" long_name="SST mod - obs interannual" unit="K"/>

         <!-- EOS -->
         <field id="alpha" long_name="thermal expansion" unit="degree_C-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"/>

         <!-- 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"/>
         <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 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" name="sothedep"/>
         <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="14d" long_name="Depth of 14C isotherm" standard_name="depth_of_isosurface_of_sea_water_potential_temperature" unit="m" name="t14d"/>
         <field id="17d" long_name="Depth of 17C isotherm" standard_name="depth_of_isosurface_of_sea_water_potential_temperature" unit="m" name="t17d"/>
         <field id="20d" long_name="Depth of 17C isotherm" standard_name="depth_of_isosurface_of_sea_water_potential_temperature" unit="m" name="t20d"/>
         <field id="26d" long_name="Depth of 26C isotherm" standard_name="depth_of_isosurface_of_sea_water_potential_temperature" unit="m" name="t26d"/>
         <field id="28d" long_name="Depth of 28C isotherm" standard_name="depth_of_isosurface_of_sea_water_potential_temperature" unit="m" name="t28d"/>
         <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" name="hc300m"/>
         <field id="sal300" long_name="Salt content 0-300m" standard_name="integral_of_sea_water_practical_salinity_wrt_depth_expressed_as_heat_content" unit="1e-3*kg/m2" name="sc300m"/>
         <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" name="hc700m"/>
         <field id="sal700" long_name="Salt content 0-700m" standard_name="integral_of_sea_water_practical_salinity_wrt_depth_expressed_as_heat_content" unit="1e-3*kg/m2" name="sc700m"/>
         <field id="hcbtm" long_name="Heat content coloum" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2"/>
         <field id="salbtm" long_name="Salt content coloum" standard_name="integral_of_sea_water_practical_salinity_wrt_depth_expressed_as_heat_content" unit="1e-3*kg/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" name="hc2000m"/>
         <field id="sal2000" long_name="Salt content 0-2000m" standard_name="integral_of_sea_water_practical_salinity_wrt_depth_expressed_as_heat_content" unit="1e-3*kg/m2" name="sc2000m"/>
         <field id="hc26c" long_name="Heat content 0-26C" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" name="hct26"/>
         <field id="hc26cd" long_name="Heat content 0-26C (relative to 26C)" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2"/>
         <field id="thetaodmp" long_name="Sea Water Potential Temperature Tendency due to Newtonian Relaxation" standard_name="sea_water_potential_temperature_tendency_due_to_newtonian_relaxation" unit="K/s" grid_ref="grid_T_3D" />
         <field id="sodmp" long_name="Sea Water Salinity Tendency due to Newtonian Relaxation" standard_name="sea_water_salinity_tendency_due_to_newtonian_relaxation" unit="PSU/s" grid_ref="grid_T_3D" />

         <!-- variables available with key_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 with key_vvl -->
         <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"/>
      </field_group>

      <!-- Tides -->

      <field_group id="Tides_T" grid_ref="grid_T_2D" operation="once">
         <!-- tidal composante -->
         <field id="M2x" long_name="M2 Elevation harmonic real part " unit="m"/>
         <field id="M2y" long_name="M2 Elevation harmonic imaginary part" unit="m"/>
         <field id="S2x" long_name="S2 Elevation harmonic real part " unit="m"/>
         <field id="S2y" long_name="S2 Elevation harmonic imaginary part" unit="m"/>
         <field id="N2x" long_name="N2 Elevation harmonic real part " unit="m"/>
         <field id="N2y" long_name="N2 Elevation harmonic imaginary part" unit="m"/>
         <field id="K1x" long_name="K1 Elevation harmonic real part " unit="m"/>
         <field id="K1y" long_name="K1 Elevation harmonic imaginary part" unit="m"/>
         <field id="O1x" long_name="O1 Elevation harmonic real part " unit="m"/>
         <field id="O1y" long_name="O1 Elevation harmonic imaginary part" unit="m"/>
         <field id="Q1x" long_name="Q1 Elevation harmonic real part " unit="m"/>
         <field id="Q1y" long_name="Q1 Elevation harmonic imaginary part" unit="m"/>
         <field id="M4x" long_name="M4 Elevation harmonic real part " unit="m"/>
         <field id="M4y" long_name="M4 Elevation harmonic imaginary part" unit="m"/>
         <field id="K2x" long_name="K2 Elevation harmonic real part " unit="m"/>
         <field id="K2y" long_name="K2 Elevation harmonic imaginary part" unit="m"/>
         <field id="P1x" long_name="P1 Elevation harmonic real part " unit="m"/>
         <field id="P1y" long_name="P1 Elevation harmonic imaginary part" unit="m"/>
         <field id="Mfx" long_name="Mf Elevation harmonic real part " unit="m"/>
         <field id="Mfy" long_name="Mf Elevation harmonic imaginary part" unit="m"/>
         <field id="Mmx" long_name="Mm Elevation harmonic real part " unit="m"/>
         <field id="Mmy" long_name="Mm Elevation harmonic imaginary part" unit="m"/>
      </field_group>
	 
      <field_group id="Tides_U" grid_ref="grid_U_2D" operation="once">
         <field id="M2x_u" long_name="M2 current barotrope along i-axis harmonic real part " unit="m/s"/>
         <field id="M2y_u" long_name="M2 current barotrope along i-axis harmonic imaginary part " unit="m/s"/>
         <field id="S2x_u" long_name="S2 current barotrope along i-axis harmonic real part " unit="m/s"/>
         <field id="S2y_u" long_name="S2 current barotrope along i-axis harmonic imaginary part " unit="m/s"/>
         <field id="N2x_u" long_name="N2 current barotrope along i-axis harmonic real part " unit="m/s"/>
         <field id="N2y_u" long_name="N2 current barotrope along i-axis harmonic imaginary part " unit="m/s"/>
         <field id="K1x_u" long_name="K1 current barotrope along i-axis harmonic real part " unit="m/s"/>
         <field id="K1y_u" long_name="K1 current barotrope along i-axis harmonic imaginary part " unit="m/s"/>
         <field id="O1x_u" long_name="O1 current barotrope along i-axis harmonic real part " unit="m/s"/>
         <field id="O1y_u" long_name="O1 current barotrope along i-axis harmonic imaginary part " unit="m/s"/>
         <field id="Q1x_u" long_name="Q1 current barotrope along i-axis harmonic real part " unit="m/s"/>
         <field id="Q1y_u" long_name="Q1 current barotrope along i-axis harmonic imaginary part " unit="m/s"/>
         <field id="M4x_u" long_name="M4 current barotrope along i-axis harmonic real part " unit="m/s"/>
         <field id="M4y_u" long_name="M4 current barotrope along i-axis harmonic imaginary part " unit="m/s"/>
         <field id="K2x_u" long_name="K2 current barotrope along i-axis harmonic real part " unit="m/s"/>
         <field id="K2y_u" long_name="K2 current barotrope along i-axis harmonic imaginary part " unit="m/s"/>
         <field id="P1x_u" long_name="P1 current barotrope along i-axis harmonic real part " unit="m/s"/>
         <field id="P1y_u" long_name="P1 current barotrope along i-axis harmonic imaginary part " unit="m/s"/>
         <field id="Mfx_u" long_name="Mf current barotrope along i-axis harmonic real part " unit="m/s"/>
         <field id="Mfy_u" long_name="Mf current barotrope along i-axis harmonic imaginary part " unit="m/s"/>
         <field id="Mmx_u" long_name="Mm current barotrope along i-axis harmonic real part " unit="m/s"/>
         <field id="Mmy_u" long_name="Mm current barotrope along i-axis harmonic imaginary part " unit="m/s"/>
      </field_group>
	 
      <field_group id="Tides_V" grid_ref="grid_V_2D" operation="once">
         <field id="M2x_v" long_name="M2 current barotrope along j-axis harmonic real part " unit="m/s"/>
         <field id="M2y_v" long_name="M2 current barotrope along j-axis harmonic imaginary part " unit="m/s"/>
         <field id="S2x_v" long_name="S2 current barotrope along j-axis harmonic real part " unit="m/s"/>
         <field id="S2y_v" long_name="S2 current barotrope along j-axis harmonic imaginary part " unit="m/s"/>
         <field id="N2x_v" long_name="N2 current barotrope along j-axis harmonic real part " unit="m/s"/>
         <field id="N2y_v" long_name="N2 current barotrope along j-axis harmonic imaginary part " unit="m/s"/>
         <field id="K1x_v" long_name="K1 current barotrope along j-axis harmonic real part " unit="m/s"/>
         <field id="K1y_v" long_name="K1 current barotrope along j-axis harmonic imaginary part " unit="m/s"/>
         <field id="O1x_v" long_name="O1 current barotrope along j-axis harmonic real part " unit="m/s"/>
         <field id="O1y_v" long_name="O1 current barotrope along j-axis harmonic imaginary part " unit="m/s"/>
         <field id="Q1x_v" long_name="Q1 current barotrope along j-axis harmonic real part " unit="m/s"/>
         <field id="Q1y_v" long_name="Q1 current barotrope along j-axis harmonic imaginary part " unit="m/s"/>
         <field id="M4x_v" long_name="M4 current barotrope along j-axis harmonic real part " unit="m/s"/>
         <field id="M4y_v" long_name="M4 current barotrope along j-axis harmonic imaginary part " unit="m/s"/>
         <field id="K2x_v" long_name="K2 current barotrope along j-axis harmonic real part " unit="m/s"/>
         <field id="K2y_v" long_name="K2 current barotrope along j-axis harmonic imaginary part " unit="m/s"/>
         <field id="P1x_v" long_name="P1 current barotrope along j-axis harmonic real part " unit="m/s"/>
         <field id="P1y_v" long_name="P1 current barotrope along j-axis harmonic imaginary part " unit="m/s"/>
         <field id="Mfx_v" long_name="Mf current barotrope along j-axis harmonic real part " unit="m/s"/>
         <field id="Mfy_v" long_name="Mf current barotrope along j-axis harmonic imaginary part " unit="m/s"/>
         <field id="Mmx_v" long_name="Mm current barotrope along j-axis harmonic real part " unit="m/s"/>
         <field id="Mmy_v" long_name="Mm current barotrope along j-axis harmonic imaginary part " unit="m/s"/>	 
      </field_group>

      <!-- SBC -->
      
      <field_group id="SBC" grid_ref="grid_T_2D"> <!-- time step automaticaly defined based on nn_fsbc -->
         <field id="empmr" long_name="Net Upward Water Flux" standard_name="water_flux_out_of_sea_ice_and_sea_water" unit="kg/m2/s" name="sowaflup"/>
         <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="Virtual Salt Flux into Sea Water" unit="0.001/m2/s" name="vsf" standard_name="virtual_salt_flux_into_sea_water"/>
         <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" name="snowpre"/>
         <field id="runoffs" long_name="Water Flux into Sea Water From Rivers" 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="qt" long_name="Downward Heat Flux at Sea Water Surface" standard_name="surface_downward_heat_flux_in_sea_water" unit="W m-2" name="hfds"/>
         <field id="qns" long_name="non solar Downward Heat Flux" unit="W/m2"/>
         <field id="qsr" long_name="Net Shortwave Surface Radiation" standard_name="surface_net_downward_shortwave_flux" unit="W m-2" name="rss"/>
         <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"/>
         <field id="qrp" long_name="Heat Flux Correction" standard_name="heat_flux_correction" unit="W m-2" name="hfcorr"/>
         <field id="erp" long_name="Virtual Salt Flux Correction" standard_name="virtual_salt_flux_correction" unit="kg m-2 s-1" name="vsfcorr"/>
         <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"/>

         <!-- diags -->
         <field id="hfxdhc" long_name="Heat content variation in snow and ice" unit="W/m2"/>
         <field id="hfxtur" long_name="turbulent heat flux at the ice base" unit="W/m2"/>
          <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="1e-3/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="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" />
          <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="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"        />

          <!-- * variable related to ice shelf forcing * -->
          <field id="fwfisf"       long_name="Ice shelf melting"                             unit="kg/m2/s"  />
          <field id="fwfisf3d"     long_name="Ice shelf melting"                             unit="kg/m2/s"  grid_ref="grid_T_3D" />
          <field id="qlatisf"      long_name="Ice shelf latent heat flux"                    unit="W/m2"     />
          <field id="qlatisf3d"    long_name="Ice shelf latent heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" />
          <field id="qhcisf"       long_name="Ice shelf heat content flux"                   unit="W/m2"     />
          <field id="qhcisf3d"     long_name="Ice shelf heat content flux"                   unit="W/m2"     grid_ref="grid_T_3D" />
          <field id="isfgammat"    long_name="transfert coefficient for isf (temperature) "  unit="m/s"      />
          <field id="isfgammas"    long_name="transfert coefficient for isf (salinity)    "  unit="m/s"      />
          <field id="stbl"         long_name="salinity in the Losh tbl                    "  unit="PSU"      />
          <field id="ttbl"         long_name="temperature in the Losh tbl                 "  unit="C"        />
          <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="thermald"     long_name="thermal driving of ice shelf melting        "  unit="C"        />
          <field id="tfrz"         long_name="top freezing point (used to compute melt)   "  unit="C"        />
          <field id="tinsitu"      long_name="top insitu temperature (used to cmpt melt)  "  unit="C"        />
          <field id="ustar"        long_name="ustar at T point used in ice shelf melting  "  unit="m/s"      />

          <!-- *_oce variables available with ln_blk_clio or ln_blk_core -->
          <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="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"  />

          <!-- 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="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_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="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"  />

          <!-- 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="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="-"    />

         <!-- * variable related to ice shelf forcing * -->
         <field id="fwfisf" long_name="Ice shelf melting" unit="kg/m2/s"/>
         <field id="fwfisf3d" long_name="Ice shelf melting" unit="kg/m2/s" grid_ref="grid_T_3D"/>
         <field id="qlatisf" long_name="Ice shelf latent heat flux" unit="W/m2"/>
         <field id="qlatisf3d" long_name="Ice shelf latent heat flux" unit="W/m2" grid_ref="grid_T_3D"/>
         <field id="qhcisf" long_name="Ice shelf heat content flux" unit="W/m2"/>
         <field id="qhcisf3d" long_name="Ice shelf heat content flux" unit="W/m2" grid_ref="grid_T_3D"/>
         <field id="isfgammat" long_name="transfert coefficient for isf (temperature)" unit="m/s"/>
         <field id="isfgammas" long_name="transfert coefficient for isf (salinity)" unit="m/s"/>
         <field id="stbl" long_name="salinity in the Losh tbl" unit="0.001"/>
         <field id="ttbl" long_name="temperature in the Losh tbl" unit="degree_C"/>

         <!-- *_oce variables available with ln_blk_clio or ln_blk_core -->
         <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="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"/>

         <!-- 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_prec_cea" long_name="heat flux due to precipitation" standard_name="temperature_flux_due_to_precipitation_expressed_as_heat_flux" 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_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="hflx_cal_cea" long_name="heat flux due to calving" standard_name="heat_flux_into_sea_water_due_to_iceberg_thermodynamics" 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"/>

	 <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 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"  />

         <!-- ice fields -->
          <field id="icemask"      long_name="Fraction of time steps with sea ice"                     standard_name="sea_ice_mask"                              unit="" />
          <field id="icemask05"    long_name="Ice mask (0 if ice conc. lower than  5%, 1 otherwise)"   standard_name="sea_ice_mask05"                            unit="" />
          <field id="icemask15"    long_name="Ice mask (0 if ice conc. lower than 15%, 1 otherwise)"   standard_name="sea_ice_mask15"                            unit="" />
	  <field id="icepres"      long_name="Fraction of time steps with sea ice"                     standard_name="sea_ice_time_fraction"                     unit="" />
          <field id="fasticepres"  long_name="Fraction of time steps with landfast ice"                standard_name="fast_ice_time_fraction"                    unit="" />
	  
	  <!-- general fields -->
          <field id="icemass"      long_name="Sea-ice mass per area"                                   standard_name="sea_ice_amount"                            unit="kg/m2"/>
          <field id="snwmass"      long_name="Snow mass per area"                                      standard_name="liquid_water_content_of_surface_snow"      unit="kg/m2"/>
          <field id="iceconc"      long_name="Sea-ice area fraction"                                   standard_name="sea_ice_area_fraction"                     unit=""    />
          <field id="icevolu"      long_name="Sea-ice volume per area"                                 standard_name="sea_ice_thickness"                         unit="m"   />
          <field id="icethic"      long_name="Sea-ice thickness per area"                              standard_name="sea_ice_thickness"                         unit="m"   />
          <field id="snwthic"      long_name="Snow thickness per area"                                 standard_name="snow_thickness"                            unit="m"   />
          <field id="icebrv"       long_name="brine volume"                                                                                                      unit="%"   />
          <field id="iceage"       long_name="ice age"                                                                                                           unit="days"/>
          <field id="icehnew"      long_name="frazil ice collection thickness"                                                                                   unit="m"   />
          <field id="snwvolu"      long_name="snow volume"                                                                                                       unit="m"   />
          <field id="icefrb"       long_name="Sea-ice freeboard"                                       standard_name="sea_ice_freeboard"                         unit="m"   />
          <field id="icealb"       long_name="Sea-ice or snow albedo"                                  standard_name="sea_ice_albedo"                            unit=""    />
	  <field id="icealbnm"       long_name="Sea-ice or snow albedo no mask"                                  standard_name="sea_ice_albedo_nomask"                            unit=""    />
	  <field id="icealbmsk"       long_name="Sea-ice or snow albedomsk"                                  standard_name="sea_ice_albedo_mask"                            unit=""    />
          <field id="tau_icebfr"   long_name="ice friction on ocean bottom for landfast ice"                                                                     unit="N/2" />
          <field id="iicevelur"    long_name="Geographical-east component of sea ice velocity"         standard_name="eastward_sea_ice_velocity"                 unit="m/s"/>
          <field id="iicevelvr"    long_name="Geographical-north component of sea ice velocity"        standard_name="northward_sea_ice_velocity"                unit="m/s"/>
          <field id="iocestrur" long_name="Wind stress u geographical" unit="Pa"/>
     	  <field id="iocestrvr" long_name="Wind stress v geographical" unit="Pa"/>
	  
	  <!-- melt ponds -->
	  <field id="iceapnd"      long_name="melt pond fraction"                                      standard_name="sea_ice_meltpond_fraction"                 unit="%" /> 
          <field id="icehpnd"      long_name="melt pond depth"                                         standard_name="sea_ice_meltpond_depth"                    unit="m" /> 
          <field id="icevpnd"      long_name="melt pond volume"                                        standard_name="sea_ice_meltpond_volume"                   unit="m" /> 
          <field id="icehlid"      long_name="melt pond lid depth"                                     standard_name="sea_ice_meltpondlid_depth"                 unit="m" /> 
          <field id="icevlid"      long_name="melt pond lid volume"                                    standard_name="sea_ice_meltpondlid_volume"                unit="m" /> 	  
	  <!-- heat -->
          <field id="icetemp"      long_name="Mean ice temperature"                                                                                              unit="degC" />
          <field id="snwtemp"      long_name="Mean snow temperature"                                                                                             unit="degC" />
          <field id="icettop"      long_name="temperature at the ice surface"                                                                                    unit="degC" />
          <field id="icetbot"      long_name="temperature at the ice bottom"                                                                                     unit="degC" />
          <field id="icetsni"      long_name="temperature at the snow-ice interface"                                                                             unit="degC" />
          <field id="icehc"        long_name="ice heat content"                                                                                                  unit="J/m2" /> 
          <field id="snwhc"        long_name="snow heat content"                                                                                                 unit="J/m2" />
	  
	  <!-- salt -->
          <field id="icesalt"      long_name="Sea ice salinity"                                                                                                  unit="g/kg"  />
          <field id="icesalm"      long_name="Mass of salt in sea ice per area"                        standard_name="sea_ice_salt_mass"                         unit="kg/m2" />
	  
	  <!-- momentum (rheology) -->
          <field id="uice"         long_name="X-component of sea ice velocity"                         standard_name="sea_ice_x_velocity"                        unit="m/s"  /> 
          <field id="vice"         long_name="Y-component of sea ice velocity"                         standard_name="sea_ice_y_velocity"                        unit="m/s"  />      
          <field id="icevel"       long_name="Sea-ice speed"                                           standard_name="sea_ice_speed"                             unit="m/s"  />
          <field id="utau_ai"      long_name="X-component of atmospheric stress on sea ice"            standard_name="surface_downward_x_stress"                 unit="N/m2" />
          <field id="vtau_ai"      long_name="Y-component of atmospheric stress on sea ice"            standard_name="surface_downward_y_stress"                 unit="N/m2" />
          <field id="utau_oi"      long_name="X-component of ocean stress on sea ice"                  standard_name="sea_ice_base_upward_x_stress"              unit="N/m2" />
          <field id="vtau_oi"      long_name="Y-component of ocean stress on sea ice"                  standard_name="sea_ice_base_upward_y_stress"              unit="N/m2" />
          <field id="utau_bi"      long_name="X-component of ocean bottom stress on sea ice -landfast" standard_name="ocean_bottom_upward_x_stress"              unit="N/m2" />
          <field id="vtau_bi"      long_name="Y-component of ocean bottom stress on sea ice -landfast" standard_name="ocean_bottom_upward_y_stress"              unit="N/m2" />
          <field id="isig1"        long_name="1st principal stress component for EVP rhg"                                                                        unit=""     />
          <field id="isig2"        long_name="2nd principal stress component for EVP rhg"                                                                        unit=""     />
          <field id="isig3"        long_name="convergence measure for EVP rheology (must be around 1)"                                                           unit=""     />
          <field id="normstr"      long_name="Average normal stress in sea ice"                        standard_name="average_normal_stress"                     unit="N/m"  />
          <field id="sheastr"      long_name="Maximum shear stress in sea ice"                         standard_name="maximum_shear_stress"                      unit="N/m"  />
          <field id="icestr"       long_name="Compressive sea ice strength"                            standard_name="compressive_strength_of_sea_ice"           unit="N/m"  />
          <field id="icediv"       long_name="Divergence of the sea-ice velocity field"                standard_name="divergence_of_sea_ice_velocity"            unit="s-1"  />
          <field id="iceshe"       long_name="Maximum shear of sea-ice velocity field"                 standard_name="maximum_shear_of_sea_ice_velocity"         unit="s-1"  />
          <field id="beta_evp"     long_name="Relaxation parameter of ice rheology (beta)"             standard_name="relaxation_parameter_of_ice_rheology"      unit=""  />

	  <!-- surface heat fluxes -->
          <field id="qt_ice"       long_name="total heat flux at ice surface"                          standard_name="surface_downward_heat_flux_in_air"         unit="W/m2" />
          <field id="qsr_ice"      long_name="solar heat flux at ice surface"                          standard_name="surface_downwelling_shortwave_flux_in_air" unit="W/m2" />
          <field id="qns_ice"      long_name="non-solar heat flux at ice surface (including E-P)"                                                                unit="W/m2" />
          <field id="qtr_ice_bot"  long_name="solar heat flux transmitted through the ice (to the ocean)"                                                        unit="W/m2" />
          <field id="qtr_ice_top"  long_name="solar heat flux transmitted through the ice surface"                                                               unit="W/m2" />
          <field id="qt_oce_ai"    long_name="total heat flux at the ocean   surface: interface oce-(ice+atm)"                                                   unit="W/m2" />
          <field id="qt_atm_oi"    long_name="total heat flux at the oce-ice surface: interface atm-(ice+oce)"                                                   unit="W/m2" />
          <field id="qemp_ice"     long_name="Downward Heat Flux from E-P over ice"                                                                              unit="W/m2" />
          <field id="albedo"       long_name="Mean albedo over sea ice and ocean"                                                                                unit=""     />
	  
	  <!-- trends -->
          <field id="afxthd"       long_name="sea-ice area fraction change from thermodynamics"  standard_name="tendency_of_sea_ice_area_fraction_due_to_dynamics" unit="s-1" />
          <field id="afxdyn"       long_name="sea-ice area fraction change from dynamics"        standard_name="tendency_of_sea_ice_area_fraction_due_to_dynamics" unit="s-1" />
	  <field id="afxtot"       long_name="area tendency (total)"                                                                                               unit="s-1" />
	  
	  <!-- momentum (advection) -->
          <field id="icemtrp"      long_name="ice  mass transport"                                          unit="kg/m2/s" />
          <field id="snwmtrp"      long_name="snw  mass transport"                                          unit="kg/m2/s" />
          <field id="salmtrp"      long_name="salt mass transport"                                          unit="kg/m2/s" />
          <field id="dihctrp"      long_name="ice heat content transport"                                   unit="W/m2"    />
          <field id="dshctrp"      long_name="snw heat content transport"                                   unit="W/m2"    />
	  
	  <!-- salt fluxes -->
          <field id="sfxice"       long_name="ice-ocean salt flux from ice growth/melt (neg = growth)"      unit="kg/m2/s" />
          <field id="sfxbog"       long_name="ice-ocean salt flux from ice bottom growth"                   unit="kg/m2/s" />
          <field id="sfxbom"       long_name="ice-ocean salt flux from ice bottom melt"                     unit="kg/m2/s" />
          <field id="sfxsum"       long_name="ice-ocean salt flux from ice surface melt"                    unit="kg/m2/s" />
          <field id="sfxlam"       long_name="ice-ocean salt flux from ice lateral melt"                    unit="kg/m2/s" />
          <field id="sfxsni"       long_name="ice-ocean salt flux from snow-ice formation"                  unit="kg/m2/s" />
          <field id="sfxopw"       long_name="ice-ocean salt flux from ice formation in open water"         unit="kg/m2/s" />
          <field id="sfxdyn"       long_name="ice-ocean salt flux from ice dynamics (ridging)"              unit="kg/m2/s" />
          <field id="sfxbri"       long_name="ice-ocean salt flux from brines"                              unit="kg/m2/s" />
          <field id="sfxres"       long_name="ice-ocean salt flux from undiagnosed processes"               unit="kg/m2/s" />
          <field id="sfxsub"       long_name="ice-ocean salt flux from ice sublimation"                     unit="kg/m2/s" />

	  <!-- mass fluxes -->
          <field id="vfxice"       long_name="ice-ocean mass flux from ice melt/growth (neg = growth)"      unit="kg/m2/s" />
          <field id="vfxbog"       long_name="ice-ocean mass flux from ice bottom growth"                   unit="kg/m2/s" />
          <field id="vfxbom"       long_name="ice-ocean mass flux from ice bottom melt"                     unit="kg/m2/s" />
          <field id="vfxsum"       long_name="ice-ocean mass flux from ice surface melt"                    unit="kg/m2/s" />
          <field id="vfxlam"       long_name="ice-ocean mass flux from ice lateral melt"                    unit="kg/m2/s" />
          <field id="vfxsni"       long_name="ice-ocean mass flux from snow-ice formation"                  unit="kg/m2/s" />
          <field id="vfxopw"       long_name="ice-ocean mass flux from ice growth in open water"            unit="kg/m2/s" />
          <field id="vfxdyn"       long_name="ice-ocean mass flux from ice dynamics (ridging)"              unit="kg/m2/s" />
          <field id="vfxres"       long_name="ice-ocean mass flux from undiagnosed processes"               unit="kg/m2/s" />
          <field id="vfxpnd"       long_name="ice-ocean mass flux from ice melt ponds"                      unit="kg/m2/s" />
          <field id="vfxsub"       long_name="ice-atm.  mass flux from ice sublimation"                     unit="kg/m2/s" />
          <field id="vfxsub_err"   long_name="ice-ocean mass flux from unbalance in sublimation (incl. in qt_oce_ai)" unit="kg/m2/s" />
          <field id="vfxthin"      long_name="ice-ocean mass flux from ice growth in open water + thin ice(20cm)"     unit="kg/m2/s" />

          <field id="vfxsnw"       long_name="ice-ocean mass flux from snw melt/growth (neg = growth)"      unit="kg/m2/s" />
          <field id="vfxsnw_sum"   long_name="ice-ocean mass flux from snw surface melt"                    unit="kg/m2/s" />
          <field id="vfxsnw_sni"   long_name="ice-ocean mass flux from snow-ice formation"                  unit="kg/m2/s" />
          <field id="vfxsnw_dyn"   long_name="ice-ocean mass flux from dynamics (ridging)"                  unit="kg/m2/s" />
          <field id="vfxsnw_sub"   long_name="ice-atm.  mass flux from snw sublimation"                     unit="kg/m2/s" />
          <field id="vfxsnw_pre"   long_name="snw precipitation on ice"                                     unit="kg/m2/s" />

	  <!-- heat fluxes -->
          <field id="hfxbog"       long_name="heat flux used for bottom ice growth (neg.)"                  unit="W/m2" />
          <field id="hfxbom"       long_name="heat flux used for bottom ice melt (pos.)"                    unit="W/m2" />
          <field id="hfxsum"       long_name="heat flux used for surface ice melt"                          unit="W/m2" />
          <field id="hfxopw"       long_name="heat flux used for open water ice formation"                  unit="W/m2" />
          <field id="hfxdif"       long_name="heat flux used for ice temperature change"                    unit="W/m2" />
          <field id="hfxsnw"       long_name="heat flux used for snow melt"                                 unit="W/m2" />
          <field id="hfxerr"       long_name="heat flux error after heat diffusion"                         unit="W/m2" />
          <!-- heat fluxes associated with mass exchange -->
          <field id="hfxthd"       long_name="heat fluxes from ice-ocean mass exchange during thermo"       unit="W/m2" />
          <field id="hfxdyn"       long_name="heat fluxes from ice-ocean mass exchange during dynamic"      unit="W/m2" />
          <field id="hfxres"       long_name="heat fluxes from ice-ocean mass exchange during resultant"    unit="W/m2" />
          <field id="hfxsub"       long_name="heat fluxes from ice-atm. mass exchange during sublimation"   unit="W/m2" />
          <field id="hfxspr"       long_name="heat fluxes from ice-atm. mass exchange during snow precip"   unit="W/m2" />

          <field id="hfxsensib"    long_name="Net sensible heat flux under sea ice (neg = ice cooling)"         standard_name="ice_ocean_heat_flux"                        unit="W/m2" />
          <field id="hfxcndtop"    long_name="Net conductive heat flux at the ice surface (neg = ice cooling)"  standard_name="conductive_heat_flux_at_sea_ice_surface"    unit="W/m2" /> 
          <field id="hfxcndbot"    long_name="Net conductive heat flux at the ice bottom (neg = ice cooling)"   standard_name="conductive_heat_flux_at_sea_ice_bottom"     unit="W/m2" />

          <!-- diags -->
          <field id="hfxdhc"       long_name="Heat content variation in snow and ice (neg = ice cooling)"   unit="W/m2" />

	  <!-- diagnostics of the negative values resulting from the advection scheme -->
	  <field id="iceneg_pres"  long_name="Fraction of time steps with negative sea ice volume"                    unit=""     />
	  <field id="iceneg_volu"  long_name="Negative sea ice volume per area arising from advection"                unit="m"    />
          <field id="iceneg_hfx"   long_name="Negative sea ice heat content (eq. heat flux) arising from advection"   unit="W/m2" />

	  <!-- categories -->
	  <field id="iceconc_cat"  long_name="Sea-ice concentration per category"                unit=""       grid_ref="grid_T_3D_ncatice" />
          <field id="icethic_cat"  long_name="Sea-ice thickness per category"                    unit="m"      grid_ref="grid_T_3D_ncatice" />
          <field id="snwthic_cat"  long_name="Snow thickness per category"                       unit="m"      grid_ref="grid_T_3D_ncatice" />
          <field id="icesalt_cat"  long_name="Sea-Ice Bulk salinity per category"                unit="g/kg"   grid_ref="grid_T_3D_ncatice" />
          <field id="icetemp_cat"  long_name="Ice temperature per category"                      unit="degC"   grid_ref="grid_T_3D_ncatice" />
          <field id="snwtemp_cat"  long_name="Snow temperature per category"                     unit="degC"   grid_ref="grid_T_3D_ncatice" />
          <field id="icettop_cat"  long_name="Ice/snow surface temperature per category"         unit="degC"   grid_ref="grid_T_3D_ncatice" />
          <field id="iceapnd_cat"  long_name="Ice melt pond concentration per category"          unit="%"      grid_ref="grid_T_3D_ncatice" /> 
          <field id="icehpnd_cat"  long_name="Ice melt pond thickness per category"              unit="m"      grid_ref="grid_T_3D_ncatice" /> 
          <field id="icehlid_cat"  long_name="Ice melt pond lid thickness per category"          unit="m"       detect_missing_value="true" /> 
          <field id="iceafpnd_cat" long_name="Ice melt pond fraction per category"               unit="m"      grid_ref="grid_T_3D_ncatice" /> 
          <field id="iceaepnd_cat" long_name="Ice melt pond effective fraction per category"     unit=""        /> 
          <field id="icemask_cat"  long_name="Fraction of time step with sea ice (per category)" unit=""       grid_ref="grid_T_3D_ncatice" />
          <field id="iceage_cat"   long_name="Ice age per category"                              unit="days"   grid_ref="grid_T_3D_ncatice" />
          <field id="icebrv_cat"   long_name="Brine volume per category"                         unit="%"      grid_ref="grid_T_3D_ncatice" />
          <field id="icealb_cat"   long_name="Sea-ice or snow albedo"                            unit=""        detect_missing_value="true" />

	  <!-- ================= -->
          <!-- Add-ons for SIMIP -->
          <!-- ================= -->
          <field id="icemass_cmip" long_name="Sea-ice mass per area"                                   standard_name="sea_ice_amount"                          detect_missing_value="true" unit="kg/m2"  > icemass * icemask   + $missval * (1.-icemask  ) </field>
          <field id="icethic_cmip" long_name="Sea-ice thickness"                                       standard_name="sea_ice_thickness"                       detect_missing_value="true" unit="m"      > icethic * icemask05 + $missval * (1.-icemask05) </field>
          <field id="snwmass_cmip" long_name="Snow mass per area"                                      standard_name="liquid_water_content_of_surface_snow"    detect_missing_value="true" unit="kg/m2"  > snwmass * icemask   + $missval * (1.-icemask  ) </field>
          <field id="snwthic_cmip" long_name="Snow thickness"                                          standard_name="surface_snow_thickness"                  detect_missing_value="true" unit="m"      > snwthic * icemask05 + $missval * (1.-icemask05) </field>
          <field id="iceconc_pct"  long_name="Sea-ice area fraction in percent"                        standard_name="sea_ice_area_fraction_in_percent"                                    unit="%"      > iceconc * 100.                                  </field>
          <field id="iceage_cmip"  long_name="Age of sea ice"                                          standard_name="age_of_sea_ice"                          detect_missing_value="true" unit="days"   > iceage  * icemask15 + $missval * (1.-icemask15) </field>
          <field id="icesalt_cmip" long_name="Sea ice salinity"                                        standard_name="sea_ice_salinity"                        detect_missing_value="true" unit="g/kg"   > icesalt * icemask   + $missval * (1.-icemask  ) </field>
          <field id="icefrb_cmip"  long_name="Sea-ice freeboard"                                       standard_name="sea_ice_freeboard"                       detect_missing_value="true" unit="m"      > icefrb  * icemask   + $missval * (1.-icemask  ) </field>

	  <!-- heat -->
          <field id="icettop_cmip" long_name="Surface temperature of sea ice"                          standard_name="sea_ice_surface_temperature"             detect_missing_value="true" unit="K"  > (icettop+273.15) * icemask + $missval * (1.-icemask) </field>
          <field id="icetsni_cmip" long_name="Temperature at snow-ice interface"                       standard_name="sea_ice_snow_interface_temperature"      detect_missing_value="true" unit="K"  > (icetsni+273.15) * icemask + $missval * (1.-icemask) </field>
          <field id="icetbot_cmip" long_name="Temperature at ice-ocean interface"                      standard_name="sea_ice_bottom_temperature"              detect_missing_value="true" unit="K"  > (icetbot+273.15) * icemask + $missval * (1.-icemask) </field>
          <field id="icehc_cmip"   long_name="Sea-ice heat content per unit area"                      standard_name="integral_of_sea_ice_temperature_wrt_depth_expressed_as_heat_content" detect_missing_value="true" unit="J/m2" > icehc * icemask + $missval * (1.-icemask) </field>
          <field id="snwhc_cmip"   long_name="Snow-heat content per unit area"                         standard_name="thermal_energy_content_of_surface_snow"                              detect_missing_value="true" unit="J/m2" > snwhc * icemask + $missval * (1.-icemask) </field>

	  <!-- fluxes (mass, heat, salt) -->
          <field id="vfxsum_cmip"    long_name="Freshwater flux from sea-ice surface"                  standard_name="freshwater_flux_from_ice_surface"        detect_missing_value="true" unit="kg/m2/s"  > vfxsum    * icemask + $missval * (1.-icemask) </field>
          <field id="vfxice_cmip"    long_name="Freshwater flux from sea ice"                          standard_name="freshwater_flux_from_ice"                detect_missing_value="true" unit="kg/m2/s"  > vfxice    * icemask + $missval * (1.-icemask) </field>
          <field id="hfxsensib_cmip" long_name="Net sensible heat flux under sea ice"                  standard_name="ice_ocean_heat_flux"                     detect_missing_value="true" unit="W/m2"     > hfxsensib * icemask + $missval * (1.-icemask) </field> 
          <field id="hfxcndtop_cmip" long_name="Net conductive heat flux at the ice surface"           standard_name="conductive_heat_flux_at_sea_ice_surface" detect_missing_value="true" unit="W/m2"     > hfxcndtop * icemask + $missval * (1.-icemask) </field> 
          <field id="hfxcndbot_cmip" long_name="Net conductive heat flux at the ice bottom"            standard_name="conductive_heat_flux_at_sea_ice_bottom"  detect_missing_value="true" unit="W/m2"     > hfxcndbot * icemask + $missval * (1.-icemask) </field>
          <field id="sfxice_cmip"    long_name="Salt flux from sea ice"                                standard_name="salt_flux_from_ice"                      detect_missing_value="true" unit="kg/m2/s"  > sfxice    * icemask + $missval * (1.-icemask) </field>

	  <!-- trends -->
          <field id="dmithd"       long_name="sea-ice mass change from thermodynamics"                 standard_name="tendency_of_sea_ice_amount_due_to_thermodynamics"                        unit="kg/m2/s" />
          <field id="dmidyn"       long_name="sea-ice mass change from dynamics"                       standard_name="tendency_of_sea_ice_amount_due_to_dynamics"                              unit="kg/m2/s" />
          <field id="dmiopw"       long_name="sea-ice mass change through growth in supercooled open water (frazil)" standard_name="tendency_of_sea_ice_amount_due_to_freezing_in_open_water"  unit="kg/m2/s" />
          <field id="dmibog"       long_name="sea-ice mass change through basal growth"                standard_name="tendency_of_sea_ice_amount_due_to_congelation_ice_accumulation"          unit="kg/m2/s" />
          <field id="dmisni"       long_name="sea-ice mass change through snow-to-ice conversion"      standard_name="tendency_of_sea_ice_amount_due_to_snow_conversion"                       unit="kg/m2/s" />
          <field id="dmtsub"       long_name="snow and sea ice mass change through sublimation"        standard_name="tendency_of_snow_and_ice_amount_due_to_sublimation"                      unit="kg/m2/s" />
          <field id="dmssub"       long_name="sea-ice mass change through evaporation and sublimation" standard_name="water_evaporation_flux"                                                  unit="kg/m2/s" />
          <field id="dmisub"       long_name="snow mass change through evaporation or sublimation"     standard_name="surface_snow_sublimation_flux"                                           unit="kg/m2/s" />
          <field id="dmisum"       long_name="sea-ice mass change through surface melting"             standard_name="tendency_of_sea_ice_amount_due_to_surface_melting"                       unit="kg/m2/s" />
          <field id="dmibom"       long_name="sea-ice mass change through bottom melting"              standard_name="tendency_of_sea_ice_amount_due_to_basal_melting"                         unit="kg/m2/s" />
          <field id="dmilam"       long_name="sea-ice mass change through lateral melting"             standard_name="tendency_of_sea_ice_amount_due_to_lateral_melting"                       unit="kg/m2/s" />
          <field id="dmsspr"       long_name="snow mass change through snow fall"                      standard_name="snowfall_flux"                                                           unit="kg/m2/s" />
          <field id="dmsmel"       long_name="snow mass change through melt"                           standard_name="surface_snow_melt_flux"                                                  unit="kg/m2/s" />
          <field id="dmsdyn"       long_name="snow mass change through advection by sea-ice dynamics"  standard_name="tendency_of_snow_mass_due_to_sea_ice_dynamics"                           unit="kg/m2/s" />
          <field id="dmsssi"       long_name="snow mass change through snow-to-ice conversion"         standard_name="tendency_of_snow_mass_due_to_snow_to_ice_conversion"                     unit="kg/m2/s" />

	  <!-- momentum (rheology) -->
          <field id="uice_cmip"    long_name="X-component of sea ice velocity"                         standard_name="sea_ice_x_velocity"                    detect_missing_value="true" unit="m/s"   > uice    * icemask + $missval * (1.-icemask) </field>
          <field id="vice_cmip"    long_name="Y-component of sea ice velocity"                         standard_name="sea_ice_y_velocity"                    detect_missing_value="true" unit="m/s"   > vice    * icemask + $missval * (1.-icemask) </field>     
          <field id="icevel_cmip"  long_name="Sea-ice speed"                                           standard_name="sea_ice_speed"                         detect_missing_value="true" unit="m/s"   > icevel  * icemask + $missval * (1.-icemask) </field> 
          <field id="utau_ai_cmip" long_name="X-component of atmospheric stress on sea ice"            standard_name="surface_downward_x_stress"             detect_missing_value="true" unit="N/m2"  > utau_ai * icemask + $missval * (1.-icemask) </field>
          <field id="vtau_ai_cmip" long_name="Y-component of atmospheric stress on sea ice"            standard_name="surface_downward_y_stress"             detect_missing_value="true" unit="N/m2"  > vtau_ai * icemask + $missval * (1.-icemask) </field>
          <field id="dssh_dx"      long_name="Sea-surface tilt term in force balance (x-component)"    standard_name="sea_surface_tilt_force_on_sea_ice_x"                               unit="N/m2" />
          <field id="dssh_dy"      long_name="Sea-surface tilt term in force balance (y-component)"    standard_name="sea_surface_tilt_force_on_sea_ice_y"                               unit="N/m2" />
          <field id="corstrx"      long_name="Coriolis force term in force balance (x-component)"      standard_name="coriolis_force_on_sea_ice_x"                                       unit="N/m2" />
          <field id="corstry"      long_name="Coriolis force term in force balance (y-component)"      standard_name="coriolis_force_on_sea_ice_y"                                       unit="N/m2" />
          <field id="intstrx"      long_name="Internal stress term in force balance (x-component)"     standard_name="internal_stress_in_sea_ice_x"                                      unit="N/m2" />
          <field id="intstry"      long_name="Internal stress term in force balance (y-component)"     standard_name="internal_stress_in_sea_ice_y"                                      unit="N/m2" />

	  <!-- momentum (advection) -->
	  <field id="xmtrpice"     long_name="X-component of ice mass transport"                       standard_name="ice_x_transport"                       unit="kg/s" />
          <field id="ymtrpice"     long_name="Y-component of ice mass transport"                       standard_name="ice_y_transport"                       unit="kg/s" />
          <field id="xmtrpsnw"     long_name="X-component of snw mass transport"                       standard_name="snw_x_transport"                       unit="kg/s" />
          <field id="ymtrpsnw"     long_name="Y-component of snw mass transport"                       standard_name="snw_y_transport"                       unit="kg/s" />
          <field id="xatrp"        long_name="X-component of ice area transport"                       standard_name="area_x_transport"                      unit="m2/s" />
          <field id="yatrp"        long_name="Y-component of ice area transport"                       standard_name="area_y_transport"                      unit="m2/s" />
          <field id="xmtrptot"     long_name="X-component of sea-ice mass transport"                   standard_name="sea_ice_x_transport"                   unit="kg/s" > xmtrpice + xmtrpsnw </field>
          <field id="ymtrptot"     long_name="Y-component of sea-ice mass transport"                   standard_name="sea_ice_y_transport"                   unit="kg/s" > ymtrpice + ymtrpsnw </field>

	  <!-- categories -->
          <field id="iceconc_cat_cmip"     long_name="Sea-ice area fractions in thickness categories"  standard_name="sea_ice_area_fraction_over_categories" detect_missing_value="true" unit=""  grid_ref="grid_T_3D_ncatice" > iceconc_cat      * icemask_cat + $missval * (1.-icemask_cat) </field>
          <field id="icethic_cat_cmip"     long_name="Sea-ice thickness in thickness categories"       standard_name="sea_ice_thickness_over_categories"     detect_missing_value="true" unit="m" grid_ref="grid_T_3D_ncatice" > icethic_cat      * icemask_cat + $missval * (1.-icemask_cat) </field>
          <field id="snwthic_cat_cmip"     long_name="Snow thickness in thickness categories"          standard_name="snow_thickness_over_categories"        detect_missing_value="true" unit="m" grid_ref="grid_T_3D_ncatice" > snwthic_cat      * icemask_cat + $missval * (1.-icemask_cat) </field>
          <field id="iceconc_cat_pct_cmip" long_name="Sea-ice area fractions in thickness categories"  standard_name="sea_ice_area_fraction_over_categories" detect_missing_value="true" unit="%" grid_ref="grid_T_3D_ncatice" > iceconc_cat*100. * icemask_cat + $missval * (1.-icemask_cat) </field>

	  <!-- categories -->
          <field id="wndi" long_name="Wind in i" standard_name="x_wind" unit="m/s"/>
          <field id="wndj" long_name="Wind in j" standard_name="y_wind" unit="m/s"/>
          <field id="tair" long_name="Air temperature" standard_name="air_temperature" unit="K"/>
          <field id="humi" long_name="Specific humidity" standard_name="specific_humidity" unit="g/g"/>
          <field id="qsri" long_name="Net Shortwave Surface Radiation" standard_name="surface_net_downward_shortwave_flux" unit="W/m2"/>
          <field id="qlwi" long_name="Net Longwave Surface Radiation"  standard_name="surface_net_downward_longwave_flux" unit="W/m2"/>
          <field id="prec" long_name="Total precipitation" standard_name="precipitation_flux" unit="kg/m2/s"/>
          <field id="snow" long_name="Snow precipitation" standard_name="snowfall_flux" unit="kg/m2/s"/>
          <field id="slp"  long_name="Mean sea level pressure" standard_name="air_pressure_at_mean_sea_level" unit="Pa"/>
          <field id="wnde" long_name="Zonal wind" standard_name="eashward_wind" unit="m/s"/>
          <field id="wndn" long_name="Meridional wind" standard_name="morthward_wind" unit="m/s"/>
	  

	</field_group> <!-- SBC_2D -->

	<!-- scalar variables -->
	<field_group id="SBC_0D"  grid_ref="grid_1point" >
          <!-- available with ln_icediaout -->
          <field id="ibgfrcvoltop"    long_name="global mean ice/snow forcing at interface ice/snow-atm (volume equivalent ocean volume)"   unit="km3"      />
          <field id="ibgfrcvolbot"    long_name="global mean ice/snow forcing at interface ice/snow-ocean (volume equivalent ocean volume)" unit="km3"      />
          <field id="ibgfrctemtop"    long_name="global mean heat on top of ice/snw/ocean-atm "                                             unit="1e20J"    />
          <field id="ibgfrctembot"    long_name="global mean heat below ice (on top of ocean) "                                             unit="1e20J"    />
          <field id="ibgfrcsal"       long_name="global mean ice/snow forcing (salt equivalent ocean volume)"                               unit="pss*km3"  />
          <field id="ibgfrchfxtop"    long_name="global mean heat flux on top of ice/snw/ocean-atm "                                        unit="W/m2"     />
          <field id="ibgfrchfxbot"    long_name="global mean heat flux below ice (on top of ocean) "                                        unit="W/m2"     />
	  
          <field id="ibgvolume"       long_name="drift in ice/snow volume (equivalent ocean volume)"                                        unit="km3"      />
          <field id="ibgsaltco"       long_name="drift in ice salt content (equivalent ocean volume)"                                       unit="pss*km3"  />
          <field id="ibgheatco"       long_name="drift in ice/snow heat content"                                                            unit="1e20J"    />
          <field id="ibgheatfx"       long_name="drift in ice/snow heat flux"                                                               unit="W/m2"     />

          <field id="ibgvol_tot"      long_name="global mean ice volume"                                                                    unit="km3"      />
          <field id="sbgvol_tot"      long_name="global mean snow volume"                                                                   unit="km3"      />
          <field id="ibgarea_tot"     long_name="global mean ice area"                                                                      unit="km2"      />
          <field id="ibgsalt_tot"     long_name="global mean ice salt content"                                                              unit="1e-3*km3" />
          <field id="ibgheat_tot"     long_name="global mean ice heat content"                                                              unit="1e20J"    />
          <field id="sbgheat_tot"     long_name="global mean snow heat content"                                                             unit="1e20J"    />

          <field id="NH_iceextt"      long_name="Sea ice extent North"                   standard_name="sea_ice_extent_n"                   unit="1e6_km2"  />
          <field id="SH_iceextt"      long_name="Sea ice extent South"                   standard_name="sea_ice_extent_s"                   unit="1e6_km2"  />
          <field id="NH_icevolu"      long_name="Sea ice volume North"                   standard_name="sea_ice_volume_n"                   unit="1e3_km3"  />
          <field id="SH_icevolu"      long_name="Sea ice volume South"                   standard_name="sea_ice_volume_s"                   unit="1e3_km3"  />
          <field id="NH_icearea"      long_name="Sea ice area North"                     standard_name="sea_ice_area_n"                     unit="1e6_km2"  />
          <field id="SH_icearea"      long_name="Sea ice area South"                     standard_name="sea_ice_area_s"                     unit="1e6_km2"  />

          <field id="strait_mifl"     long_name="Sea ice mass flux through straits"      standard_name="sea_ice_mass_transport_across_line" unit="kg/s"  grid_ref="grid_4strait" />
          <field id="strait_arfl"     long_name="Sea ice area flux through straits"      standard_name="sea_ice_area_transport_across_line" unit="m2/s"  grid_ref="grid_4strait" />  
          <field id="strait_msfl"     long_name="Sea ice snow flux through straits"      standard_name="snow_mass_transport_across_line"    unit="kg/s"  grid_ref="grid_4strait" />


	 <!-- Old style lim2 output -->
	 <field id="isnowthi" long_name="Snow thickness" unit="m" name="sisnthick" standard_name="surface_snow_thickness"/>
	 <field id="iicethic" long_name="Sea Ice Thickness" unit="m" name="sithick" standard_name="sea_ice_thickness"/>
	 <field id="iiceprod" long_name="sea-ice mass change from thermodynamics" unit="m/kt" name="sidmassth" standard_name="tendency_of_sea_ice_amount_due_to_thermodynamics"/>
	 <field id="ileadfra" long_name="Sea Ice Area Fraction (Ocean Grid)" unit="Fraction" name="siconc" standard_name="sea_ice_area_fraction"/>
         <field id="iicetemp" long_name="Surface temperature of sea ice" unit="C" name="sitemptop" standard_name="sea_ice_surface_temperature"/>
	 <field id="ioceflxb" long_name="" unit="w/m2" name="siflsensupbot" standard_name="upward_sea_ice_basal_heat_flux"/>
	 <field id="iicevelu" long_name="Ice velocity along i-axis at I-point (ice presence average)" unit="m/s" name="icevelu" standard_name="sea_ice_x_velocity"/>
	 <field id="iicevelv" long_name="Ice velocity along j-axis at I-point (ice presence average)" unit="m/s" name="icevelv" standard_name="sea_ice_y_velocity"/>
	 <field id="isstempe" long_name="Sea surface temperature" unit="C"/>
	 <field id="isssalin" long_name="Sea surface salinity" unit="PSU"/>
	 <field id="iocetflx" long_name="Total flux at ocean surface" unit="w/m2" name="iocetflx" standard_name="Total flux at ocean surface"/> 
	 <field id="iocesflx" long_name="Solar flux at ocean surface" unit="w/m2" name="iocesflx" standard_name="Solar flux at ocean surface"/>
	 <field id="iocwnsfl" long_name="Non-solar flux at ocean surface" unit="w/m2" name="iocwnsfl" standard_name="Non-solar flux at ocean surface"/>
	 <field id="iocesafl" long_name="Salt flux from sea ice" unit="kg/m2/kt" name="sfdsi" standard_name="downward_sea_ice_basal_salt_flux"/>
	 <field id="iocestru" long_name="Wind stress u" unit="Pa"/>
	 <field id="iocestrv" long_name="Wind stress v" unit="Pa"/>
	 <field id="iicesflx" long_name="Solar flux at ice/ocean surface" unit="w/m2" name="iicesflx" standard_name="Solar flux at ice/ocean surface"/>
	 <field id="iicenflx" long_name="Non-solar flux at ice/ocean surface" unit="w/m2" name="iicenflx" standard_name="Non-solar flux at ice/ocean surface"/>
	 <field id="isnowpre" long_name="Snow precipitation" unit="kg/day" name="snowpre" standard_name="snowfall_flux"/>
	 <field id="iifr1_i0" long_name="1st Qsr fraction penetrating ice" unit="-"/>
	 <field id="iifr2_i0" long_name="2nd Qsr fraction penetrating ice" unit="-"/>
	 <field id="isrftemp" long_name="Surface temperature (snow/ice/sst)" unit="C" name="ositemp" standard_name="Surface temperature (snow/ice/sst)"/>
	 <field id="iicemthi" long_name="Sea-ice volume per area" unit="m" name="sivol" standard_name="sea_ice_thickness"/>
	 <field id="iitempl1" long_name="Temp inside the ice/snow layer 1" unit="C" name="iitempl1" standard_name="Temp inside the ice/snow layer 1"/>
	 <field id="iitempl2" long_name="Temp inside the ice/snow layer 2" unit="C" name="iitempl2" standard_name="Temp inside the ice/snow layer 2"/>
	 <field id="iitempl3" long_name="Temp inside the ice/snow layer 3" unit="C" name="iitempl3" standard_name="Temp inside the ice/snow layer 3"/>

	 <!-- wave input fields -->
	 <field id="soswh" long_name="Significamt wave height" unit="m"/>
	 <field id="somwp" long_name="Mean wave period" unit="s"/>
	 <field id="sophioc" long_name="Normalized energy flux" unit="Frac."/>
	 <field id="sotauoc" long_name="Normalized stress" unit="Frac."/>
	 <field id="sowndmwam" long_name="10m neutral wind from WAM" unit="m/s"/>
	 <field id="sophifwam" long_name="Full phi from WAM" unit="kg/s3"/>
	 <field id="soicestrn" long_name="Ice mean square strain rate " unit="-"/>

	 <!-- sglexe fields -->
	 <field id="soatmsst" long_name="SST from the atmosphere" unit="degC"/>
	 <field id="soatmtsk" long_name="Skin temperature from the atmosphere" unit="degC"/>
	 <field id="socplsst" long_name="SST passed to the atmosphere" unit="degC"/>

	 <!-- IFS coupling fields -->
	 <field id="sfqns_tot" long_name="sfqns_tot" unit="-"/>
	 <field id="sfqns_ice" long_name="sfqns_ice" unit="-"/>
	 <field id="sfqsr_tot" long_name="sfqsr_tot" unit="-"/>
	 <field id="sfqsr_ice" long_name="sfqsr_ice" unit="-"/>
	 <field id="sfemp_tot" long_name="sfemp_tot" unit="-"/>
	 <field id="sfemp_ice" long_name="sfemp_ice" unit="-"/>
	 <field id="sfdqdns_ice" long_name="sfdqdns_ice" unit="-"/>
	 <field id="sfsprecip" long_name="sfsprecip" unit="-"/>
         <field id="sfqns" long_name="sfqns" unit="-"/>
         <field id="sfqsr" long_name="sfqsr" unit="-"/>
         <field id="sfemp" long_name="sfemp" unit="-"/>

	 <!-- WAM coupling fields -->
         <field id="sfswh" long_name="sfswh" unit="-"/>
         <field id="sfmwp" long_name="sfmwp" unit="-"/>
         <field id="sfphioc" long_name="sfphioc" unit="-"/>
         <field id="sftauoc" long_name="sftauoc" unit="-"/>
         <field id="sfwndmwam" long_name="sfwndmwam" unit="-"/>
         <field id="sfphifwam" long_name="sfphifwam" unit="-"/>
         <field id="sfstrn" long_name="sfstrn" unit="-"/>

      </field_group>

      <!-- U grid -->
      
      <field_group id="grid_U" grid_ref="grid_U_2D">
         <field id="e2u" long_name="U-cell width in meridional direction" standard_name="cell_width" unit="m" grid_ref="grid_U_2D"/>
         <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="Surface Downward X Stress" standard_name="surface_downward_x_stress" unit="N m-2" name="tauuo"/>
         <field id="uoce" long_name="Sea Water X Velocity" standard_name="sea_water_x_velocity" unit="m s-1" grid_ref="grid_U_3D" name="uo"/>
         <field id="uoce_e3u" long_name="ocean current along i-axis * e3u" unit="m2/s" grid_ref="grid_U_3D"> uoce * e3u </field>
         <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"/>
         <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="degree_C*m/s" grid_ref="grid_U_3D"/>
         <field id="uoces" long_name="ocean transport along i-axis times salinity (CRS)" unit="0.001*m/s" grid_ref="grid_U_3D"/>
         <field id="vostcoru" long_name="Stokes-Coriolis tendencies along i-axis" unit="m/s2" grid_ref="grid_U_3D"/>
	 <field id="soutauwam" long_name="Wind Stress along i-axis from WAM" unit="N/m2"/>

	 <!-- IFS coupling fields -->
         <field id="sfutau" long_name="sfutau" unit="-"/>
         <field id="sfutau_ice" long_name="sfutau_ice" unit="-"/>

	 <!-- WAM coupling fields -->
         <field id="sfutauwam" long_name="sfutauwam" unit="-"/>
         <field id="sfustokes" long_name="sfustokes" unit="-"/>

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

         <!-- uoce_eiv: available with key_traldf_eiv and key_diaeiv -->
         <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_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_eiv: available with key_trabbl -->
         <field id="uoce_bbl" long_name="BBL ocean current along i-axis" unit="m/s"/>
         <field id="ahu_bbl" long_name="BBL diffusive flux along i-axis" unit="m3/s"/>

        <!-- 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"/>

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

         <!-- variables available with key_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"/>
         <field id="u_heattr" long_name="ocean eulerian heat transport along i-axis" standard_name="ocean_heat_x_transport" unit="W"/>
         <field id="uadv_heattr" long_name="ocean advective heat transport along i-axis" standard_name="advectice_ocean_heat_x_transport" unit="W"/>
         <field id="u_salttr" long_name="ocean eulerian salt transport along i-axis" standard_name="ocean_salt_x_transport" unit="0.001*kg/s"/>
         <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="W"/>
         <field id="udiff_heattr" long_name="ocean diffusion heat transport along i-axis" standard_name="ocean_heat_x_transport_due_to_diffusion" unit="W"/>

        <!-- 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" />

        <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"                      />
        <field id="u_heattr"      long_name="ocean eulerian heat transport along i-axis"                      standard_name="ocean_heat_x_transport"                         unit="W"                         />
        <field id="u_salttr"      long_name="ocean eulerian salt transport along i-axis"                      standard_name="ocean_salt_x_transport"                         unit="1e-3*kg/s"                 />
        <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="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" />
        <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 meridional 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="vtau" long_name="Surface Downward Y Stress" standard_name="surface_downward_y_stress" unit="N m-2" name="tauvo"/>
         <field id="voce" long_name="Sea Water Y Velocity" standard_name="sea_water_y_velocity" unit="m s-1" grid_ref="grid_V_3D" name="vo"/>
         <field id="voce_e3v" long_name="ocean current along j-axis * e3v" unit="m2/s" grid_ref="grid_V_3D"> voce * e3v </field>
         <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"/>
         <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="degree_C*m/s" grid_ref="grid_V_3D"/>
         <field id="voces" long_name="ocean transport along j-axis times salinity (CRS)" unit="0.001*m/s" grid_ref="grid_V_3D"/>
         <field id="vostcorv" long_name="Stokes-Coriolis tendencies along j-axis" unit="m/s2" grid_ref="grid_V_3D"/>
	 <field id="sovtauwam" long_name="Wind Stress along j-axis from WAM" unit="N/m2"/>

	 <!-- IFS coupling fields -->
         <field id="sfvtau" long_name="sfvtau" unit="-"/>
         <field id="sfvtau_ice" long_name="sfvtau_ice" unit="-"/>

	 <!-- WAM coupling fields -->
         <field id="sfvtauwam" long_name="sfvtauwam" unit="-"/>
         <field id="sfvstokes" long_name="sfvstokes" 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 -->
         <field id="psiv_mle" long_name="MLE streamfunction along j-axis" unit="m3/s" grid_ref="grid_V_3D"/>

         <!-- voce_eiv: available with key_traldf_eiv and key_diaeiv -->
         <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_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_U_3D"/>

         <!-- voce_eiv: available with key_trabbl -->
         <field id="voce_bbl" long_name="BBL ocean current along j-axis" unit="m/s"/>
         <field id="ahv_bbl" long_name="BBL diffusive flux along j-axis" unit="m3/s"/>

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

         <!-- variables available with key_diaar5 -->
         <field id="v_masstr" long_name="Ocean Mass Y Transport" standard_name="ocean_mass_y_transport" unit="kg/s" grid_ref="grid_V_3D"/>
         <field id="v_heattr" long_name="Ocean Heat X Transport" standard_name="ocean_heat_y_transport" unit="W"/>
         <field id="vadv_heattr" long_name="ocean advective heat transport along j-axis" standard_name="advectice_ocean_heat_y_transport" unit="W"/>
         <field id="v_salttr" long_name="ocean eulerian salt transport along i-axis" standard_name="ocean_salt_y_transport" unit="0.001*kg/s"/>
         <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_y_transport_due_to_bolus_advection" unit="W"/>
         <field id="vdiff_heattr" long_name="ocean diffusion heat transport along j-axis" standard_name="ocean_heat_y_transport_due_to_diffusion" unit="W"/>

        <!-- 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"                         />
        <field id="v_salttr"      long_name="ocean eulerian salt transport along i-axis"    standard_name="ocean_salt_y_transport"                          unit="1e-3*kg/s"                 />
        <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="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" />
        <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="Sea Water Z Velocity" standard_name="upward_sea_water_velocity" unit="m s-1" name="wo"/>
        <field id="wocetr_eff" long_name="effective ocean vertical transport" unit="m3/s"/>

        <!-- woce_eiv: available with key_traldf_eiv and key_diaeiv -->
        <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 along k-axis" standard_name="ocean_heat_z_transport_due_to_bolus_advection" unit="W" grid_ref="grid_W_3D"/>
        <field id="weiv_salttr3d" long_name="ocean bolus salt transport along k-axis" standard_name="ocean_salt_z_transport_due_to_bolus_advection" unit="kg" grid_ref="grid_W_3D"/>

        <field id="avt" long_name="vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s"/>
        <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"/>

        <!-- avs: available with key_zdfddm -->
        <field id="avs" long_name="salt vertical eddy diffusivity" standard_name="ocean_vertical_salt_diffusivity" unit="m2/s"/>
        <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 to vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_convection" unit="m2/s"/>
        <field id="avm_evd" long_name="convective enhancement to vertical viscosity" standard_name="ocean_vertical_momentum_diffusivity_due_to_convection" unit="m2/s"/>

        <!-- avt_tide: available with key_zdftmx -->
        <field id="av_tide" long_name="tidal vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_tides" unit="m2/s"/>

       <!-- variables available with key_zdftmx_new -->
        <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="wave-induced vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_internal_waves" unit="m2/s"/>
        <field id="av_tmx" long_name="tidal-induced vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_tides" unit="m2/s"/>
        <field id="bflx_tmx" long_name="wave-induced buoyancy flux" standard_name="buoyancy_flux_due_to_internal_waves" unit="W/kg"/>
        <field id="pcmap_tmx" 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_tmx" long_name="power density available for mixing" standard_name="power_available_for_mixing_from_breaking_internal_waves" unit="W/kg"/>
        <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 key_diaar5 -->   
        <field id="w_masstr" long_name="Upward Ocean 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"/>

        <!-- aht2d and  aht2d_eiv: available with key_traldf_eiv and key_traldf_c2d -->
        <field id="aht2d" long_name="lateral eddy diffusivity" standard_name="ocean_tracer_xy_laplacian_diffusivity" unit="m2/s" grid_ref="grid_W_2D"/>
        <field id="aht2d_eiv" long_name="EIV lateral eddy diffusivity" standard_name="ocean_tracer_bolus_laplacian_diffusivity" unit="m2/s" grid_ref="grid_W_2D"/>

        <!-- variables from the TKE scheme with wave coupling  -->
        <field id="votke" long_name="Turbulent Kinetic Energy" unit="N/A"/>
        <field id="sotke" long_name="Surface Turbulent Kinetic Energy" unit="N/A" grid_ref="grid_W_2D" name="sotke" standard_name="Surface Turbulent Kinetic Energy"/>
        <field id="sotke" long_name="Surface Turbulent Kinetic Energy" unit="N/A" grid_ref="grid_W_2D"/>
        <field id="sotke0" long_name="Surface Turbulent Kinetic Energy From Waves" unit="N/A" grid_ref="grid_W_2D" name="sotke0" standard_name="Surface Turbulent Kinetic Energy From Waves"/>
        <field id="sotkedfac" long_name="Surface Turbulent Kinetic Energy Scaling From Waves" unit="N/A" grid_ref="grid_W_2D" name="sotkedfac" standard_name="Surface Turbulent Kinetic Energy Scaling From Waves"/>
	
         <!-- Variable from eosbn2 -->
         <field id="bn2" long_name="Squared buoyancy frequency measuring ocean vertical stratification" standard_name="square_of_brunt_vaisala_frequency_in_sea_water" unit="1/s^2"/>

        <!-- variables available with WAVE (ln_wave=T) -->
        <field id="wstokes"      long_name="Stokes Drift Velocity z-axis"      standard_name="StokesDrift_z_velocity"      unit="m/s"        grid_ref="grid_W_3D" />

      </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>
          
      <!-- 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_group>

      <!-- scalar variables available with key_diaar5 -->

      <field_group id="scalar" grid_ref="grid_1point">
         <field id="voltot" long_name="Sea Water 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 Average Thermosteric Sea Level Change" 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 Average Sea Water Potential Temperature" standard_name="sea_water_potential_temperature" unit="degree_C"/>
         <field id="saltot" long_name="Global Average Sea Water Salinity" standard_name="sea_water_salinity" unit="0.001"/>
         <field id="fram_trans" long_name="Sea Ice Mass Transport Through Fram Strait" standard_name="sea_ice_transport_across_line" unit="kg/s"/>
        <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="0.001"/>
    	 <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 (no vvl)" unit="degC"/>
    	 <field id="bgmissal" long_name="global mean salinity error due to free surface (no vvl)" unit="1e-3"/>
      </field_group>

      <!-- LIM3 scalar variables -->

      <field_group id="SBC_scalar" domain_ref="1point">
         <!-- available with ln_limdiaout -->
         <field id="ibgfrcvoltop" long_name="global mean ice/snow forcing at interface ice/snow-atm (volume equivalent ocean volume)" unit="km3"/>
         <field id="ibgfrcvolbot" long_name="global mean ice/snow forcing at interface ice/snow-ocean (volume equivalent ocean volume)" unit="km3"/>
         <field id="ibgfrctemtop" long_name="global mean heat on top of ice/snw/ocean-atm " unit="1e20J"/>
         <field id="ibgfrctembot" long_name="global mean heat below ice (on top of ocean) " unit="1e20J"/>
         <field id="ibgfrcsal" long_name="global mean ice/snow forcing (salt equivalent ocean volume)" unit="pss*km3"/>
         <field id="ibgfrchfxtop" long_name="global mean heat flux on top of ice/snw/ocean-atm " unit="W/m2"/>
         <field id="ibgfrchfxbot" long_name="global mean heat flux below ice (on top of ocean) " unit="W/m2"/>
 
         <field id="ibgvolume" long_name="drift in ice/snow volume (equivalent ocean volume)" unit="km3"/>
         <field id="ibgsaltco" long_name="drift in ice salt content (equivalent ocean volume)" unit="pss*km3"/>
         <field id="ibgheatco" long_name="drift in ice/snow heat content" unit="1e20J"/>
         <field id="ibgheatfx" long_name="drift in ice/snow heat flux" unit="W/m2"/>

         <field id="ibgvol_tot" long_name="global mean ice volume" unit="km3"/>
         <field id="sbgvol_tot" long_name="global mean snow volume" unit="km3"/>
         <field id="ibgarea_tot" long_name="global mean ice area" unit="km2"/>
         <field id="ibgsalt_tot" long_name="global mean ice salt content" unit="1e-3*km3"/>
         <field id="ibgheat_tot" long_name="global mean ice heat content" unit="1e20J"/>
         <field id="sbgheat_tot" long_name="global mean snow heat content" unit="1e20J"/>
      </field_group>
  
      <!-- variables available with key_float -->

      <field_group id="floatvar" domain_ref="1point" axis_ref="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="degree_C"/>
        <field id="traj_salt" long_name="floats salinity" standard_name="sea_water_practical_salinity" unit="0.001"/>
        <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="zomsfglo" long_name="Meridional Stream-Function: Global" unit="Sv" grid_ref="grid_znl_W_3D"/>
        <field id="zomsfatl" long_name="Meridional Stream-Function: Atlantic" unit="Sv" grid_ref="grid_znl_W_3D"/>
        <field id="zomsfpac" long_name="Meridional Stream-Function: Pacific" unit="Sv" grid_ref="grid_znl_W_3D"/>
        <field id="zomsfind" long_name="Meridional Stream-Function: Indian" unit="Sv" grid_ref="grid_znl_W_3D"/>
        <field id="zomsfipc" long_name="Meridional Stream-Function: Pacific+Indian" unit="Sv" grid_ref="grid_znl_W_3D"/>
        <field id="zotemglo" long_name="Zonal Mean Temperature : Global" unit="degree_C" grid_ref="grid_znl_T_3D"/>
        <field id="zotematl" long_name="Zonal Mean Temperature : Atlantic" unit="degree_C" grid_ref="grid_znl_T_3D"/>
        <field id="zotempac" long_name="Zonal Mean Temperature : Pacific" unit="degree_C" grid_ref="grid_znl_T_3D"/>
        <field id="zotemind" long_name="Zonal Mean Temperature : Indian" unit="degree_C" grid_ref="grid_znl_T_3D"/>
        <field id="zotemipc" long_name="Zonal Mean Temperature : Pacific+Indian" unit="degree_C" grid_ref="grid_znl_T_3D"/>
        <field id="zosalglo" long_name="Zonal Mean Salinity : Global" unit="0.001" grid_ref="grid_znl_T_3D"/>
        <field id="zosalatl" long_name="Zonal Mean Salinity : Atlantic" unit="0.001" grid_ref="grid_znl_T_3D"/>
        <field id="zosalpac" long_name="Zonal Mean Salinity : Pacific" unit="0.001" grid_ref="grid_znl_T_3D"/>
        <field id="zosalind" long_name="Zonal Mean Salinity : Indian" unit="0.001" grid_ref="grid_znl_T_3D"/>
        <field id="zosalipc" long_name="Zonal Mean Salinity : Pacific+Indian" unit="0.001" grid_ref="grid_znl_T_3D"/>
        <field id="zosrfglo" long_name="Zonal Mean Surface" unit="m2" grid_ref="grid_znl_T_3D"/>
        <field id="zosrfatl" long_name="Zonal Mean Surface : Atlantic" unit="m2" grid_ref="grid_znl_T_3D"/>
        <field id="zosrfpac" long_name="Zonal Mean Surface : Pacific" unit="m2" grid_ref="grid_znl_T_3D"/>
        <field id="zosrfind" long_name="Zonal Mean Surface : Indian" unit="m2" grid_ref="grid_znl_T_3D"/>
        <field id="zosrfipc" long_name="Zonal Mean Surface : Pacific+Indian" unit="m2" grid_ref="grid_znl_T_3D"/>
        <field id="sophtadv" long_name="Advective Heat Transport" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophtadv_atl" long_name="Advective Heat Transport: Atlantic" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophtadv_pac" long_name="Advective Heat Transport: Pacific" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophtadv_ind" long_name="Advective Heat Transport: Indian" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophtadv_ipc" long_name="Advective Heat Transport: Pacific+Indian" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophtldf" long_name="Diffusive Heat Transport" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophtldf_atl" long_name="Diffusive Heat Transport: Atlantic" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophtldf_pac" long_name="Diffusive Heat Transport: Pacific" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophtldf_ind" long_name="Diffusive Heat Transport: Indian" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophtldf_ipc" long_name="Diffusive Heat Transport: Pacific+Indian" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophtove" long_name="Overturning Heat Transport" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophtove_atl" long_name="Overturning Heat Transport: Atlantic" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophtove_pac" long_name="Overturning Heat Transport: Pacific" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophtove_ind" long_name="Overturning Heat Transport: Indian" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophtove_ipc" long_name="Overturning Heat Transport: Pacific+Indian" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophtbtr" long_name="Barotropic Heat Transport" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophtbtr_atl" long_name="Barotropic Heat Transport: Atlantic" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophtbtr_pac" long_name="Barotropic Heat Transport: Pacific" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophtbtr_ind" long_name="Barotropic Heat Transport: Indian" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophtbtr_ipc" long_name="Barotropic Heat Transport: Pacific+Indian" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophteiv" long_name="Heat Transport from mesoscale eddy advection" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophteiv_atl" long_name="Heat Transport from mesoscale eddy advection: Atlantic" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophteiv_pac" long_name="Heat Transport from mesoscale eddy advection: Pacific" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophteiv_ind" long_name="Heat Transport from mesoscale eddy advection: Indian" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sophteiv_ipc" long_name="Heat Transport from mesoscale eddy advection: Pacific+Indian" unit="PW" 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="sopstadv_atl" long_name="Advective Salt Transport: Atlantic" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopstadv_pac" long_name="Advective Salt Transport: Pacific" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopstadv_ind" long_name="Advective Salt Transport: Indian" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopstadv_ipc" long_name="Advective Salt Transport: Pacific+Indian" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopstove" long_name="Overturning Salt Transport" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopstove_atl" long_name="Overturning Salt Transport: Atlantic" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopstove_pac" long_name="Overturning Salt Transport: Pacific" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopstove_ind" long_name="Overturning Salt Transport: Indian" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopstove_ipc" long_name="Overturning Salt Transport: Pacific+Indian" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopstbtr" long_name="Barotropic Salt Transport" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopstbtr_atl" long_name="Barotropic Salt Transport: Atlantic" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopstbtr_pac" long_name="Barotropic Salt Transport: Pacific" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopstbtr_ind" long_name="Barotropic Salt Transport: Indian" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopstbtr_ipc" long_name="Barotropic Salt Transport: Pacific+Indian" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopstldf" long_name="Diffusive Salt Transport" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopstldf_atl" long_name="Diffusive Salt Transport: Atlantic" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopstldf_pac" long_name="Diffusive Salt Transport: Pacific" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopstldf_ind" long_name="Diffusive Salt Transport: Indian" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopstldf_ipc" long_name="Diffusive Salt Transport: Pacific+Indian" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopsteiv" long_name="Salt Transport from mesoscale eddy advection" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopsteiv_atl" long_name="Salt Transport from mesoscale eddy advection: Atlantic" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopsteiv_pac" long_name="Salt Transport from mesoscale eddy advection: Pacific" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopsteiv_ind" long_name="Salt Transport from mesoscale eddy advection: Indian" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopsteiv_ipc" long_name="Salt Transport from mesoscale eddy advection: Pacific+Indian" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>       
        <field id="sopht_vt" long_name="Heat Transport" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sopht_vt_atl" long_name="Heat Transport: Atlantic" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sopht_vt_ind" long_name="Heat Transport: Indian" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sopht_vt_pac" long_name="Heat Transport: Pacific" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sopht_vt_ipc" long_name="Heat Transport: Indo-Pacific" unit="PW" grid_ref="grid_znl_T_2D"/>
        <field id="sopst_vs" long_name="Salt Transport" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopst_vs_atl" long_name="Salt Transport: Atlantic" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopst_vs_ind" long_name="Salt Transport: Indian" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopst_vs_pac" long_name="Salt Transport: Pacific" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
        <field id="sopst_vs_ipc" long_name="Salt Transport: Indo-Pacific" unit="Giga g/s" grid_ref="grid_znl_T_2D"/>
      </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. (no-vvl)" unit="degC/s" grid_ref="grid_T_2D"/>
      <field id="strd_sad" long_name="salinity   -trend: surface adv. (no-vvl)" 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_x" long_name="temperature-trend: isopycnal diffusion x-component" unit="degC/s"/>
      <field id="strd_iso_x" long_name="salinity   -trend: isopycnal diffusion x-component" unit="1e-3/s"/>
      <field id="ttrd_iso_y" long_name="temperature-trend: isopycnal diffusion y-component" unit="degC/s"/>
      <field id="strd_iso_y" long_name="salinity   -trend: isopycnal diffusion y-component" unit="1e-3/s"/>
      <field id="ttrd_iso_z1" long_name="temperature-trend: isopycnal diffusion z-component part 1" unit="degC/s"/>
      <field id="strd_iso_z1" long_name="salinity   -trend: isopycnal diffusion z-component part 1" unit="1e-3/s"/>
      <field id="ttrd_iso_z" long_name="temperature-trend: isopycnal diffusion z-component" unit="degC/s"> ttrd_iso_z1 + ttrd_zdf - ttrd_zdfp </field>
      <field id="strd_iso_z" long_name="salinity   -trend: isopycnal diffusion z-component" unit="1e-3/s"> strd_iso_z1 + strd_zdf - strd_zdfp </field>
      <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"/>

      <!-- -->
      <field id="ttrd_dmp" long_name="temperature-trend: interior restoring" unit="degree_C/s"/>
      <field id="strd_dmp" long_name="salinity   -trend: interior restoring" unit="0.001/s"/>
      <field id="ttrd_bbl" long_name="temperature-trend: bottom boundary layer" unit="degree_C/s"/>
      <field id="strd_bbl" long_name="salinity   -trend: bottom boundary layer" unit="0.001/s"/>
      <field id="ttrd_npc" long_name="temperature-trend: non-penetrative conv." unit="degree_C/s"/>
      <field id="strd_npc" long_name="salinity   -trend: non-penetrative conv." unit="0.001/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="degree_C/s"/>
      <field id="ttrd_bbc" long_name="temperature-trend: geothermal heating" unit="degree_C/s"/>
      <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"/>
      <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_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>

      <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_x_e3t" unit="degC/s * m">  ttrd_iso_x * e3t </field>
      <field id="strd_iso_x_e3t" unit="1e-3/s * m">  strd_iso_x * e3t </field>
      <field id="ttrd_iso_y_e3t" unit="degC/s * m">  ttrd_iso_y * e3t </field>
      <field id="strd_iso_y_e3t" unit="1e-3/s * m">  strd_iso_y * e3t </field>
      <field id="ttrd_iso_z_e3t" unit="degC/s * m">  ttrd_iso_z * e3t </field>
      <field id="strd_iso_z_e3t" unit="1e-3/s * m">  strd_iso_z * e3t </field>
      <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>

      <!-- -->
      <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 * e3t_surf </field>
      <field id="strd_cdt_e3t" unit="degC/s * m">  strd_cdt * e3t_surf </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_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/m^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"/>   
      <field id="ketrd_eiv" long_name="Kinetic energy trend from mesoscale eddy advection" unit="W/s^2" grid_ref="grid_T_2D"/>   
      <!-- variables available when explicit mixing is used -->
      <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 key_traadv_eiv and key_diaeiv -->
      <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. (no-vvl)" 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="degree_C-1"/>   
      <field id="betaPE" long_name="partial deriv. of PEanom wrt S" unit="1e3"/>   
      <field id="petrd_zdfp_e3t" unit="W/m^2" grid_ref="grid_T_3D"> petrd_zdfp * e3t </field>
    </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>

    <!-- 
============================================================================================================
                                     Biogeochemistry model variables
============================================================================================================
    -->

      <!-- ptrc on T grid -->

     <field_group id="ptrc_T" grid_ref="grid_T_3D">
       <field id="DIC" long_name="Dissolved inorganic Concentration" unit="mmol/m3"/>
       <field id="DIC_E3T" long_name="DIC * E3T" unit="mmol/m2"> DIC * e3t </field>
       <field id="Alkalini" long_name="Total Alkalinity Concentration" unit="mmol/m3"/>
       <field id="Alkalini_E3T" long_name="Alkalini * E3T" unit="mmol/m2"> Alkalini * e3t </field>
       <field id="O2" long_name="Oxygen Concentration" unit="mmol/m3"/>
       <field id="O2_E3T" long_name="O2 * E3T" unit="mmol/m2"> O2 * e3t </field>
       <field id="CaCO3" long_name="Calcite Concentration" unit="mmol/m3"/>
       <field id="CaCO3_E3T" long_name="CaCO3 * E3T" unit="mmol/m2"> CaCO3 * e3t </field>
       <field id="PO4" long_name="Phosphate Concentration" unit="mmol/m3"/>
       <field id="PO4_E3T" long_name="PO4 * E3T" unit="mmol/m2"> PO4 * e3t </field>
       <field id="POC" long_name="Small organic carbon Concentration" unit="mmol/m3"/>
       <field id="POC_E3T" long_name="POC * E3T" unit="mmol/m2"> POC * e3t </field>
       <field id="Si" long_name="Silicate Concentration" unit="mmol/m3"/>
       <field id="Si_E3T" long_name="Si * E3T" unit="mmol/m2"> Si * e3t </field>
       <field id="PHY" long_name="(Nano)Phytoplankton Concentration" unit="mmol/m3"/>
       <field id="PHY_E3T" long_name="PHY * E3T" unit="mmol/m2"> PHY * e3t </field>
       <field id="ZOO" long_name="(Micro)Zooplankton Concentration" unit="mmol/m3"/>
       <field id="ZOO_E3T" long_name="ZOO2 * E3T" unit="mmol/m2"> ZOO * e3t </field>
       <field id="DOC" long_name="Dissolved organic Concentration" unit="mmol/m3"/>
       <field id="DOC_E3T" long_name="DOC * E3T" unit="mmol/m2"> DOC * e3t </field>
       <field id="PHY2" long_name="Diatoms Concentration" unit="mmol/m3"/>
       <field id="PHY2_E3T" long_name="PHY2 * E3T" unit="mmol/m2"> PHY2 * e3t </field>
       <field id="ZOO2" long_name="Mesozooplankton Concentration" unit="mmol/m3"/>
       <field id="ZOO2_E3T" long_name="ZOO2 * E3T" unit="mmol/m2"> ZOO2 * e3t </field>
       <field id="DSi" long_name="Diatoms Silicate Concentration" unit="mmol/m3"/>
       <field id="DSi_E3T" long_name="Dsi * E3T" unit="mmol/m2"> DSi * e3t </field>
       <field id="Fer" long_name="Dissolved Iron Concentration" unit="mmol/m3"/>
       <field id="Fer_E3T" long_name="Fer * E3T" unit="mmol/m2"> Fer * e3t </field>
       <field id="BFe" long_name="Big iron particles Concentration" unit="mmol/m3"/>
       <field id="BFe_E3T" long_name="BFe * E3T" unit="mmol/m2"> BFe * e3t </field>
       <field id="GOC" long_name="Big organic carbon Concentration" unit="mmol/m3"/>
       <field id="GOC_E3T" long_name="GOC * E3T" unit="mmol/m2"> GOC * e3t </field>
       <field id="SFe" long_name="Small iron particles Concentration" unit="mmol/m3"/>
       <field id="SFe_E3T" long_name="SFe * E3T" unit="mmol/m2"> SFe * e3t </field>
       <field id="DFe" long_name="Diatoms iron  Concentration" unit="mmol/m3"/>
       <field id="DFe_E3T" long_name="DFe * E3T" unit="mmol/m2"> DFe * e3t </field>
       <field id="GSi" long_name="Sinking biogenic Silicate Concentration" unit="mmol/m3"/>
       <field id="GSi_E3T" long_name="GSi * E3T" unit="mmol/m2"> GSi * e3t </field>
       <field id="NFe" long_name="Nano iron Concentration" unit="mmol/m3"/>
       <field id="NFe_E3T" long_name="NFe * E3T" unit="mmol/m2"> NFe * e3t </field>
       <field id="NCHL" long_name="Nano chlorophyl Concentration" unit="mg/m3"/>
       <field id="NCHL_E3T" long_name="NCHL * E3T" unit="mmol/m2"> NCHL * e3t </field>
       <field id="DCHL" long_name="Diatoms chlorophyl Concentration" unit="mg/m3"/>
       <field id="DCHL_E3T" long_name="DCHL * E3T" unit="mmol/m2"> DCHL * e3t </field>
       <field id="NO3" long_name="Nitrate Concentration" unit="mmol/m3"/>
       <field id="NO3_E3T" long_name="NO3 * E3T" unit="mmol/m2"> NO3 * e3t </field>
       <field id="NH4" long_name="Ammonium Concentration" unit="mmol/m3"/>
       <field id="NH4_E3T" long_name="NH4 * E3T" unit="mmol/m2"> NH4 * e3t </field>

       <!-- PISCES with Kriest parametisation : variables available with key_kriest -->
       <field id="Num" long_name="Number of organic particles" unit="1"/>
       <field id="Num_E3T" long_name="Num * E3T" unit="m"> Num * e3t </field>

       <!-- PISCES light : variables available with key_pisces_reduced -->
       <field id="DET" long_name="Detritus" unit="mmol-N/m3"/>
       <field id="DET_E3T" long_name="DET * E3T" unit="mmol-N/m2"> DET * e3t </field>
       <field id="DOM" long_name="Dissolved Organic Matter" unit="mmol-N/m3"/>
       <field id="DOM_E3T" long_name="DOM * E3T" unit="mmol-N/m2"> DOM * e3t </field>

       <!-- CFC11 : variables available with key_cfc -->
       <field id="CFC11" long_name="CFC-11 Concentration" unit="umol/m3"/>
       <field id="CFC11_E3T" long_name="CFC11 * E3T" unit="umol/m2"> CFC11 * e3t </field>
       <field id="CFC12" long_name="CFC-12 Concentration" unit="umol/m3"/>
       <field id="CFC12_E3T" long_name="CFC12 * E3T" unit="umol/m2"> CFC12 * e3t </field>
       <field id="SF6" long_name="SF6 Concentration" unit="umol/m3"/>
       <field id="SF6_E3T" long_name="SF6 * E3T" unit="umol/m2"> SF6 * e3t </field>

       <!-- Bomb C14 : variables available with key_c14b -->
       <field id="C14B" long_name="Bomb C14 Concentration" unit="1"/>
       <field id="C14B_E3T" long_name="C14B * E3T" unit="m"> C14B * e3t </field>

       <!-- AGE : variables available with key_age -->
       <field id="Age" long_name="Sea Water Age Since Surface contact" standard_name="sea_water_age_since_surface_contact" unit="yr"/>
       <field id="Age_e3t" long_name="Age * e3t" unit="yr * m"> Age * e3t </field>

     </field_group>

     <!-- PISCES additional diagnostics on T grid  -->

     <field_group id="diad_T" grid_ref="grid_T_2D">
       <field id="PH" long_name="PH" unit="1" grid_ref="grid_T_3D"/>
       <field id="CO3" long_name="Bicarbonates" unit="mol/m3" grid_ref="grid_T_3D"/>
       <field id="CO3sat" long_name="CO3 saturation" unit="mol/m3" grid_ref="grid_T_3D"/>
       <field id="PAR" long_name="Photosynthetically Available Radiation" unit="W/m2" grid_ref="grid_T_3D"/>
       <field id="PARDM" long_name="Daily mean PAR" unit="W/m2" grid_ref="grid_T_3D"/>
       <field id="PPPHY" long_name="Primary production of nanophyto" unit="mol/m3/s" grid_ref="grid_T_3D"/>
       <field id="PPPHY2" long_name="Primary production of diatoms" unit="mol/m3/s" grid_ref="grid_T_3D"/>
       <field id="PPNEWN" long_name="New Primary production of nanophyto" unit="mol/m3/s" grid_ref="grid_T_3D"/>
       <field id="PPNEWD" long_name="New Primary production of diatoms" unit="mol/m3/s" grid_ref="grid_T_3D"/>
       <field id="PBSi" long_name="Primary production of Si diatoms" unit="mol/m3/s" grid_ref="grid_T_3D"/>
       <field id="PFeN" long_name="Primary production of nano iron" unit="mol/m3/s" grid_ref="grid_T_3D"/>
       <field id="PFeD" long_name="Primary production of diatoms iron" unit="mol/m3/s" grid_ref="grid_T_3D"/>
       <field id="xfracal" long_name="Calcifying fraction" unit="1" grid_ref="grid_T_3D"/>
       <field id="PCAL" long_name="Calcite production" unit="mol/m3/s" grid_ref="grid_T_3D"/>
       <field id="DCAL" long_name="Calcite dissolution" unit="mol/m3/s" grid_ref="grid_T_3D"/>
       <field id="GRAZ1" long_name="Grazing by microzooplankton" unit="mol/m3/s" grid_ref="grid_T_3D"/>
       <field id="GRAZ2" long_name="Grazing by mesozooplankton" unit="mol/m3/s" grid_ref="grid_T_3D"/>
       <field id="REMIN" long_name="Oxic remineralization of OM" unit="mol/m3/s" grid_ref="grid_T_3D"/>
       <field id="DENIT" long_name="Anoxic remineralization of OM" unit="mol/m3/s" grid_ref="grid_T_3D"/>
       <field id="Nfix" long_name="Nitrogen fixation" unit="mol/m3/s" grid_ref="grid_T_3D"/>
       <field id="Mumax" long_name="Maximum growth rate" unit="s-1" grid_ref="grid_T_3D"/>
       <field id="MuN" long_name="Realized growth rate for nanophyto" unit="s-1" grid_ref="grid_T_3D"/>
       <field id="MuD" long_name="Realized growth rate for diatomes" unit="s-1" grid_ref="grid_T_3D"/>
       <field id="LNnut" long_name="Nutrient limitation term in Nanophyto" unit="" grid_ref="grid_T_3D"/>
       <field id="LDnut" long_name="Nutrient limitation term in Diatoms" unit="" grid_ref="grid_T_3D"/>
       <field id="LNFe" long_name="Iron limitation term in Nanophyto" unit="" grid_ref="grid_T_3D"/>
       <field id="LDFe" long_name="Iron limitation term in Diatoms" unit="" grid_ref="grid_T_3D"/>
       <field id="LNlight" long_name="Light limitation term in Nanophyto" unit="" grid_ref="grid_T_3D"/>
       <field id="LDlight" long_name="Light limitation term in Diatoms" unit="" grid_ref="grid_T_3D"/>
       <field id="Fe2" long_name="Iron II concentration" unit="nmol/m3" grid_ref="grid_T_3D"/>
       <field id="Fe3" long_name="Iron III concentration" unit="nmol/m3" grid_ref="grid_T_3D"/>
       <field id="FeL1" long_name="Complexed Iron concentration with L1" unit="nmol/m3" grid_ref="grid_T_3D"/>
       <field id="FeL2" long_name="Complexed Iron concentration with L2" unit="nmol/m3" grid_ref="grid_T_3D"/>
       <field id="FeP" long_name="Precipitated Iron III" unit="nmol/m3" grid_ref="grid_T_3D"/>
       <field id="TL1" long_name="Total L1 concentration" unit="nmol/m3" grid_ref="grid_T_3D"/>
       <field id="TL2" long_name="Total L2 concentration" unit="nmol/m3" grid_ref="grid_T_3D"/>
       <field id="pdust" long_name="dust concentration" unit="g/m3"/>
       <field id="Totlig" long_name="Total ligand concentation" unit="nmol/m3" grid_ref="grid_T_3D"/>
       <field id="Biron" long_name="Bioavailable iron" unit="nmol/m3" grid_ref="grid_T_3D"/>
       <field id="Sdenit" long_name="Nitrate reduction in the sediments" unit="molN/m2/s"/>
       <field id="SedCal" long_name="Calcite burial in the sediments" unit="molC/m2/s"/>
       <field id="SedSi" long_name="Silicon burial in the sediments" unit="molSi/m2/s"/>
       <field id="SedC" long_name="Organic C burial in the sediments" unit="molC/m2/s"/>
       <field id="Ironice" long_name="Iron input/uptake due to sea ice" unit="mol/m2/s"/>
       <field id="HYDR" long_name="Iron input from hydrothemal vents" unit="mol/m2/s" grid_ref="grid_T_3D"/>
       <field id="EPC100" long_name="Export of carbon particles at 100 m" unit="mol/m2/s"/>
       <field id="EPFE100" long_name="Export of biogenic iron at 100 m" unit="mol/m2/s"/>
       <field id="EPSI100" long_name="Export of Silicate at 100 m" unit="mol/m2/s"/>
       <field id="EPCAL100" long_name="Export of Calcite at 100 m" unit="mol/m2/s"/>
       <field id="EXPC" long_name="Export of carbon" unit="mol/m2/s" grid_ref="grid_T_3D"/>
       <field id="EXPFE" long_name="Export of biogenic iron" unit="mol/m2/s" grid_ref="grid_T_3D"/>
       <field id="EXPSI" long_name="Export of Silicate" unit="mol/m2/s" grid_ref="grid_T_3D"/>
       <field id="EXPCAL" long_name="Export of Calcite" unit="mol/m2/s" grid_ref="grid_T_3D"/>
       <field id="Cflx" long_name="DIC flux" unit="mol/m2/s"/>
       <field id="Oflx" long_name="Oxygen flux" unit="mol/m2/s"/>
       <field id="Kg" long_name="Gas transfer" unit="mol/m2/s/uatm"/>
       <field id="Dpco2" long_name="Delta CO2" unit="uatm"/>
       <field id="Dpo2" long_name="Delta O2" unit="uatm"/>
       <field id="Heup" long_name="Euphotic layer depth" unit="m"/>
       <field id="Irondep" long_name="Iron deposition from dust" unit="mol/m2/s"/>
       <field id="Ironsed" long_name="Iron deposition from sediment" unit="mol/m2/s" grid_ref="grid_T_3D"/>


       <!-- PISCES with Kriest parametisation : variables available with key_kriest -->
       <field id="EPN100" long_name="Particulate number flux at 100 m" unit="mol/m2/s"/>
       <field id="EXPN" long_name="Particulate number flux" unit="mol/m2/s" grid_ref="grid_T_3D"/>
       <field id="XNUM" long_name="Number of particles in aggregats" unit="1" grid_ref="grid_T_3D"/>
       <field id="WSC" long_name="sinking speed of mass flux" unit="m2/s" grid_ref="grid_T_3D"/>
       <field id="WSN" long_name="sinking speed of number flux" unit="m2/s" grid_ref="grid_T_3D"/>

       <!-- dbio_T on T grid : variables available with key_diaar5 -->
       <field id="TPP" long_name="Total Primary production of phyto" unit="mol/m3/s" grid_ref="grid_T_3D"/>
       <field id="TPNEW" long_name="New Primary production of phyto" unit="mol/m3/s" grid_ref="grid_T_3D"/>
       <field id="TPBFE" long_name="Total biogenic iron production" unit="mol/m3/s" grid_ref="grid_T_3D"/>
       <field id="INTDIC" long_name="DIC content" unit="kg/m2"/>
       <field id="O2MIN" long_name="Oxygen minimum concentration" unit="mol/m3"/>
       <field id="ZO2MIN" long_name="Depth of oxygen minimum concentration" unit="m"/>
       <field id="INTNFIX" long_name="Nitrogen fixation rate : vert. integrated" unit="mol/m2/s"/>
       <field id="INTPPPHY" long_name="Vertically integrated primary production by nanophy" unit="mol/m2/s"/>
       <field id="INTPPPHY2" long_name="Vertically integrated primary production by diatom" unit="mol/m2/s"/>
       <field id="INTPP" long_name="Vertically integrated primary production by phyto" unit="mol/m2/s"/>
       <field id="INTPNEW" long_name="Vertically integrated new primary production" unit="mol/m2/s"/>
       <field id="INTPBFE" long_name="Vertically integrated of biogenic iron production" unit="mol/m2/s"/>
       <field id="INTPBSI" long_name="Vertically integrated of biogenic Si production" unit="mol/m2/s"/>
       <field id="INTPCAL" long_name="Vertically integrated of calcite production" unit="mol/m2/s"/>

       <!-- PISCES light : variables available with key_pisces_reduced -->
       <field id="FNO3PHY" long_name="FNO3PHY" unit="" grid_ref="grid_T_3D"/> 
       <field id="FNH4PHY" long_name="FNH4PHY" unit="" grid_ref="grid_T_3D"/> 
       <field id="FNH4NO3" long_name="FNH4NO3" unit="" grid_ref="grid_T_3D"/> 
       <field id="TNO3PHY" long_name="TNO3PHY" unit=""/> 
       <field id="TNH4PHY" long_name="TNH4PHY" unit=""/> 
       <field id="TPHYDOM" long_name="TPHYDOM" unit=""/> 
       <field id="TPHYNH4" long_name="TPHYNH4" unit=""/> 
       <field id="TPHYZOO" long_name="TPHYZOO" unit=""/> 
       <field id="TPHYDET" long_name="TPHYDET" unit=""/> 
       <field id="TDETZOO" long_name="TDETZOO" unit=""/> 
       <field id="TZOODET" long_name="TZOODET" unit=""/> 
       <field id="TZOOBOD" long_name="TZOOBOD" unit=""/> 
       <field id="TZOONH4" long_name="TZOONH4" unit=""/> 
       <field id="TZOODOM" long_name="TZOODOM" unit=""/> 
       <field id="TNH4NO3" long_name="TNH4NO3" unit=""/> 
       <field id="TDOMNH4" long_name="TDOMNH4" unit=""/> 
       <field id="TDETNH4" long_name="TDETNH4" unit=""/> 
       <field id="TPHYTOT" long_name="TPHYTOT" unit=""/> 
       <field id="TZOOTOT" long_name="TZOOTOT" unit=""/> 
       <field id="SEDPOC" long_name="SEDPOC" unit=""/> 
       <field id="TDETSED" long_name="TDETSED" unit=""/> 

       <!-- CFC11 : variables available with key_cfc -->
       <field id="qtr_CFC11" long_name="Air-sea flux of CFC-11" unit="mol/m2/s"/>
       <field id="qint_CFC11" long_name="Cumulative air-sea flux of CFC-11" unit="mol/m2"/>
       <field id="qtr_CFC12" long_name="Air-sea flux of CFC-12" unit="mol/m2/s"/>
       <field id="qint_CFC12" long_name="Cumulative air-sea flux of CFC-12" unit="mol/m2"/>
       <field id="qtr_SF6" long_name="Air-sea flux of SF6" unit="mol/m2/s"/>
       <field id="qint_SF6" long_name="Cumulative air-sea flux of SF6" unit="mol/m2"/>

       <!-- Bomb C14 : variables available with key_c14b -->
       <field id="qtr_C14b" long_name="Air-sea flux of Bomb C14" unit="mol/m2/s"/>
       <field id="qint_C14b" long_name="Cumulative air-sea flux of Bomb C14" unit="mol/m2"/>
       <field id="fdecay" long_name="Radiactive decay of Bomb C14" unit="mol/m3" grid_ref="grid_T_3D"/>
     </field_group>

     <field_group id="PISCES_scalar" domain_ref="1point">
       <field id="pno3tot" long_name="global mean nitrate concentration" unit="mol/m3"/>
       <field id="ppo4tot" long_name="global mean phosphorus concentration" unit="mol/m3"/>
       <field id="psiltot" long_name="global mean silicate concentration" unit="mol/m3"/>
       <field id="palktot" long_name="global mean alkalinity concentration" unit="mol/m3"/>
       <field id="pfertot" long_name="global mean iron concentration" unit="mol/m3"/>
       <field id="tcflx" long_name="total Flux of Carbon out of the ocean" unit="mol/s"/>
       <field id="tcflxcum" long_name="cumulative total Flux of Carbon out of the ocean" unit="mol/s"/>
       <field id="tcexp" long_name="total Carbon export at 100m" unit="mol/s"/>
       <field id="tintpp" long_name="global total integrated primary production" unit="mol/s"/>
       <field id="tnfix" long_name="global total nitrogen fixation" unit="mol/s"/>
       <field id="tdenit" long_name="Total denitrification" unit="mol/s"/>
     </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">
      <field field_ref="toce" name="thetao" long_name="sea_water_potential_temperature"/>
      <field field_ref="soce" name="so" long_name="sea_water_salinity"/>
      <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"/>
   
      <field field_ref="sst" name="tos" long_name="sea_surface_temperature"/>
      <field field_ref="sst2" name="tossq" long_name="square_of_sea_surface_temperature"/>
      <field field_ref="sstgrad" name="tosgrad" long_name="module_of_sea_surface_temperature_gradient"/>
      <field field_ref="sss" name="sos" long_name="sea_surface_salinity"/>
      <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">
      <field field_ref="toce" name="thetao" long_name="sea_water_potential_temperature"/>
      <field field_ref="soce" name="so" long_name="sea_water_salinity"/>
      <field field_ref="sst" name="tos" long_name="sea_surface_temperature"/>
      <field field_ref="sst2" name="tossq" long_name="square_of_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"/>
      <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="ssu" name="uos" long_name="sea_surface_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="ssv" name="vos" long_name="sea_surface_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>
    

   </field_definition>
