Skip to content
Snippets Groups Projects

Compare revisions

Changes are shown as if the source revision was being merged into the target revision. Learn more about comparing revisions.

Source

Select target project
No results found

Target

Select target project
  • nemo/nemo
  • sparonuz/nemo
  • hatfield/nemo
  • extdevs/nemo
4 results
Show changes
Hide whitespace changes
Inline Side-by-side
Showing
with 974 additions and 1362 deletions
cfgs/ORCA2_OFF_TRC/EXPREF/context_nemo.xml
View file @ dd59ced6
......@@ -4,7 +4,7 @@
==============================================================================================
-->
<context id="nemo">
<!-- $id$ -->
<variable_definition>
<!-- Year/Month/Day of time origin for NetCDF files; defaults to 1800-01-01 -->
<variable id="ref_year" type="int"> 1900 </variable>
......@@ -26,7 +26,6 @@
<!-- Files definition -->
<file_definition src="./file_def_nemo-innerttrc.xml"/> <!-- NEMO ocean passive tracers -->
<!-- Axis definition -->
<axis_definition src="./axis_def_nemo.xml"/>
......
cfgs/ORCA2_OFF_TRC/EXPREF/file_def_nemo-innerttrc.xml
View file @ dd59ced6
......@@ -9,45 +9,56 @@
<file_definition type="one_file" name="@expname@_@freq@_@startdate@_@enddate@" sync_freq="10d" min_digits="4">
<file_group id="1ts" output_freq="1ts" output_level="10" enabled=".TRUE."/> <!-- 1 time step files -->
<file_group id="1h" output_freq="1h" output_level="10" enabled=".TRUE."/> <!-- 1h files -->
<file_group id="2h" output_freq="2h" output_level="10" enabled=".TRUE."/> <!-- 2h files -->
<file_group id="3h" output_freq="3h" output_level="10" enabled=".TRUE."/> <!-- 3h files -->
<file_group id="4h" output_freq="4h" output_level="10" enabled=".TRUE."/> <!-- 4h files -->
<file_group id="6h" output_freq="6h" output_level="10" enabled=".TRUE."/> <!-- 6h files -->
<file_group id="trc_1ts" output_freq="1ts" output_level="10" enabled=".TRUE."/> <!-- 1 time step files -->
<file_group id="trc_1h" output_freq="1h" output_level="10" enabled=".TRUE."/> <!-- 1h files -->
<file_group id="trc_2h" output_freq="2h" output_level="10" enabled=".TRUE."/> <!-- 2h files -->
<file_group id="trc_3h" output_freq="3h" output_level="10" enabled=".TRUE."/> <!-- 3h files -->
<file_group id="trc_4h" output_freq="4h" output_level="10" enabled=".TRUE."/> <!-- 4h files -->
<file_group id="trc_6h" output_freq="6h" output_level="10" enabled=".TRUE."/> <!-- 6h files -->
<file_group id="1d" output_freq="1d" output_level="10" enabled=".TRUE."> <!-- 1d files -->
<file_group id="trc_1d" output_freq="1d" output_level="10" enabled=".TRUE."> <!-- 1d files -->
</file_group>
<file_group id="3d" output_freq="3d" output_level="10" enabled=".TRUE."/> <!-- 3d files -->
<file_group id="5d" output_freq="5d" output_level="10" enabled=".TRUE."/> <!-- 5d files -->
<file_group id="trc_3d" output_freq="3d" output_level="10" enabled=".TRUE."/> <!-- 3d files -->
<file_group id="trc_5d" output_freq="5d" output_level="10" enabled=".TRUE."/> <!-- 5d files -->
<file_group id="1m" output_freq="1mo" output_level="10" enabled=".TRUE."> <!-- real monthly files -->
<file_group id="trc_1m" output_freq="1mo" output_level="10" enabled=".TRUE."> <!-- real monthly files -->
</file_group>
<file_group id="2m" output_freq="2mo" output_level="10" enabled=".TRUE."/> <!-- real 2m files -->
<file_group id="3m" output_freq="3mo" output_level="10" enabled=".TRUE."/> <!-- real 3m files -->
<file_group id="4m" output_freq="4mo" output_level="10" enabled=".TRUE."/> <!-- real 4m files -->
<file_group id="6m" output_freq="6mo" output_level="10" enabled=".TRUE."/> <!-- real 6m files -->
<file_group id="trc_2m" output_freq="2mo" output_level="10" enabled=".TRUE."/> <!-- real 2m files -->
<file_group id="trc_3m" output_freq="3mo" output_level="10" enabled=".TRUE."/> <!-- real 3m files -->
<file_group id="trc_4m" output_freq="4mo" output_level="10" enabled=".TRUE."/> <!-- real 4m files -->
<file_group id="trc_6m" output_freq="6mo" output_level="10" enabled=".TRUE."/> <!-- real 6m files -->
<file_group id="1y" output_freq="1y" output_level="10" enabled=".TRUE."> <!-- real yearly files -->
<file_group id="trc_1y" output_freq="1y" output_level="10" enabled=".TRUE."> <!-- real yearly files -->
<file id="file1" name_suffix="grid_T" description="passive tracers variables" >
<file id="file201" name_suffix="grid_T" description="passive tracers variables" >
<field field_ref="e3t" />
<field field_ref="ssh" name="zos" />
</file>
<file id="file1" name_suffix="_trc" description="passive tracers variables" >
<file id="file202" name_suffix="_trc" description="passive tracers variables" >
<field field_ref="Age" name="Age" operation="average" freq_op="1y" > @Age_e3t / @e3t </field>
<field field_ref="CFC11" name="CFC11" operation="average" freq_op="1y" > @CFC11_e3t / @e3t </field>
<field field_ref="CFC12" name="CFC12" operation="average" freq_op="1y" > @CFC12_e3t / @e3t </field>
<field field_ref="SF6" name="SF6" operation="average" freq_op="1y" > @SF6_e3t / @e3t </field>
<field field_ref="RC14" name="RC14" operation="average" freq_op="1y" > @RC14_e3t / @e3t </field>
<field field_ref="RC14" name="RC14" operation="average" freq_op="1y" > @RC14_e3t / @e3t </field>
<field field_ref="qtr_CFC11" />
<field field_ref="qint_CFC11" />
<field field_ref="qtr_CFC12" />
<field field_ref="qint_CFC12" />
<field field_ref="qtr_SF6" />
<field field_ref="qint_SF6" />
<field field_ref="qtr_c14" />
<field field_ref="qint_c14" />
<field field_ref="DeltaC14" />
<field field_ref="C14Age" />
<field field_ref="RAge" />
</file>
</file_group>
<file_group id="2y" output_freq="2y" output_level="10" enabled=".TRUE."/> <!-- real 2y files -->
<file_group id="5y" output_freq="5y" output_level="10" enabled=".TRUE."/> <!-- real 5y files -->
<file_group id="10y" output_freq="10y" output_level="10" enabled=".TRUE."/> <!-- real 10y files -->
<file_group id="trc_2y" output_freq="2y" output_level="10" enabled=".TRUE."/> <!-- real 2y files -->
<file_group id="trc_5y" output_freq="5y" output_level="10" enabled=".TRUE."/> <!-- real 5y files -->
<file_group id="trc_10y" output_freq="10y" output_level="10" enabled=".TRUE."/> <!-- real 10y files -->
</file_definition>
cfgs/ORCA2_OFF_TRC/EXPREF/namelist_cfg
View file @ dd59ced6
......@@ -10,10 +10,7 @@
!! namcfg parameters of the configuration (default: user defined GYRE)
!! namwad Wetting and drying (default: OFF)
!! namtsd data: temperature & salinity (default: OFF)
!! namcrs coarsened grid (for outputs and/or TOP) (ln_crs =T)
!! namc1d 1D configuration options (ln_c1d =T)
!! namc1d_dyndmp 1D newtonian damping applied on currents (ln_c1d =T)
!! namc1d_uvd 1D data (currents) (ln_c1d =T)
!!======================================================================
!
!-----------------------------------------------------------------------
......@@ -55,21 +52,9 @@
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namcrs ! coarsened grid (for outputs and/or TOP) (ln_crs =T)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namc1d ! 1D configuration options (ln_c1d =T default: PAPA station)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namc1d_dyndmp ! U & V newtonian damping (ln_c1d =T default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namc1d_uvd ! data: U & V currents (ln_c1d =T default: OFF)
!-----------------------------------------------------------------------
/
!!======================================================================
!! *** Surface Boundary Condition namelists *** !!
!! !!
......@@ -324,7 +309,7 @@
sn_mld = 'dyna_grid_T' , 120. , 'mldr10_1' , .true. , .true. , 'yearly' , '' , '' , ''
sn_emp = 'dyna_grid_T' , 120. , 'wfo' , .true. , .true. , 'yearly' , '' , '' , ''
sn_empb = 'dyna_grid_T' , 120. , 'wfob' , .true. , .true. , 'yearly' , '' , '' , ''
sn_fmf = 'dyna_grid_T' , 120. , 'fmmflx' , .true. , .true. , 'yearly' , '' , '' , ''
sn_fwf = 'dyna_grid_T' , 120. , 'iowaflup' , .true. , .true. , 'yearly' , '' , '' , ''
sn_rnf = 'dyna_grid_T' , 120. , 'runoffs' , .true. , .true. , 'yearly' , '' , '' , ''
sn_ice = 'dyna_grid_T' , 120. , 'siconc' , .true. , .true. , 'yearly' , '' , '' , ''
sn_qsr = 'dyna_grid_T' , 120. , 'rsntds' , .true. , .true. , 'yearly' , '' , '' , ''
......@@ -376,13 +361,10 @@
!! *** Diagnostics namelists *** !!
!! !!
!! namtrd dynamics and/or tracer trends (default: OFF)
!! namhsb Heat and salt budgets (default: OFF)
!! namdiu Cool skin and warm layer models (default: OFF)
!! namdiu Cool skin and warm layer models (default: OFF)
!! namflo float parameters (default: OFF)
!! namdct transports through some sections ("key_diadct")
!! nam_dia25h 25h Mean Output (default: OFF)
!! namnc4 netcdf4 chunking and compression settings ("key_netcdf4")
!! namnc4 netcdf4 chunking and compression settings
!!======================================================================
!
!-----------------------------------------------------------------------
......@@ -390,27 +372,15 @@
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namhsb ! Heat and salt budgets (default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namdiu ! Cool skin and warm layer models (default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namflo ! float parameters (default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namdct ! transports through some sections ("key_diadct")
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&nam_dia25h ! 25h Mean Output (default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namnc4 ! netcdf4 chunking and compression settings ("key_netcdf4")
&namnc4 ! netcdf4 chunking and compression settings
!-----------------------------------------------------------------------
/
!!======================================================================
......
cfgs/ORCA2_OFF_TRC/EXPREF/namelist_top_cfg
View file @ dd59ced6
!-----------------------------------------------------------------------
&namtrc_run ! run information
!-----------------------------------------------------------------------
ln_top_euler = .true.
/
!-----------------------------------------------------------------------
&namtrc ! tracers definition
......
cfgs/ORCA2_OFF_TRC/cpp_ORCA2_OFF_TRC.fcm
View file @ dd59ced6
bld::tool::fppkeys key_top key_xios key_qco
bld::tool::fppkeys key_top key_xios key_qco key_vco_1d3d
cfgs/ORCA2_SAS_ICE/EXPREF/context_nemo.xml
View file @ dd59ced6
......@@ -4,7 +4,7 @@
==============================================================================================
-->
<context id="nemo">
<!-- $id$ -->
<variable_definition>
<!-- Year/Month/Day of time origin for NetCDF files; defaults to 1800-01-01 -->
<variable id="ref_year" type="int"> 1900 </variable>
......@@ -26,7 +26,6 @@
<file_definition src="./file_def_nemo-oce.xml"/> <!-- NEMO ocean dynamics -->
<file_definition src="./file_def_nemo-ice.xml"/> <!-- NEMO sea-ice model -->
<!-- Axis definition -->
<axis_definition src="./axis_def_nemo.xml"/>
......
cfgs/ORCA2_SAS_ICE/EXPREF/file_def_nemo-ice.xml
View file @ dd59ced6
<?xml version="1.0"?>
<!-- $id$ -->
<!--
============================================================================================================
......@@ -11,8 +10,8 @@
<file_definition type="one_file" name="@expname@_@freq@_@startdate@_@enddate@" sync_freq="1mo" min_digits="4">
<file_group id="5d" output_freq="5d" output_level="10" enabled=".TRUE."> <!-- 5d files -->
<file id="file21" name_suffix="_icemod" description="ice variables" enabled=".true." >
<file_group id="ice_5d" output_freq="5d" output_level="10" enabled=".TRUE."> <!-- 5d files -->
<file id="file101" name_suffix="_icemod" description="ice variables" enabled=".true." >
<!-- ice mask -->
<field field_ref="icemask" name="simsk" />
......@@ -48,11 +47,18 @@
<!-- rheology -->
<field field_ref="icediv" name="sidive" />
<field field_ref="iceshe" name="sishea" />
<field field_ref="icedlt" name="sidelt" />
<field field_ref="icestr" name="sistre" />
<field field_ref="normstr" name="normstr" />
<field field_ref="sheastr" name="sheastr" />
<field field_ref="sig1_pnorm" name="sig1_pnorm"/>
<field field_ref="sig2_pnorm" name="sig2_pnorm"/>
<field field_ref="velo_res" name="velo_res" />
<field field_ref="velo_ero" name="velo_ero" />
<field field_ref="uice_eri" name="uice_eri" />
<field field_ref="vice_eri" name="vice_eri" />
<field field_ref="uice_cvg" name="uice_cvg" />
<!-- heat fluxes -->
<field field_ref="qt_oce_ai" name="qt_oce_ai" />
......@@ -64,7 +70,7 @@
<field field_ref="qns_ice" name="qns_ice" />
<field field_ref="qemp_ice" name="qemp_ice" />
<field field_ref="albedo" name="albedo" />
<field field_ref="hfxcndtop" name="hfxcndtop" />
<field field_ref="hfxcndbot" name="hfxcndbot" />
<field field_ref="hfxsensib" name="hfxsensib" />
......@@ -86,7 +92,7 @@
</file>
<file id="file22" name_suffix="_SBC_scalar" description="scalar variables" enabled=".true." >
<file id="file102" name_suffix="_SBC_scalar" description="scalar variables" enabled=".true." >
<!-- global contents -->
<field field_ref="ibgvol_tot" name="ibgvol_tot" />
<field field_ref="sbgvol_tot" name="sbgvol_tot" />
......@@ -113,21 +119,21 @@
</file_group>
<file_group id="1ts" output_freq="1ts" output_level="10" enabled=".TRUE."/> <!-- 1 time step files -->
<file_group id="1h" output_freq="1h" output_level="10" enabled=".TRUE."/> <!-- 1h files -->
<file_group id="2h" output_freq="2h" output_level="10" enabled=".TRUE."/> <!-- 2h files -->
<file_group id="3h" output_freq="3h" output_level="10" enabled=".TRUE."/> <!-- 3h files -->
<file_group id="4h" output_freq="4h" output_level="10" enabled=".TRUE."/> <!-- 4h files -->
<file_group id="6h" output_freq="6h" output_level="10" enabled=".TRUE."/> <!-- 6h files -->
<file_group id="1m" output_freq="1mo" output_level="10" enabled=".TRUE."/> <!-- real monthly files -->
<file_group id="2m" output_freq="2mo" output_level="10" enabled=".TRUE."/> <!-- real 2m files -->
<file_group id="3m" output_freq="3mo" output_level="10" enabled=".TRUE."/> <!-- real 3m files -->
<file_group id="4m" output_freq="4mo" output_level="10" enabled=".TRUE."/> <!-- real 4m files -->
<file_group id="6m" output_freq="6mo" output_level="10" enabled=".TRUE."/> <!-- real 6m files -->
<file_group id="1y" output_freq="1y" output_level="10" enabled=".TRUE."/> <!-- real yearly files -->
<file_group id="2y" output_freq="2y" output_level="10" enabled=".TRUE."/> <!-- real 2y files -->
<file_group id="5y" output_freq="5y" output_level="10" enabled=".TRUE."/> <!-- real 5y files -->
<file_group id="10y" output_freq="10y" output_level="10" enabled=".TRUE."/> <!-- real 10y files -->
<file_group id="ice_1ts" output_freq="1ts" output_level="10" enabled=".TRUE."/> <!-- 1 time step files -->
<file_group id="ice_1h" output_freq="1h" output_level="10" enabled=".TRUE."/> <!-- 1h files -->
<file_group id="ice_2h" output_freq="2h" output_level="10" enabled=".TRUE."/> <!-- 2h files -->
<file_group id="ice_3h" output_freq="3h" output_level="10" enabled=".TRUE."/> <!-- 3h files -->
<file_group id="ice_4h" output_freq="4h" output_level="10" enabled=".TRUE."/> <!-- 4h files -->
<file_group id="ice_6h" output_freq="6h" output_level="10" enabled=".TRUE."/> <!-- 6h files -->
<file_group id="ice_1m" output_freq="1mo" output_level="10" enabled=".TRUE."/> <!-- real monthly files -->
<file_group id="ice_2m" output_freq="2mo" output_level="10" enabled=".TRUE."/> <!-- real 2m files -->
<file_group id="ice_3m" output_freq="3mo" output_level="10" enabled=".TRUE."/> <!-- real 3m files -->
<file_group id="ice_4m" output_freq="4mo" output_level="10" enabled=".TRUE."/> <!-- real 4m files -->
<file_group id="ice_6m" output_freq="6mo" output_level="10" enabled=".TRUE."/> <!-- real 6m files -->
<file_group id="ice_1y" output_freq="1y" output_level="10" enabled=".TRUE."/> <!-- real yearly files -->
<file_group id="ice_2y" output_freq="2y" output_level="10" enabled=".TRUE."/> <!-- real 2y files -->
<file_group id="ice_5y" output_freq="5y" output_level="10" enabled=".TRUE."/> <!-- real 5y files -->
<file_group id="ice_10y" output_freq="10y" output_level="10" enabled=".TRUE."/> <!-- real 10y files -->
</file_definition>
cfgs/ORCA2_SAS_ICE/EXPREF/file_def_nemo-oce.xml
View file @ dd59ced6
......@@ -10,8 +10,8 @@
<file_definition type="one_file" name="@expname@_@freq@_@startdate@_@enddate@" sync_freq="1mo" min_digits="4">
<file_group id="5d" output_freq="5d" output_level="10" enabled=".TRUE."> <!-- 5d files -->
<file id="file11" name_suffix="_grid_T" description="ocean T grid variables" >
<file_group id="oce_5d" output_freq="5d" output_level="10" enabled=".TRUE."> <!-- 5d files -->
<file id="file1" name_suffix="_grid_T" description="ocean T grid variables" >
<field field_ref="e3t" />
<field field_ref="toce" name="thetao" operation="instant" freq_op="5d" > @toce_e3t / @e3t </field>
<field field_ref="soce" name="so" operation="instant" freq_op="5d" > @soce_e3t / @e3t </field>
......@@ -34,17 +34,18 @@
<field field_ref="qt_oce" name="qt_oce" />
<field field_ref="saltflx" name="sfx" />
<field field_ref="taum" name="taum" />
<field field_ref="utau" name="tauuo" />
<field field_ref="vtau" name="tauvo" />
<field field_ref="wspd" name="windsp" />
<field field_ref="precip" name="precip" />
<!-- ice and snow -->
<field field_ref="snowpre" />
</file>
<file id="file12" name_suffix="_grid_U" description="ocean U grid variables" >
<file id="file2" name_suffix="_grid_U" description="ocean U grid variables" >
<field field_ref="e3u" />
<field field_ref="ssu" name="uos" />
<field field_ref="uoce" name="uo" operation="instant" freq_op="5d" > @uoce_e3u / @e3u </field>
<field field_ref="utau" name="tauuo" />
<field field_ref="uocetr_eff" name="uocetr_eff" />
<!-- available with diaar5 -->
<field field_ref="u_masstr" name="vozomatr" />
......@@ -52,11 +53,10 @@
<field field_ref="u_salttr" name="sozosatr" />
</file>
<file id="file13" name_suffix="_grid_V" description="ocean V grid variables" >
<file id="file3" name_suffix="_grid_V" description="ocean V grid variables" >
<field field_ref="e3v" />
<field field_ref="ssv" name="vos" />
<field field_ref="voce" name="vo" operation="instant" freq_op="5d" > @voce_e3v / @e3v </field>
<field field_ref="vtau" name="tauvo" />
<field field_ref="vocetr_eff" name="vocetr_eff" />
<!-- available with diaar5 -->
<field field_ref="v_masstr" name="vomematr" />
......@@ -64,7 +64,7 @@
<field field_ref="v_salttr" name="somesatr" />
</file>
<file id="file14" name_suffix="_grid_W" description="ocean W grid variables" >
<file id="file4" name_suffix="_grid_W" description="ocean W grid variables" >
<field field_ref="e3w" />
<field field_ref="woce" name="wo" />
<field field_ref="avt" name="difvho" />
......@@ -78,29 +78,28 @@
<field field_ref="av_ratio" name="av_ratio" />
</file>
</file_group>
<file_group id="1m" output_freq="1mo" output_level="10" enabled=".TRUE."> <!-- real monthly files -->
<file_group id="oce_1m" output_freq="1mo" output_level="10" enabled=".TRUE."> <!-- real monthly files -->
</file_group>
<file_group id="1y" output_freq="1y" output_level="10" enabled=".TRUE."> <!-- real yearly files -->
<file_group id="oce_1y" output_freq="1y" output_level="10" enabled=".TRUE."> <!-- real yearly files -->
</file_group>
<file_group id="1ts" output_freq="1ts" output_level="10" enabled=".TRUE."/> <!-- 1 time step files -->
<file_group id="2h" output_freq="2h" output_level="10" enabled=".TRUE."/> <!-- 2h files -->
<file_group id="3h" output_freq="3h" output_level="10" enabled=".TRUE."/> <!-- 3h files -->
<file_group id="4h" output_freq="4h" output_level="10" enabled=".TRUE."/> <!-- 4h files -->
<file_group id="6h" output_freq="6h" output_level="10" enabled=".TRUE."/> <!-- 6h files -->
<file_group id="2m" output_freq="2mo" output_level="10" enabled=".TRUE."/> <!-- real 2m files -->
<file_group id="3m" output_freq="3mo" output_level="10" enabled=".TRUE."/> <!-- real 3m files -->
<file_group id="4m" output_freq="4mo" output_level="10" enabled=".TRUE."/> <!-- real 4m files -->
<file_group id="6m" output_freq="6mo" output_level="10" enabled=".TRUE."/> <!-- real 6m files -->
<file_group id="2y" output_freq="2y" output_level="10" enabled=".TRUE."/> <!-- real 2y files -->
<file_group id="5y" output_freq="5y" output_level="10" enabled=".TRUE."/> <!-- real 5y files -->
<file_group id="10y" output_freq="10y" output_level="10" enabled=".TRUE."/> <!-- real 10y files -->
<file_group id="oce_1ts" output_freq="1ts" output_level="10" enabled=".TRUE."/> <!-- 1 time step files -->
<file_group id="oce_2h" output_freq="2h" output_level="10" enabled=".TRUE."/> <!-- 2h files -->
<file_group id="oce_3h" output_freq="3h" output_level="10" enabled=".TRUE."/> <!-- 3h files -->
<file_group id="oce_4h" output_freq="4h" output_level="10" enabled=".TRUE."/> <!-- 4h files -->
<file_group id="oce_6h" output_freq="6h" output_level="10" enabled=".TRUE."/> <!-- 6h files -->
<file_group id="oce_2m" output_freq="2mo" output_level="10" enabled=".TRUE."/> <!-- real 2m files -->
<file_group id="oce_3m" output_freq="3mo" output_level="10" enabled=".TRUE."/> <!-- real 3m files -->
<file_group id="oce_4m" output_freq="4mo" output_level="10" enabled=".TRUE."/> <!-- real 4m files -->
<file_group id="oce_6m" output_freq="6mo" output_level="10" enabled=".TRUE."/> <!-- real 6m files -->
<file_group id="oce_2y" output_freq="2y" output_level="10" enabled=".TRUE."/> <!-- real 2y files -->
<file_group id="oce_5y" output_freq="5y" output_level="10" enabled=".TRUE."/> <!-- real 5y files -->
<file_group id="oce_10y" output_freq="10y" output_level="10" enabled=".TRUE."/> <!-- real 10y files -->
</file_definition>
......
cfgs/ORCA2_SAS_ICE/EXPREF/namelist_cfg
View file @ dd59ced6
......@@ -10,10 +10,7 @@
!! namcfg parameters of the configuration (default: user defined GYRE)
!! namwad Wetting and drying (default: OFF)
!! namtsd data: temperature & salinity (default: OFF)
!! namcrs coarsened grid (for outputs and/or TOP) (ln_crs =T)
!! namc1d 1D configuration options (ln_c1d =T)
!! namc1d_dyndmp 1D newtonian damping applied on currents (ln_c1d =T)
!! namc1d_uvd 1D data (currents) (ln_c1d =T)
!!======================================================================
!
!-----------------------------------------------------------------------
......@@ -27,7 +24,7 @@
!-----------------------------------------------------------------------
&namdom ! time and space domain
!-----------------------------------------------------------------------
ln_linssh = .true. ! =T linear free surface ==>> model level are fixed in time
rn_Dt = 10800. ! time step for the dynamics and tracer
/
!-----------------------------------------------------------------------
&namcfg ! parameters of the configuration (default: user defined GYRE)
......@@ -174,13 +171,10 @@
!! *** Diagnostics namelists *** !!
!! !!
!! namtrd dynamics and/or tracer trends (default: OFF)
!! namhsb Heat and salt budgets (default: OFF)
!! namdiu Cool skin and warm layer models (default: OFF)
!! namdiu Cool skin and warm layer models (default: OFF)
!! namflo float parameters (default: OFF)
!! nam_diadct transports through some sections (default: OFF)
!! nam_dia25h 25h Mean Output (default: OFF)
!! namnc4 netcdf4 chunking and compression settings ("key_netcdf4")
!! namnc4 netcdf4 chunking and compression settings
!!======================================================================
!
!!======================================================================
......
cfgs/ORCA2_SAS_ICE/cpp_ORCA2_SAS_ICE.fcm
View file @ dd59ced6
bld::tool::fppkeys key_si3 key_linssh key_xios
bld::tool::fppkeys key_si3 key_linssh key_vco_1d3d key_xios
cfgs/README.rst deleted 100644 → 0
View file @ 7b99da57
********************************
Run the Reference configurations
********************************
.. todo::
Lack of illustrations for ref. cfgs, and more generally in the guide.
NEMO is distributed with a set of reference configurations allowing both
the user to set up his own first applications and
the developer to test/validate his NEMO developments (using SETTE package).
.. contents::
:local:
:depth: 1
.. attention::
Concerning the configurations,
the NEMO System Team is only in charge of the so-called reference configurations described below.
.. hint::
Configurations developed by external research projects or initiatives that
make use of NEMO are welcome to be publicized through the website by
filling up the form :website:`to add an associated project<projects/add>`.
How to compile an experiment from a reference configuration
===========================================================
To compile the ORCA2_ICE_PISCES_ reference configuration using :file:`makenemo`,
one should use the following, by selecting among available architecture file or
providing a user defined one:
.. code-block:: console
$ ./makenemo -r 'ORCA2_ICE_PISCES' -m 'my_arch' -j '4'
A new ``EXP00`` folder will be created within the selected reference configurations,
namely ``./cfgs/ORCA2_ICE_PISCES/EXP00``.
It will be necessary to uncompress the archives listed in the above table for
the given reference configuration that includes input & forcing files.
Then it will be possible to launch the execution of the model through a runscript
(opportunely adapted to the user system).
List of Configurations
======================
All forcing files listed below in the table are available from |DOI data|_
=================== === === === === === ==================================
Configuration Component(s) Archives (input & forcing files)
------------------- ------------------- ----------------------------------
Name O S T P A
=================== === === === === === ==================================
AGRIF_DEMO_ X X X AGRIF_DEMO_v4.0.tar,
ORCA2_ICE_v4.0.tar
AMM12_ X AMM12_v4.0.tar
C1D_PAPA_ X INPUTS_C1D_PAPA_v4.0.tar
GYRE_BFM_ X X *none*
GYRE_PISCES_ X X X *none*
ORCA2_ICE_PISCES_ X X X X ORCA2_ICE_v4.0.tar,
INPUTS_PISCES_v4.0.tar
ORCA2_OFF_PISCES_ X X ORCA2_OFF_v4.0.tar,
INPUTS_PISCES_v4.0.tar
ORCA2_OFF_TRC_ X ORCA2_OFF_v4.0.tar
ORCA2_SAS_ICE_ X ORCA2_ICE_v4.0.tar,
INPUTS_SAS_v4.0.tar
SPITZ12_ X X SPITZ12_v4.0.tar
=================== === === === === === ==================================
.. admonition:: Legend for component combination
O for OCE, S for SI\ :sup:`3`, T for TOP, P for PISCES and A for AGRIF
AGRIF_DEMO
----------
``AGRIF_DEMO`` is based on the ``ORCA2_ICE_PISCES`` global configuration at 2° of resolution with
the inclusion of 3 online nested grids to demonstrate the overall capabilities of AGRIF in
a realistic context (including the nesting of sea ice models).
The configuration includes a 1:1 grid in the Pacific and two successively nested grids with
odd and even refinement ratios over the Arctic ocean,
with the finest grid spanning the whole Svalbard archipelago that is of
particular interest to test sea ice coupling.
.. image:: _static/AGRIF_DEMO_no_cap.jpg
:scale: 66%
:align: center
The 1:1 grid can be used alone as a benchmark to check that
the model solution is not corrupted by grid exchanges.
Note that since grids interact only at the baroclinic time level,
numerically exact results can not be achieved in the 1:1 case.
Perfect reproducibility is obtained only by switching to a fully explicit setup instead of
a split explicit free surface scheme.
AMM12
-----
``AMM12`` stands for *Atlantic Margin Model at 12 km* that is
a regional configuration covering the Northwest European Shelf domain on
a regular horizontal grid of ~12 km of resolution (see :cite:`ODEA2012`).
.. image:: _static/AMM_domain.png
:align: center
This configuration allows to tests several features of NEMO specifically addressed to the shelf seas.
In particular, ``AMM12`` accounts for vertical s-coordinates system, GLS turbulence scheme,
tidal lateral boundary conditions using a flather scheme (see more in ``BDY``).
Boundaries may be completely omitted by setting ``ln_bdy = .false.`` in ``nambdy``.
Sample surface fluxes, river forcing and an initial restart file are included to test a realistic model run
(``AMM12_v4.0.tar``).
Note that, the Baltic boundary is included within the river input file and is specified as a river source,
but unlike ordinary river points the Baltic inputs also include salinity and temperature data.
C1D_PAPA
--------
.. figure:: _static/Papa2015.jpg
:height: 225px
:align: left
``C1D_PAPA`` is a 1D configuration for the `PAPA station`_ located in
the northern-eastern Pacific Ocean at 50.1°N, 144.9°W.
See :gmd:`Reffray et al. (2015) <8/69/2015>` for the description of
its physical and numerical turbulent-mixing behaviour.
| The water column setup, called NEMO1D, is activated by
setting ``ln_c1d = .true.`` in ``namdom`` and
has a horizontal domain of 1x1 grid point.
| This reference configuration uses 75 vertical levels grid (1m at the surface),
GLS turbulence scheme with K-epsilon closure and the NCAR bulk formulae.
Data provided with ``INPUTS_C1D_PAPA_v4.2.tar`` file account for:
- :file:`forcing_PAPASTATION_1h_y201[0-1].nc`:
ECMWF operational analysis atmospheric forcing rescaled to 1h
(with long and short waves flux correction) for years 2010 and 2011
- :file:`init_PAPASTATION_m06d15.nc`: Initial Conditions from
observed data and Levitus 2009 climatology
- :file:`chlorophyll_PAPASTATION.nc`: surface chlorophyll file from Seawifs data
GYRE_BFM
--------
``GYRE_BFM`` shares the same physical setup of GYRE_PISCES_,
but NEMO is coupled with the `BFM`_ biogeochemical model as described in ``./cfgs/GYRE_BFM/README``.
GYRE_PISCES
-----------
``GYRE_PISCES`` is an idealized configuration representing a Northern hemisphere double gyres system,
in the Beta-plane approximation with a regular 1° horizontal resolution and 31 vertical levels,
with PISCES BGC model :cite:`gmd-8-2465-2015`.
Analytical forcing for heat, freshwater and wind-stress fields are applied.
This configuration acts also as demonstrator of the **user defined setup**
(``ln_read_cfg = .false.``) and grid setting are handled through
the ``&namusr_def`` controls in :file:`namelist_cfg`:
.. literalinclude:: ../../../cfgs/GYRE_PISCES/EXPREF/namelist_cfg
:language: fortran
:lines: 35-41
Note that, the default grid size is 30x20 grid points (with ``nn_GYRE = 1``) and
vertical levels are set by ``jpkglo``.
The specific code changes can be inspected in :file:`./src/OCE/USR`.
.. rubric:: Running GYRE as a benchmark
| This simple configuration can be used as a benchmark since it is easy to increase resolution,
with the drawback of getting results that have a very limited physical meaning.
| GYRE grid resolution can be increased at runtime by setting a different value of ``nn_GYRE``
(integer multiplier scaling factor), as described in the following table:
=========== ============ ============ ============ ===============
``nn_GYRE`` ``jpiglo`` ``jpjglo`` ``jpkglo`` Equivalent to
=========== ============ ============ ============ ===============
1 30 20 31 GYRE 1°
25 750 500 101 ORCA 1/2°
50 1500 1000 101 ORCA 1/4°
150 4500 3000 101 ORCA 1/12°
200 6000 4000 101 ORCA 1/16°
=========== ============ ============ ============ ===============
| Note that, it is necessary to set ``ln_bench = .true.`` in ``&namusr_def`` to
avoid problems in the physics computation and that
the model timestep should be adequately rescaled.
| For example if ``nn_GYRE = 150``, equivalent to an ORCA 1/12° grid,
the timestep ``rn_rdt`` should be set to 1200 seconds
Differently from previous versions of NEMO, the code uses by default the time-splitting scheme and
internally computes the number of sub-steps.
ORCA2_ICE_PISCES
----------------
``ORCA2_ICE_PISCES`` is a reference configuration for the global ocean with
a 2°x2° curvilinear horizontal mesh and 31 vertical levels,
distributed using z-coordinate system and with 10 levels in the top 100m.
ORCA is the generic name given to global ocean Mercator mesh,
(i.e. variation of meridian scale factor as cosinus of the latitude),
with two poles in the northern hemisphere so that
the ratio of anisotropy is nearly one everywhere
This configuration uses the three components
- |OCE|, the ocean dynamical core
- |ICE|, the thermodynamic-dynamic sea ice model.
- |MBG|, passive tracer transport module and PISCES BGC model :cite:`gmd-8-2465-2015`
All components share the same grid.
The model is forced with CORE-II normal year atmospheric forcing and
it uses the NCAR bulk formulae.
.. rubric:: Ocean Physics
:horizontal diffusion on momentum:
the eddy viscosity coefficient depends on the geographical position.
It is taken as 40000 m\ :sup:`2`/s, reduced in the equator regions (2000 m\ :sup:`2`/s)
excepted near the western boundaries.
:isopycnal diffusion on tracers:
the diffusion acts along the isopycnal surfaces (neutral surface) with
an eddy diffusivity coefficient of 2000 m\ :sup:`2`/s.
:Eddy induced velocity parametrization:
With a coefficient that depends on the growth rate of baroclinic instabilities
(it usually varies from 15 m\ :sup:`2`/s to 3000 m\ :sup:`2`/s).
:lateral boundary conditions:
Zero fluxes of heat and salt and no-slip conditions are applied through lateral solid boundaries.
:bottom boundary condition:
Zero fluxes of heat and salt are applied through the ocean bottom.
The Beckmann [19XX] simple bottom boundary layer parameterization is applied along
continental slopes.
A linear friction is applied on momentum.
:convection:
The vertical eddy viscosity and diffusivity coefficients are increased to 1 m\ :sup:`2`/s in
case of static instability.
:time step: is 5400sec (1h30') so that there is 16 time steps in one day.
ORCA2_OFF_PISCES
----------------
``ORCA2_OFF_PISCES`` shares the same general offline configuration of ``ORCA2_ICE_TRC``,
but only PISCES model is an active component of TOP.
ORCA2_OFF_TRC
-------------
| ``ORCA2_OFF_TRC`` is based on the ORCA2 global ocean configuration
(see ORCA2_ICE_PISCES_ for general description) along with
the tracer passive transport module (TOP),
but dynamical fields are pre-calculated and read with specific time frequency.
| This enables for an offline coupling of TOP components,
here specifically inorganic carbon compounds (CFC11, CFC12, SF6, C14) and water age module (age).
See :file:`namelist_top_cfg` to inspect the selection of
each component with the dedicated logical keys.
Pre-calculated dynamical fields are provided to NEMO using
the namelist ``&namdta_dyn`` in :file:`namelist_cfg`,
in this case with a 5 days frequency (120 hours):
.. literalinclude:: ../../namelists/namdta_dyn
:language: fortran
Input dynamical fields for this configuration (:file:`ORCA2_OFF_v4.0.tar`) comes from
a 2000 years long climatological simulation of ORCA2_ICE using ERA40 atmospheric forcing.
| Note that,
this configuration default uses linear free surface (``ln_linssh = .true.``) assuming that
model mesh is not varying in time and
it includes the bottom boundary layer parameterization (``ln_trabbl = .true.``) that
requires the provision of BBL coefficients through ``sn_ubl`` and ``sn_vbl`` fields.
| It is also possible to activate PISCES model (see ``ORCA2_OFF_PISCES``) or
a user defined set of tracers and source-sink terms with ``ln_my_trc = .true.``
(and adaptation of ``./src/TOP/MY_TRC`` routines).
In addition, the offline module (OFF) allows for the provision of further fields:
1. **River runoff** can be provided to TOP components by setting ``ln_dynrnf = .true.`` and
by including an input datastream similarly to the following:
.. code-block:: fortran
sn_rnf = 'dyna_grid_T', 120, 'sorunoff' , .true., .true., 'yearly', '', '', ''
2. **VVL dynamical fields**, in the case input data were produced by a dyamical core using
variable volume (``ln_linssh = .false.``)
it is necessary to provide also diverce and E-P at before timestep by
including input datastreams similarly to the following
.. code-block:: fortran
sn_div = 'dyna_grid_T', 120, 'e3t' , .true., .true., 'yearly', '', '', ''
sn_empb = 'dyna_grid_T', 120, 'sowaflupb', .true., .true., 'yearly', '', '', ''
More details can be found by inspecting the offline data manager in
the routine :file:`./src/OFF/dtadyn.F90`.
ORCA2_SAS_ICE
-------------
| ORCA2_SAS_ICE is a demonstrator of the Stand-Alone Surface (SAS) module and
it relies on ORCA2 global ocean configuration (see ORCA2_ICE_PISCES_ for general description).
| The standalone surface module allows surface elements such as sea-ice, iceberg drift, and
surface fluxes to be run using prescribed model state fields.
It can profitably be used to compare different bulk formulae or
adjust the parameters of a given bulk formula.
More informations about SAS can be found in :doc:`NEMO manual <cite>`.
SPITZ12
-------
``SPITZ12`` is a regional configuration around the Svalbard archipelago
at 1/12° of horizontal resolution and 75 vertical levels.
See :gmd:`Rousset et al. (2015) <8/2991/2015>` for more details.
This configuration references to year 2002,
with atmospheric forcing provided every 2 hours using NCAR bulk formulae,
while lateral boundary conditions for dynamical fields have 3 days time frequency.
.. rubric:: References
.. bibliography:: cfgs.bib
:all:
:style: unsrt
:labelprefix: C
cfgs/SHARED/README.namelists deleted 100644 → 0
View file @ 7b99da57
Simple style rules for namelists
--------------------------------
NEMO reference namelists should adhere to the following simple style rules:
1. Comments outside a namelist block start with !! in column 1
2. Each namelist block starts with 3 lines of the form:
!-----------------------------------------------------------------------
&namblockname ! short description of block
!-----------------------------------------------------------------------
with all ! and & 's starting in column 1
3. The closing / for each namelist block is in column 1
4. Comments within namelist blocks never start with !- . Use ! followed
by space or != etc.
These conventions make it possible to construct empty configuration namelists.
For example, a namelist_cfg template can be produced from namelist_ref with
the following grep command; e.g.:
grep -E '^!-|^&|^/' namelist_ref > namelist_cfg.template
head namelist_cfg.template
!-----------------------------------------------------------------------
&namrun ! parameters of the run
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namcfg ! parameters of the configuration
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namdom ! time and space domain
!-----------------------------------------------------------------------
/
.
.
If all configuration namelists are produced and maintained using this
strategy then standard, side-by-side comaparators, such as vimdiff or xxdiff,
can be used to compare and transfer lines from the reference namelist to a
configuration namelist when setting up a new configuration.
Tips and tricks
---------------
1. The following bash function is useful when checking which namelist blocks
are in active use in a configuration namelist:
function list_used_nl(){ grep -n -E '^&|^/' "$1" | sed -e 's/:/ /' \
| awk ' BEGIN { x = 0 } \
{if ( NR % 2 == 0 && $1 - x > 2 ) printf("%3d %s\n", $1 - x , n) ; \
else x = $1; n = $2}' \
| sort -k 2;}
which (assuming the namelist adheres to the conventions) will list the number
of entries in each non-empty namelist block. The list is sorted on the block
name to ease comparisons. For example:
list_used_nl ORCA2_LIM3_PISCES/EXP00/namelist_cfg
1 &nambbc
5 &nambbl
30 &namberg
10 &namcfg
4 &namctl
3 &namdom
1 &namdrg
5 &namdyn_adv
1 &namdyn_hpg
22 &namdyn_ldf
1 &namdyn_spg
5 &namdyn_vor
3 &nameos
1 &namhsb
4 &namrun
1 &namsbc
1 &namsbc_blk
3 &namtra_adv
28 &namtra_ldf
10 &namtra_ldfeiv
25 &namzdf
3 &namzdf_iwm
2. vimdiff can give garish colours in some terminals. Usually this is because
vim assumes, incorrectly, that the terminal only supports 8 colours. Try forcing
256 colours with:
:set t_Co=256
to produce more pastel shades (add this to ~/.vimrc if successful).
3. Switching between vsplit panes in vim is a multi-key sequence. The tool is
much easier to use if the sequence is mapped to a spare key. Here I use the
§ and ± key on my Mac keyboard (add to ~/.vimrc):
map § ^Wl
map ± ^Wh
4. With easy switching between panes, constructing namelists in vimdiff just
requires the following commands in addition to normal editing:
]c - Go to next block of the diff
dp - Push version of the block under cursor into the other pane
do - Pull version of the block under cursor from the other pane
cfgs/SHARED/README.rst deleted 100644 → 0
View file @ 7b99da57
***********
Diagnostics
***********
.. todo::
.. contents::
:local:
Output of diagnostics in NEMO is usually done using XIOS.
This is an efficient way of writing diagnostics because
the time averaging, file writing and even some simple arithmetic or regridding is carried out in
parallel to the NEMO model run.
This page gives a basic introduction to using XIOS with NEMO.
Much more information is available from the :xios:`XIOS homepage<>` above and from the NEMO manual.
Use of XIOS for diagnostics is activated using the pre-compiler key ``key_xios``.
Extracting and installing XIOS
==============================
1. Install the NetCDF4 library.
If you want to use single file output you will need to compile the HDF & NetCDF libraries to
allow parallel IO.
2. Download the version of XIOS that you wish to use.
The recommended version is now XIOS 2.5:
.. code-block:: console
$ svn co http://forge.ipsl.jussieu.fr/ioserver/svn/XIOS/branchs/xios-2.5
and follow the instructions in :xios:`XIOS documentation <wiki/documentation>` to compile it.
If you find problems at this stage, support can be found by subscribing to
the :xios:`XIOS mailing list <../mailman/listinfo.cgi/xios-users>` and sending a mail message to it.
XIOS Configuration files
------------------------
XIOS is controlled using XML input files that should be copied to
your model run directory before running the model.
Examples of these files can be found in the reference configurations (:file:`./cfgs`).
The XIOS executable expects to find a file called :file:`iodef.xml` in the model run directory.
In NEMO we have made the decision to use include statements in the :file:`iodef.xml` file to include:
- :file:`field_def_nemo-oce.xml` (for physics),
- :file:`field_def_nemo-ice.xml` (for ice),
- :file:`field_def_nemo-pisces.xml` (for biogeochemistry) and
- :file:`domain_def.xml` from the :file:`./cfgs/SHARED` directory.
Most users will not need to modify :file:`domain_def.xml` or :file:`field_def_nemo-???.xml` unless
they want to add new diagnostics to the NEMO code.
The definition of the output files is organized into separate :file:`file_definition.xml` files which
are included in the :file:`iodef.xml` file.
Modes
=====
Detached Mode
-------------
In detached mode the XIOS executable is executed on separate cores from the NEMO model.
This is the recommended method for using XIOS for realistic model runs.
To use this mode set ``using_server`` to ``true`` at the bottom of the :file:`iodef.xml` file:
.. code-block:: xml
<variable id="using_server" type="boolean">true</variable>
Make sure there is a copy (or link to) your XIOS executable in the working directory and
in your job submission script allocate processors to XIOS.
Attached Mode
-------------
In attached mode XIOS runs on each of the cores used by NEMO.
This method is less efficient than the detached mode but can be more convenient for testing or
with small configurations.
To activate this mode simply set ``using_server`` to false in the :file:`iodef.xml` file
.. code-block:: xml
<variable id="using_server" type="boolean">false</variable>
and don't allocate any cores to XIOS.
.. note::
Due to the different domain decompositions between XIOS and NEMO,
if the total number of cores is larger than the number of grid points in the ``j`` direction then
the model run will fail.
Adding new diagnostics
======================
If you want to add a NEMO diagnostic to the NEMO code you will need to do the following:
1. Add any necessary code to calculate you new diagnostic in NEMO
2. Send the field to XIOS using ``CALL iom_put( 'field_id', variable )`` where
``field_id`` is a unique id for your new diagnostics and
variable is the fortran variable containing the data.
This should be called at every model timestep regardless of how often you want to output the field.
No time averaging should be done in the model code.
3. If it is computationally expensive to calculate your new diagnostic
you should also use "iom_use" to determine if it is requested in the current model run.
For example,
.. code-block:: fortran
IF iom_use('field_id') THEN
!Some expensive computation
!...
!...
iom_put('field_id', variable)
ENDIF
4. Add a variable definition to the :file:`field_def_nemo-???.xml` file.
5. Add the variable to the :file:`iodef.xml` or :file:`file_definition.xml` file.
cfgs/SHARED/axis_def_nemo.xml
View file @ dd59ced6
......@@ -6,22 +6,23 @@
<axis_definition>
<axis id="deptht" long_name="Vertical T levels" unit="m" positive="down" />
<!-- Vertical zoom for a 31-levels ORCA2 grid for eORCA1 300m corresponds to n=35 -->
<!-- Vertical zoom for a 31-levels ORCA2 grid. For eORCA1 300m corresponds to n=35 -->
<axis id="deptht300" axis_ref="deptht" >
<zoom_axis begin="0" n="19" />
</axis>
<axis id="depthu" long_name="Vertical U levels" unit="m" positive="down" />
<axis id="depthv" long_name="Vertical V levels" unit="m" positive="down" />
<axis id="depthw" long_name="Vertical W levels" unit="m" positive="down" />
<axis id="depthf" long_name="Vertical F levels" unit="m" positive="down" />
<axis id="nfloat" long_name="Float number" unit="-" />
<axis id="icbcla" long_name="Iceberg class" unit="1" />
<axis id="ncatice" long_name="Ice category" unit="1" />
<axis id="iax_20C" long_name="20 degC isotherm" unit="degC" />
<axis id="iax_26C" long_name="26 degC isotherm" unit="degC" />
<axis id="iax_28C" long_name="28 degC isotherm" unit="degC" />
<axis id="depthu" long_name="Vertical U levels" unit="m" positive="down" />
<axis id="depthv" long_name="Vertical V levels" unit="m" positive="down" />
<axis id="depthw" long_name="Vertical W levels" unit="m" positive="down" />
<axis id="depthf" long_name="Vertical F levels" unit="m" positive="down" />
<axis id="icbcla" long_name="Iceberg class" unit="1" />
<axis id="icbbas" long_name="Iceberg basin" unit="1" />
<axis id="ncatice" long_name="Ice category" unit="1" />
<axis id="nlayice" long_name="Ice layer" unit="1" />
<axis id="iax_20C" long_name="20 degC isotherm" unit="degC" />
<axis id="iax_26C" long_name="26 degC isotherm" unit="degC" />
<axis id="iax_28C" long_name="28 degC isotherm" unit="degC" />
<axis id="basin" long_name="Sub-basin mask (1=Global 2=Atlantic 3=Indo-Pacific 4=Indian, 5=Pacific)" unit="1" />
<axis id="nstrait" long_name="Number of straits" unit="1" />
<axis id="nstrait" long_name="Number of straits" unit="1" />
<!-- ABL vertical axis definition -->
<axis id="ght_abl" long_name="ABL Vertical T levels" unit="m" positive="up" />
<axis id="ghw_abl" long_name="ABL Vertical W levels" unit="m" positive="up" />
......
cfgs/SHARED/field_def_nemo-ice.xml
View file @ dd59ced6
<?xml version="1.0"?>
<!-- $id$ -->
<!--
============================================================================================================
......@@ -18,7 +17,7 @@
<field_group id="SBC" > <!-- time step automaticaly defined based on nn_fsbc -->
<!-- 2D variables -->
<field_group id="SBC_2D" grid_ref="grid_T_2D" >
<field_group id="SBC_2D" grid_ref="grid_T_2D_inner" >
<!-- =================== -->
<!-- standard ice fields -->
......@@ -46,6 +45,7 @@
<!-- melt ponds -->
<field id="iceapnd" long_name="melt pond concentration" standard_name="sea_ice_meltpond_concentration" unit="" />
<field id="iceapeff" long_name="melt pond concentration effective" standard_name="sea_ice_meltpond_concentration effective" 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" />
......@@ -69,8 +69,8 @@
<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="uice" long_name="X-component of sea ice velocity" standard_name="sea_ice_x_velocity" unit="m/s" grid_ref="grid_T_2D" />
<field id="vice" long_name="Y-component of sea ice velocity" standard_name="sea_ice_y_velocity" unit="m/s" grid_ref="grid_T_2D" />
<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" />
......@@ -93,7 +93,14 @@
<field id="yield11" long_name="yield surface tensor component 11" standard_name="yield11" unit="N/m" />
<field id="yield22" long_name="yield surface tensor component 22" standard_name="yield22" unit="N/m" />
<field id="yield12" long_name="yield surface tensor component 12" standard_name="yield12" unit="N/m" />
<field id="beta_evp" long_name="Relaxation parameter of ice rheology (beta)" standard_name="relaxation_parameter_of_ice_rheology" unit="" />
<field id="beta_evp" long_name="Relaxation parameter of ice rheology (beta)" standard_name="relaxation_parameter_of_ice_rheology" unit="" grid_ref="grid_T_2D" />
<field id="lead_open" long_name="Lead area opening rate" standard_name="lead_open" unit="1/s" />
<field id="rdg_loss" long_name="Ridging ice area loss rate" standard_name="rdg_loss" unit="1/s" />
<field id="rft_loss" long_name="Rafting ice area loss rate"
standard_name="rft_loss" unit="1/s" />
<field id="rdg_gain" long_name="New ridged ice area gain rate" standard_name="rdg_gain" unit="1/s" />
<field id="rft_gain" long_name="New rafted ice area gain rate" standard_name="rft_gain" unit="1/s" />
<!-- 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" />
......@@ -111,6 +118,17 @@
<field id="Cd_ice" long_name="Drag coefficient over sea-ice" standard_name="drag_coefficient_ice" unit="" />
<field id="Ce_ice" long_name="Sublimation coefficient over sea-ice" standard_name="evap_coefficient_ice" unit="" />
<field id="Ch_ice" long_name="Sensible heat coefficient over sea-ice" standard_name="sensible_heat_coefficient_ice" unit="" />
<field id="drag_io" long_name="momentum drag coefficient ice-ocean" standard_name="sea_ice_ocean_drag" unit=" " />
<field id="drag_io_skin" long_name="momentum drag coefficient ice-ocean skin" standard_name="sea_ice_ocean_drag skin" unit=" " />
<field id="drag_io_keel" long_name="momentum drag coefficient ice-ocean keel" standard_name="sea_ice_ocean_drag keel" unit=" " />
<field id="drag_io_floe" long_name="momentum drag coefficient ice-ocean floe" standard_name="sea_ice_ocean_drag floe" unit=" " />
<field id="drag_ia" long_name="momentum drag coefficient ice-atmosphere" standard_name="sea_ice_atmosphere_drag" unit=" " />
<field id="drag_ia_skin" long_name="momentum drag coefficient ice-atmosphere skin" standard_name="sea_ice_atmosphere_drag skin" unit=" " />
<field id="drag_ia_rdg" long_name="momentum drag coefficient ice-atmosphere rdg" standard_name="sea_ice_atmosphere_drag rdg" unit=" " />
<field id="drag_ia_floe" long_name="momentum drag coefficient ice-atmosphere floe" standard_name="sea_ice_atmosphere_drag floe" unit=" " />
<field id="drag_ia_pond" long_name="momentum drag coefficient ice-atmosphere pond" standard_name="sea_ice_atmosphere_drag pond" unit=" " />
<field id="ardg_drag" long_name="ridged ice concentration" standard_name="ridged ice concentration" unit=" " />
<field id="vrdg_drag" long_name="ridged ice thickness" standard_name="ridged ice thickness" unit="m" />
<!-- 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" />
......@@ -188,18 +206,15 @@
<field id="icedrift_salt" long_name="Ice salt drift (conservation check)" unit="kg/m2/s" />
<field id="icedrift_heat" long_name="Ice heat drift (conservation check)" unit="W/m2" />
<!-- 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="-" />
<!-- rheology convergence tests -->
<field id="uice_cvg" long_name="sea ice velocity convergence" standard_name="sea_ice_velocity_convergence" unit="m/s" />
<!-- vp rheology convergence tests -->
<field id="velo_res" long_name="sea ice velocity residual" standard_name="sea_ice_velocity_residual" unit="N/m2" />
<field id="velo_ero" long_name="sea ice velocity error last outer iteration" standard_name="sea_ice_velocity_outer_error" unit="m/s" />
<field id="uice_eri" long_name="uice velocity error last inner iteration" standard_name="sea_ice_u_velocity_inner_error" unit="m/s" />
<field id="vice_eri" long_name="vice velocity error last inner iteration" standard_name="sea_ice_v_velocity_inner_error" unit="m/s" />
<!-- ================= -->
<!-- Add-ons for SIMIP -->
<!-- ================= -->
......@@ -305,7 +320,7 @@
</field_group> <!-- SBC_2D -->
<!-- categories -->
<field_group id="SBC_3D" grid_ref="grid_T_ncatice" >
<field_group id="SBC_3D" grid_ref="grid_T_ncatice_inner" >
<!-- standard ice fields -->
<field id="iceconc_cat" long_name="Sea-ice concentration per category" unit="" />
......@@ -317,6 +332,7 @@
<field id="icettop_cat" long_name="Ice/snow surface temperature per category" unit="degC" detect_missing_value="true" />
<field id="iceapnd_cat" long_name="Ice melt pond grid fraction per category" unit="" />
<field id="icevpnd_cat" long_name="Ice melt pond volume per grid area per category" unit="m" />
<field id="icevlid_cat" long_name="Ice melt pond lid volume per category" unit="m" />
<field id="icehpnd_cat" long_name="Ice melt pond thickness per category" unit="m" detect_missing_value="true" />
<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 ice fraction per category" unit="" />
......@@ -336,8 +352,23 @@
<field id="tice_cvgerr" long_name="sea ice temperature convergence error" standard_name="sea_ice_temperature_convergence_err" unit="K" />
<field id="tice_cvgstp" long_name="sea ice temperature convergence iterations" standard_name="sea_ice_temperature_convergence_stp" unit="" />
<!-- salt drainage and flushing sanity checks -->
<field id="cfl_flush" long_name="max CFL during flushing" unit="" />
<field id="cfl_drain" long_name="max CFL during drainage" unit="" />
<field id="sice_flush_neg" long_name="negative sea ice salinity during flushing" detect_missing_value="true" unit="g/kg" />
<field id="sice_drain_neg" long_name="negative sea ice salinity during drainage" detect_missing_value="true" unit="g/kg" />
</field_group> <!-- SBC_3D -->
<!-- layers -->
<field_group id="SBC_4D" grid_ref="grid_T_nlayice_inner" >
<!-- standard ice fields -->
<field id="icesalt_lay" long_name="Sea-Ice salinity per layer" unit="g/kg" detect_missing_value="true" />
<field id="icetemp_lay" long_name="Ice temperature per layer" unit="degC" detect_missing_value="true" />
</field_group> <!-- SBC_4D -->
<!-- scalar variables -->
<field_group id="SBC_scalar" grid_ref="grid_scalar" >
<field id="NH_iceextt" long_name="Sea ice extent North" standard_name="sea_ice_extent_n" unit="1e6_km2" />
......@@ -347,7 +378,6 @@
<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" />
<!-- available with ln_icediahsb -->
<!-- global forcings -->
<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" />
......@@ -380,7 +410,7 @@
-->
<!-- output variables for my configuration (example) -->
<field_group id="myvarICE" grid_ref="grid_T_2D" >
<field_group id="myvarICE" >
<!-- ice mask -->
<field field_ref="icemask" name="simsk" />
<field field_ref="icemask05" name="simsk05" />
......@@ -399,6 +429,7 @@
<field field_ref="iceconc" name="siconc" />
<field field_ref="icesalt" name="sisali" />
<field field_ref="iceapnd" name="siapnd" />
<field field_ref="iceapeff" name="siapeff" />
<field field_ref="icehpnd" name="sihpnd" />
<field field_ref="icevpnd" name="sivpnd" />
<field field_ref="iceage" name="siage" />
......@@ -421,6 +452,20 @@
<field field_ref="utau_oi" name="utau_oi" />
<field field_ref="vtau_oi" name="vtau_oi" />
<!-- drags -->
<field field_ref="Cd_ice" name="Cd_ice" />
<field field_ref="Ch_ice" name="Ch_ice" />
<field field_ref="Ce_ice" name="Ce_ice" />
<field field_ref="drag_io" name="drag_io" />
<field field_ref="drag_io_skin" name="drag_io_skin" />
<field field_ref="drag_io_keel" name="drag_io_keel" />
<field field_ref="drag_io_floe" name="drag_io_floe" />
<field field_ref="drag_ia" name="drag_ia" />
<field field_ref="drag_ia_skin" name="drag_ia_skin" />
<field field_ref="drag_ia_rdg" name="drag_ia_rdg" />
<field field_ref="drag_ia_floe" name="drag_ia_floe" />
<field field_ref="drag_ia_pond" name="drag_ia_pond" />
<!-- rheology -->
<field field_ref="icediv" name="sidive" />
<field field_ref="iceshe" name="sishea" />
......@@ -429,7 +474,7 @@
<field field_ref="sheastr" name="sheastr" />
<field field_ref="sig1_pnorm" name="sig1_pnorm"/>
<field field_ref="sig2_pnorm" name="sig2_pnorm"/>
<field field_ref="icedlt" name="sidelta" />
<field field_ref="icedlt" name="sidelt" />
<!-- heat fluxes -->
<field field_ref="qt_oce_ai" name="qt_oce_ai" />
......@@ -489,7 +534,7 @@
</field_group>
<field_group id="myvarICE_cat" grid_ref="grid_T_ncatice" >
<field_group id="myvarICE_cat" >
<!-- categories -->
<field field_ref="icemask_cat" name="simskcat"/>
......@@ -518,7 +563,7 @@
<field field_ref="ilbgvol_tot" name="ilbgvol_tot" />
</field_group>
<field_group id="ICE_budget" grid_ref="grid_T_2D" >
<field_group id="ICE_budget" >
<!-- general -->
<field field_ref="icemask" name="simsk" />
<field field_ref="iceconc" name="siconc" />
......@@ -573,7 +618,7 @@
</field_group>
<!-- SIMIP daily fields -->
<field_group id="SIday_fields" grid_ref="grid_T_2D" >
<field_group id="SIday_fields" >
<field field_ref="icepres" name="sitimefrac" />
<field field_ref="iceconc_pct" name="siconc" />
<field field_ref="icethic_cmip" name="sithick" />
......@@ -585,7 +630,7 @@
</field_group>
<!-- SIMIP monthly fields -->
<field_group id="SImon_fields" grid_ref="grid_T_2D" >
<field_group id="SImon_fields" >
<!-- Sea-ice state variables -->
<field field_ref="icepres" name="sitimefrac" />
<field field_ref="iceconc_pct" name="siconc" />
......
cfgs/SHARED/field_def_nemo-innerttrc.xml
View file @ dd59ced6
<?xml version="1.0"?>
<!-- $id$ -->
<!--
============================================================================================================
......@@ -16,7 +15,7 @@
-->
<field_group id="inerttrc" grid_ref="grid_T_2D">
<field_group id="inerttrc" grid_ref="grid_T_2D_inner">
<!-- CFC11 : variables available with ln_cfc11 -->
<field id="CFC11" long_name="Chlorofluoro carbon11 Concentration" unit="umol/m3" grid_ref="grid_T_3D" />
......@@ -39,8 +38,8 @@
<!-- C14 : variables available with ln_c14 -->
<field id="RC14" long_name="Radiocarbon ratio" unit="-" grid_ref="grid_T_3D" />
<field id="RC14_e3t" long_name="RC14 * e3t" unit="m" grid_ref="grid_T_3D" > RC14 * e3t </field >
<field id="DeltaC14" long_name="Delta C14" unit="permil" grid_ref="grid_T_3D" />
<field id="C14Age" long_name="Radiocarbon age" unit="yr" grid_ref="grid_T_3D" />
<field id="DeltaC14" long_name="Delta C14" unit="permil" grid_ref="grid_T_3D_inner" />
<field id="C14Age" long_name="Radiocarbon age" unit="yr" grid_ref="grid_T_3D_inner" />
<field id="RAge" long_name="Reservoir Age" unit="yr" />
<field id="qtr_c14" long_name="Air-sea flux of C14" unit="1/m2/s" />
<field id="qint_c14" long_name="Cumulative air-sea flux of C14" unit="1/m2" />
......@@ -52,7 +51,7 @@
<!-- AGE : variables available with ln_age -->
<field id="Age" long_name="Sea water age since surface contact" unit="yr" grid_ref="grid_T_3D" />
<field id="Age_e3t" long_name="Age * e3t" unit="yr * m" grid_ref="grid_T_3D" > Age * e3t </field >
<field id="Age_e3t" long_name="Age * e3t" unit="yr * m" grid_ref="grid_T_3D" > Age * e3t </field >
</field_group>
......
cfgs/SHARED/field_def_nemo-oce.xml
View file @ dd59ced6
<?xml version="1.0"?>
<!-- $id$ -->
<field_definition level="1" prec="4" operation="average" enabled=".TRUE." default_value="1.e20" > <!-- time step automaticaly defined -->
......@@ -25,47 +24,123 @@ that are available in the tidal-forcing implementation (see
-->
<!-- Time -->
<field id="diamlr_time" grid_ref="diamlr_grid_T_2D" prec="8" />
<!-- Regressors for tidal harmonic analysis -->
<field id="diamlr_r001" field_ref="diamlr_time" expr="sin( __TDE_M2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:M2" />
<field id="diamlr_r002" field_ref="diamlr_time" expr="cos( __TDE_M2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:M2" />
<field id="diamlr_r003" field_ref="diamlr_time" expr="sin( __TDE_N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:N2" />
<field id="diamlr_r004" field_ref="diamlr_time" expr="cos( __TDE_N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:N2" />
<field id="diamlr_r005" field_ref="diamlr_time" expr="sin( __TDE_2N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:2N2" />
<field id="diamlr_r006" field_ref="diamlr_time" expr="cos( __TDE_2N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:2N2" />
<field id="diamlr_r007" field_ref="diamlr_time" expr="sin( __TDE_S2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:S2" />
<field id="diamlr_r008" field_ref="diamlr_time" expr="cos( __TDE_S2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:S2" />
<field id="diamlr_r009" field_ref="diamlr_time" expr="sin( __TDE_K2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:K2" />
<field id="diamlr_r010" field_ref="diamlr_time" expr="cos( __TDE_K2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:K2" />
<field id="diamlr_r011" field_ref="diamlr_time" expr="sin( __TDE_K1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:K1" />
<field id="diamlr_r012" field_ref="diamlr_time" expr="cos( __TDE_K1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:K1" />
<field id="diamlr_r013" field_ref="diamlr_time" expr="sin( __TDE_O1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:O1" />
<field id="diamlr_r014" field_ref="diamlr_time" expr="cos( __TDE_O1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:O1" />
<field id="diamlr_r015" field_ref="diamlr_time" expr="sin( __TDE_Q1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Q1" />
<field id="diamlr_r016" field_ref="diamlr_time" expr="cos( __TDE_Q1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Q1" />
<field id="diamlr_r017" field_ref="diamlr_time" expr="sin( __TDE_P1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:P1" />
<field id="diamlr_r018" field_ref="diamlr_time" expr="cos( __TDE_P1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:P1" />
<field id="diamlr_r019" field_ref="diamlr_time" expr="sin( __TDE_M4_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:M4" />
<field id="diamlr_r020" field_ref="diamlr_time" expr="cos( __TDE_M4_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:M4" />
<field id="diamlr_r021" field_ref="diamlr_time" expr="sin( __TDE_Mf_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Mf" />
<field id="diamlr_r022" field_ref="diamlr_time" expr="cos( __TDE_Mf_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Mf" />
<field id="diamlr_r023" field_ref="diamlr_time" expr="sin( __TDE_Mm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Mm" />
<field id="diamlr_r024" field_ref="diamlr_time" expr="cos( __TDE_Mm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Mm" />
<field id="diamlr_r025" field_ref="diamlr_time" expr="sin( __TDE_Msqm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Msqm" />
<field id="diamlr_r026" field_ref="diamlr_time" expr="cos( __TDE_Msqm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Msqm" />
<field id="diamlr_r027" field_ref="diamlr_time" expr="sin( __TDE_Mtm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Mtm" />
<field id="diamlr_r028" field_ref="diamlr_time" expr="cos( __TDE_Mtm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Mtm" />
<field id="diamlr_r029" field_ref="diamlr_time" expr="sin( __TDE_S1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:S1" />
<field id="diamlr_r030" field_ref="diamlr_time" expr="cos( __TDE_S1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:S1" />
<field id="diamlr_r031" field_ref="diamlr_time" expr="sin( __TDE_MU2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:MU2" />
<field id="diamlr_r032" field_ref="diamlr_time" expr="cos( __TDE_MU2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:MU2" />
<field id="diamlr_r033" field_ref="diamlr_time" expr="sin( __TDE_NU2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:NU2" />
<field id="diamlr_r034" field_ref="diamlr_time" expr="cos( __TDE_NU2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:NU2" />
<field id="diamlr_r035" field_ref="diamlr_time" expr="sin( __TDE_L2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:L2" />
<field id="diamlr_r036" field_ref="diamlr_time" expr="cos( __TDE_L2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:L2" />
<field id="diamlr_r037" field_ref="diamlr_time" expr="sin( __TDE_T2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:T2" />
<field id="diamlr_r038" field_ref="diamlr_time" expr="cos( __TDE_T2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:T2" />
<field id="diamlr_time" grid_ref="diamlr_grid_T_2D_inner" prec="8" />
<!-- Regressors for tidal harmonic analysis (nn_tide_var=1, default) -->
<!-- long-period constituents -->
<field id="diamlr_r001" field_ref="diamlr_time" expr="sin( __TDE_Mf_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Mf" />
<field id="diamlr_r002" field_ref="diamlr_time" expr="cos( __TDE_Mf_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Mf" />
<field id="diamlr_r003" field_ref="diamlr_time" expr="sin( __TDE_Mm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Mm" />
<field id="diamlr_r004" field_ref="diamlr_time" expr="cos( __TDE_Mm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Mm" />
<field id="diamlr_r005" field_ref="diamlr_time" expr="sin( __TDE_Ssa_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Ssa" />
<field id="diamlr_r006" field_ref="diamlr_time" expr="cos( __TDE_Ssa_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Ssa" />
<field id="diamlr_r007" field_ref="diamlr_time" expr="sin( __TDE_Mtm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Mtm" />
<field id="diamlr_r008" field_ref="diamlr_time" expr="cos( __TDE_Mtm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Mtm" />
<field id="diamlr_r009" field_ref="diamlr_time" expr="sin( __TDE_Msf_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Msf" />
<field id="diamlr_r010" field_ref="diamlr_time" expr="cos( __TDE_Msf_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Msf" />
<field id="diamlr_r011" field_ref="diamlr_time" expr="sin( __TDE_Msqm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Msqm" />
<field id="diamlr_r012" field_ref="diamlr_time" expr="cos( __TDE_Msqm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Msqm" />
<field id="diamlr_r013" field_ref="diamlr_time" expr="sin( __TDE_Sa_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Sa" />
<field id="diamlr_r014" field_ref="diamlr_time" expr="cos( __TDE_Sa_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Sa" />
<!-- diurnal constituents -->
<field id="diamlr_r015" field_ref="diamlr_time" expr="sin( __TDE_K1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:K1" />
<field id="diamlr_r016" field_ref="diamlr_time" expr="cos( __TDE_K1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:K1" />
<field id="diamlr_r017" field_ref="diamlr_time" expr="sin( __TDE_O1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:O1" />
<field id="diamlr_r018" field_ref="diamlr_time" expr="cos( __TDE_O1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:O1" />
<field id="diamlr_r019" field_ref="diamlr_time" expr="sin( __TDE_P1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:P1" />
<field id="diamlr_r020" field_ref="diamlr_time" expr="cos( __TDE_P1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:P1" />
<field id="diamlr_r021" field_ref="diamlr_time" expr="sin( __TDE_Q1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Q1" />
<field id="diamlr_r022" field_ref="diamlr_time" expr="cos( __TDE_Q1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Q1" />
<field id="diamlr_r023" field_ref="diamlr_time" expr="sin( __TDE_J1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:J1" />
<field id="diamlr_r024" field_ref="diamlr_time" expr="cos( __TDE_J1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:J1" />
<field id="diamlr_r025" field_ref="diamlr_time" expr="sin( __TDE_S1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:S1" />
<field id="diamlr_r026" field_ref="diamlr_time" expr="cos( __TDE_S1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:S1" />
<!-- semidiurnal constituents -->
<field id="diamlr_r027" field_ref="diamlr_time" expr="sin( __TDE_M2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:M2" />
<field id="diamlr_r028" field_ref="diamlr_time" expr="cos( __TDE_M2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:M2" />
<field id="diamlr_r029" field_ref="diamlr_time" expr="sin( __TDE_S2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:S2" />
<field id="diamlr_r030" field_ref="diamlr_time" expr="cos( __TDE_S2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:S2" />
<field id="diamlr_r031" field_ref="diamlr_time" expr="sin( __TDE_N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:N2" />
<field id="diamlr_r032" field_ref="diamlr_time" expr="cos( __TDE_N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:N2" />
<field id="diamlr_r033" field_ref="diamlr_time" expr="sin( __TDE_K2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:K2" />
<field id="diamlr_r034" field_ref="diamlr_time" expr="cos( __TDE_K2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:K2" />
<field id="diamlr_r035" field_ref="diamlr_time" expr="sin( __TDE_nu2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:nu2" />
<field id="diamlr_r036" field_ref="diamlr_time" expr="cos( __TDE_nu2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:nu2" />
<field id="diamlr_r037" field_ref="diamlr_time" expr="sin( __TDE_mu2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:mu2" />
<field id="diamlr_r038" field_ref="diamlr_time" expr="cos( __TDE_mu2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:mu2" />
<field id="diamlr_r039" field_ref="diamlr_time" expr="sin( __TDE_2N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:2N2" />
<field id="diamlr_r040" field_ref="diamlr_time" expr="cos( __TDE_2N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:2N2" />
<field id="diamlr_r041" field_ref="diamlr_time" expr="sin( __TDE_L2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:L2" />
<field id="diamlr_r042" field_ref="diamlr_time" expr="cos( __TDE_L2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:L2" />
<field id="diamlr_r043" field_ref="diamlr_time" expr="sin( __TDE_T2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:T2" />
<field id="diamlr_r044" field_ref="diamlr_time" expr="cos( __TDE_T2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:T2" />
<field id="diamlr_r045" field_ref="diamlr_time" expr="sin( __TDE_eps2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:eps2" />
<field id="diamlr_r046" field_ref="diamlr_time" expr="cos( __TDE_eps2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:eps2" />
<field id="diamlr_r047" field_ref="diamlr_time" expr="sin( __TDE_lam2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:lam2" />
<field id="diamlr_r048" field_ref="diamlr_time" expr="cos( __TDE_lam2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:lam2" />
<field id="diamlr_r049" field_ref="diamlr_time" expr="sin( __TDE_R2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:R2" />
<field id="diamlr_r050" field_ref="diamlr_time" expr="cos( __TDE_R2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:R2" />
<!-- terdiurnal constituent -->
<field id="diamlr_r051" field_ref="diamlr_time" expr="sin( __TDE_M3_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:M3" />
<field id="diamlr_r052" field_ref="diamlr_time" expr="cos( __TDE_M3_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:M3" />
<!-- compound tides -->
<field id="diamlr_r053" field_ref="diamlr_time" expr="sin( __TDE_MKS2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:MKS2" />
<field id="diamlr_r054" field_ref="diamlr_time" expr="cos( __TDE_MKS2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:MKS2" />
<field id="diamlr_r055" field_ref="diamlr_time" expr="sin( __TDE_MN4_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:MN4" />
<field id="diamlr_r056" field_ref="diamlr_time" expr="cos( __TDE_MN4_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:MN4" />
<field id="diamlr_r057" field_ref="diamlr_time" expr="sin( __TDE_MS4_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:MS4" />
<field id="diamlr_r058" field_ref="diamlr_time" expr="cos( __TDE_MS4_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:MS4" />
<!-- overtides -->
<field id="diamlr_r059" field_ref="diamlr_time" expr="sin( __TDE_M4_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:M4" />
<field id="diamlr_r060" field_ref="diamlr_time" expr="cos( __TDE_M4_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:M4" />
<field id="diamlr_r061" field_ref="diamlr_time" expr="sin( __TDE_N4_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:N4" />
<field id="diamlr_r062" field_ref="diamlr_time" expr="cos( __TDE_N4_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:N4" />
<field id="diamlr_r063" field_ref="diamlr_time" expr="sin( __TDE_S4_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:S4" />
<field id="diamlr_r064" field_ref="diamlr_time" expr="cos( __TDE_S4_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:S4" />
<field id="diamlr_r065" field_ref="diamlr_time" expr="sin( __TDE_M6_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:M6" />
<field id="diamlr_r066" field_ref="diamlr_time" expr="cos( __TDE_M6_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:M6" />
<field id="diamlr_r067" field_ref="diamlr_time" expr="sin( __TDE_M8_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:M8" />
<field id="diamlr_r068" field_ref="diamlr_time" expr="cos( __TDE_M8_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:M8" />
<!-- legacy option (nn_tide_var=0)
<field id="diamlr_r001" field_ref="diamlr_time" expr="sin( __TDE_M2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:M2" />
<field id="diamlr_r002" field_ref="diamlr_time" expr="cos( __TDE_M2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:M2" />
<field id="diamlr_r003" field_ref="diamlr_time" expr="sin( __TDE_N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:N2" />
<field id="diamlr_r004" field_ref="diamlr_time" expr="cos( __TDE_N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:N2" />
<field id="diamlr_r005" field_ref="diamlr_time" expr="sin( __TDE_2N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:2N2" />
<field id="diamlr_r006" field_ref="diamlr_time" expr="cos( __TDE_2N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:2N2" />
<field id="diamlr_r007" field_ref="diamlr_time" expr="sin( __TDE_S2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:S2" />
<field id="diamlr_r008" field_ref="diamlr_time" expr="cos( __TDE_S2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:S2" />
<field id="diamlr_r009" field_ref="diamlr_time" expr="sin( __TDE_K2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:K2" />
<field id="diamlr_r010" field_ref="diamlr_time" expr="cos( __TDE_K2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:K2" />
<field id="diamlr_r011" field_ref="diamlr_time" expr="sin( __TDE_K1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:K1" />
<field id="diamlr_r012" field_ref="diamlr_time" expr="cos( __TDE_K1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:K1" />
<field id="diamlr_r013" field_ref="diamlr_time" expr="sin( __TDE_O1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:O1" />
<field id="diamlr_r014" field_ref="diamlr_time" expr="cos( __TDE_O1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:O1" />
<field id="diamlr_r015" field_ref="diamlr_time" expr="sin( __TDE_Q1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Q1" />
<field id="diamlr_r016" field_ref="diamlr_time" expr="cos( __TDE_Q1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Q1" />
<field id="diamlr_r017" field_ref="diamlr_time" expr="sin( __TDE_P1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:P1" />
<field id="diamlr_r018" field_ref="diamlr_time" expr="cos( __TDE_P1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:P1" />
<field id="diamlr_r019" field_ref="diamlr_time" expr="sin( __TDE_M4_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:M4" />
<field id="diamlr_r020" field_ref="diamlr_time" expr="cos( __TDE_M4_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:M4" />
<field id="diamlr_r021" field_ref="diamlr_time" expr="sin( __TDE_Mf_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Mf" />
<field id="diamlr_r022" field_ref="diamlr_time" expr="cos( __TDE_Mf_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Mf" />
<field id="diamlr_r023" field_ref="diamlr_time" expr="sin( __TDE_Mm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Mm" />
<field id="diamlr_r024" field_ref="diamlr_time" expr="cos( __TDE_Mm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Mm" />
<field id="diamlr_r025" field_ref="diamlr_time" expr="sin( __TDE_Msqm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Msqm" />
<field id="diamlr_r026" field_ref="diamlr_time" expr="cos( __TDE_Msqm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Msqm" />
<field id="diamlr_r027" field_ref="diamlr_time" expr="sin( __TDE_Mtm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Mtm" />
<field id="diamlr_r028" field_ref="diamlr_time" expr="cos( __TDE_Mtm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Mtm" />
<field id="diamlr_r029" field_ref="diamlr_time" expr="sin( __TDE_S1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:S1" />
<field id="diamlr_r030" field_ref="diamlr_time" expr="cos( __TDE_S1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:S1" />
<field id="diamlr_r031" field_ref="diamlr_time" expr="sin( __TDE_MU2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:MU2" />
<field id="diamlr_r032" field_ref="diamlr_time" expr="cos( __TDE_MU2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:MU2" />
<field id="diamlr_r033" field_ref="diamlr_time" expr="sin( __TDE_NU2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:NU2" />
<field id="diamlr_r034" field_ref="diamlr_time" expr="cos( __TDE_NU2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:NU2" />
<field id="diamlr_r035" field_ref="diamlr_time" expr="sin( __TDE_L2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:L2" />
<field id="diamlr_r036" field_ref="diamlr_time" expr="cos( __TDE_L2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:L2" />
<field id="diamlr_r037" field_ref="diamlr_time" expr="sin( __TDE_T2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:T2" />
<field id="diamlr_r038" field_ref="diamlr_time" expr="cos( __TDE_T2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:T2" />
-->
<field id="diamlr_r101" field_ref="diamlr_time" expr="diamlr_time^0.0" enabled=".TRUE." comment="mean" />
<!-- Fields selected for regression analysis -->
......@@ -109,76 +184,74 @@ that are available in the tidal-forcing implementation (see
<!-- 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="e3ts" long_name="T-cell thickness" field_ref="e3t" standard_name="cell_thickness" unit="m" grid_ref="grid_T_SFC"/>
<field id="e3t_0" long_name="Initial T-cell thickness" standard_name="ref_cell_thickness" unit="m" grid_ref="grid_T_3D" />
<field id="e3tb" long_name="bottom T-cell thickness" standard_name="bottom_cell_thickness" unit="m" grid_ref="grid_T_2D"/>
<field id="e3t_300" field_ref="e3t" grid_ref="grid_T_zoom_300" detect_missing_value="true" />
<field id="e3t_vsum300" field_ref="e3t_300" grid_ref="grid_T_vsum" detect_missing_value="true" />
<field id="masscello" long_name="Sea Water Mass per unit area" standard_name="sea_water_mass_per_unit_area" unit="kg/m2" grid_ref="grid_T_3D"/>
<field id="volcello" long_name="Ocean Volume" standard_name="ocean_volume" unit="m3" grid_ref="grid_T_3D"/>
<field id="toce" long_name="temperature" standard_name="sea_water_potential_temperature" unit="degC" grid_ref="grid_T_3D"/>
<field id="toce_e3t" long_name="temperature (thickness weighted)" unit="degC" grid_ref="grid_T_3D" > toce * e3t </field >
<field id="soce" long_name="salinity" standard_name="sea_water_practical_salinity" unit="1e-3" grid_ref="grid_T_3D"/>
<field id="soce_e3t" long_name="salinity (thickness weighted)" unit="1e-3" grid_ref="grid_T_3D" > soce * e3t </field >
<field id="toce_e3t_300" field_ref="toce_e3t" unit="degree_C" grid_ref="grid_T_zoom_300" detect_missing_value="true" />
<field id="toce_e3t_vsum300" field_ref="toce_e3t_300" unit="degress_C*m" grid_ref="grid_T_vsum" detect_missing_value="true" />
<field id="toce_vmean300" field_ref="toce_e3t_vsum300" unit="degree_C" grid_ref="grid_T_vsum" detect_missing_value="true" > toce_e3t_vsum300/e3t_vsum300 </field>
<field id="e3t" long_name="T-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_T_3D_inner" />
<field id="e3ts" long_name="T-cell thickness" field_ref="e3t" standard_name="cell_thickness" unit="m" grid_ref="grid_T_SFC_inner" />
<field id="e3t_0" long_name="Initial T-cell thickness" standard_name="ref_cell_thickness" unit="m" grid_ref="grid_T_3D_inner" />
<field id="e3tb" long_name="bottom T-cell thickness" standard_name="bottom_cell_thickness" unit="m" grid_ref="grid_T_2D_inner" />
<field id="e3t_300" field_ref="e3t" grid_ref="grid_T_zoom_300_inner" detect_missing_value="true" />
<field id="e3t_vsum300" field_ref="e3t_300" grid_ref="grid_T_vsum_inner" 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_inner"/>
<field id="volcello" long_name="Ocean Volume" standard_name="ocean_volume" unit="m3" grid_ref="grid_T_3D_inner"/>
<field id="toce" long_name="temperature" standard_name="sea_water_potential_temperature" unit="degC" grid_ref="grid_T_3D"/>
<field id="toce_e3t" long_name="temperature (thickness weighted)" unit="degC" grid_ref="grid_T_3D" > toce * e3t </field >
<field id="soce" long_name="salinity" standard_name="sea_water_practical_salinity" unit="1e-3" grid_ref="grid_T_3D"/>
<field id="soce_e3t" long_name="salinity (thickness weighted)" unit="1e-3" grid_ref="grid_T_3D" > soce * e3t </field >
<field id="toce_e3t_300" field_ref="toce_e3t" unit="degree_C" grid_ref="grid_T_zoom_300" detect_missing_value="true" />
<field id="toce_e3t_vsum300" field_ref="toce_e3t_300" unit="degress_C*m" grid_ref="grid_T_vsum" detect_missing_value="true" />
<field id="toce_vmean300" field_ref="toce_e3t_vsum300" unit="degree_C" grid_ref="grid_T_vsum" detect_missing_value="true" > toce_e3t_vsum300/e3t_vsum300 </field>
<!-- AGRIF sponge -->
<field id="agrif_spt" long_name=" AGRIF t-sponge coefficient" unit=" " />
<!-- additions to diawri.F90 -->
<field id="sssgrad" long_name="module of surface salinity gradient" unit="1e-3/m" grid_ref="grid_T_2D_inner"/>
<field id="sssgrad2" long_name="square of module of surface salinity gradient" unit="1e-6/m2" grid_ref="grid_T_2D_inner"/>
<field id="ke" long_name="kinetic energy" standard_name="specific_kinetic_energy_of_sea_water" unit="m2/s2" grid_ref="grid_T_3D" />
<field id="ke_int" long_name="vertical integration of kinetic energy" unit="m3/s2" grid_ref="grid_T_2D_inner" />
<field id="sssgrad" long_name="module of surface salinity gradient" unit="1e-3/m" grid_ref="grid_T_2D_inner" />
<field id="sssgrad2" long_name="square of module of surface salinity gradient" unit="1e-6/m2" grid_ref="grid_T_2D_inner" />
<field id="ke" long_name="kinetic energy" standard_name="specific_kinetic_energy_of_sea_water" unit="m2/s2" grid_ref="grid_T_3D_inner" />
<field id="ke_int" long_name="vertical integration of kinetic energy" unit="m3/s2" grid_ref="grid_T_2D_inner" />
<field id="taubot" long_name="bottom stress module" unit="N/m2" grid_ref="grid_T_2D_inner" />
<!-- t-eddy viscosity coefficients (ldfdyn) -->
<field id="ahmt_2d" long_name=" surface t-eddy viscosity coefficient" unit="m2/s or m4/s" />
<field id="ahmt_3d" long_name=" 3D t-eddy viscosity coefficient" unit="m2/s or m4/s" grid_ref="grid_T_3D"/>
<field id="sst" long_name="Bulk sea surface temperature" standard_name="bulk_sea_surface_temperature" unit="degC" />
<field id="t_skin" long_name="Skin temperature aka SSST" standard_name="skin_temperature" unit="degC" />
<field id="sss" long_name="sea surface salinity" standard_name="sea_surface_salinity" unit="1e-3" />
<field id="sst2" long_name="square of sea surface temperature" standard_name="square_of_sea_surface_temperature" unit="degC2" > sst * sst </field >
<field id="sss2" long_name="square of sea surface salinity" unit="1e-6" > sss * sss </field >
<field id="sstmax" long_name="max of sea surface temperature" field_ref="sst" operation="maximum" />
<field id="sssmax" long_name="max of sea surface salinity" field_ref="sss" operation="maximum" />
<field id="sstmin" long_name="min of sea surface temperature" field_ref="sst" operation="minimum" />
<field id="sssmin" long_name="min of sea surface salinity" field_ref="sss" operation="minimum" />
<field id="sstgrad" long_name="module of sst gradient" unit="degC/m" grid_ref="grid_T_2D_inner" />
<field id="sstgrad2" long_name="square of module of sst gradient" unit="degC2/m2" grid_ref="grid_T_2D_inner" />
<field id="sbt" long_name="sea bottom temperature" unit="degC" grid_ref="grid_T_2D_inner" />
<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" grid_ref="grid_T_2D_inner" />
<field id="sst_wl" long_name="Delta SST of warm layer" unit="degC" />
<field id="sst_cs" long_name="Delta SST of cool skin" unit="degC" />
<field id="temp_3m" long_name="temperature at 3m" unit="degC" />
<field id="sss" long_name="sea surface salinity" standard_name="sea_surface_salinity" unit="1e-3" />
<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" grid_ref="grid_T_2D_inner" />
<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)" grid_ref="grid_T_2D_inner" />
<field id="sbs" long_name="sea bottom salinity" unit="0.001" grid_ref="grid_T_2D_inner" />
<field id="sst_wl" long_name="Delta SST of warm layer" unit="degC" grid_ref="grid_T_2D_inner" />
<field id="sst_cs" long_name="Delta SST of cool skin" unit="degC" grid_ref="grid_T_2D_inner" />
<field id="taubot" long_name="bottom stress module" unit="N/m2" grid_ref="grid_T_2D_inner" />
<field id="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" grid_ref="grid_T_2D_inner" />
<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)" grid_ref="grid_T_2D_inner" />
<!-- Case EOS = TEOS-10 : output potential temperature -->
<field id="toce_pot" long_name="Sea Water Potential Temperature" standard_name="sea_water_potential_temperature" unit="degC" grid_ref="grid_T_3D"/>
<field id="sst_pot" long_name="potential sea surface temperature" standard_name="sea_surface_temperature" unit="degC" />
<field id="tosmint_pot" long_name="vertical integral of potential temperature times density" standard_name="integral_wrt_depth_of_product_of_density_and_potential_temperature" unit="(kg m2) degree_C" />
<field id="toce_pot" long_name="Sea Water Potential Temperature" standard_name="sea_water_potential_temperature" unit="degC" grid_ref="grid_T_3D_inner"/>
<field id="sst_pot" long_name="potential sea surface temperature" standard_name="sea_surface_temperature" unit="degC" grid_ref="grid_T_2D_inner"/>
<field id="tosmint_pot" long_name="vertical integral of potential temperature times density" standard_name="integral_wrt_depth_of_product_of_density_and_potential_temperature" unit="(kg m2) degree_C" grid_ref="grid_T_2D_inner"/>
<field id="ht" long_name="water column height at T point" standard_name="water_column_height_T" unit="m" />
<field id="ssh" long_name="sea surface height" standard_name="sea_surface_height_above_geoid" unit="m" />
<field id="ssh2" long_name="square of sea surface height" standard_name="square_of_sea_surface_height_above_geoid" unit="m2" > ssh * ssh </field >
<field id="wetdep" long_name="wet depth" standard_name="wet_depth" unit="m" />
<field id="sshmax" long_name="max of sea surface height" field_ref="ssh" operation="maximum" />
<field id="sshmax" long_name="max of sea surface height" field_ref="ssh" operation="maximum" grid_ref="grid_T_2D_inner"/>
<field id="mldkz5" long_name="Turbocline depth (Kz = 5e-4)" standard_name="ocean_mixed_layer_thickness_defined_by_vertical_tracer_diffusivity" unit="m" />
<field id="mldr10_1" long_name="Mixed Layer Depth (dsigma = 0.01 wrt 10m)" standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta" unit="m" />
<field id="mldr10_1max" long_name="Max of Mixed Layer Depth (dsigma = 0.01 wrt 10m)" field_ref="mldr10_1" operation="maximum" />
<field id="mldr10_1min" long_name="Min of Mixed Layer Depth (dsigma = 0.01 wrt 10m)" field_ref="mldr10_1" operation="minimum" />
<field id="heatc" long_name="Heat content vertically integrated" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" grid_ref="grid_T_2D_inner" />
<field id="saltc" long_name="Salt content vertically integrated" unit="PSU*kg/m2" grid_ref="grid_T_2D_inner" />
<field id="salt2c" long_name="square of Salt content vertically integrated" unit="PSU2*kg/m2" grid_ref="grid_T_2D_inner" />
<field id="mldkz5" long_name="Turbocline depth (Kz = 5e-4)" standard_name="ocean_mixed_layer_thickness_defined_by_vertical_tracer_diffusivity" unit="m" grid_ref="grid_T_2D_inner" />
<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" grid_ref="grid_T_2D_inner" />
<field id="mldr10_1max" long_name="Max of Mixed Layer Depth (dsigma = 0.01 wrt 10m)" field_ref="mldr10_1" operation="maximum" grid_ref="grid_T_2D_inner" />
<field id="mldr10_1min" long_name="Min of Mixed Layer Depth (dsigma = 0.01 wrt 10m)" field_ref="mldr10_1" operation="minimum" grid_ref="grid_T_2D_inner" />
<field id="heatc" long_name="Heat content vertically integrated" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" grid_ref="grid_T_2D_inner" />
<field id="saltc" long_name="Salt content vertically integrated" unit="PSU*kg/m2" grid_ref="grid_T_2D_inner" />
<field id="salt2c" long_name="square of Salt content vertically integrated" unit="PSU2*kg/m2" grid_ref="grid_T_2D_inner" />
<!-- EOS -->
<field id="alpha" long_name="thermal expansion" unit="degC-1" grid_ref="grid_T_3D" />
......@@ -187,98 +260,113 @@ that are available in the tidal-forcing implementation (see
<!-- Energy - horizontal divergence -->
<field id="sKE" long_name="surface kinetic energy" standard_name="specific_kinetic_energy_of_sea_water" unit="m2/s2" grid_ref="grid_T_2D_inner" />
<field id="hdiv" long_name="horizontal divergence" unit="s-1" grid_ref="grid_T_3D" />
<field id="hdiv" long_name="horizontal divergence" unit="s-1" grid_ref="grid_T_3D_inner" />
<!-- variables available with MLE -->
<field id="Lf_NHpf" long_name="MLE: Lf = N H / f" unit="m" />
<field id="Lf_NHpf" long_name="MLE: Lf = N H / f" unit="m" grid_ref="grid_T_2D_inner" />
<!-- next variables available with ln_zad_Aimp=.true. -->
<field id="Courant" long_name="Courant number" unit="#" grid_ref="grid_T_3D" />
<field id="wimp" long_name="Implicit vertical velocity" unit="m/s" grid_ref="grid_T_3D" />
<field id="Courant" long_name="Courant number" unit="#" grid_ref="grid_T_3D_inner" />
<field id="wimp" long_name="Implicit vertical velocity" unit="m/s" grid_ref="grid_T_3D_inner" />
<field id="wexp" long_name="Explicit vertical velocity" unit="m/s" grid_ref="grid_T_3D" />
<field id="wi_cff" long_name="Fraction of implicit vertical velocity" unit="#" grid_ref="grid_T_3D" />
<field id="wi_cff" long_name="Fraction of implicit vertical velocity" unit="#" grid_ref="grid_T_3D_inner" />
<field id="Aimp_loc" long_name="Locations of columns with active adaptive-implicit vertical advection" unit="%ts" grid_ref="grid_T_2D_inner" />
<field id="Aimp_Cmx" long_name="Column maximum of Courant number with active adaptive-implicit vertical advection" operation="maximum" unit="#" grid_ref="grid_T_2D_inner" />
<field id="Aimp_Cmx_h" long_name="Column maximum of horizontal Courant number with active adaptive-implicit vertical advection" operation="maximum" unit="#" grid_ref="grid_T_2D_inner" />
<field id="Aimp_Cmx_v" long_name="Column maximum of vertical Courant number with active adaptive-implicit vertical advection" operation="maximum" unit="#" grid_ref="grid_T_2D" />
<!-- next variables available with key_diahth -->
<field id="mlddzt" long_name="Thermocline Depth (depth of max dT/dz)" standard_name="depth_at_maximum_upward_derivative_of_sea_water_potential_temperature" unit="m" />
<field id="mldr10_3" long_name="Mixed Layer Depth (dsigma = 0.03 wrt 10m)" standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta" unit="m" />
<field id="mldr0_1" long_name="Mixed Layer Depth (dsigma = 0.01 wrt sfc)" standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta" unit="m" />
<field id="mldr0_3" long_name="Mixed Layer Depth (dsigma = 0.03 wrt sfc)" standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta" unit="m" />
<field id="mld_dt02" long_name="Mixed Layer Depth (|dT| = 0.2 wrt 10m)" standard_name="ocean_mixed_layer_thickness_defined_by_temperature" unit="m" />
<field id="topthdep" long_name="Top of Thermocline Depth (dT = -0.2 wrt 10m)" standard_name="ocean_mixed_layer_thickness_defined_by_temperature" unit="m" />
<field id="pycndep" long_name="Pycnocline Depth (dsigma[dT=-0.2] wrt 10m)" standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta" unit="m" />
<field id="BLT" long_name="Barrier Layer Thickness" unit="m" > topthdep - pycndep </field>
<field id="tinv" long_name="Max of vertical invertion of temperature" unit="degC" />
<field id="depti" long_name="Depth of max. vert. inv. of temperature" unit="m" />
<field id="20d" long_name="Depth of 20C isotherm" standard_name="depth_of_isosurface_of_sea_water_potential_temperature" unit="m" axis_ref="iax_20C" />
<field id="26d" long_name="Depth of 26C isotherm" standard_name="depth_of_isosurface_of_sea_water_potential_temperature" unit="m" axis_ref="iax_26C" />
<field id="28d" long_name="Depth of 28C isotherm" standard_name="depth_of_isosurface_of_sea_water_potential_temperature" unit="m" axis_ref="iax_28C" />
<field id="hc300" long_name="Heat content 0-300m" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" />
<field id="hc700" long_name="Heat content 0-700m" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" />
<field id="hc2000" long_name="Heat content 0-2000m" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" />
<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" grid_ref="grid_T_2D_inner" />
<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" grid_ref="grid_T_2D_inner" />
<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" grid_ref="grid_T_2D_inner" />
<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" grid_ref="grid_T_2D_inner" />
<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" grid_ref="grid_T_2D_inner" />
<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" grid_ref="grid_T_2D_inner" />
<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" grid_ref="grid_T_2D_inner" />
<field id="BLT" long_name="Barrier Layer Thickness" unit="m" grid_ref="grid_T_2D_inner" > topthdep - pycndep </field>
<field id="tinv" long_name="Max of vertical invertion of temperature" unit="degC" grid_ref="grid_T_2D_inner" />
<field id="depti" long_name="Depth of max. vert. inv. of temperature" unit="m" grid_ref="grid_T_2D_inner" />
<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" grid_ref="grid_T_2D_inner" />
<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" grid_ref="grid_T_2D_inner" />
<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" grid_ref="grid_T_2D_inner" />
<!-- variables available with diaar5 -->
<field id="botpres" long_name="Sea Water Pressure at Sea Floor" standard_name="sea_water_pressure_at_sea_floor" unit="dbar" />
<field id="botpres" long_name="Sea Water Pressure at Sea Floor" standard_name="sea_water_pressure_at_sea_floor" unit="dbar" grid_ref="grid_T_2D_inner" />
<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" />
<field id="tnpeo" long_name="Tendency of ocean potential energy content" unit="W/m2" grid_ref="grid_T_2D_inner" />
<!-- variables available ln_linssh=.FALSE. -->
<field id="tpt_dep" long_name="T-point depth" standard_name="depth_below_geoid" unit="m" grid_ref="grid_T_3D" />
<field id="e3tdef" long_name="T-cell thickness deformation" unit="%" grid_ref="grid_T_3D" />
<field id="tpt_dep" long_name="T-point depth" standard_name="depth_below_geoid" unit="m" grid_ref="grid_T_3D_inner" />
<field id="e3tdef" long_name="T-cell thickness deformation" unit="%" grid_ref="grid_T_3D_inner" />
<!-- variables available with ln_diacfl=.true. -->
<field id="cfl_cu" long_name="u-courant number" unit="#" />
<field id="cfl_cv" long_name="v-courant number" unit="#" />
<field id="cfl_cw" long_name="w-courant number" unit="#" />
<field id="cfl_cu" long_name="u-courant number" unit="#" grid_ref="grid_T_2D_inner" />
<field id="cfl_cv" long_name="v-courant number" unit="#" grid_ref="grid_T_2D_inner" />
<field id="cfl_cw" long_name="w-courant number" unit="#" grid_ref="grid_T_2D_inner" />
<!-- GEOMETRIC fields (requires nn_aei_ijk_t = 32) -->
<field id="eke" long_name="total EKE (EKE+EPE)" unit="m3/s2" />
<field id="trd_eke_adv_ubt" long_name="ubt advective trend of EKE (LHS)" unit="m3/s3" grid_ref="grid_T_2D_inner" />
<field id="trd_eke_adv_wav" long_name="wav advective trend of EKE (LHS)" unit="m3/s3" grid_ref="grid_T_2D_inner" />
<field id="trd_eke_lap" long_name="diffusive trend of EKE (RHS)" unit="m3/s3" grid_ref="grid_T_2D_inner" />
<field id="trd_eke_peS" long_name="PE to EKE source trend (RHS)" unit="m3/s3" grid_ref="grid_T_2D_inner" />
<field id="trd_eke_keS" long_name="KE to EKE source trend (RHS)" unit="m3/s3" grid_ref="grid_T_2D_inner" />
<field id="trd_eke_dis" long_name="dissipation trend of EKE (RHS)" unit="m3/s3" grid_ref="grid_T_2D_inner" />
<!-- variables available with ln_zdfmfc=.true. -->
<field id="mf_Tp" long_name="plume_temperature" standard_name="plume_temperature" unit="degC" grid_ref="grid_T_3D" />
<field id="mf_Sp" long_name="plume_salinity" standard_name="plume_salinity" unit="1e-3" grid_ref="grid_T_3D" />
<field id="mf_mf" long_name="mass flux" standard_name="mf_mass_flux" unit="m" grid_ref="grid_T_3D" />
<field id="mf_Tp" long_name="plume_temperature" standard_name="plume_temperature" unit="degC" grid_ref="grid_T_3D_inner" />
<field id="mf_Sp" long_name="plume_salinity" standard_name="plume_salinity" unit="1e-3" grid_ref="grid_T_3D_inner" />
<field id="mf_mf" long_name="mass flux" standard_name="mf_mass_flux" unit="m" grid_ref="grid_T_3D_inner" />
<!-- fluxes from damping -->
<field id="sflx_dmp_cea" long_name="salt flux due to damping" standard_name="salt_flux_due_to_damping" unit="g/m2/s" />
<field id="hflx_dmp_cea" long_name="heat flux due to damping" standard_name="heat_flux_due_to_damping" unit="W/m2" />
<field id="sflx_dmp_cea" long_name="salt flux due to damping" standard_name="salt_flux_due_to_damping" unit="g/m2/s" grid_ref="grid_T_2D_inner" />
<field id="hflx_dmp_cea" long_name="heat flux due to damping" standard_name="heat_flux_due_to_damping" unit="W/m2" grid_ref="grid_T_2D_inner" />
<!-- * variable related to ice shelf forcing * -->
<!-- * fwf * -->
<field id="fwfisf_cav" long_name="Ice shelf fresh water flux ( from isf to oce )" unit="kg/m2/s" />
<field id="fwfisf_par" long_name="Ice shelf fresh water flux ( from isf to oce )" unit="kg/m2/s" />
<field id="fwfisf3d_cav" long_name="3d Ice shelf fresh water flux ( from isf to oce )" unit="kg/m2/s" grid_ref="grid_T_3D" />
<field id="fwfisf3d_par" long_name="3d Ice shelf fresh water flux ( from isf to oce )" unit="kg/m2/s" grid_ref="grid_T_3D" />
<field id="fwfisf_cav" long_name="Ice shelf fresh water flux ( from isf to oce )" unit="kg/m2/s" grid_ref="grid_T_2D_inner" />
<field id="fwfisf_par" long_name="Ice shelf fresh water flux ( from isf to oce )" unit="kg/m2/s" grid_ref="grid_T_2D_inner" />
<field id="fwfisf3d_cav" long_name="3d Ice shelf fresh water flux ( from isf to oce )" unit="kg/m2/s" grid_ref="grid_T_3D_inner" />
<field id="fwfisf3d_par" long_name="3d Ice shelf fresh water flux ( from isf to oce )" unit="kg/m2/s" grid_ref="grid_T_3D_inner" />
<!-- * heat fluxes * -->
<field id="qoceisf_cav" long_name="Ice shelf ocean heat flux ( from isf to oce )" unit="W/m2" />
<field id="qoceisf_par" long_name="Ice shelf ocean heat flux ( from isf to oce )" unit="W/m2" />
<field id="qlatisf_cav" long_name="Ice shelf latent heat flux ( from isf to oce )" unit="W/m2" />
<field id="qlatisf_par" long_name="Ice shelf latent heat flux ( from isf to oce )" unit="W/m2" />
<field id="qhcisf_cav" long_name="Ice shelf heat content flux of injected water ( from isf to oce )" unit="W/m2" />
<field id="qhcisf_par" long_name="Ice shelf heat content flux of injected water ( from isf to oce )" unit="W/m2" />
<field id="qoceisf3d_cav" long_name="Ice shelf ocean heat flux ( from isf to oce )" unit="W/m2" grid_ref="grid_T_3D" />
<field id="qoceisf3d_par" long_name="Ice shelf ocean heat flux ( from isf to oce )" unit="W/m2" grid_ref="grid_T_3D" />
<field id="qlatisf3d_cav" long_name="Ice shelf latent heat flux ( from isf to oce )" unit="W/m2" grid_ref="grid_T_3D" />
<field id="qlatisf3d_par" long_name="Ice shelf latent heat flux ( from isf to oce )" unit="W/m2" grid_ref="grid_T_3D" />
<field id="qhcisf3d_cav" long_name="Ice shelf heat content flux of injected water ( from isf to oce )" unit="W/m2" grid_ref="grid_T_3D" />
<field id="qhcisf3d_par" long_name="Ice shelf heat content flux of injected water ( from isf to oce )" unit="W/m2" grid_ref="grid_T_3D" />
<field id="qconisf" long_name="Conductive heat flux through the ice shelf ( from isf to oce )" unit="W/m2" />
<field id="qoceisf_cav" long_name="Ice shelf ocean heat flux ( from isf to oce )" unit="W/m2" grid_ref="grid_T_2D_inner" />
<field id="qoceisf_par" long_name="Ice shelf ocean heat flux ( from isf to oce )" unit="W/m2" grid_ref="grid_T_2D_inner" />
<field id="qlatisf_cav" long_name="Ice shelf latent heat flux ( from isf to oce )" unit="W/m2" grid_ref="grid_T_2D_inner" />
<field id="qlatisf_par" long_name="Ice shelf latent heat flux ( from isf to oce )" unit="W/m2" grid_ref="grid_T_2D_inner" />
<field id="qhcisf_cav" long_name="Ice shelf heat content flux of injected water ( from isf to oce )" unit="W/m2" grid_ref="grid_T_2D_inner" />
<field id="qhcisf_par" long_name="Ice shelf heat content flux of injected water ( from isf to oce )" unit="W/m2" grid_ref="grid_T_2D_inner" />
<field id="qoceisf3d_cav" long_name="Ice shelf ocean heat flux ( from isf to oce )" unit="W/m2" grid_ref="grid_T_3D_inner" />
<field id="qoceisf3d_par" long_name="Ice shelf ocean heat flux ( from isf to oce )" unit="W/m2" grid_ref="grid_T_3D_inner" />
<field id="qlatisf3d_cav" long_name="Ice shelf latent heat flux ( from isf to oce )" unit="W/m2" grid_ref="grid_T_3D_inner" />
<field id="qlatisf3d_par" long_name="Ice shelf latent heat flux ( from isf to oce )" unit="W/m2" grid_ref="grid_T_3D_inner" />
<field id="qhcisf3d_cav" long_name="Ice shelf heat content flux of injected water ( from isf to oce )" unit="W/m2" grid_ref="grid_T_3D_inner" />
<field id="qhcisf3d_par" long_name="Ice shelf heat content flux of injected water ( from isf to oce )" unit="W/m2" grid_ref="grid_T_3D_inner" />
<field id="qconisf" long_name="Conductive heat flux through the ice shelf ( from isf to oce )" unit="W/m2" grid_ref="grid_T_2D_inner" />
<!-- top boundary layer properties -->
<field id="isftfrz_cav" long_name="freezing point temperature at ocean/isf interface" unit="degC" />
<field id="isftfrz_par" long_name="freezing point temperature in the parametrization boundary layer" unit="degC" />
<field id="isfthermald_cav" long_name="thermal driving of ice shelf melting" unit="degC" />
<field id="isfthermald_par" long_name="thermal driving of ice shelf melting" unit="degC" />
<field id="isfgammat" long_name="Ice shelf heat-transfert velocity" unit="m/s" />
<field id="isfgammas" long_name="Ice shelf salt-transfert velocity" unit="m/s" />
<field id="ttbl_cav" long_name="temperature in Losch tbl" unit="degC" />
<field id="ttbl_par" long_name="temperature in the parametrisation boundary layer" unit="degC" />
<field id="stbl" long_name="salinity in the Losh tbl" unit="1e-3" />
<field id="utbl" long_name="zonal current in the Losh tbl at T point" unit="m/s" />
<field id="vtbl" long_name="merid current in the Losh tbl at T point" unit="m/s" />
<field id="isfustar" long_name="ustar at T point used in ice shelf melting" unit="m/s" />
<field id="isftfrz_cav" long_name="freezing point temperature at ocean/isf interface" unit="degC" grid_ref="grid_T_2D_inner" />
<field id="isftfrz_par" long_name="freezing point temperature in the parametrization boundary layer" unit="degC" grid_ref="grid_T_2D_inner" />
<field id="isfthermald_cav" long_name="thermal driving of ice shelf melting" unit="degC" grid_ref="grid_T_2D_inner" />
<field id="isfthermald_par" long_name="thermal driving of ice shelf melting" unit="degC" grid_ref="grid_T_2D_inner" />
<field id="isfgammat" long_name="Ice shelf heat-transfert velocity" unit="m/s" grid_ref="grid_T_2D_inner" />
<field id="isfgammas" long_name="Ice shelf salt-transfert velocity" unit="m/s" grid_ref="grid_T_2D_inner" />
<field id="ttbl_cav" long_name="temperature in Losch tbl" unit="degC" grid_ref="grid_T_2D_inner" />
<field id="ttbl_par" long_name="temperature in the parametrisation boundary layer" unit="degC" grid_ref="grid_T_2D_inner" />
<field id="stbl" long_name="salinity in the Losh tbl" unit="1e-3" grid_ref="grid_T_2D_inner" />
<field id="utbl" long_name="zonal current in the Losh tbl" unit="m/s" grid_ref="grid_T_2D_inner" />
<field id="vtbl" long_name="meridional current in the Losh tbl" unit="m/s" grid_ref="grid_T_2D_inner" />
<field id="isfustar" long_name="ustar at T point used in ice shelf melting" unit="m/s" grid_ref="grid_T_2D_inner" />
</field_group> <!-- grid_T -->
<field_group id="grid_T_iso" domain_ref="grid_T_inner" >
<field id="20d" long_name="Depth of 20C isotherm" standard_name="depth_of_isosurface_of_sea_water_potential_temperature" unit="m" axis_ref="iax_20C" />
<field id="26d" long_name="Depth of 26C isotherm" standard_name="depth_of_isosurface_of_sea_water_potential_temperature" unit="m" axis_ref="iax_26C" />
<field id="28d" long_name="Depth of 28C isotherm" standard_name="depth_of_isosurface_of_sea_water_potential_temperature" unit="m" axis_ref="iax_28C" />
</field_group>
<!-- Tides -->
<field_group id="Tides_T" grid_ref="grid_T_2D" operation="once" >
......@@ -330,7 +418,12 @@ that are available in the tidal-forcing implementation (see
<field id="zustar" long_name="friction velocity" unit="m/s" />
<field id="zwstrl" long_name="langmuir velocity scale" unit="m/s" />
<field id="zvstr" long_name="mixed velocity scale" unit="m/s" />
<field id="zla" long_name="langmuir number" unit="m/s" />
<field id="zla" long_name="langmuir number" unit="#" />
<field id="wmp" long_name="wave mean period" unit="s" />
<field id="hsw" long_name="significant wave height" unit="m" />
<field id="wmp_NP" long_name="wave mean period from NP spectrum" unit="s" />
<field id="hsw_NP" long_name="significant wave height from NP spectrum" unit="m" />
<field id="wndm" long_name="U_10" unit="m/s" />
<field id="wind_wave_abs_power" long_name="\rho |U_s| x u*^2" unit="mW" />
<field id="wind_wave_power" long_name="U_s \dot tau" unit="mW" />
<field id="wind_power" long_name="\rho u*^3" unit="mW" />
......@@ -339,8 +432,6 @@ that are available in the tidal-forcing implementation (see
<field id="zwbav" long_name="av turb flux of buoyancy in ml" unit="m^2/s^3" />
<field id="zt_ml" long_name="av T in ml" unit="deg" />
<field id="zhol" long_name="Hoenekker number" unit="#" />
<field id="zws_ent" long_name="entrainment turb flux of S" unit="10^-3 m/s" />
<field id="zwth_ent" long_name="entrainment turb flux of T" unit="deg m/s" />
<field id="zwb_ent" long_name="entrainment turb flux of buoyancy" unit="m^2/s^-3" />
<field id="zdt_bl" long_name="temperature jump at base of BL" unit="deg" />
......@@ -359,6 +450,15 @@ that are available in the tidal-forcing implementation (see
<field id="zsc_uw_2_f" long_name="2nd zsc u-momentum flux on T after Coriolis" unit="m^2/s^2" />
<field id="zsc_vw_2_f" long_name="2nd zsc v-momentum flux on T after Coriolis" unit="m^2/s^2" />
<!-- Conditional OSMOSIS diagnostics (selectable via ln_osm_cdiags) -->
<field id="cmask" long_name="Convection mask" unit="" />
<field id="wb_lang" long_name="Entrainment buoyancy flux from Langmuire where convective" unit="m^2/s^3" />
<field id="wb_ustr" long_name="Entrainment buoyancy flux from surface shear where convective" unit="m^2/s^3" />
<field id="wb_conv" long_name="Convective entrainment buoyancy flux" unit="m^2/s^3" />
<field id="wb_shr" long_name="Entrainment buoyancy flux from basal shear where convective" unit="m^2/s^3" />
<field id="shr_prd" long_name="Shear production where convective" unit="m^2/s^3" />
<field id="h_bl" long_name="hbl where convective" unit="m" />
<!-- FK_OSM OSMOSIS diagnostics (require also ln_osm_mle=.true.-->
<field id="hmle" long_name="OBL FK-layer thickness" unit="m" />
<field id="mld_prof" long_name="FK-layer depth index" unit="#" />
......@@ -410,61 +510,63 @@ that are available in the tidal-forcing implementation (see
<!-- SBC -->
<field_group id="SBC" > <!-- time step automaticaly defined based on nn_fsbc -->
<field_group id="SBC_2D" grid_ref="grid_T_2D" >
<field_group id="SBC_2D" grid_ref="grid_T_2D_inner" >
<field id="empmr" long_name="Net Upward Water Flux" standard_name="water_flux_out_of_sea_ice_and_sea_water" unit="kg/m2/s" />
<field id="empbmr" long_name="Net Upward Water Flux at pre. tstep" standard_name="water_flux_out_of_sea_ice_and_sea_water" unit="kg/m2/s" />
<field id="emp_oce" long_name="Evap minus Precip over ocean" standard_name="evap_minus_precip_over_sea_water" unit="kg/m2/s" />
<field id="emp_ice" long_name="Evap minus Precip over ice" standard_name="evap_minus_precip_over_sea_ice" unit="kg/m2/s" />
<field id="saltflx" long_name="Downward salt flux" unit="g/m2/s" />
<field id="fmmflx" long_name="Water flux due to freezing/melting" unit="kg/m2/s" />
<field id="fwfice" long_name="Ice-Ocean Freshwater Flux (>0 to the ocean)" 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="runoffs" long_name="River Runoffs" standard_name="water_flux_into_sea_water_from_rivers" unit="kg/m2/s" grid_ref="grid_T_2D" />
<field id="precip" long_name="Total precipitation" standard_name="precipitation_flux" unit="kg/m2/s" />
<field id="wclosea" long_name="closed sea empmr correction" standard_name="closea_empmr" unit="kg/m2/s" />
<field id="qt" long_name="Net Downward Heat Flux" standard_name="surface_downward_heat_flux_in_sea_water" unit="W/m2" />
<field id="qns" long_name="non solar Downward Heat Flux" unit="W/m2" />
<field id="qsr" long_name="Shortwave Radiation" standard_name="net_downward_shortwave_flux_at_sea_water_surface" unit="W/m2" />
<field id="qsr3d" long_name="Shortwave Radiation 3D distribution" standard_name="downwelling_shortwave_flux_in_sea_water" unit="W/m2" grid_ref="grid_T_3D" />
<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="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="qrp" long_name="Surface Heat Flux: Damping" standard_name="heat_flux_into_sea_water_due_to_newtonian_relaxation" unit="W/m2" />
<field id="qclosea" long_name="closed sea heat content flux" standard_name="closea_heat_content_downward_flux" unit="W/m2" />
<field id="erp" long_name="Surface Water Flux: Damping" standard_name="water_flux_out_of_sea_water_due_to_newtonian_relaxation" unit="kg/m2/s" />
<field id="taum" long_name="wind stress module" standard_name="magnitude_of_surface_downward_stress" unit="N/m2" />
<field id="wspd" long_name="wind speed module" standard_name="wind_speed" unit="m/s" />
<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" />
<field id="utau" long_name="Wind Stress along i-axis" standard_name="surface_downward_x_stress" unit="N/m2" grid_ref="grid_T_2D" />
<field id="vtau" long_name="Wind Stress along j-axis" standard_name="surface_downward_y_stress" unit="N/m2" grid_ref="grid_T_2D" />
<!-- * 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" />
<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" grid_ref="grid_T_2D" />
<!-- *_oce variables available with ln_blk_clio or ln_blk_core -->
<field id="rho_air" long_name="Air density at 10m above sea surface" standard_name="rho_air_10m" unit="kg/m3" />
<field id="dt_skin" long_name="SSST-SST temperature difference" standard_name="SSST-SST" unit="K" />
<field id="qlw_oce" long_name="Longwave Downward Heat Flux over open ocean" standard_name="surface_net_downward_longwave_flux" unit="W/m2" />
<field id="qsb_oce" long_name="Sensible Downward Heat Flux over open ocean" standard_name="surface_downward_sensible_heat_flux" unit="W/m2" />
<field id="qla_oce" long_name="Latent Downward Heat Flux over open ocean" standard_name="surface_downward_latent_heat_flux" unit="W/m2" />
<field id="evap_oce" long_name="Evaporation over open ocean" standard_name="evaporation" unit="kg/m2/s" />
<field id="qt_oce" long_name="total flux at ocean surface" standard_name="surface_downward_heat_flux_in_sea_water" unit="W/m2" />
<field id="qsr_oce" long_name="solar heat flux at ocean surface" standard_name="net_downward_shortwave_flux_at_sea_water_surface" unit="W/m2" />
<field id="qns_oce" long_name="non-solar heat flux at ocean surface (including E-P)" unit="W/m2" />
<field id="qemp_oce" long_name="Downward Heat Flux from E-P over open ocean" unit="W/m2" />
<field id="taum_oce" long_name="wind stress module over open ocean" standard_name="magnitude_of_surface_downward_stress" unit="N/m2" />
<field id="utau_oce" long_name="Wind Stress along i-axis over open ocean (T-points)" standard_name="surf_down_x_stress_open_oce_Tpoints" unit="N/m2" />
<field id="vtau_oce" long_name="Wind Stress along j-axis over open ocean (T-points)" standard_name="surf_down_y_stress_open_oce_Tpoints" unit="N/m2" />
<field id="rho_air" long_name="Air density at 10m above sea surface" standard_name="rho_air_10m" unit="kg/m3" />
<field id="t_skin" long_name="Skin temperature aka SSST" standard_name="skin_temperature" unit="degC" />
<field id="dt_skin" long_name="SSST-SST temperature difference" standard_name="SSST-SST" unit="K" />
<field id="qlw_oce" long_name="Longwave Downward Heat Flux over open ocean" standard_name="surface_net_downward_longwave_flux" unit="W/m2" />
<field id="qsb_oce" long_name="Sensible Downward Heat Flux over open ocean" standard_name="surface_downward_sensible_heat_flux" unit="W/m2" />
<field id="qla_oce" long_name="Latent Downward Heat Flux over open ocean" standard_name="surface_downward_latent_heat_flux" unit="W/m2" />
<field id="evap_oce" long_name="Evaporation over open ocean" standard_name="evaporation" unit="kg/m2/s"/>
<field id="qt_oce" long_name="total flux at ocean surface" standard_name="surface_downward_heat_flux_in_sea_water" unit="W/m2" />
<field id="qsr_oce" long_name="solar heat flux at ocean surface" standard_name="net_downward_shortwave_flux_at_sea_water_surface" unit="W/m2" />
<field id="qns_oce" long_name="non-solar heat flux at ocean surface (including E-P)" unit="W/m2" />
<field id="qemp_oce" long_name="Downward Heat Flux from E-P over open ocean" unit="W/m2" />
<field id="taum_oce" long_name="wind stress module over open ocean" standard_name="magnitude_of_surface_downward_stress" unit="N/m2" />
<field id="utau_oce" long_name="Wind Stress along i-axis over open ocean (T-points)" standard_name="surf_down_x_stress_open_oce_Tpoints" unit="N/m2" grid_ref="grid_T_2D" />
<field id="vtau_oce" long_name="Wind Stress along j-axis over open ocean (T-points)" standard_name="surf_down_y_stress_open_oce_Tpoints" unit="N/m2" grid_ref="grid_T_2D" />
<!-- variables computed by the bulk parameterization algorithms (ln_blk) -->
<field id="Cd_oce" long_name="Drag coefficient over open ocean" standard_name="drag_coefficient_water" unit="" />
<field id="Ce_oce" long_name="Evaporaion coefficient over open ocean" standard_name="evap_coefficient_water" unit="" />
<field id="Ch_oce" long_name="Sensible heat coefficient over open ocean" standard_name="sensible_heat_coefficient_water" unit="" />
<field id="Cd_oce" long_name="Drag coefficient over open ocean" standard_name="drag_coefficient_water" unit="" />
<field id="Ce_oce" long_name="Evaporaion coefficient over open ocean" standard_name="evap_coefficient_water" unit="" />
<field id="Ch_oce" long_name="Sensible heat coefficient over open ocean" standard_name="sensible_heat_coefficient_water" unit="" />
<field id="theta_zt" long_name="Potential air temperature at z=zt" standard_name="potential_air_temperature_at_zt" unit="degC" />
<field id="q_zt" long_name="Specific air humidity at z=zt" standard_name="specific_air_humidity_at_zt" unit="kg/kg" />
<field id="q_zt" long_name="Specific air humidity at z=zt" standard_name="specific_air_humidity_at_zt" unit="kg/kg"/>
<field id="theta_zu" long_name="Potential air temperature at z=zu" standard_name="potential_air_temperature_at_zu" unit="degC" />
<field id="q_zu" long_name="Specific air humidity at z=zu" standard_name="specific_air_humidity_at_zu" unit="kg/kg" />
<field id="ssq" long_name="Saturation specific humidity of air at z=0" standard_name="surface_air_saturation_spec_humidity" unit="kg/kg" />
<field id="wspd_blk" long_name="Bulk wind speed at z=zu" standard_name="bulk_wind_speed_at_zu" unit="m/s" />
<field id="q_zu" long_name="Specific air humidity at z=zu" standard_name="specific_air_humidity_at_zu" unit="kg/kg"/>
<field id="ssq" long_name="Saturation specific humidity of air at z=0" standard_name="surface_air_saturation_spec_humidity" unit="kg/kg"/>
<field id="wspd_blk" long_name="Bulk wind speed at z=zu" standard_name="bulk_wind_speed_at_zu" unit="m/s" />
<!-- ln_blk + key_si3 -->
<field id="Cd_ice" long_name="Drag coefficient over ice" standard_name="drag_coefficient_ice" unit="" />
<field id="Ce_ice" long_name="Evaporaion coefficient over ice" standard_name="evap_coefficient_ice" unit="" />
<field id="Ch_ice" long_name="Sensible heat coefficient over ice" standard_name="sensible_heat_coefficient_ice" unit="" />
<field id="Cd_ice" long_name="Drag coefficient over ice" standard_name="drag_coefficient_ice" unit="" />
<field id="Ce_ice" long_name="Evaporaion coefficient over ice" standard_name="evap_coefficient_ice" unit="" />
<field id="Ch_ice" long_name="Sensible heat coefficient over ice" standard_name="sensible_heat_coefficient_ice" unit="" />
<!-- available key_oasis3 -->
<field id="snow_ao_cea" long_name="Snow over ice-free ocean (cell average)" standard_name="snowfall_flux" unit="kg/m2/s" />
......@@ -510,22 +612,22 @@ that are available in the tidal-forcing implementation (see
<field id="vflx_fwb_cea" long_name="volume flux due to fwb" standard_name="volume_flux_due_to_fwb" unit="kg/m2/s" />
<!-- ice field (nn_ice=1) -->
<field id="ice_cover" long_name="Ice fraction" standard_name="sea_ice_area_fraction" unit="1" />
<field id="ice_cover" long_name="Ice fraction" standard_name="sea_ice_area_fraction" unit="1" grid_ref="grid_T_2D" />
<!-- 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" />
<field id="emp_x_sst" long_name="Concentration/Dilution term on SST" unit="kg*degC/m2/s" grid_ref="grid_T_2D" />
<field id="emp_x_sss" long_name="Concentration/Dilution term on SSS" unit="kg*1e-3/m2/s" grid_ref="grid_T_2D" />
<field id="rnf_x_sst" long_name="Runoff term on SST" unit="kg*degC/m2/s" grid_ref="grid_T_2D" />
<field id="rnf_x_sss" long_name="Runoff term on SSS" unit="kg*1e-3/m2/s" grid_ref="grid_T_2D" />
<!-- 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="-" />
<field id="sst_m" unit="degC" grid_ref="grid_T_2D" />
<field id="sss_m" unit="psu" grid_ref="grid_T_2D" />
<field id="ssu_m" unit="m/s" grid_ref="grid_T_2D" />
<field id="ssv_m" unit="m/s" grid_ref="grid_T_2D" />
<field id="ssh_m" unit="m" grid_ref="grid_T_2D" />
<field id="e3t_m" unit="m" grid_ref="grid_T_2D" />
<field id="frq_m" unit="-" grid_ref="grid_T_2D" />
</field_group>
......@@ -580,46 +682,43 @@ that are available in the tidal-forcing implementation (see
<field_group id="grid_U" grid_ref="grid_U_2D">
<field id="hu" long_name="water column height at U point" standard_name="water_column_height_U" unit="m" />
<field id="e2u" long_name="U-cell width in meridional direction" standard_name="cell_width" unit="m" />
<field id="e3u" long_name="U-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_U_3D" />
<field id="e3u_0" long_name="Initial U-cell thickness" standard_name="ref_cell_thickness" unit="m" grid_ref="grid_U_3D"/>
<field id="utau" long_name="Wind Stress along i-axis" standard_name="surface_downward_x_stress" unit="N/m2" />
<field id="uoce" long_name="ocean current along i-axis" standard_name="sea_water_x_velocity" unit="m/s" grid_ref="grid_U_3D" />
<field id="uoce_e3u" long_name="ocean current along i-axis (thickness weighted)" unit="m/s" grid_ref="grid_U_3D" > uoce * e3u </field>
<field id="uoce_e3u_vsum" long_name="ocean current along i-axis * e3u summed on the vertical" field_ref="uoce_e3u" unit="m3/s" grid_ref="grid_U_vsum"/>
<field id="uocetr_vsum" long_name="ocean transport along i-axis summed on the vertical" field_ref="e2u" unit="m3/s"> this * uoce_e3u_vsum </field>
<field id="uocetr_vsum_op" long_name="ocean current along i-axis * e3u * e2u summed on the vertical" read_access="true" freq_op="1mo" field_ref="e2u" unit="m3/s"> @uocetr_vsum </field>
<field id="e2u" long_name="U-cell width in meridional direction" standard_name="cell_width" unit="m" />
<field id="e3u" long_name="U-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_U_3D_inner" />
<field id="e3u_0" long_name="Initial U-cell thickness" standard_name="ref_cell_thickness" unit="m" grid_ref="grid_U_3D_inner" />
<field id="uoce" long_name="ocean current along i-axis" standard_name="sea_water_x_velocity" unit="m/s" grid_ref="grid_U_3D" />
<field id="uoce_e3u" long_name="ocean current along i-axis (thickness weighted)" unit="m/s" grid_ref="grid_U_3D" > uoce * e3u </field>
<field id="uoce_e3u_vsum" long_name="ocean current along i-axis * e3u summed on the vertical" field_ref="uoce_e3u" unit="m3/s" grid_ref="grid_U_vsum" />
<field id="uocetr_vsum" long_name="ocean transport along i-axis summed on the vertical" field_ref="e2u" unit="m3/s" > this * uoce_e3u_vsum </field>
<field id="uocetr_vsum_op" long_name="ocean current along i-axis * e3u * e2u summed on the vertical" read_access="true" freq_op="1mo" field_ref="e2u" unit="m3/s" > @uocetr_vsum </field>
<field id="uocetr_vsum_cumul" long_name="ocean current along i-axis * e3u * e2u cumulated from southwest point" freq_offset="_reset_" operation="instant" freq_op="1mo" unit="m3/s" />
<field id="msftbarot" long_name="ocean_barotropic_mass_streamfunction" unit="kg s-1" > uocetr_vsum_cumul * $rho0 </field>
<field id="msftbarot" long_name="ocean_barotropic_mass_streamfunction" unit="kg s-1" > uocetr_vsum_cumul * $rho0 </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" grid_ref="grid_U_2D_inner" />
<field id="ubar" long_name="ocean barotropic current along i-axis" unit="m/s" />
<field id="uocetr_eff" long_name="Effective ocean transport along i-axis" standard_name="ocean_volume_x_transport" unit="m3/s" grid_ref="grid_U_3D" />
<field id="uocet" long_name="ocean transport along i-axis times temperature (CRS)" unit="degC*m/s" grid_ref="grid_U_3D" />
<field id="uoces" long_name="ocean transport along i-axis times salinity (CRS)" unit="1e-3*m/s" grid_ref="grid_U_3D" />
<field id="ssuww" long_name="ocean surface wind work along i-axis" standard_name="surface_x_wind_work" unit="N/m*s" > utau * ssu </field>
<!-- AGRIF sponge -->
<field id="agrif_spu" long_name=" AGRIF u-sponge coefficient" unit=" " />
<!-- u-eddy diffusivity coefficients (available if ln_traldf_OFF=F) -->
<field id="ahtu_2d" long_name=" surface u-eddy diffusivity coefficient" unit="m2/s or m4/s" />
<field id="ahtu_3d" long_name=" 3D u-EIV coefficient" unit="m2/s or m4/s" grid_ref="grid_U_3D"/>
<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"/>
<field id="aeiu_2d" long_name=" surface u-EIV coefficient" unit="m2/s" />
<field id="aeiu_3d" long_name=" 3D u-EIV coefficient" unit="m2/s" grid_ref="grid_U_3D" />
<!-- variables available with MLE (ln_mle=T) -->
<field id="psiu_mle" long_name="MLE streamfunction along i-axis" unit="m3/s" grid_ref="grid_U_3D" />
<field id="psiu_mle" long_name="MLE streamfunction along i-axis" unit="m3/s" grid_ref="grid_U_3D_inner" />
<!-- uoce_eiv: available EIV (ln_ldfeiv=T and ln_ldfeiv_dia=T) -->
<field id="uoce_eiv" long_name="EIV ocean current along i-axis" standard_name="bolus_sea_water_x_velocity" unit="m/s" grid_ref="grid_U_3D" />
<field id="ueiv_masstr" long_name="EIV Ocean Mass X Transport" standard_name="bolus_ocean_mass_x_transport" unit="kg/s" grid_ref="grid_U_3D" />
<field id="ueiv_heattr" long_name="ocean bolus heat transport along i-axis" standard_name="ocean_heat_x_transport_due_to_bolus_advection" unit="W" />
<field id="ueiv_salttr" long_name="ocean bolus salt transport along i-axis" standard_name="ocean_salt_x_transport_due_to_bolus_advection" unit="Kg" />
<field id="ueiv_heattr3d" long_name="ocean bolus heat transport along i-axis" standard_name="ocean_heat_x_transport_due_to_bolus_advection" unit="W" grid_ref="grid_U_3D" />
<field id="ueiv_salttr3d" long_name="ocean bolus salt transport along i-axis" standard_name="ocean_salt_x_transport_due_to_bolus_advection" unit="kg" grid_ref="grid_U_3D" />
<!-- uoce_eiv: available EIV (ln_ldfeiv=T) -->
<field id="uoce_eiv" long_name="EIV ocean current along i-axis" standard_name="bolus_sea_water_x_velocity" unit="m/s" grid_ref="grid_U_3D_inner" />
<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_inner" />
<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" grid_ref="grid_U_2D_inner" />
<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" grid_ref="grid_U_2D_inner" />
<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_inner" />
<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_inner" />
<!-- uoce_bbl: available with ln_trabbl=T and nn_bbl_adv=1 -->
<field id="uoce_bbl" long_name="BBL ocean current along i-axis" unit="m/s" />
......@@ -630,56 +729,50 @@ that are available in the tidal-forcing implementation (see
<field id="ustokes" long_name="Stokes Drift Velocity i-axis" standard_name="StokesDrift_x_velocity" unit="m/s" grid_ref="grid_U_3D" />
<field id="ustokes_e3u" long_name="Stokes Drift Velocity i-axis (thickness weighted)" unit="m/s" grid_ref="grid_U_3D" > ustokes * e3u </field>
<!-- variable for ice shelves -->
<field id="utbl" long_name="zonal current in the Losh tbl" unit="m/s" />
<!-- variables available with diaar5 -->
<field id="u_masstr" long_name="Ocean Mass X Transport" standard_name="ocean_mass_x_transport" unit="kg/s" grid_ref="grid_U_3D" />
<field id="u_masstr_vint" long_name="vertical integral of ocean eulerian mass transport along i-axis" standard_name="vertical_integral_of_ocean_mass_x_transport" unit="kg/s" grid_ref="grid_U_2D_inner" />
<field id="u_heattr" long_name="ocean eulerian heat transport along i-axis" standard_name="ocean_heat_x_transport" unit="W" grid_ref="grid_U_2D_inner" />
<field id="u_masstr" long_name="Ocean Mass X Transport" standard_name="ocean_mass_x_transport" unit="kg/s" grid_ref="grid_U_3D_inner" />
<field id="u_masstr_vint" long_name="vertical integral of ocean eulerian mass transport along i-axis" standard_name="vertical_integral_of_ocean_mass_x_transport" unit="kg/s" grid_ref="grid_U_2D_inner" />
<field id="u_heattr" long_name="ocean eulerian heat transport along i-axis" standard_name="ocean_heat_x_transport" unit="W" grid_ref="grid_U_2D_inner" />
<field id="u_salttr" long_name="ocean eulerian salt transport along i-axis" standard_name="ocean_salt_x_transport" unit="1e-3*kg/s" grid_ref="grid_U_2D_inner" />
<field id="uadv_heattr" long_name="ocean advective heat transport along i-axis" standard_name="advectice_ocean_heat_x_transport" unit="W" />
<field id="uadv_salttr" long_name="ocean advective salt transport along i-axis" standard_name="advectice_ocean_salt_x_transport" unit="1e-3*kg/s" />
<field id="udiff_heattr" long_name="ocean diffusion heat transport along i-axis" standard_name="ocean_heat_x_transport_due_to_diffusion" unit="W" />
<field id="udiff_salttr" long_name="ocean diffusion salt transport along i-axis" standard_name="ocean_salt_x_transport_due_to_diffusion" unit="1e-3*kg/s" />
<field id="uadv_heattr" long_name="ocean advective heat transport along i-axis" standard_name="advectice_ocean_heat_x_transport" unit="W" grid_ref="grid_U_2D_inner" />
<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" grid_ref="grid_U_2D_inner" />
<field id="udiff_heattr" long_name="ocean diffusion heat transport along i-axis" standard_name="ocean_heat_x_transport_due_to_diffusion" unit="W" grid_ref="grid_U_2D_inner" />
<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" grid_ref="grid_U_2D_inner" />
</field_group>
<!-- V grid -->
<field_group id="grid_V" grid_ref="grid_V_2D">
<field id="e1v" long_name="V-cell width in longitudinal direction" standard_name="cell_width" unit="m" />
<field id="e3v" long_name="V-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_V_3D" />
<field id="e3v_0" long_name="Initial V-cell thickness" standard_name="ref_cell_thickness" unit="m" grid_ref="grid_V_3D" />
<field id="e1v" long_name="V-cell width in longitudinal direction" standard_name="cell_width" unit="m" />
<field id="e3v" long_name="V-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_V_3D_inner" />
<field id="e3v_0" long_name="Initial V-cell thickness" standard_name="ref_cell_thickness" unit="m" grid_ref="grid_V_3D_inner" />
<field id="hv" long_name="water column height at V point" standard_name="water_column_height_V" unit="m" />
<field id="vtau" long_name="Wind Stress along j-axis" standard_name="surface_downward_y_stress" unit="N/m2" />
<field id="voce" long_name="ocean current along j-axis" standard_name="sea_water_y_velocity" unit="m/s" grid_ref="grid_V_3D" />
<field id="voce_e3v" long_name="ocean current along j-axis (thickness weighted)" unit="m/s" grid_ref="grid_V_3D" > voce * e3v </field>
<field id="ssv" long_name="ocean surface current along j-axis" unit="m/s" />
<field id="sbv" long_name="ocean bottom current along j-axis" unit="m/s" grid_ref="grid_V_2D_inner" />
<field id="vbar" long_name="ocean barotropic current along j-axis" unit="m/s" />
<field id="vocetr_eff" long_name="Effective ocean transport along j-axis" standard_name="ocean_volume_y_transport" unit="m3/s" grid_ref="grid_V_3D" />
<field id="vocet" long_name="ocean transport along j-axis times temperature (CRS)" unit="degC*m/s" grid_ref="grid_V_3D" />
<field id="voces" long_name="ocean transport along j-axis times salinity (CRS)" unit="1e-3*m/s" grid_ref="grid_V_3D" />
<field id="ssvww" long_name="ocean surface wind work along j-axis" standard_name="surface_y_wind_work" unit="N/m*s" > vtau * ssv </field>
<!-- AGRIF sponge -->
<field id="agrif_spv" long_name=" AGRIF v-sponge coefficient" unit=" " />
<!-- v-eddy diffusivity coefficients (available if ln_traldf_OFF=F) -->
<field id="ahtv_2d" long_name=" surface v-eddy diffusivity coefficient" unit="m2/s or (m4/s)^1/2" />
<field id="ahtv_3d" long_name=" 3D v-eddy diffusivity coefficient" unit="m2/s or (m4/s)^1/2" grid_ref="grid_V_3D"/>
<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" />
<field id="aeiv_2d" long_name=" surface v-EIV coefficient" unit="m2/s" />
<field id="aeiv_3d" long_name=" 3D v-EIV coefficient" unit="m2/s" grid_ref="grid_V_3D" />
<!-- variables available with MLE (ln_mle=T) -->
<field id="psiv_mle" long_name="MLE streamfunction along j-axis" unit="m3/s" grid_ref="grid_V_3D" />
<field id="psiv_mle" long_name="MLE streamfunction along j-axis" unit="m3/s" grid_ref="grid_V_3D_inner" />
<!-- voce_eiv: available EIV (ln_ldfeiv=T and ln_ldfeiv_dia=T) -->
<field id="voce_eiv" long_name="EIV ocean current along j-axis" standard_name="bolus_sea_water_y_velocity" unit="m/s" grid_ref="grid_V_3D" />
<field id="veiv_masstr" long_name="EIV Ocean Mass Y Transport" standard_name="bolus_ocean_mass_y_transport" unit="kg/s" grid_ref="grid_V_3D" />
<field id="veiv_heattr" long_name="ocean bolus heat transport along j-axis" standard_name="ocean_heat_y_transport_due_to_bolus_advection" unit="W" />
<field id="veiv_salttr" long_name="ocean bolus salt transport along j-axis" standard_name="ocean_salt_x_transport_due_to_bolus_advection" unit="Kg" />
<field id="veiv_heattr3d" long_name="ocean bolus heat transport along j-axis" standard_name="ocean_heat_y_transport_due_to_bolus_advection" unit="W" grid_ref="grid_V_3D" />
<field id="veiv_salttr3d" long_name="ocean bolus salt transport along j-axis" standard_name="ocean_salt_y_transport_due_to_bolus_advection" unit="kg" grid_ref="grid_V_3D" />
<!-- voce_eiv: available EIV (ln_ldfeiv=T) -->
<field id="voce_eiv" long_name="EIV ocean current along j-axis" standard_name="bolus_sea_water_y_velocity" unit="m/s" grid_ref="grid_V_3D_inner" />
<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_inner" />
<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" grid_ref="grid_V_2D_inner" />
<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" grid_ref="grid_V_2D_inner" />
<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_inner" />
<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_inner" />
<!-- voce_bbl: available with ln_trabbl=T and nn_bbl_adv=1 -->
......@@ -691,91 +784,97 @@ that are available in the tidal-forcing implementation (see
<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_masstr" long_name="ocean eulerian mass transport along j-axis" standard_name="ocean_mass_y_transport" unit="kg/s" grid_ref="grid_V_3D_inner" />
<field id="v_heattr" long_name="ocean eulerian heat transport along j-axis" standard_name="ocean_heat_y_transport" unit="W" grid_ref="grid_V_2D_inner" />
<field id="v_salttr" long_name="ocean eulerian salt transport along i-axis" standard_name="ocean_salt_y_transport" unit="1e-3*kg/s" grid_ref="grid_V_2D_inner" />
<field id="vadv_heattr" long_name="ocean advective heat transport along j-axis" standard_name="advectice_ocean_heat_y_transport" unit="W" />
<field id="vadv_salttr" long_name="ocean advective salt transport along j-axis" standard_name="advectice_ocean_salt_y_transport" unit="1e-3*kg/s" />
<field id="vdiff_heattr" long_name="ocean diffusion heat transport along j-axis" standard_name="ocean_heat_y_transport_due_to_diffusion" unit="W" />
<field id="vdiff_salttr" long_name="ocean diffusion salt transport along j-axis" standard_name="ocean_salt_y_transport_due_to_diffusion" unit="1e-3*kg/s" />
<field id="vadv_heattr" long_name="ocean advective heat transport along j-axis" standard_name="advectice_ocean_heat_y_transport" unit="W" grid_ref="grid_V_2D_inner" />
<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" grid_ref="grid_V_2D_inner" />
<field id="vdiff_heattr" long_name="ocean diffusion heat transport along j-axis" standard_name="ocean_heat_y_transport_due_to_diffusion" unit="W" grid_ref="grid_V_2D_inner" />
<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" grid_ref="grid_V_2D_inner" />
</field_group>
<!-- W grid -->
<field_group id="grid_W" grid_ref="grid_W_3D">
<field id="e3w" long_name="W-cell thickness" standard_name="cell_thickness" unit="m" />
<field id="woce" long_name="ocean vertical velocity" standard_name="upward_sea_water_velocity" unit="m/s" />
<field id="e3w" long_name="W-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_W_3D_inner" />
<field id="woce" long_name="ocean vertical velocity" standard_name="upward_sea_water_velocity" unit="m/s" grid_ref="grid_W_3D_inner"/>
<field id="woce_e3w" long_name="ocean vertical velocity * e3w" unit="m2/s" > woce * e3w </field>
<field id="wocetr_eff" long_name="effective ocean vertical transport" unit="m3/s" />
<!-- woce_eiv: available with EIV (ln_ldfeiv=T and ln_ldfeiv_dia=T) -->
<field id="woce_eiv" long_name="EIV ocean vertical velocity" standard_name="bolus_upward_sea_water_velocity" unit="m/s" />
<field id="weiv_masstr" long_name="EIV Upward Ocean Mass Transport" standard_name="bolus_upward_ocean_mass_transport" unit="kg/s" />
<field id="weiv_heattr3d" long_name="ocean bolus heat transport" standard_name="ocean_heat_z_transport_due_to_bolus_advection" unit="W" />
<field id="weiv_salttr3d" long_name="ocean bolus salt transport" standard_name="ocean_salt_z_transport_due_to_bolus_advection" unit="kg" />
<!-- variables available with WAVE (ln_wave=T) -->
<field id="wstokes" long_name="Stokes Drift vertical velocity" standard_name="upward_StokesDrift_velocity" unit="m/s" />
<!-- woce_eiv: available with EIV (ln_ldfeiv=T) -->
<field id="woce_eiv" long_name="EIV ocean vertical velocity" standard_name="bolus_upward_sea_water_velocity" unit="m/s" grid_ref="grid_W_3D_inner" />
<field id="weiv_masstr" long_name="EIV Upward Ocean Mass Transport" standard_name="bolus_upward_ocean_mass_transport" unit="kg/s" grid_ref="grid_W_3D_inner" />
<!-- currently not coded
<field id="weiv_heattr3d" long_name="ocean bolus heat transport" standard_name="ocean_heat_z_transport_due_to_bolus_advection" unit="W" grid_ref="grid_W_3D_inner" />
<field id="weiv_salttr3d" long_name="ocean bolus salt transport" standard_name="ocean_salt_z_transport_due_to_bolus_advection" unit="kg" grid_ref="grid_W_3D_inner" />
-->
<field id="avt" long_name="vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" />
<field id="avt_e3w" long_name="vertical heat diffusivity * e3w" unit="m3/s" > avt * e3w </field>
<field id="logavt" long_name="logarithm of vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" />
<field id="avm" long_name="vertical eddy viscosity" standard_name="ocean_vertical_momentum_diffusivity" unit="m2/s" />
<field id="avm_e3w" long_name="vertical eddy viscosity * e3w" unit="m3/s" > avm * e3w </field>
<!-- avt, avm -->
<field id="avt" long_name="vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" grid_ref="grid_W_3D_inner" />
<field id="avt_e3w" long_name="vertical heat diffusivity * e3w" unit="m3/s" > avt * e3w </field>
<field id="logavt" long_name="logarithm of vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" grid_ref="grid_W_3D_inner" />
<field id="avm" long_name="vertical eddy viscosity" standard_name="ocean_vertical_momentum_diffusivity" unit="m2/s" />
<field id="avm_e3w" long_name="vertical eddy viscosity * e3w" unit="m3/s" > avm * e3w </field>
<!-- avs: /= avt with ln_zdfddm=T -->
<field id="avs" long_name="salt vertical eddy diffusivity" standard_name="ocean_vertical_salt_diffusivity" unit="m2/s" />
<field id="avs_e3w" long_name="vertical salt diffusivity * e3w" unit="m3/s" > avs * e3w </field>
<field id="logavs" long_name="logarithm of salt vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" />
<field id="avs" long_name="salt vertical eddy diffusivity" standard_name="ocean_vertical_salt_diffusivity" unit="m2/s" grid_ref="grid_W_3D_inner" />
<field id="avs_e3w" long_name="vertical salt diffusivity * e3w" unit="m3/s" > avs * e3w </field>
<field id="logavs" long_name="logarithm of salt vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" grid_ref="grid_W_3D_inner" />
<!-- avt_evd and avm_evd: available with ln_zdfevd -->
<field id="avt_evd" long_name="convective enhancement of vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_convection" unit="m2/s" />
<field id="avt_evd_e3w" long_name="convective enhancement to vertical diffusivity * e3w " unit="m3/s" > avt_evd * e3w </field>
<field id="avm_evd" long_name="convective enhancement of vertical viscosity" standard_name="ocean_vertical_momentum_diffusivity_due_to_convection" unit="m2/s" />
<field id="avt_evd" long_name="convective enhancement of vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_convection" unit="m2/s" grid_ref="grid_W_3D_inner" />
<field id="avt_evd_e3w" long_name="convective enhancement to vertical diffusivity * e3w " unit="m3/s" > avt_evd * e3w </field>
<field id="avm_evd" long_name="convective enhancement of vertical viscosity" standard_name="ocean_vertical_momentum_diffusivity_due_to_convection" unit="m2/s" grid_ref="grid_W_3D_inner" />
<!-- mf_app and mf_wp: available with ln_zdfmfc -->
<field id="mf_app" long_name="convective area" standard_name="mf_convective_area" unit="%" grid_ref="grid_W_3D" />
<field id="mf_wp" long_name="convective velocity" standard_name="mf_convective_velo" unit="m/s" grid_ref="grid_W_3D" />
<field id="mf_app" long_name="convective area" standard_name="mf_convective_area" unit="%" grid_ref="grid_W_3D_inner" />
<field id="mf_wp" long_name="convective velocity" standard_name="mf_convective_velo" unit="m/s" grid_ref="grid_W_3D_inner" />
<!-- avt_tide: available with ln_zdfiwm=T -->
<field id="av_ratio" long_name="S over T diffusivity ratio" standard_name="salinity_over_temperature_diffusivity_ratio" unit="1" />
<field id="av_wave" long_name="internal wave-induced vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_internal_waves" unit="m2/s" />
<field id="bflx_iwm" long_name="internal wave-induced buoyancy flux" standard_name="buoyancy_flux_due_to_internal_waves" unit="W/kg" />
<field id="pcmap_iwm" long_name="power consumed by wave-driven mixing" standard_name="vertically_integrated_power_consumption_by_wave_driven_mixing" unit="W/m2" grid_ref="grid_W_2D" />
<field id="emix_iwm" long_name="power density available for mixing" standard_name="power_available_for_mixing_from_breaking_internal_waves" unit="W/kg" />
<!-- variables available with WAVE (ln_wave=T) -->
<field id="wstokes" long_name="Stokes Drift vertical velocity" standard_name="upward_StokesDrift_velocity" unit="m/s" />
<field id="av_ratio" long_name="S over T diffusivity ratio" standard_name="salinity_over_temperature_diffusivity_ratio" unit="1" grid_ref="grid_W_3D_inner" />
<field id="av_wave" long_name="internal wave-induced vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_internal_waves" unit="m2/s" grid_ref="grid_W_3D_inner" />
<field id="bflx_iwm" long_name="internal wave-induced buoyancy flux" standard_name="buoyancy_flux_due_to_internal_waves" unit="W/kg" grid_ref="grid_W_3D_inner" />
<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_inner" />
<field id="emix_iwm" long_name="power density available for mixing" standard_name="power_available_for_mixing_from_breaking_internal_waves" unit="W/kg" grid_ref="grid_W_3D_inner" />
<!-- variables available with diaar5 -->
<field id="w_masstr" long_name="vertical mass transport" standard_name="upward_ocean_mass_transport" unit="kg/s" />
<field id="w_masstr2" long_name="square of vertical mass transport" standard_name="square_of_upward_ocean_mass_transport" unit="kg2/s2" />
<field id="w_masstr" long_name="vertical mass transport" standard_name="upward_ocean_mass_transport" unit="kg/s" grid_ref="grid_W_3D_inner" />
<field id="w_masstr2" long_name="square of vertical mass transport" standard_name="square_of_upward_ocean_mass_transport" unit="kg2/s2" grid_ref="grid_W_3D_inner" />
<!-- EOS -->
<field id="bn2" long_name="squared Brunt-Vaisala frequency" unit="s-2" />
<!-- GEOMETRIC fields (requires nn_aei_ijk_t = 32) -->
<field id="aeiv_geom" long_name="3D w-EIV coefficient from GEOMETRIC param." unit="m2/s" grid_ref="grid_W_3D_inner"/>
<field id="rossby_rad" long_name="internal Rossby deformation radius" unit="m" grid_ref="grid_W_2D_inner"/>
<field id="bn2" long_name="squared Brunt-Vaisala frequency" unit="s-1" />
<field id="c1_vert" long_name="1st baroclinic mode phase speed" unit="m/s" grid_ref="grid_W_2D_inner"/>
<field id="c_ros" long_name="long Rossby phase speed" unit="m/s" grid_ref="grid_W_2D_inner"/>
<!-- dissipation diagnostics (note: ediss_k is only available with tke scheme) -->
<field id="avt_k" long_name="vertical eddy diffusivity from closure schemes" standard_name="ocean_vertical_eddy_diffusivity" unit="m2/s" />
<field id="avm_k" long_name="vertical eddy viscosity from closure schemes" standard_name="ocean_vertical_eddy_viscosity" unit="m2/s" />
<field id="ediss_k" long_name="Kolmogorov energy dissipation (tke scheme)" standard_name="Kolmogorov_energy_dissipation" unit="W/kg" />
<field id="eshear_k" long_name="energy source from vertical shear" standard_name="energy_source_from_shear" unit="W/kg" />
<field id="estrat_k" long_name="energy sink from stratification" standard_name="energy_sink_from_stratification" unit="W/kg" />
<field id="avt_k" long_name="vertical eddy diffusivity from closure schemes" standard_name="ocean_vertical_eddy_diffusivity" unit="m2/s" grid_ref="grid_W_3D_inner" />
<field id="avm_k" long_name="vertical eddy viscosity from closure schemes" standard_name="ocean_vertical_eddy_viscosity" unit="m2/s" grid_ref="grid_W_3D_inner" />
<field id="ediss_k" long_name="Kolmogorov energy dissipation (tke scheme)" standard_name="Kolmogorov_energy_dissipation" unit="W/kg" grid_ref="grid_W_3D_inner" />
<field id="eshear_k" long_name="energy source from vertical shear" standard_name="energy_source_from_shear" unit="W/kg" grid_ref="grid_W_3D_inner" />
<field id="estrat_k" long_name="energy sink from stratification" standard_name="energy_sink_from_stratification" unit="W/kg" grid_ref="grid_W_3D_inner" />
</field_group>
<!-- F grid -->
<field_group id="grid_F" grid_ref="grid_F_2D">
<field id="e3f" long_name="F-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_F_3D" />
<field id="e3f_0" long_name="F-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_F_3D" />
<field id="hf" long_name="water column height at F point" standard_name="water_column_height_F" unit="m" />
<field id="ssKEf" long_name="surface kinetic energy at F point" standard_name="specific_kinetic_energy_of_sea_water" unit="m2/s2" grid_ref="grid_F_2D_inner" />
<field id="ssrelvor" long_name="surface relative vorticity" standard_name="relative_vorticity" unit="1/s" grid_ref="grid_F_2D_inner" />
<field id="ssplavor" long_name="surface planetary vorticity" standard_name="planetary_vorticity" unit="1/s" />
<field id="ssrelpotvor" long_name="surface relative potential vorticity" standard_name="relpot_vorticity" unit="1/m.s" grid_ref="grid_F_2D_inner" />
<field id="ssabspotvor" long_name="surface absolute potential vorticity" standard_name="abspot_vorticity" unit="1/m.s" grid_ref="grid_F_2D_inner" />
<field id="ssEns" long_name="surface enstrophy" standard_name="enstrophy" unit="1/m2.s2" grid_ref="grid_F_2D_inner" />
<field id="e3f" long_name="F-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_F_3D_inner" />
<field id="e3f_0" long_name="F-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_F_3D_inner" />
<field id="hf" long_name="water column height at F point" standard_name="water_column_height_F" unit="m" />
<field id="ssKEf" long_name="surface kinetic energy at F point" standard_name="specific_kinetic_energy_of_sea_water" unit="m2/s2" grid_ref="grid_F_2D_inner" />
<field id="ssrelvor" long_name="surface relative vorticity" standard_name="relative_vorticity" unit="1/s" grid_ref="grid_F_2D_inner" />
<field id="ssplavor" long_name="surface planetary vorticity" standard_name="planetary_vorticity" unit="1/s" grid_ref="grid_F_2D_inner" />
<field id="ssrelpotvor" long_name="surface relative potential vorticity" standard_name="relpot_vorticity" unit="1/m.s" grid_ref="grid_F_2D_inner" />
<field id="ssabspotvor" long_name="surface absolute potential vorticity" standard_name="abspot_vorticity" unit="1/m.s" grid_ref="grid_F_2D_inner" />
<field id="ssEns" long_name="surface enstrophy" standard_name="enstrophy" unit="1/m2.s2" grid_ref="grid_F_2D_inner" />
</field_group>
<!-- AGRIF sponge -->
......@@ -790,6 +889,7 @@ that are available in the tidal-forcing implementation (see
<field id="sshtot" long_name="global mean ssh" standard_name="global_average_sea_level_change" unit="m" />
<field id="sshsteric" long_name="global mean ssh steric" standard_name="global_average_steric_sea_level_change" unit="m" />
<field id="sshthster" long_name="global mean ssh thermosteric" standard_name="global_average_thermosteric_sea_level_change" unit="m" />
<field id="sshice" long_name="global mean ssh equivalent for ice and snow" standard_name="global_average_ice_snow_sea_level_change" unit="m" />
<field id="masstot" long_name="global total mass" standard_name="sea_water_mass" unit="kg" />
<field id="temptot" long_name="global mean temperature" standard_name="sea_water_potential_temperature" unit="degC" />
<field id="saltot" long_name="global mean salinity" standard_name="sea_water_salinity" unit="1e-3" />
......@@ -797,7 +897,7 @@ that are available in the tidal-forcing implementation (see
<!-- 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 -->
<!-- global diags -->
<field id="bgtemper" long_name="drift in global mean temperature wrt timestep 1" standard_name="change_over_time_in_sea_water_potential_temperature" unit="degC" />
<field id="bgsaline" long_name="drift in global mean salinity wrt timestep 1" standard_name="change_over_time_in_sea_water_practical_salinity" unit="1e-3" />
<field id="bgheatco" long_name="drift in global mean heat content wrt timestep 1" unit="1.e20J" />
......@@ -816,34 +916,27 @@ that are available in the tidal-forcing implementation (see
<!-- transects -->
<field_group id="oce_straits">
<field id="uoce_e3u_ave" long_name="Monthly average of u*e3u" field_ref="uoce_e3u" freq_op="1mo" freq_offset="_reset_" > @uoce_e3u </field>
<field id="uoce_e3u_ave_vsum" long_name="Vertical sum of u*e3u" field_ref="uoce_e3u_ave" grid_ref="grid_U_vsum" />
<field id="uoce_e3u_ave" long_name="Monthly average of u*e3u" field_ref="uoce_e3u" freq_op="1mo" freq_offset="_reset_" > @uoce_e3u </field>
<field id="uoce_e3u_ave_vsum" long_name="Vertical sum of u*e3u" field_ref="uoce_e3u_ave" grid_ref="grid_U_vsum" />
<field id="uocetr_vsum_section" long_name="Total 2D transport in i-direction" field_ref="uoce_e3u_ave_vsum" grid_ref="grid_U_scalar" detect_missing_value="true"> this * e2u </field>
<field id="uocetr_strait" long_name="Total transport across lines in i-direction" field_ref="uocetr_vsum_section" grid_ref="grid_U_4strait" />
<field id="u_masstr_strait" long_name="Sea water transport across line in i-direction" field_ref="uocetr_strait" grid_ref="grid_U_4strait_hsum" unit="kg/s"> this * maskMFO_u * $rho0 </field>
<field id="voce_e3v_ave" long_name="Monthly average of v*e3v" field_ref="voce_e3v" freq_op="1mo" freq_offset="_reset_" > @voce_e3v </field>
<field id="voce_e3v_ave_vsum" long_name="Vertical sum of v*e3v" field_ref="voce_e3v_ave" grid_ref="grid_V_vsum" />
<field id="voce_e3v_ave" long_name="Monthly average of v*e3v" field_ref="voce_e3v" freq_op="1mo" freq_offset="_reset_" > @voce_e3v </field>
<field id="voce_e3v_ave_vsum" long_name="Vertical sum of v*e3v" field_ref="voce_e3v_ave" grid_ref="grid_V_vsum" />
<field id="vocetr_vsum_section" long_name="Total 2D transport of in j-direction" field_ref="voce_e3v_ave_vsum" grid_ref="grid_V_scalar" detect_missing_value="true"> this * e1v </field>
<field id="vocetr_strait" long_name="Total transport across lines in j-direction" field_ref="vocetr_vsum_section" grid_ref="grid_V_4strait" />
<field id="vocetr_strait" long_name="Total transport across lines in j-direction" field_ref="vocetr_vsum_section" grid_ref="grid_V_4strait" />
<field id="v_masstr_strait" long_name="Sea water transport across line in j-direction" field_ref="vocetr_strait" grid_ref="grid_V_4strait_hsum" unit="kg/s"> this * maskMFO_v * $rho0 </field>
<field id="masstr_strait" long_name="Sea water transport across line" grid_ref="grid_4strait" > u_masstr_strait + v_masstr_strait </field>
</field_group>
<!-- variables available with ln_floats -->
<field_group id="floatvar" grid_ref="grid_T_nfloat" operation="instant" >
<field id="traj_lon" long_name="floats longitude" unit="degrees_east" />
<field id="traj_lat" long_name="floats latitude" unit="degrees_north" />
<field id="traj_dep" long_name="floats depth" unit="m" />
<field id="traj_temp" long_name="floats temperature" standard_name="sea_water_potential_temperature" unit="degC" />
<field id="traj_salt" long_name="floats salinity" standard_name="sea_water_practical_salinity" unit="1e-3" />
<field id="traj_dens" long_name="floats in-situ density" standard_name="sea_water_density" unit="kg/m3" />
<field id="traj_group" long_name="floats group" unit="1" />
</field_group>
<!-- variables available with iceberg trajectories -->
<field_group id="icbvar_bas" grid_ref="grid_T_icbbas" >
<field id="berg_melt_basins" long_name="icb melt rate of icebergs per basin" unit="kg/m2/s" />
<field id="berg_sumthic_basins" long_name="sum of icb thickness per cell and per basin" unit="m" operation="accumulate" prec="8" />
<field id="berg_numb_basins" long_name="sum of icb number per cell and per basin" unit="" operation="accumulate" prec="8" />
</field_group>
<field_group id="icbvar" domain_ref="grid_T" >
<field id="berg_melt" long_name="icb melt rate of icebergs" unit="kg/m2/s" />
......@@ -901,20 +994,9 @@ that are available in the tidal-forcing implementation (see
-->
<!-- 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_group id="trendT_even" grid_ref="grid_T_3D_inner">
<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" />
......@@ -925,8 +1007,8 @@ that are available in the tidal-forcing implementation (see
<field id="strd_ad" long_name="salinity -trend: advection" standard_name="tendency_of_sea_water_salinity_due_to_advection" unit="1e-3/s" > sqrt( strd_xad^2 + strd_yad^2 + strd_zad^2 ) </field>
<field id="ttrd_totad" long_name="temperature-trend: total advection" standard_name="tendency_of_sea_water_salinity_due_to_advection" unit="degC/s" />
<field id="strd_totad" long_name="salinity -trend: total advection" standard_name="tendency_of_sea_water_salinity_due_to_advection" unit="1e-3/s" />
<field id="ttrd_sad" long_name="temperature-trend: surface adv. (linssh true)" unit="degC/s" grid_ref="grid_T_2D" />
<field id="strd_sad" long_name="salinity -trend: surface adv. (linssh true)" unit="1e-3/s" grid_ref="grid_T_2D" />
<field id="ttrd_sad" long_name="temperature-trend: surface adv. (linssh true)" unit="degC/s" grid_ref="grid_T_2D_inner" />
<field id="strd_sad" long_name="salinity -trend: surface adv. (linssh true)" unit="1e-3/s" grid_ref="grid_T_2D_inner" />
<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" />
......@@ -952,8 +1034,8 @@ that are available in the tidal-forcing implementation (see
<field id="strd_bbl" long_name="salinity -trend: bottom boundary layer" unit="1e-3/s" />
<field id="ttrd_npc" long_name="temperature-trend: non-penetrative conv." unit="degC/s" />
<field id="strd_npc" long_name="salinity -trend: non-penetrative conv." unit="1e-3/s" />
<field id="ttrd_qns" long_name="temperature-trend: non-solar flux + runoff" unit="degC/s" grid_ref="grid_T_2D" />
<field id="strd_cdt" long_name="salinity -trend: C/D term + runoff" unit="degC/s" grid_ref="grid_T_2D" />
<field id="ttrd_qns" long_name="temperature-trend: non-solar flux + runoff" unit="degC/s" grid_ref="grid_T_2D_inner" />
<field id="strd_cdt" long_name="salinity -trend: C/D term + runoff" unit="degC/s" grid_ref="grid_T_2D_inner" />
<field id="ttrd_qsr" long_name="temperature-trend: solar penetr. heating" unit="degC/s" />
<field id="ttrd_bbc" long_name="temperature-trend: geothermal heating" unit="degC/s" />
......@@ -1008,17 +1090,15 @@ that are available in the tidal-forcing implementation (see
<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_qns_li" long_name="layer integrated heat-trend: non-solar flux + runoff" unit="W/m^2" grid_ref="grid_T_2D_inner"> 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_inner"> 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_group id="trendT" grid_ref="grid_T_3D_inner">
<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: -->
......@@ -1034,10 +1114,10 @@ that are available in the tidal-forcing implementation (see
<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="ssh_flt" long_name="filtered contribution to ssh (dynspg_flt)" unit="m" grid_ref="grid_T_2D_inner" />
<field id="w0" long_name="surface vertical velocity" unit="m/s" grid_ref="grid_T_2D_inner" />
<field id="pw0_exp" long_name="surface pressure flux due to ssh" unit="W/s^2" grid_ref="grid_T_2D_inner" />
<field id="pw0_flt" long_name="surface pressure flux due to filtered ssh" unit="W/s^2" grid_ref="grid_T_2D_inner" />
<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" />
......@@ -1045,24 +1125,17 @@ that are available in the tidal-forcing implementation (see
<field id="ketrd_udx" long_name="ke-trend: U.dx[U]" unit="W/s^3" />
<field id="ketrd_ldf" long_name="ke-trend: lateral diffusion" unit="W/s^3" />
<field id="ketrd_zdf" long_name="ke-trend: vertical diffusion" unit="W/s^3" />
<field id="ketrd_tau" long_name="ke-trend: wind stress " unit="W/s^3" grid_ref="grid_T_2D" />
<field id="ketrd_tau" long_name="ke-trend: wind stress " unit="W/s^3" grid_ref="grid_T_2D_inner" />
<field id="ketrd_bfr" long_name="ke-trend: bottom friction (explicit)" unit="W/s^3" />
<field id="ketrd_bfri" long_name="ke-trend: bottom friction (implicit)" unit="W/s^3" />
<field id="ketrd_atf" long_name="ke-trend: asselin time filter trend" unit="W/s^3" />
<field id="ketrd_convP2K" long_name="ke-trend: conversion (potential to kinetic)" unit="W/s^3" />
<field id="KE" long_name="kinetic energy: u(n)*u(n+1)/2" unit="W/s^2" />
<!-- variables available when explicit lateral mixing is used (ln_dynldf_OFF=F) -->
<field id="dispkexyfo" long_name="KE-trend: lateral mixing induced dissipation" standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_xy_friction" unit="W/m^2" grid_ref="grid_T_2D" />
<field id="dispkevfo" long_name="KE-trend: vertical mixing induced dissipation" standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_vertical_friction" unit="W/m^2" grid_ref="grid_T_2D" />
<!-- variables available with ln_traadv_eiv=T and ln_diaeiv=T -->
<field id="eketrd_eiv" long_name="EKE-trend due to parameterized eddy advection" standard_name="tendency_of_ocean_eddy_kinetic_energy_content_due_to_parameterized_eddy_advection" unit="W/m^2" grid_ref="grid_T_2D" />
<!-- variables available with ln_PE_trd -->
<field id="petrd_xad" long_name="pe-trend: i-advection" unit="W/m^3" />
<field id="petrd_yad" long_name="pe-trend: j-advection" unit="W/m^3" />
<field id="petrd_zad" long_name="pe-trend: k-advection" unit="W/m^3" />
<field id="petrd_sad" long_name="pe-trend: surface adv. (linssh true)" unit="W/m^3" grid_ref="grid_T_2D" />
<field id="petrd_sad" long_name="pe-trend: surface adv. (linssh true)" unit="W/m^3" grid_ref="grid_T_2D_inner" />
<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" />
......@@ -1072,50 +1145,47 @@ that are available in the tidal-forcing implementation (see
<field id="petrd_nsr" long_name="pe-trend: surface forcing + runoff" unit="W/m^3" />
<field id="petrd_qsr" long_name="pe-trend: solar penetr. heating" unit="W/m^3" />
<field id="petrd_bbc" long_name="pe-trend: geothermal heating" unit="W/m^3" />
<field id="petrd_atf" long_name="pe-trend: asselin time filter" unit="W/m^3" />
<field id="PEanom" long_name="potential energy anomaly" unit="1" />
<field id="alphaPE" long_name="partial deriv. of PEanom wrt T" unit="degC-1" />
<field id="betaPE" long_name="partial deriv. of PEanom wrt S" unit="1e3" />
</field_group>
<field_group id="trendU" grid_ref="grid_U_3D">
<field_group id="trendU" grid_ref="grid_U_3D_inner">
<!-- 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 id="utrd_hpg" long_name="i-trend: hydrostatic pressure gradient" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="utrd_spg" long_name="i-trend: surface pressure gradient" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="utrd_spgexp" long_name="i-trend: surface pressure gradient (explicit)" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="utrd_spgflt" long_name="i-trend: surface pressure gradient (filtered)" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="utrd_keg" long_name="i-trend: KE gradient or hor. adv." unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="utrd_rvo" long_name="i-trend: relative vorticity or metric term" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="utrd_pvo" long_name="i-trend: planetary vorticity" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="utrd_zad" long_name="i-trend: vertical advection" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="utrd_udx" long_name="i-trend: U.dx[U]" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="utrd_ldf" long_name="i-trend: lateral diffusion" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="utrd_zdf" long_name="i-trend: vertical diffusion" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="utrd_tau" long_name="i-trend: wind stress " unit="m/s^2" grid_ref="grid_U_2D_inner" />
<field id="utrd_bfr" long_name="i-trend: bottom friction (explicit)" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="utrd_bfri" long_name="i-trend: bottom friction (implicit)" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="utrd_tot" long_name="i-trend: total momentum trend" unit="m/s^2" grid_ref="grid_T_3D_inner" />
</field_group>
<field_group id="trendV" grid_ref="grid_V_3D">
<field_group id="trendV" grid_ref="grid_V_3D_inner">
<!-- 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 id="vtrd_hpg" long_name="j-trend: hydrostatic pressure gradient" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="vtrd_spg" long_name="j-trend: surface pressure gradient" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="vtrd_spgexp" long_name="j-trend: surface pressure gradient (explicit)" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="vtrd_spgflt" long_name="j-trend: surface pressure gradient (filtered)" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="vtrd_keg" long_name="j-trend: KE gradient or hor. adv." unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="vtrd_rvo" long_name="j-trend: relative vorticity or metric term" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="vtrd_pvo" long_name="j-trend: planetary vorticity" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="vtrd_zad" long_name="j-trend: vertical advection" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="vtrd_vdy" long_name="i-trend: V.dx[V]" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="vtrd_ldf" long_name="j-trend: lateral diffusion" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="vtrd_zdf" long_name="j-trend: vertical diffusion" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="vtrd_tau" long_name="j-trend: wind stress " unit="m/s^2" grid_ref="grid_V_2D_inner" />
<field id="vtrd_bfr" long_name="j-trend: bottom friction (explicit)" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="vtrd_bfri" long_name="j-trend: bottom friction (implicit)" unit="m/s^2" grid_ref="grid_T_3D_inner" />
<field id="vtrd_tot" long_name="j-trend: total momentum trend" unit="m/s^2" grid_ref="grid_T_3D_inner" />
</field_group>
<!-- shared variables available with TOP interface -->
......@@ -1189,6 +1259,8 @@ that are available in the tidal-forcing implementation (see
<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="utau" name="tauuo" long_name="surface_downward_x_stress" />
<field field_ref="vtau" name="tauvo" long_name="surface_downward_y_stress" />
<field field_ref="20d" />
<field field_ref="mldkz5" />
<field field_ref="mldr10_1" />
......@@ -1203,12 +1275,10 @@ that are available in the tidal-forcing implementation (see
<field_group id="groupU" >
<field field_ref="uoce" name="uo" long_name="sea_water_x_velocity" />
<field field_ref="utau" name="tauuo" long_name="surface_downward_x_stress" />
</field_group>
<field_group id="groupV" >
<field field_ref="voce" name="vo" long_name="sea_water_y_velocity" />
<field field_ref="vtau" name="tauvo" long_name="surface_downward_y_stress" />
</field_group>
<field_group id="groupW" >
......@@ -1242,26 +1312,26 @@ that are available in the tidal-forcing implementation (see
<!-- 25h diagnostic output -->
<field_group id="25h_grid_T" grid_ref="grid_T_3D_inner" operation="instant">
<field id="temper25h" name="potential temperature 25h mean" unit="degC" />
<field id="tempis25h" name="insitu temperature 25h mean" unit="degC" />
<field id="salin25h" name="salinity 25h mean" unit="psu" />
<field id="ssh25h" name="sea surface height 25h mean" grid_ref="grid_T_2D_inner" unit="m" />
<field id="temper25h" long_name="potential temperature 25h mean" unit="degC" />
<field id="tempis25h" long_name="insitu temperature 25h mean" unit="degC" />
<field id="salin25h" long_name="salinity 25h mean" unit="psu" />
<field id="ssh25h" long_name="sea surface height 25h mean" grid_ref="grid_T_2D_inner" unit="m" />
</field_group>
<field_group id="25h_grid_U" grid_ref="grid_U_3D_inner" operation="instant" >
<field id="vozocrtx25h" name="i current 25h mean" unit="m/s" />
<field id="vozocrtx25h" long_name="i current 25h mean" unit="m/s" />
</field_group>
<field_group id="25h_grid_V" grid_ref="grid_V_3D_inner" operation="instant">
<field id="vomecrty25h" name="j current 25h mean" unit="m/s" />
<field id="vomecrty25h" long_name="j current 25h mean" unit="m/s" />
</field_group>
<field_group id="25h_grid_W" grid_ref="grid_W_3D_inner" operation="instant">
<field id="vovecrtz25h" name="k current 25h mean" unit="m/s" />
<field id="avt25h" name="vertical diffusivity25h mean" unit="m2/s" />
<field id="avm25h" name="vertical viscosity 25h mean" unit="m2/s" />
<field id="tke25h" name="turbulent kinetic energy 25h mean" />
<field id="mxln25h" name="mixing length 25h mean" unit="m" />
<field id="vovecrtz25h" long_name="k current 25h mean" unit="m/s" />
<field id="avt25h" long_name="vertical diffusivity25h mean" unit="m2/s" />
<field id="avm25h" long_name="vertical viscosity 25h mean" unit="m2/s" />
<field id="tke25h" long_name="turbulent kinetic energy 25h mean" />
<field id="mxln25h" long_name="mixing length 25h mean" unit="m" />
</field_group>
<!--
......
cfgs/SHARED/field_def_nemo-pisces.xml
View file @ dd59ced6
<?xml version="1.0"?>
<!-- $id$ -->
<!--
============================================================================================================
......@@ -105,12 +104,6 @@
<field id="LGW" long_name="Weak ligands concentration" unit="mmol/m3" />
<field id="LGW_e3t" long_name="LGW * e3t" unit="mmol/m2" > LGW * e3t </field >
<!-- PISCES light : variables available with ln_p2z -->
<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 >
</field_group>
<!-- SEDIMENT variables on T sediment grid -->
......@@ -118,201 +111,220 @@
<field id="SedDIC" long_name="Dissolved inorganic Concentration" unit="mmol/m3" />
<field id="SedAlkalini" long_name="Total Alkalinity Concentration" unit="mmol/m3" />
<field id="SedO2" long_name="Oxygen Concentration" unit="mmol/m3" />
<field id="SedCaCO3" long_name="Calcite Concentration" unit="%" />
<field id="SedPOS" long_name="Semi-ref POC Concentration" unit="%" />
<field id="SedPOR" long_name="Refractory POC Concentration" unit="%" />
<field id="SedCaCO3" long_name="Calcite Concentration" unit="g/cm3" />
<field id="SedPOC1" long_name="POC1 Concentration" unit="g/cm3" />
<field id="SedPO4" long_name="Phosphate Concentration" unit="mmol/m3" />
<field id="SedPOC" long_name="Labile POC Concentration" unit="%" />
<field id="SedPOC2" long_name="POC2 Concentration" unit="g/cm3" />
<field id="SedPOC3" long_name="POC3 Concentration" unit="g/cm3" />
<field id="SedPOC4" long_name="POC4 Concentration" unit="g/cm3" />
<field id="SedPOC5" long_name="POC5 Concentration" unit="g/cm3" />
<field id="SedPOC6" long_name="POC6 Concentration" unit="g/cm3" />
<field id="SedSil" long_name="Silicate Concentration" unit="mmol/m3" />
<field id="SedFe2" long_name="Fe2+ Concentration" unit="mmol/m3" />
<field id="SedBSi" long_name="Biogenic Silicate Concentration" unit="%" />
<field id="SedBSi" long_name="Biogenic Silicate Concentration" unit="g/cm3" />
<field id="SedNO3" long_name="Nitrate Concentration" unit="mmol/m3" />
<field id="SedNH4" long_name="Ammonium Concentration" unit="mmol/m3" />
<field id="SedH2S" long_name="H2S Concentration" unit="mmol/m3" />
<field id="SedSO4" long_name="SO4 Concentration" unit="mmol/m3" />
<field id="SedClay" long_name="Clay Concentration" unit="%" />
<field id="SedFeO" long_name="Fe(OH)3 Concentration" unit="%" />
<field id="SedFeS" long_name="FeS Concentration" unit="%" />
<field id="SedpH" long_name="PH" unit="1" />
<field id="SedCO3por" long_name="Bicarbonates" unit="mol/m3" />
<field id="Sedligand" long_name="Ligands" unit="mol/m3" />
<field id="SaturCO3" long_name="CO3 Saturation" unit="-" />
<field id="SedClay" long_name="Clay Concentration" unit="g/cm3" />
<field id="SedFeO" long_name="Fe(OH)3 Concentration" unit="g/cm3" />
<field id="SedFeS" long_name="FeS Concentration" unit="g/cm3" />
<field id="Sedligand" long_name="Ligands" unit="mol/m3" />
</field_group>
<!-- SEDIMENT additional variables on T sediment grid -->
<field_group id="Diag_S" grid_ref="grid_T_2D">
<field id="FlxSi" long_name="Si sediment flux" unit="mol/cm2/s" />
<field id="FlxO2" long_name="O2 sediment flux" unit="mol/cm2/s" />
<field id="FlxDIC" long_name="DIC sediment flux" unit="mol/cm2/s" />
<field id="FlxNO3" long_name="NO3 sediment flux" unit="mol/cm2/s" />
<field id="FlxPO4" long_name="PO4 sediment flux" unit="mol/cm2/s" />
<field id="FlxAlkalini" long_name="Alkalinity sediment flux" unit="mol/cm2/s" />
<field id="FlxNH4" long_name="Ammonium sediment flux" unit="mol/cm2/s" />
<field id="FlxH2S" long_name="H2S sediment flux" unit="mol/cm2/s" />
<field id="FlxSO4" long_name="SO4 sediment flux" unit="mol/cm2/s" />
<field id="FlxFe2" long_name="Fe2+ sediment flux" unit="mol/cm2/s" />
<field id="FlxLig" long_name="Ligand sediment flux" unit="mol/cm2/s" />
<field id="dzdep" long_name="Sedimentation rate" unit="cm/s" />
<field id="sflxclay" long_name="Clay sedimentation rate" unit="g/cm2/s" />
<field id="sflxbsi" long_name="BSi sedimentation rate" unit="g/cm2/s" />
<field id="sflxpoc" long_name="POC sedimentation rate" unit="g/cm2/s" />
<field id="sflxcal" long_name="Calcite sedimentation rate" unit="g/cm2/s" />
<field id="sflxfeo" long_name="Fe(OH)3 Sedimentation rate" unit="mol/cm2/s" />
<field id="FlxClay" long_name="Clay burial rate" unit="g/cm2/s" />
<field id="FlxCaCO3" long_name="Calcite burial rate" unit="g/cm2/s" />
<field id="FlxBSi" long_name="BSi burial rate" unit="g/cm2/s" />
<field id="FlxPOC" long_name="POC burial rate" unit="g/cm2/s" />
<field id="FlxFeO" long_name="Fe(OH)3 burial rate" unit="g/cm2/s" />
<field id="FlxFeS" long_name="FeS burial rate" unit="g/cm2/s" />
<field id="FlxPOR" long_name="POR burial rate" unit="g/cm2/s" />
<field id="FlxPOS" long_name="POS burial rate" unit="g/cm2/s" />
<field id="Flxtot" long_name="total burial flux" unit="g/cm2/s" />
<field id="Rstepros" long_name="Number of iterations" unit="-" />
<field id="SaturCO3" long_name="CO3 Saturation" unit="-" grid_ref="grid_T_3DS" />
<field id="SedpH" long_name="PH" unit="1" grid_ref="grid_T_3DS" />
<field id="SedCO3por" long_name="Bicarbonates" unit="mol/m3" grid_ref="grid_T_3DS" />
<field id="FlxSi" long_name="Si sediment flux" unit="mol/cm2/s" />
<field id="FlxO2" long_name="O2 sediment flux" unit="mol/cm2/s" />
<field id="FlxDIC" long_name="DIC sediment flux" unit="mol/cm2/s" />
<field id="FlxNO3" long_name="NO3 sediment flux" unit="mol/cm2/s" />
<field id="FlxPO4" long_name="PO4 sediment flux" unit="mol/cm2/s" />
<field id="FlxAlkalini" long_name="Alkalinity sediment flux" unit="mol/cm2/s" />
<field id="FlxNH4" long_name="Ammonium sediment flux" unit="mol/cm2/s" />
<field id="FlxH2S" long_name="H2S sediment flux" unit="mol/cm2/s" />
<field id="FlxSO4" long_name="SO4 sediment flux" unit="mol/cm2/s" />
<field id="FlxFe2" long_name="Fe2+ sediment flux" unit="mol/cm2/s" />
<field id="FlxLig" long_name="Ligand sediment flux" unit="mol/cm2/s" />
<field id="dzdep" long_name="Sedimentation rate" unit="cm/s" />
<field id="FlxClay" long_name="Clay burial rate" unit="g/cm2/s" />
<field id="FlxCaCO3" long_name="Calcite burial rate" unit="g/cm2/s" />
<field id="FlxBSi" long_name="BSi burial rate" unit="g/cm2/s" />
<field id="FlxPOC1" long_name="POC1 burial rate" unit="g/cm2/s" />
<field id="FlxFeO" long_name="Fe(OH)3 burial rate" unit="g/cm2/s" />
<field id="FlxFeS" long_name="FeS burial rate" unit="g/cm2/s" />
<field id="FlxPOC2" long_name="POC2 burial rate" unit="g/cm2/s" />
<field id="FlxPOC3" long_name="POC3 burial rate" unit="g/cm2/s" />
<field id="FlxPOC4" long_name="POC4 burial rate" unit="g/cm2/s" />
<field id="FlxPOC5" long_name="POC5 burial rate" unit="g/cm2/s" />
<field id="FlxPOC6" long_name="POC6 burial rate" unit="g/cm2/s" />
<field id="Flxtot" long_name="total burial flux" unit="g/cm2/s" />
<field id="Rstepros" long_name="Number of iterations" unit="-" />
<field id="sflxclay" long_name="Clay sedimentation rate" unit="g/cm2/s" />
<field id="sflxbsi" long_name="BSi sedimentation rate" unit="g/cm2/s" />
<field id="sflxpoc" long_name="POC sedimentation rate" unit="g/cm2/s" />
<field id="sflxcal" long_name="Calcite sedimentation rate" unit="mol/cm2/s" />
<field id="OceDIC" long_name="DIC at ocean bottom " unit="mmol/m3" />
<field id="OceALK" long_name="ALK at ocean bottom " unit="mmol/m3" />
<field id="OceOXY" long_name="OXY at ocean bottom " unit="mmol/m3" />
<field id="OceSIL" long_name="SIL at ocean bottom " unit="mmol/m3" />
<field id="OceNO3" long_name="NO3 at ocean bottom " unit="mmol/m3" />
<field id="OcePO4" long_name="PO4 at ocean bottom " unit="mmol/m3" />
<field id="OceNH4" long_name="NH4 at ocean bottom " unit="mmol/m3" />
<field id="OceSO4" long_name="SO4 at ocean bottom " unit="mmol/m3" />
<field id="OceFE2" long_name="Dissolved Iron at ocean bottom " unit="mmol/m3" />
<field id="OceTEM" long_name="Temperature at ocean bottom " unit="degC" />
<field id="OceSAL" long_name="Salinity at oceaan bottom " unit="1e-3" />
<field id="OceSFE" long_name="Sinking Iron at ocean bottom " unit="-" />
</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="PPPHYN" long_name="Primary production of nanophyto" unit="molC/m3/s" grid_ref="grid_T_3D" />
<field id="PPPHYP" long_name="Primary production of picophyto" unit="molC/m3/s" grid_ref="grid_T_3D" />
<field id="PPPHYD" long_name="Primary production of diatoms" unit="molC/m3/s" grid_ref="grid_T_3D" />
<field id="PPNEWN" long_name="New Primary production of nanophyto" unit="molC/m3/s" grid_ref="grid_T_3D" />
<field id="PPNEWP" long_name="New Primary production of picophyto" unit="molC/m3/s" grid_ref="grid_T_3D" />
<field id="PPNEWD" long_name="New Primary production of diatoms" unit="molC/m3/s" grid_ref="grid_T_3D" />
<field id="PBSi" long_name="Primary production of Si diatoms" unit="molC/m3/s" grid_ref="grid_T_3D" />
<field id="PFeN" long_name="Primary production of nano iron" unit="molC/m3/s" grid_ref="grid_T_3D" />
<field id="PFeP" long_name="Primary production of pico iron" unit="molC/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="REMINP" long_name="Oxic remineralization rate of POC" unit="d-1" grid_ref="grid_T_3D" />
<field id="REMING" long_name="Oxic remineralization rate of GOC" unit="d-1" 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="MuP" long_name="Realized growth rate for picophyto" 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="MunetN" long_name="Net growth rate for nanophyto" unit="s-1" grid_ref="grid_T_3D" />
<field id="MunetP" long_name="Net growth rate for picophyto" unit="s-1" grid_ref="grid_T_3D" />
<field id="MunetD" long_name="Net 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="LPnut" long_name="Nutrient limitation term in Picophyto" 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="LPFe" long_name="Iron limitation term in Picophyto" 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="LPlight" long_name="Light limitation term in Picophyto" unit="-" grid_ref="grid_T_3D" />
<field id="LDlight" long_name="Light limitation term in Diatoms" unit="" grid_ref="grid_T_3D" />
<field id="SIZEN" long_name="Mean relative size of nanophyto." unit="-" grid_ref="grid_T_3D" />
<field id="SIZEP" long_name="Mean relative size of picophyto." unit="-" grid_ref="grid_T_3D" />
<field id="SIZED" long_name="Mean relative size of diatoms" unit="-" grid_ref="grid_T_3D" />
<field id="RASSD" long_name="Size of the protein machinery (Diat.)" unit="-" grid_ref="grid_T_3D" />
<field id="RASSN" long_name="Size of the protein machinery (Nano.)" unit="-" grid_ref="grid_T_3D" />
<field id="RASSP" long_name="Size of the protein machinery (Pico.)" unit="-" 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="TL1" long_name="Total L1 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="mol/m2/s" />
<field id="Ironice" long_name="Iron input/uptake due to sea ice" unit="mol/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="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="pCO2sea" long_name="surface ocean pCO2" unit="uatm" />
<field id="Dpo2" long_name="Delta O2" unit="uatm" />
<field id="Heup" long_name="Euphotic layer depth" unit="m" />
<field id="AtmCo2" long_name="Atmospheric CO2 concentration" unit="ppm" />
<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" />
<field id="FESCAV" long_name="Scavenging of Iron" unit="mmol-Fe/m3/s" grid_ref="grid_T_3D" />
<field id="FECOLL" long_name="Colloidal Pumping of FeL" unit="mmol-FeL/m3/s" grid_ref="grid_T_3D" />
<field id="LGWCOLL" long_name="Coagulation loss of ligands" unit="mmol-L/m3/s" grid_ref="grid_T_3D" />
<field id="REMINF" long_name="Oxic remineralization suppy of Fe" unit="mmol-Fe/m3/s" grid_ref="grid_T_3D" />
<field id="BACT" long_name="Bacterial Biomass" unit="mmol/m3" grid_ref="grid_T_3D" />
<field id="FEBACT" long_name="Bacterial uptake of Fe" unit="molFe/m3/s" grid_ref="grid_T_3D" />
<field id="FEPREC" long_name="Precipitation of Fe" unit="molFe/m3/s" grid_ref="grid_T_3D" />
<field id="LPRODR" long_name="OM remineralisation ligand production rate" unit="nmol-L/m3/s" grid_ref="grid_T_3D" />
<field id="LPRODP" long_name="phytoplankton ligand production rate" unit="nmol-L/m3/s" grid_ref="grid_T_3D" />
<field id="LIGREM" long_name="Remineralisation loss of ligands" unit="nmol-L/m3/s" grid_ref="grid_T_3D" />
<field id="LIGPR" long_name="Photochemical loss of ligands" unit="nmol-L/m3/s" grid_ref="grid_T_3D" />
<field id="LDETP" long_name="Ligand destruction during phytoplankton uptake" unit="nmol-L/m3/s" grid_ref="grid_T_3D" />
<field id="LPRODZ2" long_name="mesozooplankton ligand production rate" unit="nmol-L/m3/s" grid_ref="grid_T_3D" />
<field id="LPRODZ" long_name="microzooplankton ligand production rate" unit="nmol-L/m3/s" grid_ref="grid_T_3D" />
<field id="FEZOO" long_name="microzooplankton iron recycling rate" unit="nmol-FeL/m3/s" grid_ref="grid_T_3D" />
<field id="FEZOO2" long_name="mesozooplankton iron recycling rate" unit="nmol-FeL/m3/s" grid_ref="grid_T_3D" />
<field_group id="diad_T" grid_ref="grid_T_2D_inner" >
<field id="PH" long_name="PH" unit="1" grid_ref="grid_T_3D_inner" />
<field id="CO3" long_name="Bicarbonates" unit="mol/m3" grid_ref="grid_T_3D_inner" />
<field id="CO3sat" long_name="CO3 saturation" unit="mol/m3" grid_ref="grid_T_3D_inner" />
<field id="DCAL" long_name="Calcite dissolution" unit="mol/m3/s" grid_ref="grid_T_3D_inner" />
<field id="PAR" long_name="Photosynthetically Available Radiation" unit="W/m2" grid_ref="grid_T_3D_inner" />
<field id="PPPHYN" long_name="Primary production of nanophyto" unit="molC/m3/s" grid_ref="grid_T_3D_inner" />
<field id="PPPHYP" long_name="Primary production of picophyto" unit="molC/m3/s" grid_ref="grid_T_3D_inner" />
<field id="PPPHYD" long_name="Primary production of diatoms" unit="molC/m3/s" grid_ref="grid_T_3D_inner" />
<field id="PPNEWN" long_name="New Primary production of nanophyto" unit="molC/m3/s" grid_ref="grid_T_3D_inner" />
<field id="PPNEWP" long_name="New Primary production of picophyto" unit="molC/m3/s" grid_ref="grid_T_3D_inner" />
<field id="PPNEWD" long_name="New Primary production of diatoms" unit="molC/m3/s" grid_ref="grid_T_3D_inner" />
<field id="GPPHYN" long_name="Gross Primary production of nanophyto" unit="molC/m3/s" grid_ref="grid_T_3D_inner" />
<field id="GPPHYP" long_name="Gross Primary production of picophyto" unit="molC/m3/s" grid_ref="grid_T_3D_inner" />
<field id="GPPHYD" long_name="Gross Primary production of diatoms" unit="molC/m3/s" grid_ref="grid_T_3D_inner" />
<field id="PBSi" long_name="Primary production of Si diatoms" unit="molC/m3/s" grid_ref="grid_T_3D_inner" />
<field id="PFeN" long_name="Primary production of nano iron" unit="molC/m3/s" grid_ref="grid_T_3D_inner" />
<field id="PFeP" long_name="Primary production of pico iron" unit="molC/m3/s" grid_ref="grid_T_3D_inner" />
<field id="PFeD" long_name="Primary production of diatoms iron" unit="mol/m3/s" grid_ref="grid_T_3D_inner" />
<field id="xfracal" long_name="Calcifying fraction" unit="1" grid_ref="grid_T_3D_inner" />
<field id="PCAL" long_name="Calcite production" unit="mol/m3/s" grid_ref="grid_T_3D_inner" />
<field id="GRAZ1" long_name="Grazing by microzooplankton" unit="mol/m3/s" grid_ref="grid_T_3D_inner" />
<field id="GRAZ2" long_name="Grazing by mesozooplankton" unit="mol/m3/s" grid_ref="grid_T_3D_inner" />
<field id="REMIN" long_name="Oxic remineralization of OM" unit="mol/m3/s" grid_ref="grid_T_3D_inner" />
<field id="DENIT" long_name="Anoxic remineralization of OM" unit="mol/m3/s" grid_ref="grid_T_3D_inner" />
<field id="REMINP" long_name="Oxic remineralization rate of POC" unit="d-1" grid_ref="grid_T_3D_inner" />
<field id="REMING" long_name="Oxic remineralization rate of GOC" unit="d-1" grid_ref="grid_T_3D_inner" />
<field id="Nfix" long_name="Nitrogen fixation" unit="mol/m3/s" grid_ref="grid_T_3D_inner" />
<field id="Mumax" long_name="Maximum growth rate" unit="s-1" grid_ref="grid_T_3D_inner" />
<field id="MuN" long_name="Realized growth rate for nanophyto" unit="s-1" grid_ref="grid_T_3D_inner" />
<field id="MuP" long_name="Realized growth rate for picophyto" unit="s-1" grid_ref="grid_T_3D_inner" />
<field id="MuD" long_name="Realized growth rate for diatomes" unit="s-1" grid_ref="grid_T_3D_inner" />
<field id="MunetN" long_name="Net growth rate for nanophyto" unit="s-1" grid_ref="grid_T_3D_inner" />
<field id="MunetP" long_name="Net growth rate for picophyto" unit="s-1" grid_ref="grid_T_3D_inner" />
<field id="MunetD" long_name="Net growth rate for diatomes" unit="s-1" grid_ref="grid_T_3D_inner" />
<field id="LNnut" long_name="Nutrient limitation term in Nanophyto" unit="" grid_ref="grid_T_3D_inner" />
<field id="LPnut" long_name="Nutrient limitation term in Picophyto" unit="-" grid_ref="grid_T_3D_inner" />
<field id="LDnut" long_name="Nutrient limitation term in Diatoms" unit="" grid_ref="grid_T_3D_inner" />
<field id="LNFe" long_name="Iron limitation term in Nanophyto" unit="" grid_ref="grid_T_3D_inner" />
<field id="LPFe" long_name="Iron limitation term in Picophyto" unit="-" grid_ref="grid_T_3D_inner" />
<field id="LDFe" long_name="Iron limitation term in Diatoms" unit="" grid_ref="grid_T_3D_inner" />
<field id="LNlight" long_name="Light limitation term in Nanophyto" unit="" grid_ref="grid_T_3D_inner" />
<field id="LPlight" long_name="Light limitation term in Picophyto" unit="-" grid_ref="grid_T_3D_inner" />
<field id="LDlight" long_name="Light limitation term in Diatoms" unit="" grid_ref="grid_T_3D_inner" />
<field id="SIZEN" long_name="Mean relative size of nanophyto." unit="-" grid_ref="grid_T_3D_inner" />
<field id="SIZEP" long_name="Mean relative size of picophyto." unit="-" grid_ref="grid_T_3D_inner" />
<field id="SIZED" long_name="Mean relative size of diatoms" unit="-" grid_ref="grid_T_3D_inner" />
<field id="THETANANO" long_name="Diagnostic Chl:C ratio" unit="mgChl/mgC" grid_ref="grid_T_3D_inner" />
<field id="CHL" long_name="Total Chl" unit="gChl/m3" grid_ref="grid_T_3D_inner"/>
<field id="RASSD" long_name="Size of the protein machinery (Diat.)" unit="-" grid_ref="grid_T_3D_inner" />
<field id="RASSN" long_name="Size of the protein machinery (Nano.)" unit="-" grid_ref="grid_T_3D_inner" />
<field id="RASSP" long_name="Size of the protein machinery (Pico.)" unit="-" grid_ref="grid_T_3D_inner" />
<field id="Fe3" long_name="Iron III concentration" unit="nmol/m3" grid_ref="grid_T_3D_inner" />
<field id="FeL1" long_name="Complexed Iron concentration with L1" unit="nmol/m3" grid_ref="grid_T_3D_inner" />
<field id="TL1" long_name="Total L1 concentration" unit="nmol/m3" grid_ref="grid_T_3D_inner" />
<field id="pdust" long_name="dust concentration" unit="kg/m3" />
<field id="Totlig" long_name="Total ligand concentation" unit="nmol/m3" grid_ref="grid_T_3D_inner" />
<field id="Biron" long_name="Bioavailable iron" unit="nmol/m3" grid_ref="grid_T_3D_inner" />
<field id="Sdenit" long_name="Nitrate reduction in the sediments" unit="mol/m2/s" />
<field id="Ironice" long_name="Iron input/uptake due to sea ice" unit="mol/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="HYDR" long_name="Iron input from hydrothemal vents" unit="mol/m2/s" grid_ref="grid_T_3D_inner" />
<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_inner" />
<field id="EXPFE" long_name="Export of biogenic iron" unit="mol/m2/s" grid_ref="grid_T_3D_inner" />
<field id="EXPSI" long_name="Export of Silicate" unit="mol/m2/s" grid_ref="grid_T_3D_inner" />
<field id="EXPCAL" long_name="Export of Calcite" unit="mol/m2/s" grid_ref="grid_T_3D_inner" />
<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="pCO2sea" long_name="surface ocean pCO2" unit="uatm" />
<field id="Dpo2" long_name="Delta O2" unit="uatm" />
<field id="Heup" long_name="Euphotic layer depth" unit="m" grid_ref="grid_T_2D_inner" />
<field id="AtmCo2" long_name="Atmospheric CO2 concentration" unit="ppm" />
<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_inner" />
<field id="FESCAV" long_name="Scavenging of Iron" unit="mmol-Fe/m3/s" grid_ref="grid_T_3D_inner" />
<field id="FECOLL" long_name="Colloidal Pumping of FeL" unit="mmol-FeL/m3/s" grid_ref="grid_T_3D_inner" />
<field id="LGWCOLL" long_name="Coagulation loss of ligands" unit="mmol-L/m3/s" grid_ref="grid_T_3D_inner" />
<field id="REMINF" long_name="Oxic remineralization suppy of Fe" unit="mmol-Fe/m3/s" grid_ref="grid_T_3D_inner" />
<field id="BACT" long_name="Bacterial Biomass" unit="mmol/m3" grid_ref="grid_T_3D_inner" />
<field id="FEBACT" long_name="Bacterial uptake of Fe" unit="molFe/m3/s" grid_ref="grid_T_3D_inner" />
<field id="FEPREC" long_name="Precipitation of Fe" unit="molFe/m3/s" grid_ref="grid_T_3D_inner" />
<field id="LPRODR" long_name="OM remineralisation ligand production rate" unit="nmol-L/m3/s" grid_ref="grid_T_3D_inner" />
<field id="LPRODP" long_name="phytoplankton ligand production rate" unit="nmol-L/m3/s" grid_ref="grid_T_3D_inner" />
<field id="LIGREM" long_name="Remineralisation loss of ligands" unit="nmol-L/m3/s" grid_ref="grid_T_3D_inner" />
<field id="LIGPR" long_name="Photochemical loss of ligands" unit="nmol-L/m3/s" grid_ref="grid_T_3D_inner" />
<field id="LDETP" long_name="Ligand destruction during phytoplankton uptake" unit="nmol-L/m3/s" grid_ref="grid_T_3D_inner" />
<field id="LPRODZ2" long_name="mesozooplankton ligand production rate" unit="nmol-L/m3/s" grid_ref="grid_T_3D_inner" />
<field id="LPRODZ" long_name="microzooplankton ligand production rate" unit="nmol-L/m3/s" grid_ref="grid_T_3D_inner" />
<field id="FEZOO" long_name="microzooplankton iron recycling rate" unit="nmol-FeL/m3/s" grid_ref="grid_T_3D_inner" />
<field id="FEZOO2" long_name="mesozooplankton iron recycling rate" unit="nmol-FeL/m3/s" grid_ref="grid_T_3D_inner" />
<field id="MicroZo2" long_name="O2 consumption by microzooplankton" unit="molC/m2/s" grid_ref="grid_T_3D_inner" />
<field id="Remino2" long_name="O2 consumption by remineralisation" unit="molC/m2/s" grid_ref="grid_T_3D_inner" />
<field id="Nfixo2" long_name="O2 production by Nfix" unit="molN/m2/s" grid_ref="grid_T_3D_inner" />
<field id="PPNEWo2" long_name="Oxygen production by the New Produc." unit="molC/m3/s" grid_ref="grid_T_3D_inner" />
<field id="PPRego2" long_name="Oxygen production by the Regen Produc." unit="molC/m3/s" grid_ref="grid_T_3D_inner" />
<field id="MesoZo2" long_name="O2 consumption by mesozooplankton" unit="molC/m2/s" grid_ref="grid_T_3D_inner" />
<field id="Nitrifo2" long_name="O2 consumption by nitrification" unit="molC/m2/s" grid_ref="grid_T_3D_inner" />
<field id="No3dep" long_name="Nitrate deposition at surface" unit="molN/m2/s" />
<field id="Sildep" long_name="silicate deposition at surface" unit="molC/m2/s" />
<field id="Po4dep" long_name="Phosphorus deposition at surface" unit="molP/m2/s" />
<field id="Nh4dep" long_name="Ammonium deposition at surface" unit="molN/m2/s" />
<!-- PISCES tracers trends -->
<field id="INTdtAlk" long_name="Vertically int. of change of alkalinity" unit="mol/m2/s" />
<field id="INTdtDIC" long_name="Vertically int. of change of dissic " unit="mol/m2/s" />
<field id="INTdtFer" long_name="Vertically int. of change of iron " unit="mol/m2/s" />
<field id="INTdtDIN" long_name="Vertically int. of change of nitrogen " unit="mol/m2/s" />
<field id="INTdtDIP" long_name="Vertically int. of change of phophate " unit="mol/m2/s" />
<field id="INTdtSil" long_name="Vertically int. of change of silicon " unit="mol/m2/s" />
<field id="INTdtAlk" long_name="Vertically int. of change of alkalinity" unit="mol/m2/s" />
<field id="INTdtDIC" long_name="Vertically int. of change of dissic " unit="mol/m2/s" />
<field id="INTdtFer" long_name="Vertically int. of change of iron " unit="mol/m2/s" />
<field id="INTdtDIN" long_name="Vertically int. of change of nitrogen " unit="mol/m2/s" />
<field id="INTdtDIP" long_name="Vertically int. of change of phophate " unit="mol/m2/s" />
<field id="INTdtSil" long_name="Vertically int. of change of silicon " unit="mol/m2/s" />
<field id="TPP" long_name="Total Primary production of phyto" unit="mol/m3/s" grid_ref="grid_T_3D_inner" />
<field id="TPNEW" long_name="New Primary production of phyto" unit="mol/m3/s" grid_ref="grid_T_3D_inner" />
<field id="TPBFE" long_name="Total biogenic iron production" unit="mol/m3/s" grid_ref="grid_T_3D_inner" />
<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="Nfix_e3t" long_name="Nfix * e3t" unit="molN/m2/s" grid_ref="grid_T_3D_inner" > Nfix * e3t </field >
<field id="PPPHYN_e3t" long_name="PPPHYN * e3t" unit="mol/m2/s" grid_ref="grid_T_3D_inner" > PPPHYN * e3t </field >
<field id="PPPHYD_e3t" long_name="PPPHYD * e3t" unit="mol/m2/s" grid_ref="grid_T_3D_inner" > PPPHYD * e3t </field >
<field id="PPPHYP_e3t" long_name="PPPHYP * e3t" unit="mol/m2/s" grid_ref="grid_T_3D_inner" > PPPHYP * e3t </field >
<field id="PBSi_e3t" long_name="PBSi * e3t" unit="mol/m2/s" grid_ref="grid_T_3D_inner" > PBSi * e3t </field >
<field id="TPP_e3t" long_name="TPP * e3t" unit="mol/m2/s" grid_ref="grid_T_3D_inner" > TPP * e3t </field >
<field id="TPNEW_e3t" long_name="TPNEW * e3t" unit="mol/m2/s" grid_ref="grid_T_3D_inner" > TPNEW * e3t </field >
<field id="TPBFE_e3t" long_name="TPBFE * e3t" unit="mol/m2/s" grid_ref="grid_T_3D_inner" > TPBFE * e3t </field >
<!-- dbio_T on T grid : variables available with 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" grid_ref="grid_T_vsum" detect_missing_value="true" > Nfix * e3t </field >
<field id="INTPPPHYN" long_name="Vertically integrated primary production by nanophy" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" > PPPHYN * e3t </field >
<field id="INTPPPHYD" long_name="Vertically integrated primary production by diatom" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" > PPPHYD * e3t </field >
<field id="INTPPPHYP" long_name="Vertically integrated primary production by picophy" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" > PPPHYP * e3t </field >
<field id="INTPP" long_name="Vertically integrated primary production by phyto" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" > TPP * e3t </field >
<field id="INTPNEW" long_name="Vertically integrated new primary production" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" > TPNEW * e3t </field >
<field id="INTPBFE" long_name="Vertically integrated of biogenic iron production" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" > TPBFE * e3t </field >
<field id="INTPBSI" long_name="Vertically integrated of biogenic Si production" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" > PBSi * e3t </field >
<field id="INTNFIX" long_name="Nitrogen fixation rate : vert. integrated" field_ref="Nfix_e3t" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" />
<field id="INTPPPHYN" long_name="Vertically integrated primary production by nanophy" field_ref="PPPHYN_e3t" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" />
<field id="INTPPPHYD" long_name="Vertically integrated primary production by diatom" field_ref="PPPHYD_e3t" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" />
<field id="INTPPPHYP" long_name="Vertically integrated primary production by picophy" field_ref="PPPHYP_e3t" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" />
<field id="INTPP" long_name="Vertically integrated primary production by phyto" field_ref="TPP_e3t" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" />
<field id="INTPNEW" long_name="Vertically integrated new primary production" field_ref="TPNEW_e3t" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" />
<field id="INTPBFE" long_name="Vertically integrated of biogenic iron production" field_ref="TPBFE_e3t" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" />
<field id="INTPBSI" long_name="Vertically integrated of biogenic Si production" field_ref="PBSi_e3t" unit="mol/m2/s" grid_ref="grid_T_vsum" detect_missing_value="true" />
<!-- 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="" />
</field_group>
<field_group id="tracer_scalar" grid_ref="grid_scalar" >
......
cfgs/SHARED/grid_def_nemo.xml
View file @ dd59ced6
......@@ -12,7 +12,7 @@
<domain domain_ref="grid_T" />
</grid>
<grid id="grid_T_2D_inner" >
<domain domain_ref="grid_T_inner" />
<domain domain_ref="grid_T_inner" name="grid_T" /> <!-- use name="grid_T" so we don't duplicate x, y, dimensions -->
</grid>
<!-- -->
<grid id="grid_T_ncatice" >
......@@ -20,16 +20,36 @@
<axis axis_ref="ncatice" />
</grid>
<grid id="grid_T_ncatice_inner" >
<domain domain_ref="grid_T_inner" />
<domain domain_ref="grid_T_inner" name="grid_T" /> <!-- use name="grid_T" so we don't duplicate x, y, dimensions -->
<axis axis_ref="ncatice" />
</grid>
<!-- -->
<grid id="grid_T_nlayice" >
<domain domain_ref="grid_T" />
<axis axis_ref="nlayice" />
<axis axis_ref="ncatice" />
</grid>
<grid id="grid_T_nlayice_inner" >
<domain domain_ref="grid_T_inner" name="grid_T" /> <!-- use name="grid_T" so we don't duplicate x, y, dimensions -->
<axis axis_ref="nlayice" />
<axis axis_ref="ncatice" />
</grid>
<!-- -->
<grid id="grid_T_icbbas" >
<domain domain_ref="grid_T" />
<axis axis_ref="icbbas" />
</grid>
<grid id="grid_T_icbbas_inner" >
<domain domain_ref="grid_T_inner" name="grid_T" />
<axis axis_ref="icbbas" />
</grid>
<!-- -->
<grid id="grid_T_3D" >
<domain domain_ref="grid_T" />
<axis axis_ref="deptht" />
</grid>
<grid id="grid_T_3D_inner" >
<domain domain_ref="grid_T_inner" />
<domain domain_ref="grid_T_inner" name="grid_T" /> <!-- use name="grid_T" so we don't duplicate x, y, dimensions -->
<axis axis_ref="deptht" />
</grid>
<!-- -->
......@@ -38,7 +58,7 @@
<axis axis_ref="profsed" />
</grid>
<grid id="grid_T_3DS_inner" >
<domain domain_ref="grid_T_inner" />
<domain domain_ref="grid_T_inner" name="grid_T" /> <!-- use name="grid_T" so we don't duplicate x, y, dimensions -->
<axis axis_ref="profsed" />
</grid>
<!-- -->
......@@ -46,7 +66,7 @@
<domain domain_ref="grid_U" />
</grid>
<grid id="grid_U_2D_inner" >
<domain domain_ref="grid_U_inner" />
<domain domain_ref="grid_U_inner" name="grid_U" /> <!-- use name="grid_U" so we don't duplicate x, y, dimensions -->
</grid>
<!-- -->
<grid id="grid_U_3D" >
......@@ -54,7 +74,7 @@
<axis axis_ref="depthu" />
</grid>
<grid id="grid_U_3D_inner" >
<domain domain_ref="grid_U_inner" />
<domain domain_ref="grid_U_inner" name="grid_U" /> <!-- use name="grid_U" so we don't duplicate x, y, dimensions -->
<axis axis_ref="depthu" />
</grid>
<!-- -->
......@@ -62,7 +82,7 @@
<domain domain_ref="grid_V" />
</grid>
<grid id="grid_V_2D_inner" >
<domain domain_ref="grid_V_inner" />
<domain domain_ref="grid_V_inner" name="grid_V" /> <!-- use name="grid_V" so we don't duplicate x, y, dimensions -->
</grid>
<!-- -->
<grid id="grid_V_3D" >
......@@ -70,7 +90,7 @@
<axis axis_ref="depthv" />
</grid>
<grid id="grid_V_3D_inner" >
<domain domain_ref="grid_V_inner" />
<domain domain_ref="grid_V_inner" name="grid_V" /> <!-- use name="grid_V" so we don't duplicate x, y, dimensions -->
<axis axis_ref="depthv" />
</grid>
<!-- -->
......@@ -78,7 +98,7 @@
<domain domain_ref="grid_W" />
</grid>
<grid id="grid_W_2D_inner" >
<domain domain_ref="grid_W_inner" />
<domain domain_ref="grid_W_inner" name="grid_W" /> <!-- use name="grid_W" so we don't duplicate x, y, dimensions -->
</grid>
<!-- -->
<grid id="grid_W_3D" >
......@@ -86,7 +106,7 @@
<axis axis_ref="depthw" />
</grid>
<grid id="grid_W_3D_inner" >
<domain domain_ref="grid_W_inner" />
<domain domain_ref="grid_W_inner" name="grid_W" /> <!-- use name="grid_W" so we don't duplicate x, y, dimensions -->
<axis axis_ref="depthw" />
</grid>
<!-- -->
......@@ -94,7 +114,7 @@
<domain domain_ref="grid_F" />
</grid>
<grid id="grid_F_2D_inner" >
<domain domain_ref="grid_F_inner" />
<domain domain_ref="grid_F_inner" name="grid_F" /> <!-- use name="grid_F" so we don't duplicate x, y, dimensions -->
</grid>
<!-- -->
<grid id="grid_F_3D" >
......@@ -102,7 +122,7 @@
<axis axis_ref="depthf" />
</grid>
<grid id="grid_F_3D_inner" >
<domain domain_ref="grid_F_inner" />
<domain domain_ref="grid_F_inner" name="grid_F" /> <!-- use name="grid_F" so we don't duplicate x, y, dimensions -->
<axis axis_ref="depthf" />
</grid>
<!-- -->
......@@ -110,11 +130,6 @@
<domain domain_ref="1point"/>
</grid>
<!-- -->
<grid id="grid_T_nfloat" >
<domain domain_ref="grid_T" />
<axis axis_ref="nfloat" />
</grid>
<!-- -->
<grid id="grid_EqT" >
<domain domain_ref="EqT" />
</grid>
......@@ -184,6 +199,12 @@
<extract_axis position="0" />
</scalar>
</grid>
<grid id="grid_T_SFC_inner">
<domain domain_ref="grid_T_inner" name="grid_T" /> <!-- use name="grid_T" so we don't duplicate x, y, dimensions -->
<scalar>
<extract_axis position="0" />
</scalar>
</grid>
<grid id="grid_T_vsum">
<domain domain_ref="grid_T"/>
......@@ -191,6 +212,12 @@
<reduce_axis operation="sum" />
</scalar>
</grid>
<grid id="grid_T_vsum_inner">
<domain domain_ref="grid_T_inner" name="grid_T"/> <!-- use name="grid_T" so we don't duplicate x, y, dimensions -->
<scalar>
<reduce_axis operation="sum" />
</scalar>
</grid>
<grid id="grid_U_vsum">
<domain domain_ref="grid_U"/>
......@@ -232,6 +259,10 @@
<domain domain_ref="grid_T" />
<axis axis_ref="deptht300" />
</grid>
<grid id="grid_T_zoom_300_inner">
<domain domain_ref="grid_T_inner" name="grid_T" /> <!-- use name="grid_T" so we don't duplicate x, y, dimensions -->
<axis axis_ref="deptht300" />
</grid>
<grid id="grid_U_scalar" >
<domain domain_ref="grid_U" />
......@@ -318,14 +349,14 @@
<!-- ABL grid definition -->
<grid id="grid_TA_2D">
<domain domain_ref="grid_T" />
<domain domain_ref="grid_T_inner" name="grid_T" /> <!-- use name="grid_T" so we don't duplicate x, y, dimensions -->
</grid>
<grid id="grid_TA_3D">
<domain domain_ref="grid_T" />
<domain domain_ref="grid_T_inner" name="grid_T" /> <!-- use name="grid_T" so we don't duplicate x, y, dimensions -->
<axis id="ght_abl" />
</grid>
<grid id="grid_WA_3D">
<domain domain_ref="grid_T" />
<domain domain_ref="grid_T_inner" name="grid_T" /> <!-- use name="grid_T" so we don't duplicate x, y, dimensions -->
<axis id="ghw_abl" />
</grid>
<!-- -->
......@@ -336,84 +367,87 @@
<scalar />
</grid>
<grid id="diamlr_grid_T_2D" >
<domain domain_ref="grid_T" />
<domain domain_ref="grid_T_inner" name="grid_T" /> <!-- use name="grid_T" so we don't duplicate x, y, dimensions -->
<scalar />
</grid>
<grid id="diamlr_grid_T_2D_inner" >
<domain domain_ref="grid_T_inner" name="grid_T" /> <!-- use name="grid_T" so we don't duplicate x, y, dimensions -->
</grid>
<grid id="diamlr_grid_U_2D" >
<domain domain_ref="grid_U" />
<domain domain_ref="grid_U_inner" name="grid_U" /> <!-- use name="grid_U" so we don't duplicate x, y, dimensions -->
<scalar />
</grid>
<grid id="diamlr_grid_V_2D" >
<domain domain_ref="grid_V" />
<domain domain_ref="grid_V_inner" name="grid_V" /> <!-- use name="grid_V" so we don't duplicate x, y, dimensions -->
<scalar />
</grid>
<grid id="diamlr_grid_W_2D" >
<domain domain_ref="grid_W" />
<domain domain_ref="grid_W_inner" name="grid_W" /> <!-- use name="grid_W" so we don't duplicate x, y, dimensions -->
<scalar />
</grid>
<grid id="diamlr_grid_2D_to_grid_T_3D" >
<domain domain_ref="grid_T" />
<domain domain_ref="grid_T_inner" name="grid_T" /> <!-- use name="grid_T" so we don't duplicate x, y, dimensions -->
<axis axis_ref="deptht">
<duplicate_scalar />
</axis>
</grid>
<grid id="diamlr_grid_2D_to_grid_U_3D" >
<domain domain_ref="grid_U" />
<domain domain_ref="grid_U_inner" name="grid_U" /> <!-- use name="grid_U" so we don't duplicate x, y, dimensions -->
<axis axis_ref="depthu">
<duplicate_scalar />
</axis>
</grid>
<grid id="diamlr_grid_2D_to_grid_V_3D" >
<domain domain_ref="grid_V" />
<domain domain_ref="grid_V_inner" name="grid_V" /> <!-- use name="grid_V" so we don't duplicate x, y, dimensions -->
<axis axis_ref="depthv">
<duplicate_scalar />
</axis>
</grid>
<grid id="diamlr_grid_2D_to_grid_W_3D" >
<domain domain_ref="grid_W" />
<domain domain_ref="grid_W_inner" name="grid_W" /> <!-- use name="grid_W" so we don't duplicate x, y, dimensions -->
<axis axis_ref="depthw">
<duplicate_scalar />
</axis>
</grid>
<grid id="diamlr_grid_2D_to_scalar" >
<scalar>
<reduce_domain operation="average" />
<reduce_domain operation="average" local="true" />
</scalar>
<scalar />
</grid>
<!-- grid definitions for the computation of daily detided model diagnostics (diadetide) -->
<grid id="diadetide_grid_T_2D" >
<domain domain_ref="grid_T" />
<domain domain_ref="grid_T_inner" name="grid_T" /> <!-- use name="grid_T" so we don't duplicate x, y, dimensions -->
<scalar />
</grid>
<grid id="diadetide_grid_U_2D" >
<domain domain_ref="grid_U" />
<domain domain_ref="grid_U_inner" name="grid_U" /> <!-- use name="grid_U" so we don't duplicate x, y, dimensions -->
<scalar />
</grid>
<grid id="diadetide_grid_V_2D" >
<domain domain_ref="grid_V" />
<domain domain_ref="grid_V_inner" name="grid_V" /> <!-- use name="grid_V" so we don't duplicate x, y, dimensions -->
<scalar />
</grid>
<grid id="diadetide_grid_2D_to_grid_T_3D" >
<domain domain_ref="grid_T" />
<domain domain_ref="grid_T_inner" name="grid_T" /> <!-- use name="grid_T" so we don't duplicate x, y, dimensions -->
<axis axis_ref="deptht">
<duplicate_scalar />
</axis>
</grid>
<grid id="diadetide_grid_2D_to_grid_U_3D" >
<domain domain_ref="grid_U" />
<domain domain_ref="grid_U_inner" name="grid_U" /> <!-- use name="grid_U" so we don't duplicate x, y, dimensions -->
<axis axis_ref="depthu">
<duplicate_scalar />
</axis>
</grid>
<grid id="diadetide_grid_2D_to_grid_V_3D" >
<domain domain_ref="grid_V" />
<domain domain_ref="grid_V_inner" name="grid_V" /> <!-- use name="grid_V" so we don't duplicate x, y, dimensions -->
<axis axis_ref="depthv">
<duplicate_scalar />
</axis>
</grid>
<grid id="diadetide_grid_2D_to_grid_W_3D" >
<domain domain_ref="grid_W" />
<domain domain_ref="grid_W_inner" name="grid_W" /> <!-- use name="grid_W" so we don't duplicate x, y, dimensions -->
<axis axis_ref="depthw">
<duplicate_scalar />
</axis>
......
cfgs/SHARED/namelist_ice_ref
View file @ dd59ced6
......@@ -22,8 +22,8 @@
&nampar ! Generic parameters
!------------------------------------------------------------------------------
jpl = 5 ! number of ice categories
nlay_i = 2 ! number of ice layers
nlay_s = 2 ! number of snow layers
nlay_i = 10 ! number of ice layers
nlay_s = 5 ! number of snow layers
ln_virtual_itd = .false. ! virtual ITD mono-category parameterization (jpl=1 only)
! i.e. enhanced thermal conductivity & virtual thin ice melting
ln_icedyn = .true. ! ice dynamics (T) or not (F)
......@@ -62,11 +62,12 @@
rn_lf_relax = 1.e-5 ! relaxation time scale to reach static friction [s-1]
rn_lf_tensile = 0.05 ! isotropic tensile strength [0-0.5??]
cn_dir = './' ! root directory for the grounded icebergs mask data location
cn_dir = './' ! root directory for the grounded icebergs and landfast mask data location
!___________!________________!___________________!___________!_____________!________!___________!__________!__________!_______________!
! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename !
sn_icbmsk = 'NOT USED' , -12. , 'icb_mask', .false. , .true. , 'yearly' , '' , '' , ''
sn_icbmsk = 'NOT USED' , -12. , 'icb_mask' , .false. , .true. , 'yearly' , '' , '' , ''
sn_fastmsk = 'NOT USED' , -12. , 'fast_mask', .false. , .true. , 'yearly' , '' , '' , ''
/
!------------------------------------------------------------------------------
&namdyn_rdgrft ! Ice ridging/rafting
......@@ -82,7 +83,7 @@
ln_str_smooth = .true. ! spatial smoothing of the ice strength
! -- ice_rdgrft -- !
ln_distf_lin = .true. ! redistribution function of ridged ice: linear (Hibler 1980)
ln_distf_exp = .false. ! redistribution function of ridged ice: exponential => not coded yet
ln_distf_exp = .false. ! redistribution function of ridged ice: exponential (Lipscomb et al 2017)
rn_murdg = 3.0 ! e-folding scale of ridged ice (m**.5)
rn_csrdg = 0.5 ! fraction of shearing energy contributing to ridging
! -- ice_rdgrft_prep -- !
......@@ -132,7 +133,7 @@
!------------------------------------------------------------------------------
&namsbc ! Ice surface boundary conditions
!------------------------------------------------------------------------------
rn_cio = 5.0e-03 ! ice-ocean drag coefficient (-)
rn_Cd_io = 5.0e-03 ! ice-ocean drag coefficient (-)
nn_snwfra = 2 ! calculate the fraction of ice covered by snow (for zdf and albedo)
! = 0 fraction = 1 (if snow) or 0 (if no snow)
! = 1 fraction = 1-exp(-0.2*rhos*hsnw) [MetO formulation]
......@@ -156,8 +157,6 @@
ln_icedH = .true. ! activate ice thickness change from growing/melting (T) or not (F)
ln_icedA = .true. ! activate lateral melting param. (T) or not (F)
ln_icedO = .true. ! activate ice growth in open-water (T) or not (F)
ln_icedS = .true. ! activate brine drainage (T) or not (F)
!
ln_leadhfx = .true. ! heat in the leads is used to melt sea-ice before warming the ocean
/
!------------------------------------------------------------------------------
......@@ -197,32 +196,66 @@
!------------------------------------------------------------------------------
&namthd_sal ! Ice salinity
!------------------------------------------------------------------------------
nn_icesal = 2 ! ice salinity option
nn_icesal = 4 ! ice salinity option
! 1: constant ice salinity (S=rn_icesal)
! 2: varying salinity parameterization S(z,t)
! 3: prescribed salinity profile S(z) (Schwarzacher 1959)
rn_icesal = 4. ! (nn_icesal=1) ice salinity (g/kg)
rn_sal_gd = 5. ! (nn_icesal=2) restoring ice salinity, gravity drainage (g/kg)
rn_time_gd = 1.73e+6 ! (nn_icesal=2) restoring time scale, gravity drainage (s)
rn_sal_fl = 2. ! (nn_icesal=2) restoring ice salinity, flushing (g/kg)
rn_time_fl = 8.64e+5 ! (nn_icesal=2) restoring time scale, flushing (s)
rn_simax = 20. ! maximum tolerated ice salinity (g/kg)
! 4: Gravity Drainage and Flushing
ln_flushing = .true. ! activate ice salt flushing
ln_drainage = .true. ! activate ice salt gravity drainage
rn_simin = 0.1 ! minimum tolerated ice salinity (g/kg)
rn_sinew = 0.75 ! fraction of sss that is entrapped in new ice
! (must be set to 0.25 if nn_icesal=2)
! -- nn_icesal=1 -- !
rn_icesal = 4. ! ice salinity (g/kg)
! -- nn_icesal=2 -- !
rn_sal_gd = 5. ! restoring ice salinity, gravity drainage (g/kg)
rn_time_gd = 1.73e+6 ! restoring time scale, gravity drainage (s)
rn_sal_fl = 2. ! restoring ice salinity, flushing (g/kg)
rn_time_fl = 8.64e+5 ! restoring time scale, flushing (s)
! -- nn_icesal=4 -- !
rn_sal_himin = 0.1 ! min ice thickness for gravity drainage and flushing calculation
nn_liquidus = 2 ! formulation of liquidus (also used for outputs)
! 1 = linear liquidus
! 2 = Vancopenolle et al (2019) formulation
! 3 = Weast 71 (used in RJW2014)
nn_drainage = 50 ! maximum number of subcycles for gravity drainage (10 is usually enough)
nn_flushing = 5 ! maximum number of subcycles for flushing ( 1 is usually enough)
rn_flushrate = 0.3 ! rate of flushing (fraction of melt water used for flushing)
rn_vbrc = 0.05 ! critical brines volume above which flushing can occur
! ** drainage convection scheme **
nn_sal_scheme = 1 ! 1 = Rees Jones and Worster (2014) => RJW2014
! 2 = Griewank and Notz (2013) => GN2013
! 3 = Cox and Weeks (1988) => CW1988
! ** parameters for each scheme **
rn_alpha_RJW = 0.03 ! 1: Intensity of the Brine flow ==> minmax tested [0.03 ; 0.037]?
rn_Rc_RJW = 9.0 ! 1: Critical Rayleigh number ==> [2.9 ; 9]?
rn_alpha_GN = 5.1e-4 ! 2: Intensity of the Brine flow ==> [0.510e-3 ; 0.681e-3] from Griewank and Notz 2015
rn_Rc_GN = 10.0 ! 2: Critical Rayleigh number ==> [7.10 ; 3.23] from Griewank and Notz 2015
rn_alpha_CW = 4.8e-7 ! 3: Intensity of the Brine flow ==> [4e-7 ; 7e-7]?
!
ln_sal_chk = .FALSE. ! sanity checks for drainage and flushing
/
!------------------------------------------------------------------------------
&namthd_pnd ! Melt ponds
!------------------------------------------------------------------------------
ln_pnd = .true. ! activate melt ponds or not
ln_pnd_TOPO = .false. ! topographic melt ponds
ln_pnd_LEV = .true. ! level ice melt ponds
rn_apnd_min = 0.15 ! minimum meltwater fraction contributing to pond growth (TOPO and LEV)
rn_apnd_max = 0.85 ! maximum meltwater fraction contributing to pond growth (TOPO and LEV)
rn_pnd_flush= 0.1 ! pond flushing efficiency (tuning parameter) (LEV)
ln_pnd_CST = .false. ! constant melt ponds
rn_apnd = 0.2 ! prescribed pond fraction, at Tsu=0 degC
rn_hpnd = 0.05 ! prescribed pond depth, at Tsu=0 degC
ln_pnd_lids = .true. ! frozen lids on top of the ponds (only for ln_pnd_LEV)
ln_pnd_alb = .true. ! effect of melt ponds on ice albedo
!
rn_apnd_min = 0.15 ! minimum meltwater fraction contributing to pond growth (TOPO and LEV)
rn_apnd_max = 0.80 ! maximum meltwater fraction contributing to pond growth (TOPO and LEV)
!
ln_pnd_TOPO = .false. ! topographic melt ponds
ln_pnd_rain = .true. ! rain added to melt ponds
ln_pnd_LEV = .true. ! level ice melt ponds
rn_pnd_flush = 0.1 ! pond flushing efficiency (tuning parameter)
ln_pnd_lids = .true. ! frozen lids on top of the ponds
ln_pnd_CST = .false. ! constant melt ponds
rn_apnd = 0.2 ! prescribed pond fraction, at Tsu=0 degC
rn_hpnd = 0.05 ! prescribed pond depth, at Tsu=0 degC
!
ln_pnd_alb = .true. ! effect of melt ponds on ice albedo
rn_pnd_hl_min = 0.005 ! pond lid thickness below which full pond area used in albedo calculation
rn_pnd_hl_max = 0.015 ! pond lid thickness above which ponds disappear from albedo calculation
/
!------------------------------------------------------------------------------
&namini ! Ice initialization
......@@ -232,10 +265,10 @@
! 1 = Initialise sea ice from single category netcdf file
! 2 = Initialise sea ice from multi category restart file
rn_thres_sst = 2.0 ! max temp. above Tfreeze with initial ice = (sst - tfreeze)
rn_hti_ini_n = 3.0 ! initial ice thickness (m), North
rn_hti_ini_n = 2.0 ! initial ice thickness (m), North
rn_hti_ini_s = 1.0 ! " " South
rn_hts_ini_n = 0.3 ! initial snow thickness (m), North
rn_hts_ini_s = 0.3 ! " " South
rn_hts_ini_n = 0.2 ! initial snow thickness (m), North
rn_hts_ini_s = 0.2 ! " " South
rn_ati_ini_n = 0.9 ! initial ice concentration (-), North
rn_ati_ini_s = 0.9 ! " " South
rn_smi_ini_n = 6.3 ! initial ice salinity (g/kg), North
......@@ -272,9 +305,10 @@
! ! ! obs range (cloud-sky)
rn_alb_sdry = 0.85 ! dry snow albedo : 0.85 -- 0.87
rn_alb_smlt = 0.75 ! melting snow albedo : 0.72 -- 0.82
rn_alb_idry = 0.60 ! dry ice albedo : 0.54 -- 0.65
rn_alb_imlt = 0.50 ! bare puddled ice albedo : 0.49 -- 0.58
rn_alb_dpnd = 0.27 ! ponded ice albedo : 0.10 -- 0.30
rn_alb_idry = 0.64 ! dry ice albedo : 0.54 -- 0.65
rn_alb_imlt = 0.53 ! bare puddled ice albedo : 0.49 -- 0.58
rn_alb_dpnd = 0.18 ! dark pond albedo : 0.12 -- 0.25
rn_alb_lpnd = 0.30 ! lidded pond albedo : 0.2 -- 0.4
rn_alb_hpiv = 1.00 ! pivotal ice thickness in m (above which albedo is constant)
/
!------------------------------------------------------------------------------
......@@ -284,7 +318,6 @@
! ! rate of ice spuriously gained/lost at each time step => rn_icechk=1 <=> 1.e-6 m/hour
rn_icechk_cel = 1. ! check at each gridcell (1.e-06m/h)=> stops the code if violated (and writes a file)
rn_icechk_glo = 1.e-04 ! check over the entire ice cover (1.e-10m/h)=> only prints warnings
ln_icediahsb = .false. ! output the heat, mass & salt budgets (T) or not (F)
ln_icectl = .false. ! ice points output for debug (T or F)
iiceprt = 10 ! i-index for debug
jiceprt = 10 ! j-index for debug
......