diff --git a/EXP_Apr19/context_nemo.xml b/EXP_Apr19/context_nemo.xml
new file mode 100755
index 0000000000000000000000000000000000000000..4ab125fed1f6cf8e790a1832dc7359ca6165b2f9
--- /dev/null
+++ b/EXP_Apr19/context_nemo.xml
@@ -0,0 +1,120 @@
+<!--
+ ==============================================================================================
+ NEMO context
+==============================================================================================
+-->
+<context id="nemo">
+<!-- $id$ -->
+<!-- Fields definition -->
+ <field_definition src="./field_def_nemo-opa.xml"/> <!-- NEMO ocean dynamics -->
+
+<!-- Files definition -->
+ <file_definition src="./file_def_nemo.xml"/> <!-- NEMO ocean dynamics -->
+ <!--
+============================================================================================================
+= grid definition = = DO NOT CHANGE =
+============================================================================================================
+ -->
+
+ <axis_definition>
+ <axis id="deptht" long_name="Vertical T levels" unit="m" positive="down" />
+ <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="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_28C" long_name="28 degC isotherm" unit="degC" />
+ </axis_definition>
+
+ <domain_definition src="./domain_def_nemo.xml"/>
+
+ <grid_definition>
+
+ <!-- -->
+ <grid id="grid_T_2D" >
+ <domain id="grid_T" />
+ </grid>
+ <!-- -->
+ <grid id="grid_T_3D_ncatice" >
+ <domain id="grid_T" />
+ <axis id="ncatice" />
+ </grid>
+ <!-- -->
+ <grid id="grid_T_3D" >
+ <domain id="grid_T" />
+ <axis id="deptht" />
+ </grid>
+ <!-- -->
+ <grid id="grid_U_2D" >
+ <domain id="grid_U" />
+ </grid>
+ <!-- -->
+ <grid id="grid_U_3D" >
+ <domain id="grid_U" />
+ <axis id="depthu" />
+ </grid>
+ <!-- -->
+ <grid id="grid_V_2D" >
+ <domain id="grid_V" />
+ </grid>
+ <!-- -->
+ <grid id="grid_V_3D" >
+ <domain id="grid_V" />
+ <axis id="depthv" />
+ </grid>
+ <!-- -->
+ <grid id="grid_W_2D" >
+ <domain id="grid_W" />
+ </grid>
+ <!-- -->
+ <grid id="grid_W_3D" >
+ <domain id="grid_W" />
+ <axis id="depthw" />
+ </grid>
+ <!-- -->
+ <grid id="grid_1point" >
+ <domain domain_ref="1point"/>
+ </grid>
+ <!-- -->
+ <grid id="grid_T_nfloat" >
+ <domain id="grid_T" />
+ <axis id="nfloat" />
+ </grid>
+ <!-- -->
+ <grid id="grid_EqT" >
+ <domain domain_ref="EqT" />
+ </grid>
+ <!-- -->
+ <grid id="grid_znl_T_2D">
+ <domain id="gznl" />
+ </grid>
+ <!-- -->
+ <grid id="grid_znl_T_3D">
+ <domain id="gznl" />
+ <axis id="deptht" />
+ </grid>
+ <!-- -->
+ <grid id="grid_znl_W_3D">
+ <domain id="gznl" />
+ <axis id="depthw" />
+ </grid>
+
+ <grid id="grid_ptr_T_2D">
+ <domain id="ptr" />
+ </grid>
+
+ <grid id="grid_ptr_T_3D">
+ <domain id="ptr" />
+ <axis id="deptht" />
+ </grid>
+
+ <grid id="grid_ptr_W_3D">
+ <domain id="ptr" />
+ <axis id="depthw" />
+ </grid>
+
+ </grid_definition>
+
+</context>
diff --git a/EXP_Apr19/domain_def_nemo.xml b/EXP_Apr19/domain_def_nemo.xml
new file mode 100755
index 0000000000000000000000000000000000000000..a8f71f897b9aeb83825acb937f58ed47099806ce
--- /dev/null
+++ b/EXP_Apr19/domain_def_nemo.xml
@@ -0,0 +1,258 @@
+
+ <domain_definition>
+ <domain_group id="grid_T">
+ <domain id="grid_T" long_name="grid T"/>
+ <!-- My zoom: example of hand defined zoom -->
+
+ <domain id="moorT" ibegin="39" jbegin="82" ni="1" nj="1" />
+
+ <domain id="sec01T" ibegin="33" jbegin="62" ni="28" nj="1" />
+ <domain id="sec02T" ibegin="34" jbegin="68" ni="27" nj="1" />
+ <domain id="sec03T" ibegin="29" jbegin="75" ni="18" nj="1" />
+ <domain id="sec04T" ibegin="29" jbegin="80" ni="16" nj="1" />
+ <domain id="sec05T" ibegin="25" jbegin="88" ni="12" nj="1" />
+ <domain id="sec06T" ibegin="25" jbegin="91" ni="7" nj="1" />
+ <domain id="sec07T" ibegin="25" jbegin="92" ni="3" nj="1" />
+ <domain id="sec08T" ibegin="57" jbegin="77" ni="4" nj="1" />
+ <domain id="sec09T" ibegin="61" jbegin="62" ni="1" nj="15" />
+
+ <domain id="sec00ST" ibegin="2" jbegin="2" ni="115" nj="1" />
+ <domain id="sec00NT" ibegin="2" jbegin="104" ni="115" nj="1" />
+ <domain id="sec00WT" ibegin="2" jbegin="2" ni="1" nj="103" />
+ <domain id="sec00ET" ibegin="116" jbegin="2" ni="1" nj="103" />
+
+ <domain id="1point" domain_ref="grid_T" >
+ <zoom_domain id="1point" ibegin="139" jbegin="119" ni="1" nj="1" />
+ </domain>
+ <!-- Eq section -->
+ <domain id="EqT" domain_ref="grid_T" >
+ <zoom_domain id="EqT" ibegin="0" jbegin="0000" ni="0000" nj="1" />
+ </domain>
+ <!-- TAO : see example above -->
+ <!-- 137e -->
+ <!-- <domain id="2n137eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5n137eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="8n137eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- 147e -\-> -->
+ <!-- <domain id="0n147eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="2n147eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5n147eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- 156e -\-> -->
+ <!-- <domain id="5s156eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="2s156eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="0n156eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="2n156eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5n156eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="8n156eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- 165e -\-> -->
+ <!-- <domain id="8s165eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5s165eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="2s165eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="0n165eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="2n165eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5n165eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="8n165eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- 180w -\-> -->
+ <!-- <domain id="8s180wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5s180wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="2s180wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="0n180wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="2n180wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5n180wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="8n180wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- 170w -\-> -->
+ <!-- <domain id="8s170wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5s170wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="2s170wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="0n170wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="2n170wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5n170wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="8n170wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- 155w -\-> -->
+ <!-- <domain id="8s155wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5s155wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="2s155wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="0n155wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="2n155wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5n155wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="8n155wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- 140w -\-> -->
+ <!-- <domain id="8s140wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5s140wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="2s140wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="0n140wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="2n140wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5n140wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="8n140wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- 125w -\-> -->
+ <!-- <domain id="8s125wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5s125wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="2s125wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="0n125wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="2n125wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5n125wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="8n125wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- 110w -\-> -->
+ <!-- <domain id="8s110wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5s110wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="2s110wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="0n110wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="2n110wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5n110wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="8n110wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- 95w -\-> -->
+ <!-- <domain id="8s95wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5s95wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="2s95wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="0n95wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="2n95wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5n95wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="8n95wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- RAMA -\-> -->
+ <!-- <\!-- 55e -\-> -->
+ <!-- <domain id="16s55eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="12s55eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="8s55eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="4s55eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="1.5s55eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="0n55eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="1.5n55eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="4n55eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- 65e -\-> -->
+ <!-- <domain id="15n65eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- 67e -\-> -->
+ <!-- <domain id="16s67eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="12s67eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="8s67eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="4s67eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="1.5s67eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="0n67eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="1.5n67eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="4n67eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="8n67eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- 80.5e -\-> -->
+ <!-- <domain id="16s80.5eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="12s80.5eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="8s80.5eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="4s80.5eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="1.5s80.5eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="0n80.5eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="1.5n80.5eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="4n80.5eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- 90e -\-> -->
+ <!-- <domain id="1.5s90eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="0n90eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="1.5n90eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="4n90eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="8n90eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="12n90eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="15n90eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- 95e -\-> -->
+ <!-- <domain id="16s95eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="12s95eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="8s95eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="5s95eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- PIRATA -\-> -->
+ <!-- <\!-- 38w-30w -\-> -->
+ <!-- <domain id="19s34wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="14s32wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="8s30wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="0n35wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="4n38wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="8n38wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="12n38wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="15n38wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="20n38wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- 23w -\-> -->
+ <!-- <domain id="0n23wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="4n23wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="12n23wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="21n23wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- 10w -\-> -->
+ <!-- <domain id="10s10wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="6s10wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <domain id="0n10wT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+ <!-- <\!-- 0e -\-> -->
+ <!-- <domain id="0n0eT" domain_ref="grid_T" ibegin="0000" jbegin="0000" ni="1" nj="1" /> -->
+
+ </domain_group>
+
+ <domain_group id="grid_U">
+ <domain id="grid_U" long_name="grid U"/>
+ <!-- My zoom: example of hand defined zoom -->
+ <domain id="moorU" ibegin="39" jbegin="82" ni="1" nj="1" />
+ <domain id="sec01U" ibegin="33" jbegin="62" ni="28" nj="1" />
+ <domain id="sec02U" ibegin="34" jbegin="68" ni="27" nj="1" />
+ <domain id="sec03U" ibegin="29" jbegin="75" ni="18" nj="1" />
+ <domain id="sec04U" ibegin="29" jbegin="80" ni="16" nj="1" />
+ <domain id="sec05U" ibegin="25" jbegin="88" ni="12" nj="1" />
+ <domain id="sec06U" ibegin="25" jbegin="91" ni="7" nj="1" />
+ <domain id="sec07U" ibegin="25" jbegin="92" ni="3" nj="1" />
+ <domain id="sec08U" ibegin="57" jbegin="77" ni="4" nj="1" />
+ <domain id="sec09U" ibegin="61" jbegin="62" ni="1" nj="15" />
+
+ <domain id="sec00SU" ibegin="2" jbegin="2" ni="115" nj="1" />
+ <domain id="sec00NU" ibegin="2" jbegin="104" ni="115" nj="1" />
+ <domain id="sec00WU" ibegin="2" jbegin="2" ni="1" nj="103" />
+ <domain id="sec00EU" ibegin="116" jbegin="2" ni="1" nj="103" />
+ <!-- Eq section -->
+ <domain id="EqU" ibegin="1" jbegin="0000" ni="0000" nj="1" />
+
+ </domain_group>
+
+ <domain_group id="grid_V">
+ <domain id="grid_V" long_name="grid V"/>
+ <!-- My zoom: example of hand defined zoom -->
+ <domain id="moorV" ibegin="39" jbegin="82" ni="1" nj="1" />
+ <domain id="sec01V" ibegin="33" jbegin="62" ni="28" nj="1" />
+ <domain id="sec02V" ibegin="34" jbegin="68" ni="27" nj="1" />
+ <domain id="sec03V" ibegin="29" jbegin="75" ni="18" nj="1" />
+ <domain id="sec04V" ibegin="29" jbegin="80" ni="16" nj="1" />
+ <domain id="sec05V" ibegin="25" jbegin="88" ni="12" nj="1" />
+ <domain id="sec06V" ibegin="25" jbegin="91" ni="7" nj="1" />
+ <domain id="sec07V" ibegin="25" jbegin="92" ni="3" nj="1" />
+ <domain id="sec08V" ibegin="57" jbegin="77" ni="4" nj="1" />
+ <domain id="sec09V" ibegin="61" jbegin="62" ni="1" nj="15" />
+
+ <domain id="sec00SV" ibegin="2" jbegin="2" ni="115" nj="1" />
+ <domain id="sec00NV" ibegin="2" jbegin="104" ni="115" nj="1" />
+ <domain id="sec00WV" ibegin="2" jbegin="2" ni="1" nj="103" />
+ <domain id="sec00EV" ibegin="116" jbegin="2" ni="1" nj="103" />
+ </domain_group>
+
+ <domain_group id="grid_W">
+ <domain id="grid_W" long_name="grid W"/>
+ <!-- My zoom: example of hand defined zoom -->
+ <domain id="moorW" ibegin="39" jbegin="82" ni="1" nj="1" />
+ <domain id="sec01W" ibegin="33" jbegin="62" ni="28" nj="1" />
+ <domain id="sec02W" ibegin="34" jbegin="68" ni="27" nj="1" />
+ <domain id="sec03W" ibegin="29" jbegin="75" ni="18" nj="1" />
+ <domain id="sec04W" ibegin="29" jbegin="80" ni="16" nj="1" />
+ <domain id="sec05W" ibegin="25" jbegin="88" ni="12" nj="1" />
+ <domain id="sec06W" ibegin="25" jbegin="91" ni="7" nj="1" />
+ <domain id="sec07W" ibegin="25" jbegin="92" ni="3" nj="1" />
+ <domain id="sec08W" ibegin="57" jbegin="77" ni="4" nj="1" />
+ <domain id="sec09W" ibegin="61" jbegin="62" ni="1" nj="15" />
+
+ <domain id="sec00SW" ibegin="2" jbegin="2" ni="115" nj="1" />
+ <domain id="sec00NW" ibegin="2" jbegin="104" ni="115" nj="1" />
+ <domain id="sec00WW" ibegin="2" jbegin="2" ni="1" nj="103" />
+ <domain id="sec00EW" ibegin="116" jbegin="2" ni="1" nj="103" />
+ <!-- Eq section -->
+ <domain id="EqW" ibegin="1" jbegin="0000" ni="0000" nj="1" />
+ </domain_group>
+
+ <domain_group id="scalarpoint">
+ <domain id="scalarpoint" long_name="scalar"/>
+ </domain_group>
+
+ <domain_group id="gznl">
+ <domain id="gznl" long_name="gznl"/>
+ <domain id="ptr" domain_ref="gznl" >
+ <zoom_domain id="ptr" ibegin="0000" jbegin="0" ni="1" nj="0000" />
+ </domain>
+ </domain_group>
+
+
+ </domain_definition>
+
diff --git a/EXP_Apr19/field_def_nemo-lim.xml b/EXP_Apr19/field_def_nemo-lim.xml
new file mode 100755
index 0000000000000000000000000000000000000000..556fe89df6b7955c3eb682ea7d7a7d195592dabf
--- /dev/null
+++ b/EXP_Apr19/field_def_nemo-lim.xml
@@ -0,0 +1,182 @@
+<?xml version="1.0"?>
+ <!-- $id$ -->
+
+ <!--
+============================================================================================================
+= definition of all existing variables =
+= DO NOT CHANGE =
+============================================================================================================
+ -->
+ <field_definition level="1" prec="4" operation="average" enabled=".TRUE." default_value="1.e20" > <!-- time step automaticaly defined -->
+
+ <field_group id="SBC" grid_ref="grid_T_2D" > <!-- time step automaticaly defined based on nn_fsbc -->
+
+ <!-- LIM2(only) fields -->
+ <field id="qsr_ai_cea" long_name="Air-Ice downward solar heat flux (cell average)" standard_name="surface_downwelling_shortwave_flux_in_air" unit="W/m2" />
+ <field id="qns_ai_cea" long_name="Air-Ice downward non-solar heat flux (cell average)" unit="W/m2" />
+ <field id="qla_ai_cea" long_name="Air-Ice downward Latent heat flux (cell average)" standard_name="surface_downward_latent_heat_flux" unit="W/m2" />
+ <field id="qsr_io_cea" long_name="Ice-Oce downward solar heat flux (cell average)" standard_name="net_downward_shortwave_flux_at_sea_water_surface" unit="W/m2" />
+ <field id="qns_io_cea" long_name="Ice-Oce downward non-solar heat flux (cell average)" unit="W/m2" />
+ <field id="iceprod_cea" long_name="Ice production (cell average)" unit="m/s" />
+ <field id="iiceconc" long_name="Ice concentration" standard_name="sea_ice_area_fraction" unit="1" />
+ <field id="ice_pres" long_name="Ice presence" unit="" />
+ <field id="ist_cea" long_name="Ice surface temperature (cell average)" standard_name="surface_temperature" unit="degC" />
+ <field id="ist_ipa" long_name="Ice surface temperature (ice presence average)" standard_name="surface_temperature" unit="degC" />
+ <field id="u_imasstr" long_name="Sea-ice mass transport along i-axis" standard_name="sea_ice_x_transport" unit="kg/s" />
+ <field id="v_imasstr" long_name="Sea-ice mass transport along j-axis" standard_name="sea_ice_y_transport" unit="kg/s" />
+
+ <!-- LIM3 fields -->
+ <field id="ice_cover" long_name="Ice fraction" standard_name="sea_ice_area_fraction" unit="1" />
+
+ <field id="snowthic_cea" long_name="Snow thickness (cell average)" standard_name="surface_snow_thickness" unit="m" />
+ <field id="icethic_cea" long_name="Ice thickness (cell average)" standard_name="sea_ice_thickness" unit="m" />
+
+ <field id="uice_ipa" long_name="Ice velocity along i-axis at I-point (ice presence average)" standard_name="sea_ice_x_velocity" unit="m/s" />
+ <field id="vice_ipa" long_name="Ice velocity along j-axis at I-point (ice presence average)" standard_name="sea_ice_y_velocity" unit="m/s" />
+ <field id="utau_ice" long_name="Wind stress along i-axis over the ice at i-point" standard_name="surface_downward_x_stress" unit="N/m2" />
+ <field id="vtau_ice" long_name="Wind stress along j-axis over the ice at i-point" standard_name="surface_downward_y_stress" unit="N/m2" />
+
+ <field id="iceconc" long_name="ice concentration" standard_name="sea_ice_area_fraction" unit="%" />
+ <field id="isst" long_name="sea surface temperature" standard_name="sea_surface_temperature" unit="degC" />
+ <field id="isss" long_name="sea surface salinity" standard_name="sea_surface_salinity" unit="1e-3" />
+ <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="qt_ice" long_name="total heat flux at ice surface: sum over categories" standard_name="surface_downward_heat_flux_in_air" unit="W/m2" />
+ <field id="qsr_ice" long_name="solar heat flux at ice surface: sum over categories" standard_name="surface_downwelling_shortwave_flux_in_air" unit="W/m2" />
+ <field id="qns_ice" long_name="non-solar heat flux at ice surface: sum over categories" unit="W/m2" />
+ <field id="qtr_ice" long_name="solar heat flux transmitted through ice: sum over categories" unit="W/m2" />
+ <field id="qemp_ice" long_name="Downward Heat Flux from E-P over ice" unit="W/m2" />
+ <field id="micesalt" long_name="Mean ice salinity" unit="1e-3" />
+ <field id="miceage" long_name="Mean ice age" unit="years" />
+ <field id="alb_ice" long_name="Mean albedo over sea ice" unit="" />
+ <field id="albedo" long_name="Mean albedo over sea ice and ocean" unit="" />
+
+ <field id="iceage_cat" long_name="Ice age for categories" unit="days" grid_ref="grid_T_3D_ncatice" />
+ <field id="iceconc_cat" long_name="Ice concentration for categories" unit="%" grid_ref="grid_T_3D_ncatice" />
+ <field id="icethic_cat" long_name="Ice thickness for categories" unit="m" grid_ref="grid_T_3D_ncatice" />
+ <field id="snowthic_cat" long_name="Snow thicknessi for categories" unit="m" grid_ref="grid_T_3D_ncatice" />
+ <field id="salinity_cat" long_name="Sea-Ice Bulk salinity for categories" unit="g/kg" grid_ref="grid_T_3D_ncatice" />
+ <field id="brinevol_cat" long_name="Brine volume for categories" unit="%" grid_ref="grid_T_3D_ncatice" />
+ <field id="icetemp_cat" long_name="Ice temperature for categories" unit="degC" grid_ref="grid_T_3D_ncatice" />
+ <field id="snwtemp_cat" long_name="Snow temperature for categories" unit="degC" grid_ref="grid_T_3D_ncatice" />
+
+ <field id="micet" long_name="Mean ice temperature" unit="degC" />
+ <field id="icehc" long_name="ice total heat content" unit="10^9J" />
+ <field id="isnowhc" long_name="snow total heat content" unit="10^9J" />
+ <field id="icest" long_name="ice surface temperature" unit="degC" />
+ <field id="ibrinv" long_name="brine volume" unit="%" />
+ <field id="icecolf" long_name="frazil ice collection thickness" unit="m" />
+ <field id="icestr" long_name="ice strength" unit="N/m" />
+ <field id="icevel" long_name="ice velocity" unit="m/s" />
+ <field id="idive" long_name="divergence" unit="1e-8s-1" />
+ <field id="ishear" long_name="shear" unit="1e-8s-1" />
+ <field id="icevolu" long_name="ice volume" unit="m" />
+ <field id="snowvol" long_name="snow volume" unit="m" />
+ <field id="tau_icebfr" long_name="ice friction on ocean bottom for landfast ice" unit="N/2" />
+
+ <field id="icetrp" long_name="ice volume transport" unit="m/day" />
+ <field id="snwtrp" long_name="snw volume transport" unit="m/day" />
+ <field id="saltrp" long_name="salt content transport" unit="1e-3*kg/m2/day" />
+ <field id="deitrp" long_name="advected ice enthalpy" unit="W/m2" />
+ <field id="destrp" long_name="advected snw enthalpy" unit="W/m2" />
+
+ <field id="sfxbri" long_name="brine salt flux" unit="1e-3*kg/m2/day" />
+ <field id="sfxdyn" long_name="salt flux from ridging rafting" unit="1e-3*kg/m2/day" />
+ <field id="sfxres" long_name="salt flux from lipupdate (resultant)" unit="1e-3*kg/m2/day" />
+ <field id="sfxbog" long_name="salt flux from bot growth" unit="1e-3*kg/m2/day" />
+ <field id="sfxbom" long_name="salt flux from bot melt" unit="1e-3*kg/m2/day" />
+ <field id="sfxsum" long_name="salt flux from surf melt" unit="1e-3*kg/m2/day" />
+ <field id="sfxlam" long_name="salt flux from lateral melt" unit="1e-3*kg/m2/day" />
+ <field id="sfxsni" long_name="salt flux from snow-ice formation" unit="1e-3*kg/m2/day" />
+ <field id="sfxopw" long_name="salt flux from open water ice formation" unit="1e-3*kg/m2/day" />
+ <field id="sfxsub" long_name="salt flux from sublimation" unit="1e-3*kg/m2/day" />
+ <field id="sfx" long_name="salt flux total" unit="1e-3*kg/m2/day" />
+
+ <field id="vfxbog" long_name="daily bottom thermo ice prod." unit="m/day" />
+ <field id="vfxdyn" long_name="daily dynamic ice prod." unit="m/day" />
+ <field id="vfxopw" long_name="daily lateral thermo ice prod." unit="m/day" />
+ <field id="vfxsni" long_name="daily snowice ice prod." unit="m/day" />
+ <field id="vfxsum" long_name="surface melt" unit="m/day" />
+ <field id="vfxlam" long_name="lateral melt" unit="m/day" />
+ <field id="vfxbom" long_name="bottom melt" unit="m/day" />
+ <field id="vfxres" long_name="daily resultant ice prod./melting from limupdate" unit="m/day" />
+ <field id="vfxice" long_name="ice melt/growth" unit="m/day" />
+ <field id="vfxsnw" long_name="snw melt/growth" unit="m/day" />
+ <field id="vfxsub" long_name="snw sublimation" unit="m/day" />
+ <field id="vfxsub_err" long_name="excess of snw sublimation sent to ocean" unit="m/day" />
+ <field id="vfxspr" long_name="snw precipitation on ice" unit="m/day" />
+ <field id="vfxthin" long_name="daily thermo ice prod. for thin ice(20cm) + open water" unit="m/day" />
+
+ <field id="afxtot" long_name="area tendency (total)" unit="day-1" />
+ <field id="afxdyn" long_name="area tendency (dynamics)" unit="day-1" />
+ <field id="afxthd" long_name="area tendency (thermo)" unit="day-1" />
+
+ <field id="hfxsum" long_name="heat fluxes causing surface ice melt" unit="W/m2" />
+ <field id="hfxbom" long_name="heat fluxes causing bottom ice melt" unit="W/m2" />
+ <field id="hfxbog" long_name="heat fluxes causing bottom ice growth" unit="W/m2" />
+ <field id="hfxdif" long_name="heat fluxes causing ice temperature change" unit="W/m2" />
+ <field id="hfxopw" long_name="heat fluxes causing open water ice formation" unit="W/m2" />
+ <field id="hfxsnw" long_name="heat fluxes causing snow melt" unit="W/m2" />
+ <field id="hfxerr" long_name="heat fluxes error after heat diffusion" unit="W/m2" />
+ <field id="hfxerr_rem" long_name="heat fluxes error after remapping" unit="W/m2" />
+ <field id="hfxout" long_name="total heat fluxes received by the ocean" unit="W/m2" />
+ <field id="hfxin" long_name="total heat fluxes at the ice/ocean surface" unit="W/m2" />
+
+ <!-- heat flux associated with mass exchange -->
+ <field id="hfxthd" long_name="heat fluxes from ice-ocean mass exchange during thermo" unit="W/m2" />
+ <field id="hfxdyn" long_name="heat fluxes from ice-ocean mass exchange during dynamic" unit="W/m2" />
+ <field id="hfxres" long_name="heat fluxes from ice-ocean mass exchange during resultant" unit="W/m2" />
+ <field id="hfxsub" long_name="heat fluxes from ice-atm. mass exchange during sublimation" unit="W/m2" />
+ <field id="hfxspr" long_name="heat fluxes from ice-atm. mass exchange during snow precip" unit="W/m2" />
+
+ <!-- diags -->
+ <field id="hfxdhc" long_name="Heat content variation in snow and ice" unit="W/m2" />
+ <field id="hfxtur" long_name="turbulent heat flux at the ice base" unit="W/m2" />
+
+ <!-- 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_group>
+
+ <!-- LIM3 scalar variables -->
+ <field_group id="SBC_scalar" grid_ref="grid_T_2D" >
+ <!-- available with ln_limdiaout -->
+ <field id="ibgfrcvoltop" long_name="global mean ice/snow forcing at interface ice/snow-atm (volume equivalent ocean volume)" unit="km3" />
+ <field id="ibgfrcvolbot" long_name="global mean ice/snow forcing at interface ice/snow-ocean (volume equivalent ocean volume)" unit="km3" />
+ <field id="ibgfrctemtop" long_name="global mean heat on top of ice/snw/ocean-atm " unit="1e20J" />
+ <field id="ibgfrctembot" long_name="global mean heat below ice (on top of ocean) " unit="1e20J" />
+ <field id="ibgfrcsal" long_name="global mean ice/snow forcing (salt equivalent ocean volume)" unit="pss*km3" />
+ <field id="ibgfrchfxtop" long_name="global mean heat flux on top of ice/snw/ocean-atm " unit="W/m2" />
+ <field id="ibgfrchfxbot" long_name="global mean heat flux below ice (on top of ocean) " unit="W/m2" />
+
+ <field id="ibgvolume" long_name="drift in ice/snow volume (equivalent ocean volume)" unit="km3" />
+ <field id="ibgsaltco" long_name="drift in ice salt content (equivalent ocean volume)" unit="pss*km3" />
+ <field id="ibgheatco" long_name="drift in ice/snow heat content" unit="1e20J" />
+ <field id="ibgheatfx" long_name="drift in ice/snow heat flux" unit="W/m2" />
+
+ <field id="ibgvol_tot" long_name="global mean ice volume" unit="km3" />
+ <field id="sbgvol_tot" long_name="global mean snow volume" unit="km3" />
+ <field id="ibgarea_tot" long_name="global mean ice area" unit="km2" />
+ <field id="ibgsalt_tot" long_name="global mean ice salt content" unit="1e-3*km3" />
+ <field id="ibgheat_tot" long_name="global mean ice heat content" unit="1e20J" />
+ <field id="sbgheat_tot" long_name="global mean snow heat content" unit="1e20J" />
+ </field_group>
+
+ <!--
+============================================================================================================
+ -->
+ <!-- output variables for my configuration (example) -->
+
+ <field_group id="myvarICE" >
+ <field field_ref="icethic_cea" name="sithic" long_name="sea_ice_thickness" />
+ <field field_ref="icevolu" name="sivolu" />
+ <field field_ref="iceconc" name="siconc" />
+ </field_group>
+
+ </field_definition>
diff --git a/EXP_Apr19/field_def_nemo-opa.xml b/EXP_Apr19/field_def_nemo-opa.xml
new file mode 100755
index 0000000000000000000000000000000000000000..65a6e780c4ff417ff5b2ee3348971feaf80e18bb
--- /dev/null
+++ b/EXP_Apr19/field_def_nemo-opa.xml
@@ -0,0 +1,1107 @@
+<?xml version="1.0"?>
+ <!-- $id$ -->
+
+ <!--
+============================================================================================================
+= definition of all existing variables =
+= DO NOT CHANGE =
+============================================================================================================
+ -->
+ <field_definition level="1" prec="4" operation="average" enabled=".TRUE." default_value="1.e20" > <!-- time step automaticaly defined -->
+
+ <!--
+============================================================================================================
+ Physical ocean model variables
+============================================================================================================
+ -->
+
+ <!-- T grid -->
+
+ <field_group id="grid_T" grid_ref="grid_T_2D" >
+ <field id="e3t" long_name="T-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_T_3D"/>
+ <field id="e3t_surf" long_name="T-cell thickness" field_ref="e3t" standard_name="cell_thickness" unit="m" grid_ref="grid_T_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="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 >
+
+ <!-- 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="sea surface temperature" standard_name="sea_surface_temperature" unit="degC" />
+ <field id="sst2" long_name="square of sea surface temperature" standard_name="square_of_sea_surface_temperature" unit="degC2" > sst * sst </field >
+ <field id="sstmax" long_name="max of sea surface temperature" field_ref="sst" operation="maximum" />
+ <field id="sstmin" long_name="min of sea surface temperature" field_ref="sst" operation="minimum" />
+ <field id="sstgrad" long_name="module of sst gradient" unit="degC/m" />
+ <field id="sstgrad2" long_name="square of module of sst gradient" unit="degC2/m2" />
+ <field id="sbt" long_name="sea bottom temperature" unit="degC" />
+ <field id="tosmint" long_name="vertical integral of temperature times density" standard_name="integral_wrt_depth_of_product_of_density_and_potential_temperature" unit="(kg m2) degree_C" />
+ <field id="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" />
+ <field id="somint" long_name="vertical integral of salinity times density" standard_name="integral_wrt_depth_of_product_of_density_and_salinity" unit="(kg m2) x (1e-3)" />
+
+ <field id="taubot" long_name="bottom stress module" unit="N/m2" />
+
+ <field id="ssh" long_name="sea surface height" standard_name="sea_surface_height_above_geoid" unit="m" />
+ <field id="ssh2" long_name="square of sea surface height" standard_name="square_of_sea_surface_height_above_geoid" unit="m2" > ssh * ssh </field >
+ <field id="wetdep" long_name="wet depth" standard_name="wet_depth" unit="m" />
+ <field id="sshmax" long_name="max of sea surface height" field_ref="ssh" operation="maximum" />
+
+ <field id="mldkz5" long_name="Turbocline depth (Kz = 5e-4)" standard_name="ocean_mixed_layer_thickness_defined_by_vertical_tracer_diffusivity" unit="m" />
+ <field id="mldr10_1" long_name="Mixed Layer Depth (dsigma = 0.01 wrt 10m)" standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta" unit="m" />
+ <field id="mldr10_1max" long_name="Max of Mixed Layer Depth (dsigma = 0.01 wrt 10m)" field_ref="mldr10_1" operation="maximum" />
+ <field id="mldr10_1min" long_name="Min of Mixed Layer Depth (dsigma = 0.01 wrt 10m)" field_ref="mldr10_1" operation="minimum" />
+ <field id="heatc" long_name="Heat content vertically integrated" standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content" unit="J/m2" />
+ <field id="saltc" long_name="Salt content vertically integrated" unit="1e-3*kg/m2" />
+
+ <!-- EOS -->
+ <field id="alpha" long_name="thermal expansion" unit="degC-1" grid_ref="grid_T_3D" />
+ <field id="beta" long_name="haline contraction" unit="1e3" grid_ref="grid_T_3D" />
+ <field id="bn2" long_name="squared Brunt-Vaisala frequency" unit="s-1" grid_ref="grid_T_3D" />
+ <field id="rhop" long_name="potential density (sigma0)" standard_name="sea_water_sigma_theta" unit="kg/m3" grid_ref="grid_T_3D" />
+
+ <!-- Energy - horizontal divergence -->
+ <field id="eken" long_name="kinetic energy" standard_name="specific_kinetic_energy_of_sea_water" unit="m2/s2" grid_ref="grid_T_3D" />
+ <field id="hdiv" long_name="horizontal divergence" unit="s-1" grid_ref="grid_T_3D" />
+
+ <!-- variables available with MLE -->
+ <field id="Lf_NHpf" long_name="MLE: Lf = N H / f" unit="m" />
+
+ <!-- next variables available with key_diahth -->
+ <field id="mlddzt" long_name="Thermocline Depth (depth of max dT/dz)" standard_name="depth_at_maximum_upward_derivative_of_sea_water_potential_temperature" unit="m" />
+ <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="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" />
+
+ <!-- variables available with diaar5 -->
+ <field id="botpres" long_name="Sea Water Pressure at Sea Floor" standard_name="sea_water_pressure_at_sea_floor" unit="dbar" />
+ <field id="sshdyn" long_name="dynamic sea surface height" standard_name="dynamic_sea_surface_height_above_geoid" unit="m" />
+ <field id="sshdyn2" long_name="square of dynamic sea surface height" standard_name="dynamic_sea_surface_height_above_geoid_squared" unit="m2" > sshdyn * sshdyn </field>
+ <field id="tnpeo" long_name="Tendency of ocean potential energy content" unit="W/m2" />
+
+ <!-- variables available ln_linssh=.FALSE. -->
+ <field id="tpt_dep" long_name="T-point depth" standard_name="depth_below_geoid" unit="m" grid_ref="grid_T_3D" />
+ <field id="e3tdef" long_name="T-cell thickness deformation" unit="%" grid_ref="grid_T_3D" />
+ </field_group>
+
+ <!-- Tides -->
+
+ <field_group id="Tides_T" grid_ref="grid_T_2D" operation="once" >
+ <!-- tidal composante -->
+ <field id="M2x" long_name="M2 Elevation harmonic real part " unit="m" />
+ <field id="M2y" long_name="M2 Elevation harmonic imaginary part" unit="m" />
+ <field id="S2x" long_name="S2 Elevation harmonic real part " unit="m" />
+ <field id="S2y" long_name="S2 Elevation harmonic imaginary part" unit="m" />
+ <field id="N2x" long_name="N2 Elevation harmonic real part " unit="m" />
+ <field id="N2y" long_name="N2 Elevation harmonic imaginary part" unit="m" />
+ <field id="K1x" long_name="K1 Elevation harmonic real part " unit="m" />
+ <field id="K1y" long_name="K1 Elevation harmonic imaginary part" unit="m" />
+ <field id="O1x" long_name="O1 Elevation harmonic real part " unit="m" />
+ <field id="O1y" long_name="O1 Elevation harmonic imaginary part" unit="m" />
+ <field id="Q1x" long_name="Q1 Elevation harmonic real part " unit="m" />
+ <field id="Q1y" long_name="Q1 Elevation harmonic imaginary part" unit="m" />
+ <field id="M4x" long_name="M4 Elevation harmonic real part " unit="m" />
+ <field id="M4y" long_name="M4 Elevation harmonic imaginary part" unit="m" />
+ <field id="K2x" long_name="K2 Elevation harmonic real part " unit="m" />
+ <field id="K2y" long_name="K2 Elevation harmonic imaginary part" unit="m" />
+ <field id="P1x" long_name="P1 Elevation harmonic real part " unit="m" />
+ <field id="P1y" long_name="P1 Elevation harmonic imaginary part" unit="m" />
+ <field id="Mfx" long_name="Mf Elevation harmonic real part " unit="m" />
+ <field id="Mfy" long_name="Mf Elevation harmonic imaginary part" unit="m" />
+ <field id="Mmx" long_name="Mm Elevation harmonic real part " unit="m" />
+ <field id="Mmy" long_name="Mm Elevation harmonic imaginary part" unit="m" />
+
+ <field id="M2amp" long_name="M2 Elevation harmonic Amplitude" unit="m" />
+ <field id="M2phase" long_name="M2 Elevation harmonic Phase" unit="degree" />
+
+ <field id="S2amp" long_name="S2 Elevation harmonic Amplitude" unit="m" />
+ <field id="S2phase" long_name="S2 Elevation harmonic Phase" unit="degree" />
+
+ <field id="N2amp" long_name="N2 Elevation harmonic Amplitude" unit="m" />
+ <field id="N2phase" long_name="N2 Elevation harmonic Phase" unit="degree" />
+
+ <field id="K1amp" long_name="K1 Elevation harmonic Amplitude" unit="m" />
+ <field id="K1phase" long_name="K1 Elevation harmonic Phase" unit="degree" />
+
+ <field id="O1amp" long_name="O1 Elevation harmonic Amplitude" unit="m" />
+ <field id="O1phase" long_name="O1 Elevation harmonic Phase" unit="degree" />
+
+ <field id="Q1amp" long_name="Q1 Elevation harmonic Amplitude" unit="m" />
+ <field id="Q1phase" long_name="Q1 Elevation harmonic Phase" unit="degree" />
+
+ <field id="M4amp" long_name="M4 Elevation harmonic Amplitude" unit="m" />
+ <field id="M4phase" long_name="M4 Elevation harmonic Phase" unit="degree" />
+
+ <field id="MS4amp" long_name="MS4 Elevation harmonic Amplitude" unit="m" />
+ <field id="MS4phase" long_name="MS4 Elevation harmonic Phase" unit="degree" />
+
+ <field id="MN4amp" long_name="MN4 Elevation harmonic Amplitude" unit="m" />
+ <field id="MN4phase" long_name="MN4 Elevation harmonic Phase" unit="degree" />
+
+ <field id="K2amp" long_name="K2 Elevation harmonic Amplitude" unit="m" />
+ <field id="K2phase" long_name="K2 Elevation harmonic Phase" unit="degree" />
+
+ <field id="P1amp" long_name="P1 Elevation harmonic Amplitude" unit="m" />
+ <field id="P1phase" long_name="P1 Elevation harmonic Phase" unit="degree" />
+
+ <field id="Mfamp" long_name="Mf Elevation harmonic Amplitude" unit="m" />
+ <field id="Mfphase" long_name="Mf Elevation harmonic Phase" unit="degree" />
+
+ <field id="Mmamp" long_name="Mm Elevation harmonic Amplitude" unit="m" />
+ <field id="Mmphase" long_name="Mm Elevation harmonic Phase" unit="degree" />
+
+ <field id="T2amp" long_name="T2 Elevation harmonic Amplitude" unit="m" />
+ <field id="T2phase" long_name="T2 Elevation harmonic Phase" unit="degree" />
+
+ <field id="L2amp" long_name="L2 Elevation harmonic Amplitude" unit="m" />
+ <field id="L2phase" long_name="L2 Elevation harmonic Phase" unit="degree" />
+
+ <field id="S1amp" long_name="S1 Elevation harmonic Amplitude" unit="m" />
+ <field id="S1phase" long_name="S1 Elevation harmonic Phase" unit="degree" />
+
+ <field id="2N2amp" long_name="2N2 Elevation harmonic Amplitude" unit="m" />
+ <field id="2N2phase" long_name="2N2 Elevation harmonic Phase" unit="degree" />
+
+ <field id="MU2amp" long_name="MU2 Elevation harmonic Amplitude" unit="m" />
+ <field id="MU2phase" long_name="MU2 Elevation harmonic Phase" unit="degree" />
+
+ <field id="NU2amp" long_name="NU2 Elevation harmonic Amplitude" unit="m" />
+ <field id="NU2phase" long_name="NU2 Elevation harmonic Phase" unit="degree" />
+
+ </field_group>
+
+ <field_group id="Tides_U" grid_ref="grid_U_2D" operation="once" >
+ <field id="M2x_u" long_name="M2 current barotrope along i-axis harmonic real part " unit="m/s" />
+ <field id="M2y_u" long_name="M2 current barotrope along i-axis harmonic imaginary part " unit="m/s" />
+ <field id="S2x_u" long_name="S2 current barotrope along i-axis harmonic real part " unit="m/s" />
+ <field id="S2y_u" long_name="S2 current barotrope along i-axis harmonic imaginary part " unit="m/s" />
+ <field id="N2x_u" long_name="N2 current barotrope along i-axis harmonic real part " unit="m/s" />
+ <field id="N2y_u" long_name="N2 current barotrope along i-axis harmonic imaginary part " unit="m/s" />
+ <field id="K1x_u" long_name="K1 current barotrope along i-axis harmonic real part " unit="m/s" />
+ <field id="K1y_u" long_name="K1 current barotrope along i-axis harmonic imaginary part " unit="m/s" />
+ <field id="O1x_u" long_name="O1 current barotrope along i-axis harmonic real part " unit="m/s" />
+ <field id="O1y_u" long_name="O1 current barotrope along i-axis harmonic imaginary part " unit="m/s" />
+ <field id="Q1x_u" long_name="Q1 current barotrope along i-axis harmonic real part " unit="m/s" />
+ <field id="Q1y_u" long_name="Q1 current barotrope along i-axis harmonic imaginary part " unit="m/s" />
+ <field id="M4x_u" long_name="M4 current barotrope along i-axis harmonic real part " unit="m/s" />
+ <field id="M4y_u" long_name="M4 current barotrope along i-axis harmonic imaginary part " unit="m/s" />
+ <field id="K2x_u" long_name="K2 current barotrope along i-axis harmonic real part " unit="m/s" />
+ <field id="K2y_u" long_name="K2 current barotrope along i-axis harmonic imaginary part " unit="m/s" />
+ <field id="P1x_u" long_name="P1 current barotrope along i-axis harmonic real part " unit="m/s" />
+ <field id="P1y_u" long_name="P1 current barotrope along i-axis harmonic imaginary part " unit="m/s" />
+ <field id="Mfx_u" long_name="Mf current barotrope along i-axis harmonic real part " unit="m/s" />
+ <field id="Mfy_u" long_name="Mf current barotrope along i-axis harmonic imaginary part " unit="m/s" />
+ <field id="Mmx_u" long_name="Mm current barotrope along i-axis harmonic real part " unit="m/s" />
+ <field id="Mmy_u" long_name="Mm current barotrope along i-axis harmonic imaginary part " unit="m/s" />
+
+ <field id="M2amp_u2D" long_name="M2 U barotropic harmonic Amplitude" unit="m/s" />
+ <field id="M2phase_u2D" long_name="M2 U barotropic harmonic Phase" unit="degree" />
+
+ <field id="S2amp_u2D" long_name="S2 U barotropic harmonic Amplitude" unit="m/s" />
+ <field id="S2phase_u2D" long_name="S2 U barotropic harmonic Phase" unit="degree" />
+
+ <field id="N2amp_u2D" long_name="N2 U barotropic harmonic Amplitude" unit="m/s" />
+ <field id="N2phase_u2D" long_name="N2 U barotropic harmonic Phase" unit="degree" />
+
+ <field id="K1amp_u2D" long_name="K1 U barotropic harmonic Amplitude" unit="m/s" />
+ <field id="K1phase_u2D" long_name="K1 U barotropic harmonic Phase" unit="degree" />
+
+ <field id="O1amp_u2D" long_name="O1 U barotropic harmonic Amplitude" unit="m/s" />
+ <field id="O1phase_u2D" long_name="O1 U barotropic harmonic Phase" unit="degree" />
+
+ <field id="Q1amp_u2D" long_name="Q1 U barotropic harmonic Amplitude" unit="m/s" />
+ <field id="Q1phase_u2D" long_name="Q1 U barotropic harmonic Phase" unit="degree" />
+
+ <field id="M4amp_u2D" long_name="M4 U barotropic harmonic Amplitude" unit="m/s" />
+ <field id="M4phase_u2D" long_name="M4 U barotropic harmonic Phase" unit="degree" />
+
+ <field id="MS4amp_u2D" long_name="MS4 U barotropic harmonic Amplitude" unit="m/s" />
+ <field id="MS4phase_u2D" long_name="MS4 U barotropic harmonic Phase" unit="degree" />
+
+ <field id="MN4amp_u2D" long_name="MN4 U barotropic harmonic Amplitude" unit="m/s" />
+ <field id="MN4phase_u2D" long_name="MN4 U barotropic harmonic Phase" unit="degree" />
+
+ <field id="K2amp_u2D" long_name="K2 U barotropic harmonic Amplitude" unit="m/s" />
+ <field id="K2phase_u2D" long_name="K2 U barotropic harmonic Phase" unit="degree" />
+
+ <field id="P1amp_u2D" long_name="P1 U barotropic harmonic Amplitude" unit="m/s" />
+ <field id="P1phase_u2D" long_name="P1 U barotropic harmonic Phase" unit="degree" />
+
+ <field id="Mfamp_u2D" long_name="Mf U barotropic harmonic Amplitude" unit="m/s" />
+ <field id="Mfphase_u2D" long_name="Mf U barotropic harmonic Phase" unit="degree" />
+
+ <field id="Mmamp_u2D" long_name="Mm U barotropic harmonic Amplitude" unit="m/s" />
+ <field id="Mmphase_u2D" long_name="Mm U barotropic harmonic Phase" unit="degree" />
+ <field id="T2amp_u2D" long_name="T2 U barotropic harmonic Amplitude" unit="m/s" />
+
+ <field id="T2phase_u2D" long_name="T2 U barotropic harmonic Phase" unit="degree" />
+
+ <field id="L2amp_u2D" long_name="L2 U barotropic harmonic Amplitude" unit="m/s" />
+ <field id="L2phase_u2D" long_name="L2 U barotropic harmonic Phase" unit="degree" />
+
+ <field id="S1amp_u2D" long_name="S1 U barotropic harmonic Amplitude" unit="m/s" />
+ <field id="S1phase_u2D" long_name="S1 U barotropic harmonic Phase" unit="degree" />
+
+ <field id="2N2amp_u2D" long_name="2N2 U barotropic harmonic Amplitude" unit="m/s" />
+ <field id="2N2phase_u2D" long_name="2N2 U barotropic harmonic Phase" unit="degree" />
+
+ <field id="MU2amp_u2D" long_name="MU2 U barotropic harmonic Amplitude" unit="m/s" />
+ <field id="MU2phase_u2D" long_name="MU2 U barotropic harmonic Phase" unit="degree" />
+
+ <field id="NU2amp_u2D" long_name="NU2 U barotropic harmonic Amplitude" unit="m/s" />
+ <field id="NU2phase_u2D" long_name="NU2 U barotropic harmonic Phase" unit="degree" />
+
+ </field_group>
+
+ <field_group id="Tides_V" grid_ref="grid_V_2D" operation="once" >
+ <field id="M2x_v" long_name="M2 current barotrope along j-axis harmonic real part " unit="m/s" />
+ <field id="M2y_v" long_name="M2 current barotrope along j-axis harmonic imaginary part " unit="m/s" />
+ <field id="S2x_v" long_name="S2 current barotrope along j-axis harmonic real part " unit="m/s" />
+ <field id="S2y_v" long_name="S2 current barotrope along j-axis harmonic imaginary part " unit="m/s" />
+ <field id="N2x_v" long_name="N2 current barotrope along j-axis harmonic real part " unit="m/s" />
+ <field id="N2y_v" long_name="N2 current barotrope along j-axis harmonic imaginary part " unit="m/s" />
+ <field id="K1x_v" long_name="K1 current barotrope along j-axis harmonic real part " unit="m/s" />
+ <field id="K1y_v" long_name="K1 current barotrope along j-axis harmonic imaginary part " unit="m/s" />
+ <field id="O1x_v" long_name="O1 current barotrope along j-axis harmonic real part " unit="m/s" />
+ <field id="O1y_v" long_name="O1 current barotrope along j-axis harmonic imaginary part " unit="m/s" />
+ <field id="Q1x_v" long_name="Q1 current barotrope along j-axis harmonic real part " unit="m/s" />
+ <field id="Q1y_v" long_name="Q1 current barotrope along j-axis harmonic imaginary part " unit="m/s" />
+ <field id="M4x_v" long_name="M4 current barotrope along j-axis harmonic real part " unit="m/s" />
+ <field id="M4y_v" long_name="M4 current barotrope along j-axis harmonic imaginary part " unit="m/s" />
+ <field id="K2x_v" long_name="K2 current barotrope along j-axis harmonic real part " unit="m/s" />
+ <field id="K2y_v" long_name="K2 current barotrope along j-axis harmonic imaginary part " unit="m/s" />
+ <field id="P1x_v" long_name="P1 current barotrope along j-axis harmonic real part " unit="m/s" />
+ <field id="P1y_v" long_name="P1 current barotrope along j-axis harmonic imaginary part " unit="m/s" />
+ <field id="Mfx_v" long_name="Mf current barotrope along j-axis harmonic real part " unit="m/s" />
+ <field id="Mfy_v" long_name="Mf current barotrope along j-axis harmonic imaginary part " unit="m/s" />
+ <field id="Mmx_v" long_name="Mm current barotrope along j-axis harmonic real part " unit="m/s" />
+ <field id="Mmy_v" long_name="Mm current barotrope along j-axis harmonic imaginary part " unit="m/s" />
+
+ <field id="M2amp_v2D" long_name="M2 V barotropic harmonic Amplitude" unit="m/s" />
+ <field id="M2phase_v2D" long_name="M2 V barotropic harmonic Phase" unit="degree" />
+
+ <field id="S2amp_v2D" long_name="S2 V barotropic harmonic Amplitude" unit="m/s" />
+ <field id="S2phase_v2D" long_name="S2 V barotropic harmonic Phase" unit="degree" />
+
+ <field id="N2amp_v2D" long_name="N2 V barotropic harmonic Amplitude" unit="m/s" />
+ <field id="N2phase_v2D" long_name="N2 V barotropic harmonic Phase" unit="degree" />
+
+ <field id="K1amp_v2D" long_name="K1 V barotropic harmonic Amplitude" unit="m/s" />
+ <field id="K1phase_v2D" long_name="K1 V barotropic harmonic Phase" unit="degree" />
+
+ <field id="O1amp_v2D" long_name="O1 V barotropic harmonic Amplitude" unit="m/s" />
+ <field id="O1phase_v2D" long_name="O1 V barotropic harmonic Phase" unit="degree" />
+
+ <field id="Q1amp_v2D" long_name="Q1 V barotropic harmonic Amplitude" unit="m/s" />
+ <field id="Q1phase_v2D" long_name="Q1 V barotropic harmonic Phase" unit="degree" />
+
+ <field id="M4amp_v2D" long_name="M4 V barotropic harmonic Amplitude" unit="m/s" />
+ <field id="M4phase_v2D" long_name="M4 V barotropic harmonic Phase" unit="degree" />
+
+ <field id="MS4amp_v2D" long_name="MS4 V barotropic harmonic Amplitude" unit="m/s" />
+ <field id="MS4phase_v2D" long_name="MS4 V barotropic harmonic Phase" unit="degree" />
+
+ <field id="MN4amp_v2D" long_name="MN4 V barotropic harmonic Amplitude" unit="m/s" />
+ <field id="MN4phase_v2D" long_name="MN4 V barotropic harmonic Phase" unit="degree" />
+
+ <field id="K2amp_v2D" long_name="K2 V barotropic harmonic Amplitude" unit="m/s" />
+ <field id="K2phase_v2D" long_name="K2 V barotropic harmonic Phase" unit="degree" />
+
+ <field id="P1amp_v2D" long_name="P1 V barotropic harmonic Amplitude" unit="m/s" />
+ <field id="P1phase_v2D" long_name="P1 V barotropic harmonic Phase" unit="degree" />
+
+ <field id="Mfamp_v2D" long_name="Mf V barotropic harmonic Amplitude" unit="m/s" />
+ <field id="Mfphase_v2D" long_name="Mf V barotropic harmonic Phase" unit="degree" />
+
+ <field id="Mmamp_v2D" long_name="Mm V barotropic harmonic Amplitude" unit="m/s" />
+ <field id="Mmphase_v2D" long_name="Mm V barotropic harmonic Phase" unit="degree" />
+
+ <field id="T2amp_v2D" long_name="T2 V barotropic harmonic Amplitude" unit="m/s" />
+ <field id="T2phase_v2D" long_name="T2 V barotropic harmonic Phase" unit="degree" />
+
+ <field id="L2amp_v2D" long_name="L2 V barotropic harmonic Amplitude" unit="m/s" />
+ <field id="L2phase_v2D" long_name="L2 V barotropic harmonic Phase" unit="degree" />
+
+ <field id="S1amp_v2D" long_name="S1 V barotropic harmonic Amplitude" unit="m/s" />
+ <field id="S1phase_v2D" long_name="S1 V barotropic harmonic Phase" unit="degree" />
+
+ <field id="2N2amp_v2D" long_name="2N2 V barotropic harmonic Amplitude" unit="m/s" />
+ <field id="2N2phase_v2D" long_name="2N2 V barotropic harmonic Phase" unit="degree" />
+
+ <field id="MU2amp_v2D" long_name="MU2 V barotropic harmonic Amplitude" unit="m/s" />
+ <field id="MU2phase_v2D" long_name="MU2 V barotropic harmonic Phase" unit="degree" />
+
+ <field id="NU2amp_v2D" long_name="NU2 V barotropic harmonic Amplitude" unit="m/s" />
+ <field id="NU2phase_v2D" long_name="NU2 V barotropic harmonic Phase" unit="degree" />
+
+ </field_group>
+
+ <!-- SBC -->
+
+ <field_group id="SBC" grid_ref="grid_T_2D" > <!-- time step automaticaly defined based on nn_fsbc -->
+ <field id="empmr" long_name="Net Upward Water Flux" standard_name="water_flux_out_of_sea_ice_and_sea_water" unit="kg/m2/s" />
+ <field id="empbmr" long_name="Net Upward Water Flux at pre. tstep" standard_name="water_flux_out_of_sea_ice_and_sea_water" unit="kg/m2/s" />
+ <field id="emp_oce" long_name="Evap minus Precip over ocean" standard_name="evap_minus_precip_over_sea_water" unit="kg/m2/s" />
+ <field id="emp_ice" long_name="Evap minus Precip over ice" standard_name="evap_minus_precip_over_sea_ice" unit="kg/m2/s" />
+ <field id="saltflx" long_name="Downward salt flux" unit="1e-3/m2/s" />
+ <field id="fmmflx" long_name="Water flux due to freezing/melting" unit="kg/m2/s" />
+ <field id="snowpre" long_name="Snow precipitation" standard_name="snowfall_flux" unit="kg/m2/s" />
+ <field id="runoffs" long_name="River Runoffs" standard_name="water_flux_into_sea_water_from_rivers" unit="kg/m2/s" />
+ <field id="precip" long_name="Total precipitation" standard_name="precipitation_flux" unit="kg/m2/s" />
+
+ <field id="qt" long_name="Net Downward Heat Flux" standard_name="surface_downward_heat_flux_in_sea_water" unit="W/m2" />
+ <field id="qns" long_name="non solar Downward Heat Flux" unit="W/m2" />
+ <field id="qsr" long_name="Shortwave Radiation" standard_name="net_downward_shortwave_flux_at_sea_water_surface" unit="W/m2" />
+ <field id="qsr3d" long_name="Shortwave Radiation 3D distribution" standard_name="downwelling_shortwave_flux_in_sea_water" unit="W/m2" grid_ref="grid_T_3D" />
+ <field id="qrp" long_name="Surface Heat Flux: Damping" standard_name="heat_flux_into_sea_water_due_to_newtonian_relaxation" unit="W/m2" />
+ <field id="erp" long_name="Surface Water Flux: Damping" standard_name="water_flux_out_of_sea_water_due_to_newtonian_relaxation" unit="kg/m2/s" />
+ <field id="taum" long_name="wind stress module" standard_name="magnitude_of_surface_downward_stress" unit="N/m2" />
+ <field id="wspd" long_name="wind speed module" standard_name="wind_speed" unit="m/s" />
+ <field id="uwnd" long_name="u component of wind" unit="m/s" grid_ref="grid_U_2D" />
+ <field id="vwnd" long_name="v component of wind" unit="m/s" grid_ref="grid_V_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" />
+
+ <!-- * variable related to ice shelf forcing * -->
+ <field id="fwfisf" long_name="Ice shelf melting" unit="kg/m2/s" />
+ <field id="fwfisf3d" long_name="Ice shelf melting" unit="kg/m2/s" grid_ref="grid_T_3D" />
+ <field id="qlatisf" long_name="Ice shelf latent heat flux" unit="W/m2" />
+ <field id="qlatisf3d" long_name="Ice shelf latent heat flux" unit="W/m2" grid_ref="grid_T_3D" />
+ <field id="qhcisf" long_name="Ice shelf heat content flux" unit="W/m2" />
+ <field id="qhcisf3d" long_name="Ice shelf heat content flux" unit="W/m2" grid_ref="grid_T_3D" />
+ <field id="isfgammat" long_name="transfert coefficient for isf (temperature) " unit="m/s" />
+ <field id="isfgammas" long_name="transfert coefficient for isf (salinity) " unit="m/s" />
+ <field id="stbl" long_name="salinity in the Losh tbl " unit="PSU" />
+ <field id="ttbl" long_name="temperature in the Losh tbl " unit="C" />
+ <field id="utbl" long_name="zonal current in the Losh tbl at T point " unit="m/s" />
+ <field id="vtbl" long_name="merid current in the Losh tbl at T point " unit="m/s" />
+ <field id="thermald" long_name="thermal driving of ice shelf melting " unit="C" />
+ <field id="tfrz" long_name="top freezing point (used to compute melt) " unit="C" />
+ <field id="tinsitu" long_name="top insitu temperature (used to cmpt melt) " unit="C" />
+ <field id="ustar" long_name="ustar at T point used in ice shelf melting " unit="m/s" />
+
+ <!-- *_oce variables available with ln_blk_clio or ln_blk_core -->
+ <field id="qlw_oce" long_name="Longwave Downward Heat Flux over open ocean" standard_name="surface_net_downward_longwave_flux" unit="W/m2" />
+ <field id="qsb_oce" long_name="Sensible Downward Heat Flux over open ocean" standard_name="surface_downward_sensible_heat_flux" unit="W/m2" />
+ <field id="qla_oce" long_name="Latent Downward Heat Flux over open ocean" standard_name="surface_downward_latent_heat_flux" unit="W/m2" />
+ <field id="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="qns_oce" long_name="Downward non-solar heat flux" unit="W/m2" />
+ <field id="qsr_oce" long_name="Downward solar heat flux" unit="W/m2" />
+ <field id="qt_oce" long_name="total flux at ocean surface" standard_name="surface_downward_heat_flux_in_sea_water" unit="W/m2" />
+ <!-- available key_oasis3 -->
+ <field id="snow_ao_cea" long_name="Snow over ice-free ocean (cell average)" standard_name="snowfall_flux" unit="kg/m2/s" />
+ <field id="snow_ai_cea" long_name="Snow over sea-ice (cell average)" standard_name="snowfall_flux" unit="kg/m2/s" />
+ <field id="subl_ai_cea" long_name="Sublimation over sea-ice (cell average)" standard_name="surface_snow_and_ice_sublimation_flux" unit="kg/m2/s" />
+ <field id="icealb_cea" long_name="Ice albedo (cell average)" standard_name="sea_ice_albedo" unit="1" />
+ <field id="calving_cea" long_name="Calving" standard_name="water_flux_into_sea_water_from_icebergs" unit="kg/m2/s" />
+ <field id="iceberg_cea" long_name="Iceberg" standard_name="water_flux_into_sea_water_from_icebergs" unit="kg/m2/s" />
+ <field id="iceshelf_cea" long_name="Iceshelf" standard_name="water_flux_into_sea_water_from_iceshelf" unit="kg/m2/s" />
+
+
+ <!-- available if key_oasis3 + conservative method -->
+ <field id="rain" long_name="Liquid precipitation" standard_name="rainfall_flux" unit="kg/m2/s" />
+ <field id="evap_ao_cea" long_name="Evaporation over ice-free ocean (cell average)" standard_name="water_evaporation_flux" unit="kg/m2/s" />
+ <field id="isnwmlt_cea" long_name="Snow over Ice melting (cell average)" standard_name="surface_snow_melt_flux" unit="kg/m2/s" />
+ <field id="fsal_virt_cea" long_name="Virtual salt flux due to ice formation (cell average)" standard_name="virtual_salt_flux_into_sea_water_due_to_sea_ice_thermodynamics" unit="kg/m2/s" />
+ <field id="fsal_real_cea" long_name="Real salt flux due to ice formation (cell average)" standard_name="downward_sea_ice_basal_salt_flux" unit="kg/m2/s" />
+ <field id="hflx_rain_cea" long_name="heat flux due to rainfall" standard_name="temperature_flux_due_to_rainfall_expressed_as_heat_flux_into_sea_water" unit="W/m2" />
+ <field id="hflx_evap_cea" long_name="heat flux due to evaporation" standard_name="temperature_flux_due_to_evaporation_expressed_as_heat_flux_out_of_sea_water" unit="W/m2" />
+ <field id="hflx_snow_cea" long_name="heat flux due to snow falling" standard_name="heat_flux_onto_ocean_and_ice_due_to_snow_thermodynamics" unit="W/m2" />
+ <field id="hflx_snow_ai_cea" long_name="heat flux due to snow falling over ice" standard_name="heat_flux_onto_ice_due_to_snow_thermodynamics" unit="W/m2" />
+ <field id="hflx_snow_ao_cea" long_name="heat flux due to snow falling over ice-free ocean" standard_name="heat_flux_onto_sea_water_due_to_snow_thermodynamics" unit="W/m2" />
+ <field id="hflx_ice_cea" long_name="heat flux due to ice thermodynamics" standard_name="heat_flux_into_sea_water_due_to_sea_ice_thermodynamics" unit="W/m2" />
+ <field id="hflx_rnf_cea" long_name="heat flux due to runoffs" standard_name="temperature_flux_due_to_runoff_expressed_as_heat_flux_into_sea_water" unit="W/m2" />
+ <field id="hflx_cal_cea" long_name="heat flux due to calving" standard_name="heat_flux_into_sea_water_due_to_calving" unit="W/m2" />
+ <field id="hflx_icb_cea" long_name="heat flux due to iceberg" standard_name="heat_flux_into_sea_water_due_to_icebergs" unit="W/m2" />
+ <field id="hflx_isf_cea" long_name="heat flux due to iceshelf" standard_name="heat_flux_into_sea_water_due_to_iceshelf" unit="W/m2" />
+ <field id="bicemel_cea" long_name="Rate of Melt at Sea Ice Base (cell average)" standard_name="tendency_of_sea_ice_amount_due_to_basal_melting" unit="kg/m2/s" />
+ <field id="licepro_cea" long_name="Lateral Sea Ice Growth Rate (cell average)" standard_name="tendency_of_sea_ice_amount_due_to_lateral_growth_of_ice_floes" unit="kg/m2/s" />
+ <field id="snowmel_cea" long_name="Snow Melt Rate (cell average)" standard_name="surface_snow_melt_flux" unit="kg/m2/s" />
+ <field id="sntoice_cea" long_name="Snow-Ice Formation Rate (cell average)" standard_name="tendency_of_sea_ice_amount_due_to_snow_conversion" unit="kg/m2/s" />
+ <field id="ticemel_cea" long_name="Rate of Melt at Upper Surface of Sea Ice (cell average)" standard_name="tendency_of_sea_ice_amount_due_to_surface_melting" unit="kg/m2/s" />
+
+ <!-- ice field (nn_ice=1) -->
+ <field id="ice_cover" long_name="Ice fraction" standard_name="sea_ice_area_fraction" unit="1" />
+
+ <!-- dilution -->
+ <field id="emp_x_sst" long_name="Concentration/Dilution term on SST" unit="kg*degC/m2/s" />
+ <field id="emp_x_sss" long_name="Concentration/Dilution term on SSS" unit="kg*1e-3/m2/s" />
+ <field id="rnf_x_sst" long_name="Runoff term on SST" unit="kg*degC/m2/s" />
+ <field id="rnf_x_sss" long_name="Runoff term on SSS" unit="kg*1e-3/m2/s" />
+
+ <!-- sbcssm variables -->
+ <field id="sst_m" unit="degC" />
+ <field id="sss_m" unit="psu" />
+ <field id="ssu_m" unit="m/s" />
+ <field id="ssv_m" unit="m/s" />
+ <field id="ssh_m" unit="m" />
+ <field id="e3t_m" unit="m" />
+ <field id="frq_m" unit="-" />
+
+ </field_group>
+
+ <!-- U grid -->
+
+ <field_group id="grid_U" grid_ref="grid_U_2D">
+ <field id="e3u" long_name="U-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_U_3D" />
+ <field id="e3u_0" long_name="Initial U-cell thickness" standard_name="ref_cell_thickness" unit="m" grid_ref="grid_U_3D"/>
+ <field id="utau" long_name="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="ssu" long_name="ocean surface current along i-axis" unit="m/s" />
+ <field id="sbu" long_name="ocean bottom current along i-axis" unit="m/s" />
+ <field id="ubar" long_name="ocean barotropic current along i-axis" unit="m/s" />
+ <field id="uocetr_eff" long_name="Effective ocean transport along i-axis" standard_name="ocean_volume_x_transport" unit="m3/s" grid_ref="grid_U_3D" />
+ <field id="uocet" long_name="ocean transport along i-axis times temperature (CRS)" unit="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" />
+
+ <!-- u-eddy coefficients (ldftra) -->
+ <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="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 -->
+ <field id="psiu_mle" long_name="MLE streamfunction along i-axis" unit="m3/s" grid_ref="grid_U_3D" />
+
+ <!-- uoce_eiv: available EIV -->
+ <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" />
+
+ <!-- uoce_eiv: available with key_trabbl -->
+ <field id="uoce_bbl" long_name="BBL ocean current along i-axis" unit="m/s" />
+ <field id="ahu_bbl" long_name="BBL diffusive flux along i-axis" unit="m3/s" />
+
+ <!-- variable for ice shelves -->
+ <field id="utbl" long_name="zonal current in the Losh tbl" unit="m/s" />
+
+ <field id="u_masstr" long_name="Ocean Mass X Transport" standard_name="ocean_mass_x_transport" unit="kg/s" grid_ref="grid_U_3D" />
+ <field id="u_masstr_vint" long_name="vertical integral of ocean eulerian mass transport along i-axis" standard_name="vertical_integral_of_ocean_mass_x_transport" unit="kg/s" />
+ <field id="u_heattr" long_name="ocean eulerian heat transport along i-axis" standard_name="ocean_heat_x_transport" unit="W" />
+ <field id="u_salttr" long_name="ocean eulerian salt transport along i-axis" standard_name="ocean_salt_x_transport" unit="1e-3*kg/s" />
+ <field id="uadv_heattr" long_name="ocean advective heat transport along i-axis" standard_name="advectice_ocean_heat_x_transport" unit="W" />
+ <field id="uadv_salttr" long_name="ocean advective salt transport along i-axis" standard_name="advectice_ocean_salt_x_transport" unit="1e-3*kg/s" />
+ <field id="ueiv_heattr" long_name="ocean bolus heat transport along i-axis" standard_name="ocean_heat_x_transport_due_to_bolus_advection" unit="W" />
+ <field id="ueiv_salttr" long_name="ocean bolus salt transport along i-axis" standard_name="ocean_salt_x_transport_due_to_bolus_advection" unit="Kg" />
+ <field id="ueiv_heattr3d" long_name="ocean bolus heat transport along i-axis" standard_name="ocean_heat_x_transport_due_to_bolus_advection" unit="W" grid_ref="grid_U_3D" />
+ <field id="ueiv_salttr3d" long_name="ocean bolus salt transport along i-axis" standard_name="ocean_salt_x_transport_due_to_bolus_advection" unit="kg" grid_ref="grid_U_3D" />
+ <field id="udiff_heattr" long_name="ocean diffusion heat transport along i-axis" standard_name="ocean_heat_x_transport_due_to_diffusion" unit="W" />
+ <field id="udiff_salttr" long_name="ocean diffusion salt transport along i-axis" standard_name="ocean_salt_x_transport_due_to_diffusion" unit="1e-3*kg/s" />
+ </field_group>
+
+ <!-- V grid -->
+
+ <field_group id="grid_V" grid_ref="grid_V_2D">
+ <field id="e3v" long_name="V-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_V_3D" />
+ <field id="e3v_0" long_name="Initial V-cell thickness" standard_name="ref_cell_thickness" unit="m" grid_ref="grid_V_3D"/>
+ <field id="vtau" long_name="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" />
+ <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" />
+
+ <!-- v-eddy coefficients (ldftra, ldfdyn) -->
+ <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="aeiv_2d" long_name=" surface v-EIV coefficient" unit="m2/s" />
+ <field id="aeiv_3d" long_name=" 3D v-EIV coefficient" unit="m2/s" grid_ref="grid_V_3D" />
+
+ <!-- variables available with MLE -->
+ <field id="psiv_mle" long_name="MLE streamfunction along j-axis" unit="m3/s" grid_ref="grid_V_3D" />
+
+ <!-- voce_eiv: available with EIV -->
+ <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" />
+
+ <!-- voce_eiv: available with key_trabbl -->
+ <field id="voce_bbl" long_name="BBL ocean current along j-axis" unit="m/s" />
+ <field id="ahv_bbl" long_name="BBL diffusive flux along j-axis" unit="m3/s" />
+
+ <!-- variable for ice shelves -->
+ <field id="vtbl" long_name="meridional current in the Losh tbl" unit="m/s" />
+
+ <!-- variables available with diaar5 -->
+ <field id="v_masstr" long_name="ocean eulerian mass transport along j-axis" standard_name="ocean_mass_y_transport" unit="kg/s" grid_ref="grid_V_3D" />
+ <field id="v_heattr" long_name="ocean eulerian heat transport along j-axis" standard_name="ocean_heat_y_transport" unit="W" />
+ <field id="v_salttr" long_name="ocean eulerian salt transport along i-axis" standard_name="ocean_salt_y_transport" unit="1e-3*kg/s" />
+ <field id="vadv_heattr" long_name="ocean advective heat transport along j-axis" standard_name="advectice_ocean_heat_y_transport" unit="W" />
+ <field id="vadv_salttr" long_name="ocean advective salt transport along j-axis" standard_name="advectice_ocean_salt_y_transport" unit="1e-3*kg/s" />
+ <field id="veiv_heattr" long_name="ocean bolus heat transport along j-axis" standard_name="ocean_heat_y_transport_due_to_bolus_advection" unit="W" />
+ <field id="veiv_salttr" long_name="ocean bolus salt transport along j-axis" standard_name="ocean_salt_x_transport_due_to_bolus_advection" unit="Kg" />
+ <field id="veiv_heattr3d" long_name="ocean bolus heat transport along j-axis" standard_name="ocean_heat_y_transport_due_to_bolus_advection" unit="W" grid_ref="grid_V_3D" />
+ <field id="veiv_salttr3d" long_name="ocean bolus salt transport along j-axis" standard_name="ocean_salt_y_transport_due_to_bolus_advection" unit="kg" grid_ref="grid_V_3D" />
+ <field id="vdiff_heattr" long_name="ocean diffusion heat transport along j-axis" standard_name="ocean_heat_y_transport_due_to_diffusion" unit="W" />
+ <field id="vdiff_salttr" long_name="ocean diffusion salt transport along j-axis" standard_name="ocean_salt_y_transport_due_to_diffusion" unit="1e-3*kg/s" />
+ </field_group>
+
+ <!-- W grid -->
+
+ <field_group id="grid_W" grid_ref="grid_W_3D">
+ <field id="e3w" long_name="W-cell thickness" standard_name="cell_thickness" unit="m" />
+ <field id="woce" long_name="ocean vertical velocity" standard_name="upward_sea_water_velocity" unit="m/s" />
+ <field id="wocetr_eff" long_name="effective ocean vertical transport" unit="m3/s" />
+
+ <!-- woce_eiv: available with EIV -->
+ <field id="woce_eiv" long_name="EIV ocean vertical velocity" standard_name="bolus_upward_sea_water_velocity" unit="m/s" />
+
+ <field id="avt" long_name="vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" />
+ <field id="logavt" long_name="logarithm of vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" />
+ <field id="avm" long_name="vertical eddy viscosity" standard_name="ocean_vertical_momentum_diffusivity" unit="m2/s" />
+
+ <!-- avs: available with key_zdfddm -->
+ <field id="avs" long_name="salt vertical eddy diffusivity" standard_name="ocean_vertical_salt_diffusivity" unit="m2/s" />
+ <field id="logavs" long_name="logarithm of salt vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" />
+
+ <!-- avt_evd and avm_evd: available with ln_zdfevd -->
+ <field id="avt_evd" long_name="convective enhancement of vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_convection" unit="m2/s" />
+ <field id="avm_evd" long_name="convective enhancement of vertical viscosity" standard_name="ocean_vertical_momentum_diffusivity_due_to_convection" unit="m2/s" />
+
+ <!-- avt_tide: available with key_zdftmx -->
+ <field id="av_tide" long_name="tidal vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_tides" unit="m2/s" />
+
+ <!-- variables available with key_zdftmx_new -->
+ <field id="av_ratio" long_name="S over T diffusivity ratio" standard_name="salinity_over_temperature_diffusivity_ratio" unit="1" />
+ <field id="av_wave" long_name="wave-induced vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_internal_waves" unit="m2/s" />
+ <field id="bflx_tmx" long_name="wave-induced buoyancy flux" standard_name="buoyancy_flux_due_to_internal_waves" unit="W/kg" />
+ <field id="pcmap_tmx" long_name="power consumed by wave-driven mixing" standard_name="vertically_integrated_power_consumption_by_wave_driven_mixing" unit="W/m2" grid_ref="grid_W_2D" />
+ <field id="emix_tmx" long_name="power density available for mixing" standard_name="power_available_for_mixing_from_breaking_internal_waves" unit="W/kg" />
+
+ <!-- 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" />
+
+ <!-- aht2d and aht2d_eiv -->
+ <field id="aht2d" long_name="lateral eddy diffusivity" standard_name="ocean_tracer_xy_laplacian_diffusivity" unit="m2/s" grid_ref="grid_W_2D" />
+ <field id="aht2d_eiv" long_name="EIV lateral eddy diffusivity" standard_name="ocean_tracer_bolus_laplacian_diffusivity" unit="m2/s" grid_ref="grid_W_2D" />
+
+ </field_group>
+
+ <!-- F grid -->
+ <!-- f-eddy viscosity coefficients (ldfdyn) -->
+ <field id="ahmf_2d" long_name=" surface f-eddy viscosity coefficient" unit="m2/s or m4/s" />
+ <field id="ahmf_3d" long_name=" 3D f-eddy viscosity coefficient" unit="m2/s or m4/s" grid_ref="grid_T_3D"/>
+
+ <field_group id="scalar" grid_ref="grid_T_2D" >
+ <field id="voltot" long_name="global total volume" standard_name="sea_water_volume" unit="m3" />
+ <field id="sshtot" long_name="global mean ssh" standard_name="global_average_sea_level_change" unit="m" />
+ <field id="sshsteric" long_name="global mean ssh steric" standard_name="global_average_steric_sea_level_change" unit="m" />
+ <field id="sshthster" long_name="global mean ssh thermosteric" standard_name="global_average_thermosteric_sea_level_change" unit="m" />
+ <field id="masstot" long_name="global total mass" standard_name="sea_water_mass" unit="kg" />
+ <field id="temptot" long_name="global mean temperature" standard_name="sea_water_potential_temperature" unit="degC" />
+ <field id="saltot" long_name="global mean salinity" standard_name="sea_water_salinity" unit="1e-3" />
+ <field id="fram_trans" long_name="Sea Ice Mass Transport Through Fram Strait" standard_name="sea_ice_transport_across_line" unit="kg/s" />
+
+ <!-- available with ln_diahsb -->
+ <field id="bgtemper" long_name="drift in global mean temperature wrt timestep 1" standard_name="change_over_time_in_sea_water_potential_temperature" unit="degC" />
+ <field id="bgsaline" long_name="drift in global mean salinity wrt timestep 1" standard_name="change_over_time_in_sea_water_practical_salinity" unit="1e-3" />
+ <field id="bgheatco" long_name="drift in global mean heat content wrt timestep 1" unit="1.e20J" />
+ <field id="bgheatfx" long_name="drift in global mean heat flux wrt timestep 1" unit="W/m2" />
+ <field id="bgsaltco" long_name="drift in global mean salt content wrt timestep 1" unit="1e-3*km3" />
+ <field id="bgvolssh" long_name="drift in global mean ssh volume wrt timestep 1" unit="km3" />
+ <field id="bgvole3t" long_name="drift in global mean volume variation (e3t) wrt timestep 1" unit="km3" />
+ <field id="bgfrcvol" long_name="global mean volume from forcing" unit="km3" />
+ <field id="bgfrctem" long_name="global mean heat content from forcing" unit="1.e20J" />
+ <field id="bgfrchfx" long_name="global mean heat flux from forcing" unit="W/m2" />
+ <field id="bgfrcsal" long_name="global mean salt content from forcing" unit="1e-3*km3" />
+ <field id="bgmistem" long_name="global mean temperature error due to free surface (linssh true)" unit="degC" />
+ <field id="bgmissal" long_name="global mean salinity error due to free surface (linssh true)" unit="1e-3" />
+ </field_group>
+
+ <!-- variables available with key_float -->
+
+ <field_group id="floatvar" grid_ref="grid_T_nfloat" operation="instant" >
+ <field id="traj_lon" long_name="floats longitude" unit="degrees_east" />
+ <field id="traj_lat" long_name="floats latitude" unit="degrees_north" />
+ <field id="traj_dep" long_name="floats depth" unit="m" />
+ <field id="traj_temp" long_name="floats temperature" standard_name="sea_water_potential_temperature" unit="degC" />
+ <field id="traj_salt" long_name="floats salinity" standard_name="sea_water_practical_salinity" unit="1e-3" />
+ <field id="traj_dens" long_name="floats in-situ density" standard_name="sea_water_density" unit="kg/m3" />
+ <field id="traj_group" long_name="floats group" unit="1" />
+ </field_group>
+
+ <!-- variables available with iceberg trajectories -->
+
+ <field_group id="icbvar" domain_ref="grid_T" >
+ <field id="berg_melt" long_name="icb melt rate of icebergs" unit="kg/m2/s" />
+ <field id="berg_buoy_melt" long_name="icb buoyancy component of iceberg melt rate" unit="kg/m2/s" />
+ <field id="berg_eros_melt" long_name="icb erosion component of iceberg melt rate" unit="kg/m2/s" />
+ <field id="berg_conv_melt" long_name="icb convective component of iceberg melt rate" unit="kg/m2/s" />
+ <field id="berg_virtual_area" long_name="icb virtual coverage by icebergs" unit="m2" />
+ <field id="bits_src" long_name="icb mass source of bergy bits" unit="kg/m2/s" />
+ <field id="bits_melt" long_name="icb melt rate of bergy bits" unit="kg/m2/s" />
+ <field id="bits_mass" long_name="icb bergy bit density field" unit="kg/m2" />
+ <field id="berg_mass" long_name="icb iceberg density field" unit="kg/m2" />
+ <field id="calving" long_name="icb calving mass input" unit="kg/s" />
+ <field id="berg_floating_melt" long_name="icb melt rate of icebergs + bits" unit="kg/m2/s" />
+ <field id="berg_real_calving" long_name="icb calving into iceberg class" unit="kg/s" axis_ref="icbcla" />
+ <field id="berg_stored_ice" long_name="icb accumulated ice mass by class" unit="kg" axis_ref="icbcla" />
+ </field_group>
+
+ <!-- Poleward transport : ptr -->
+ <field_group id="diaptr" >
+ <field id="zomsfglo" long_name="Meridional Stream-Function: Global" unit="Sv" grid_ref="gznl_W_3D" />
+ <field id="zomsfatl" long_name="Meridional Stream-Function: Atlantic" unit="Sv" grid_ref="gznl_W_3D" />
+ <field id="zomsfpac" long_name="Meridional Stream-Function: Pacific" unit="Sv" grid_ref="gznl_W_3D" />
+ <field id="zomsfind" long_name="Meridional Stream-Function: Indian" unit="Sv" grid_ref="gznl_W_3D" />
+ <field id="zomsfipc" long_name="Meridional Stream-Function: Pacific+Indian" unit="Sv" grid_ref="gznl_W_3D" />
+ <field id="zotemglo" long_name="Zonal Mean Temperature : Global" unit="degree_C" grid_ref="gznl_T_3D" />
+ <field id="zotematl" long_name="Zonal Mean Temperature : Atlantic" unit="degree_C" grid_ref="gznl_T_3D" />
+ <field id="zotempac" long_name="Zonal Mean Temperature : Pacific" unit="degree_C" grid_ref="gznl_T_3D" />
+ <field id="zotemind" long_name="Zonal Mean Temperature : Indian" unit="degree_C" grid_ref="gznl_T_3D" />
+ <field id="zotemipc" long_name="Zonal Mean Temperature : Pacific+Indian" unit="degree_C" grid_ref="gznl_T_3D" />
+ <field id="zosalglo" long_name="Zonal Mean Salinity : Global" unit="0.001" grid_ref="gznl_T_3D" />
+ <field id="zosalatl" long_name="Zonal Mean Salinity : Atlantic" unit="0.001" grid_ref="gznl_T_3D" />
+ <field id="zosalpac" long_name="Zonal Mean Salinity : Pacific" unit="0.001" grid_ref="gznl_T_3D" />
+ <field id="zosalind" long_name="Zonal Mean Salinity : Indian" unit="0.001" grid_ref="gznl_T_3D" />
+ <field id="zosalipc" long_name="Zonal Mean Salinity : Pacific+Indian" unit="0.001" grid_ref="gznl_T_3D" />
+ <field id="zosrfglo" long_name="Zonal Mean Surface" unit="m2" grid_ref="gznl_T_3D" />
+ <field id="zosrfatl" long_name="Zonal Mean Surface : Atlantic" unit="m2" grid_ref="gznl_T_3D" />
+ <field id="zosrfpac" long_name="Zonal Mean Surface : Pacific" unit="m2" grid_ref="gznl_T_3D" />
+ <field id="zosrfind" long_name="Zonal Mean Surface : Indian" unit="m2" grid_ref="gznl_T_3D" />
+ <field id="zosrfipc" long_name="Zonal Mean Surface : Pacific+Indian" unit="m2" grid_ref="gznl_T_3D" />
+ <field id="sophtadv" long_name="Advective Heat Transport" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophtadv_atl" long_name="Advective Heat Transport: Atlantic" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophtadv_pac" long_name="Advective Heat Transport: Pacific" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophtadv_ind" long_name="Advective Heat Transport: Indian" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophtadv_ipc" long_name="Advective Heat Transport: Pacific+Indian" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophtldf" long_name="Diffusive Heat Transport" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophtldf_atl" long_name="Diffusive Heat Transport: Atlantic" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophtldf_pac" long_name="Diffusive Heat Transport: Pacific" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophtldf_ind" long_name="Diffusive Heat Transport: Indian" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophtldf_ipc" long_name="Diffusive Heat Transport: Pacific+Indian" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophtove" long_name="Overturning Heat Transport" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophtove_atl" long_name="Overturning Heat Transport: Atlantic" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophtove_pac" long_name="Overturning Heat Transport: Pacific" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophtove_ind" long_name="Overturning Heat Transport: Indian" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophtove_ipc" long_name="Overturning Heat Transport: Pacific+Indian" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophtbtr" long_name="Barotropic Heat Transport" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophtbtr_atl" long_name="Barotropic Heat Transport: Atlantic" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophtbtr_pac" long_name="Barotropic Heat Transport: Pacific" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophtbtr_ind" long_name="Barotropic Heat Transport: Indian" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophtbtr_ipc" long_name="Barotropic Heat Transport: Pacific+Indian" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophteiv" long_name="Heat Transport from mesoscale eddy advection" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophteiv_atl" long_name="Heat Transport from mesoscale eddy advection: Atlantic" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophteiv_pac" long_name="Heat Transport from mesoscale eddy advection: Pacific" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophteiv_ind" long_name="Heat Transport from mesoscale eddy advection: Indian" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sophteiv_ipc" long_name="Heat Transport from mesoscale eddy advection: Pacific+Indian" unit="PW" grid_ref="gznl_T_2D" />
+ <field id="sopstadv" long_name="Advective Salt Transport" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopstadv_atl" long_name="Advective Salt Transport: Atlantic" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopstadv_pac" long_name="Advective Salt Transport: Pacific" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopstadv_ind" long_name="Advective Salt Transport: Indian" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopstadv_ipc" long_name="Advective Salt Transport: Pacific+Indian" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopstove" long_name="Overturning Salt Transport" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopstove_atl" long_name="Overturning Salt Transport: Atlantic" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopstove_pac" long_name="Overturning Salt Transport: Pacific" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopstove_ind" long_name="Overturning Salt Transport: Indian" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopstove_ipc" long_name="Overturning Salt Transport: Pacific+Indian" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopstbtr" long_name="Barotropic Salt Transport" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopstbtr_atl" long_name="Barotropic Salt Transport: Atlantic" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopstbtr_pac" long_name="Barotropic Salt Transport: Pacific" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopstbtr_ind" long_name="Barotropic Salt Transport: Indian" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopstbtr_ipc" long_name="Barotropic Salt Transport: Pacific+Indian" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopstldf" long_name="Diffusive Salt Transport" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopstldf_atl" long_name="Diffusive Salt Transport: Atlantic" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopstldf_pac" long_name="Diffusive Salt Transport: Pacific" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopstldf_ind" long_name="Diffusive Salt Transport: Indian" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopstldf_ipc" long_name="Diffusive Salt Transport: Pacific+Indian" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopsteiv" long_name="Salt Transport from mesoscale eddy advection" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopsteiv_atl" long_name="Salt Transport from mesoscale eddy advection: Atlantic" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopsteiv_pac" long_name="Salt Transport from mesoscale eddy advection: Pacific" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopsteiv_ind" long_name="Salt Transport from mesoscale eddy advection: Indian" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ <field id="sopsteiv_ipc" long_name="Salt Transport from mesoscale eddy advection: Pacific+Indian" unit="Giga g/s" grid_ref="gznl_T_2D" />
+ </field_group>
+
+ <!--
+============================================================================================================
+ Physical ocean model trend diagnostics : temperature, KE, PE, momentum
+============================================================================================================
+ -->
+
+ <!-- variables available with ln_tra_trd -->
+ <!-- Asselin trends calculated on odd time steps-->
+ <field_group id="trendT_odd" grid_ref="grid_T_3D">
+ <field id="ttrd_atf" long_name="temperature-trend: asselin time filter" unit="degree_C/s" />
+ <field id="strd_atf" long_name="salinity -trend: asselin time filter" unit="0.001/s" />
+ <!-- Thickness weighted versions: -->
+ <field id="ttrd_atf_e3t" unit="degC/s * m" > ttrd_atf * e3t </field>
+ <field id="strd_atf_e3t" unit="1e-3/s * m" > strd_atf * e3t </field>
+ <!-- OMIP layer-integrated trends -->
+ <field id="ttrd_atf_li" long_name="layer integrated heat-trend: asselin time filter " unit="W/m^2" > ttrd_atf_e3t * 1026.0 * 3991.86795711963 </field>
+ <field id="strd_atf_li" long_name="layer integrated salt -trend: asselin time filter " unit="kg/(m^2 s)" > strd_atf_e3t * 1026.0 * 0.001 </field>
+ </field_group>
+
+ <!-- Other trends calculated on even time steps-->
+ <field_group id="trendT_even" grid_ref="grid_T_3D">
+ <field id="ttrd_xad" long_name="temperature-trend: i-advection" unit="degC/s" />
+ <field id="strd_xad" long_name="salinity -trend: i-advection" unit="1e-3/s" />
+ <field id="ttrd_yad" long_name="temperature-trend: j-advection" unit="degC/s" />
+ <field id="strd_yad" long_name="salinity -trend: j-advection" unit="1e-3/s" />
+ <field id="ttrd_zad" long_name="temperature-trend: k-advection" unit="degC/s" />
+ <field id="strd_zad" long_name="salinity -trend: k-advection" unit="1e-3/s" />
+ <field id="ttrd_ad" long_name="temperature-trend: advection" standard_name="tendency_of_sea_water_temperature_due_to_advection" unit="degC/s" > sqrt( ttrd_xad^2 + ttrd_yad^2 + ttrd_zad^2 ) </field>
+ <field id="strd_ad" long_name="salinity -trend: advection" standard_name="tendency_of_sea_water_salinity_due_to_advection" unit="1e-3/s" > sqrt( strd_xad^2 + strd_yad^2 + strd_zad^2 ) </field>
+ <field id="ttrd_totad" long_name="temperature-trend: total advection" standard_name="tendency_of_sea_water_salinity_due_to_advection" unit="degC/s" />
+ <field id="strd_totad" long_name="salinity -trend: total advection" standard_name="tendency_of_sea_water_salinity_due_to_advection" unit="1e-3/s" />
+ <field id="ttrd_sad" long_name="temperature-trend: surface adv. (linssh true)" unit="degC/s" grid_ref="grid_T_2D" />
+ <field id="strd_sad" long_name="salinity -trend: surface adv. (linssh true)" unit="1e-3/s" grid_ref="grid_T_2D" />
+ <field id="ttrd_ldf" long_name="temperature-trend: lateral diffusion" standard_name="tendency_of_sea_water_temperature_due_to_horizontal_mixing" unit="degC/s" />
+ <field id="strd_ldf" long_name="salinity -trend: lateral diffusion" standard_name="tendency_of_sea_water_salinity_due_to_horizontal_mixing" unit="1e-3/s" />
+ <field id="ttrd_zdf" long_name="temperature-trend: vertical diffusion" standard_name="tendency_of_sea_water_temperature_due_to_vertical_mixing" unit="degC/s" />
+ <field id="strd_zdf" long_name="salinity -trend: vertical diffusion" standard_name="tendency_of_sea_water_salinity_due_to_vertical_mixing" unit="1e-3/s" />
+ <field id="ttrd_evd" long_name="temperature-trend: EVD convection" unit="degC/s" />
+ <field id="strd_evd" long_name="salinity -trend: EVD convection" unit="1e-3/s" />
+
+ <!-- ln_traldf_iso=T only (iso-neutral diffusion) -->
+ <field id="ttrd_iso" long_name="temperature-trend: isopycnal diffusion" unit="degC/s" > ttrd_ldf + ttrd_zdf - ttrd_zdfp </field>
+ <field id="strd_iso" long_name="salinity -trend: isopycnal diffusion" unit="1e-3/s" > strd_ldf + strd_zdf - strd_zdfp </field>
+ <field id="ttrd_zdfp" long_name="temperature-trend: pure vert. diffusion" unit="degC/s" />
+ <field id="strd_zdfp" long_name="salinity -trend: pure vert. diffusion" unit="1e-3/s" />
+
+ <!-- -->
+ <field id="ttrd_dmp" long_name="temperature-trend: interior restoring" unit="degC/s" />
+ <field id="strd_dmp" long_name="salinity -trend: interior restoring" unit="1e-3/s" />
+ <field id="ttrd_bbl" long_name="temperature-trend: bottom boundary layer" unit="degC/s" />
+ <field id="strd_bbl" long_name="salinity -trend: bottom boundary layer" unit="1e-3/s" />
+ <field id="ttrd_npc" long_name="temperature-trend: non-penetrative conv." unit="degC/s" />
+ <field id="strd_npc" long_name="salinity -trend: non-penetrative conv." unit="1e-3/s" />
+ <field id="ttrd_qns" long_name="temperature-trend: non-solar flux + runoff" unit="degC/s" grid_ref="grid_T_2D" />
+ <field id="strd_cdt" long_name="salinity -trend: C/D term + runoff" unit="degC/s" grid_ref="grid_T_2D" />
+ <field id="ttrd_qsr" long_name="temperature-trend: solar penetr. heating" unit="degC/s" />
+ <field id="ttrd_bbc" long_name="temperature-trend: geothermal heating" unit="degC/s" />
+
+ <!-- Thickness weighted versions: -->
+ <field id="ttrd_xad_e3t" unit="degC/s * m" > ttrd_xad * e3t </field>
+ <field id="strd_xad_e3t" unit="1e-3/s * m" > strd_xad * e3t </field>
+ <field id="ttrd_yad_e3t" unit="degC/s * m" > ttrd_yad * e3t </field>
+ <field id="strd_yad_e3t" unit="1e-3/s * m" > strd_yad * e3t </field>
+ <field id="ttrd_zad_e3t" unit="degC/s * m" > ttrd_zad * e3t </field>
+ <field id="strd_zad_e3t" unit="1e-3/s * m" > strd_zad * e3t </field>
+ <field id="ttrd_ad_e3t" unit="degC/s * m" > ttrd_ad * e3t </field>
+ <field id="strd_ad_e3t" unit="1e-3/s * m" > strd_ad * e3t </field>
+ <field id="ttrd_totad_e3t" unit="degC/s * m" > ttrd_totad * e3t </field>
+ <field id="strd_totad_e3t" unit="1e-3/s * m" > strd_totad * e3t </field>
+ <field id="ttrd_ldf_e3t" unit="degC/s * m" > ttrd_ldf * e3t </field>
+ <field id="strd_ldf_e3t" unit="1e-3/s * m" > strd_ldf * e3t </field>
+ <field id="ttrd_zdf_e3t" unit="degC/s * m" > ttrd_zdf * e3t </field>
+ <field id="strd_zdf_e3t" unit="1e-3/s * m" > strd_zdf * e3t </field>
+ <field id="ttrd_evd_e3t" unit="degC/s * m" > ttrd_evd * e3t </field>
+ <field id="strd_evd_e3t" unit="1e-3/s * m" > strd_evd * e3t </field>
+
+ <!-- ln_traldf_iso=T only (iso-neutral diffusion) -->
+ <field id="ttrd_iso_e3t" unit="degC/s * m" > ttrd_iso * e3t </field>
+ <field id="strd_iso_e3t" unit="1e-3/s * m" > strd_iso * e3t </field>
+ <field id="ttrd_zdfp_e3t" unit="degC/s * m" > ttrd_zdfp * e3t </field>
+ <field id="strd_zdfp_e3t" unit="1e-3/s * m" > strd_zdfp * e3t </field>
+
+ <!-- -->
+ <field id="ttrd_dmp_e3t" unit="degC/s * m" > ttrd_dmp * e3t </field>
+ <field id="strd_dmp_e3t" unit="1e-3/s * m" > strd_dmp * e3t </field>
+ <field id="ttrd_bbl_e3t" unit="degC/s * m" > ttrd_bbl * e3t </field>
+ <field id="strd_bbl_e3t" unit="1e-3/s * m" > strd_bbl * e3t </field>
+ <field id="ttrd_npc_e3t" unit="degC/s * m" > ttrd_npc * e3t </field>
+ <field id="strd_npc_e3t" unit="1e-3/s * m" > strd_npc * e3t </field>
+ <field id="ttrd_qns_e3t" unit="degC/s * m" > ttrd_qns * e3t_surf </field>
+ <field id="strd_cdt_e3t" unit="degC/s * m" > strd_cdt * e3t_surf </field>
+ <field id="ttrd_qsr_e3t" unit="degC/s * m" > ttrd_qsr * e3t </field>
+ <field id="ttrd_bbc_e3t" unit="degC/s * m" > ttrd_bbc * e3t </field>
+
+ <!-- OMIP layer-integrated trends -->
+ <field id="ttrd_totad_li" long_name="layer integrated heat-trend : total advection" unit="W/m^2" > ttrd_totad_e3t * 1026.0 * 3991.86795711963 </field>
+ <field id="strd_totad_li" long_name="layer integrated salt -trend : total advection" unit="kg/(m^2 s)" > strd_totad_e3t * 1026.0 * 0.001 </field>
+ <field id="ttrd_evd_li" long_name="layer integrated heat-trend : EVD convection" unit="W/m^2" > ttrd_evd_e3t * 1026.0 * 3991.86795711963 </field>
+ <field id="strd_evd_li" long_name="layer integrated salt -trend : EVD convection" unit="kg/(m^2 s)" > strd_evd_e3t * 1026.0 * 0.001 </field>
+ <field id="ttrd_iso_li" long_name="layer integrated heat-trend : isopycnal diffusion" unit="W/m^2" > ttrd_iso_e3t * 1026.0 * 3991.86795711963 </field>
+ <field id="strd_iso_li" long_name="layer integrated salt -trend : isopycnal diffusion" unit="kg/(m^2 s)" > strd_iso_e3t * 1026.0 * 0.001 </field>
+ <field id="ttrd_zdfp_li" long_name="layer integrated heat-trend : pure vert. diffusion" unit="W/m^2" > ttrd_zdfp_e3t * 1026.0 * 3991.86795711963 </field>
+ <field id="strd_zdfp_li" long_name="layer integrated salt -trend : pure vert. diffusion" unit="kg/(m^2 s)" > strd_zdfp_e3t * 1026.0 * 0.001 </field>
+ <field id="ttrd_qns_li" long_name="layer integrated heat-trend : non-solar flux + runoff" unit="W/m^2" grid_ref="grid_T_2D"> ttrd_qns_e3t * 1026.0 * 3991.86795711963 </field>
+ <field id="ttrd_qsr_li" long_name="layer integrated heat-trend : solar flux" unit="W/m^2" grid_ref="grid_T_3D"> ttrd_qsr_e3t * 1026.0 * 3991.86795711963 </field>
+ <field id="ttrd_bbl_li" long_name="layer integrated heat-trend: bottom boundary layer " unit="W/m^2" > ttrd_bbl_e3t * 1026.0 * 3991.86795711963 </field>
+ <field id="strd_bbl_li" long_name="layer integrated salt -trend: bottom boundary layer " unit="kg/(m^2 s)" > strd_bbl_e3t * 1026.0 * 0.001 </field>
+ <field id="ttrd_evd_li" long_name="layer integrated heat -trend: evd convection " unit="W/m^2" >ttrd_evd_e3t * 1026.0 * 3991.86795711963 </field>
+ <field id="strd_evd_li" long_name="layer integrated salt -trend: evd convection " unit="kg/(m^2 s)" > strd_evd_e3t * 1026.0 * 0.001 </field>
+
+ </field_group>
+
+ <!-- Total trends calculated every time step-->
+ <field_group id="trendT" grid_ref="grid_T_3D">
+ <field id="ttrd_tot" long_name="temperature-trend: total model trend" unit="degC/s" />
+ <field id="strd_tot" long_name="salinity -trend: total model trend" unit="1e-3/s" />
+ <!-- Thickness weighted versions: -->
+ <field id="ttrd_tot_e3t" unit="degC/s * m" > ttrd_tot * e3t </field>
+ <field id="strd_tot_e3t" unit="1e-3/s * m" > strd_tot * e3t </field>
+ <!-- OMIP layer-integrated total trends -->
+ <field id="ttrd_tot_li" long_name="layer integrated heat-trend: total model trend :" unit="W/m^2" > ttrd_tot_e3t * 1026.0 * 3991.86795711963 </field>
+ <field id="strd_tot_li" long_name="layer integrated salt -trend: total model trend :" unit="kg/(m^2 s)" > strd_tot_e3t * 1026.0 * 0.001 </field>
+
+ <!-- **** these trends have not been apportioned to all/even/odd ts yet **** -->
+ <!-- variables available with ln_KE_trd -->
+ <field id="ketrd_hpg" long_name="ke-trend: hydrostatic pressure gradient" unit="W/s^3" />
+ <field id="ketrd_spg" long_name="ke-trend: surface pressure gradient" unit="W/s^3" />
+ <field id="ketrd_spgexp" long_name="ke-trend: surface pressure gradient (explicit)" unit="W/s^3" />
+ <field id="ketrd_spgflt" long_name="ke-trend: surface pressure gradient (filter)" unit="W/s^3" />
+ <field id="ssh_flt" long_name="filtered contribution to ssh (dynspg_flt)" unit="m" grid_ref="grid_T_2D" />
+ <field id="w0" long_name="surface vertical velocity" unit="m/s" grid_ref="grid_T_2D" />
+ <field id="pw0_exp" long_name="surface pressure flux due to ssh" unit="W/s^2" grid_ref="grid_T_2D" />
+ <field id="pw0_flt" long_name="surface pressure flux due to filtered ssh" unit="W/s^2" grid_ref="grid_T_2D" />
+ <field id="ketrd_keg" long_name="ke-trend: KE gradient or hor. adv." unit="W/s^3" />
+ <field id="ketrd_rvo" long_name="ke-trend: relative vorticity or metric term" unit="W/s^3" />
+ <field id="ketrd_pvo" long_name="ke-trend: planetary vorticity" unit="W/s^3" />
+ <field id="ketrd_zad" long_name="ke-trend: vertical advection" unit="W/s^3" />
+ <field id="ketrd_udx" long_name="ke-trend: U.dx[U]" unit="W/s^3" />
+ <field id="ketrd_ldf" long_name="ke-trend: lateral diffusion" unit="W/s^3" />
+ <field id="ketrd_zdf" long_name="ke-trend: vertical diffusion" unit="W/s^3" />
+ <field id="ketrd_tau" long_name="ke-trend: wind stress " unit="W/s^3" grid_ref="grid_T_2D" />
+ <field id="ketrd_bfr" long_name="ke-trend: bottom friction (explicit)" unit="W/s^3" />
+ <field id="ketrd_bfri" long_name="ke-trend: bottom friction (implicit)" unit="W/s^3" />
+ <field id="ketrd_atf" long_name="ke-trend: asselin time filter trend" unit="W/s^3" />
+ <field id="ketrd_convP2K" long_name="ke-trend: conversion (potential to kinetic)" unit="W/s^3" />
+ <field id="KE" long_name="kinetic energy: u(n)*u(n+1)/2" unit="W/s^2" />
+
+ <!-- variables available with ln_PE_trd -->
+ <field id="petrd_xad" long_name="pe-trend: i-advection" unit="W/m^3" />
+ <field id="petrd_yad" long_name="pe-trend: j-advection" unit="W/m^3" />
+ <field id="petrd_zad" long_name="pe-trend: k-advection" unit="W/m^3" />
+ <field id="petrd_sad" long_name="pe-trend: surface adv. (linssh true)" unit="W/m^3" grid_ref="grid_T_2D" />
+ <field id="petrd_ldf" long_name="pe-trend: lateral diffusion" unit="W/m^3" />
+ <field id="petrd_zdf" long_name="pe-trend: vertical diffusion" unit="W/m^3" />
+ <field id="petrd_zdfp" long_name="pe-trend: pure vert. diffusion" unit="W/m^3" />
+ <field id="petrd_dmp" long_name="pe-trend: interior restoring" unit="W/m^3" />
+ <field id="petrd_bbl" long_name="pe-trend: bottom boundary layer" unit="W/m^3" />
+ <field id="petrd_npc" long_name="pe-trend: non-penetrative conv." unit="W/m^3" />
+ <field id="petrd_nsr" long_name="pe-trend: surface forcing + runoff" unit="W/m^3" />
+ <field id="petrd_qsr" long_name="pe-trend: solar penetr. heating" unit="W/m^3" />
+ <field id="petrd_bbc" long_name="pe-trend: geothermal heating" unit="W/m^3" />
+ <field id="petrd_atf" long_name="pe-trend: asselin time filter" unit="W/m^3" />
+ <field id="PEanom" long_name="potential energy anomaly" unit="1" />
+ <field id="alphaPE" long_name="partial deriv. of PEanom wrt T" unit="degC-1" />
+ <field id="betaPE" long_name="partial deriv. of PEanom wrt S" unit="1e3" />
+ </field_group>
+
+ <field_group id="trendU" grid_ref="grid_U_3D">
+ <!-- variables available with ln_dyn_trd -->
+ <field id="utrd_hpg" long_name="i-trend: hydrostatic pressure gradient" unit="m/s^2" />
+ <field id="utrd_spg" long_name="i-trend: surface pressure gradient" unit="m/s^2" />
+ <field id="utrd_spgexp" long_name="i-trend: surface pressure gradient (explicit)" unit="m/s^2" />
+ <field id="utrd_spgflt" long_name="i-trend: surface pressure gradient (filtered)" unit="m/s^2" />
+ <field id="utrd_keg" long_name="i-trend: KE gradient or hor. adv." unit="m/s^2" />
+ <field id="utrd_rvo" long_name="i-trend: relative vorticity or metric term" unit="m/s^2" />
+ <field id="utrd_pvo" long_name="i-trend: planetary vorticity" unit="m/s^2" />
+ <field id="utrd_zad" long_name="i-trend: vertical advection" unit="m/s^2" />
+ <field id="utrd_udx" long_name="i-trend: U.dx[U]" unit="m/s^2" />
+ <field id="utrd_ldf" long_name="i-trend: lateral diffusion" unit="m/s^2" />
+ <field id="utrd_zdf" long_name="i-trend: vertical diffusion" unit="m/s^2" />
+ <field id="utrd_tau" long_name="i-trend: wind stress " unit="m/s^2" grid_ref="grid_U_2D" />
+ <field id="utrd_bfr" long_name="i-trend: bottom friction (explicit)" unit="m/s^2" />
+ <field id="utrd_bfri" long_name="i-trend: bottom friction (implicit)" unit="m/s^2" />
+ <field id="utrd_tot" long_name="i-trend: total momentum trend before atf" unit="m/s^2" />
+ <field id="utrd_atf" long_name="i-trend: asselin time filter trend" unit="m/s^2" />
+ </field_group>
+
+ <field_group id="trendV" grid_ref="grid_V_3D">
+ <!-- variables available with ln_dyn_trd -->
+ <field id="vtrd_hpg" long_name="j-trend: hydrostatic pressure gradient" unit="m/s^2" />
+ <field id="vtrd_spg" long_name="j-trend: surface pressure gradient" unit="m/s^2" />
+ <field id="vtrd_spgexp" long_name="j-trend: surface pressure gradient (explicit)" unit="m/s^2" />
+ <field id="vtrd_spgflt" long_name="j-trend: surface pressure gradient (filtered)" unit="m/s^2" />
+ <field id="vtrd_keg" long_name="j-trend: KE gradient or hor. adv." unit="m/s^2" />
+ <field id="vtrd_rvo" long_name="j-trend: relative vorticity or metric term" unit="m/s^2" />
+ <field id="vtrd_pvo" long_name="j-trend: planetary vorticity" unit="m/s^2" />
+ <field id="vtrd_zad" long_name="j-trend: vertical advection" unit="m/s^2" />
+ <field id="vtrd_vdy" long_name="i-trend: V.dx[V]" unit="m/s^2" />
+ <field id="vtrd_ldf" long_name="j-trend: lateral diffusion" unit="m/s^2" />
+ <field id="vtrd_zdf" long_name="j-trend: vertical diffusion" unit="m/s^2" />
+ <field id="vtrd_tau" long_name="j-trend: wind stress " unit="m/s^2" grid_ref="grid_V_2D" />
+ <field id="vtrd_bfr" long_name="j-trend: bottom friction (explicit)" unit="m/s^2" />
+ <field id="vtrd_bfri" long_name="j-trend: bottom friction (implicit)" unit="m/s^2" />
+ <field id="vtrd_tot" long_name="j-trend: total momentum trend before atf" unit="m/s^2" />
+ <field id="vtrd_atf" long_name="j-trend: asselin time filter trend" unit="m/s^2" />
+ </field_group>
+
+
+ <!--
+============================================================================================================
+ Definitions for iodef_demo.xml
+============================================================================================================
+ -->
+
+ <field_group id="TRD" >
+ <field field_ref="ttrd_totad_li" name="opottempadvect" />
+ <field field_ref="ttrd_iso_li" name="opottemppmdiff" />
+ <field field_ref="ttrd_zdfp_li" name="opottempdiff" />
+ <field field_ref="ttrd_evd_li" name="opottempevd" />
+ <field field_ref="strd_evd_li" name="osaltevd" />
+ <field field_ref="ttrd_qns_li" name="opottempqns" />
+ <field field_ref="ttrd_qsr_li" name="rsdoabsorb" operation="accumulate" />
+ <field field_ref="strd_totad_li" name="osaltadvect" />
+ <field field_ref="strd_iso_li" name="osaltpmdiff" />
+ <field field_ref="strd_zdfp_li" name="osaltdiff" />
+ </field_group>
+
+ <field_group id="mooring" >
+ <field field_ref="toce" name="thetao" long_name="sea_water_potential_temperature" />
+ <field field_ref="soce" name="so" long_name="sea_water_salinity" />
+ <field field_ref="uoce" name="uo" long_name="sea_water_x_velocity" />
+ <field field_ref="voce" name="vo" long_name="sea_water_y_velocity" />
+ <field field_ref="woce" name="wo" long_name="sea_water_z_velocity" />
+ <field field_ref="avt" name="difvho" long_name="ocean_vertical_heat_diffusivity" />
+ <field field_ref="avm" name="difvmo" long_name="ocean_vertical_momentum_diffusivity" />
+
+ <field field_ref="sst" name="tos" long_name="sea_surface_temperature" />
+ <field field_ref="sst2" name="tossq" long_name="square_of_sea_surface_temperature" />
+ <field field_ref="sstgrad" name="tosgrad" long_name="module_of_sea_surface_temperature_gradient" />
+ <field field_ref="sss" name="sos" long_name="sea_surface_salinity" />
+ <field field_ref="ssh" name="zos" long_name="sea_surface_height_above_geoid" />
+ <field field_ref="empmr" name="wfo" long_name="water_flux_into_sea_water" />
+ <field field_ref="qsr" name="rsntds" long_name="surface_net_downward_shortwave_flux" />
+ <field field_ref="qt" name="tohfls" long_name="surface_net_downward_total_heat_flux" />
+ <field field_ref="taum" />
+ <field field_ref="20d" />
+ <field field_ref="mldkz5" />
+ <field field_ref="mldr10_1" />
+ <field field_ref="mldr10_3" />
+ <field field_ref="mldr0_1" />
+ <field field_ref="mldr0_3" />
+ <field field_ref="mld_dt02" />
+ <field field_ref="topthdep" />
+ <field field_ref="pycndep" />
+ <field field_ref="tinv" />
+ <field field_ref="depti" />
+ <field field_ref="BLT" name="blt" long_name="barrier_layer_thickness" />
+ <field field_ref="utau" name="tauuo" long_name="surface_downward_x_stress" />
+ <field field_ref="vtau" name="tauvo" long_name="surface_downward_y_stress" />
+ </field_group>
+
+ <field_group id="groupT" >
+ <field field_ref="toce" name="thetao" long_name="sea_water_potential_temperature" />
+ <field field_ref="soce" name="so" long_name="sea_water_salinity" />
+ <field field_ref="sst" name="tos" long_name="sea_surface_temperature" />
+ <field field_ref="sst2" name="tossq" long_name="square_of_sea_surface_temperature" />
+ <field field_ref="sss" name="sos" long_name="sea_surface_salinity" />
+ <field field_ref="ssh" name="zos" long_name="sea_surface_height_above_geoid" />
+ <field field_ref="empmr" name="wfo" long_name="water_flux_into_sea_water" />
+ <field field_ref="qsr" name="rsntds" long_name="surface_net_downward_shortwave_flux" />
+ <field field_ref="qt" name="tohfls" long_name="surface_net_downward_total_heat_flux" />
+ <field field_ref="taum" />
+ <field field_ref="20d" />
+ <field field_ref="mldkz5" />
+ <field field_ref="mldr10_1" />
+ <field field_ref="mldr10_3" />
+ <field field_ref="mld_dt02" />
+ <field field_ref="topthdep" />
+ <field field_ref="pycndep" />
+ <field field_ref="tinv" />
+ <field field_ref="depti" />
+ <field field_ref="BLT" name="blt" long_name="Barrier Layer Thickness" />
+ </field_group>
+
+ <field_group id="groupU" >
+ <field field_ref="uoce" name="uo" long_name="sea_water_x_velocity" />
+ <field field_ref="ssu" name="uos" long_name="sea_surface_x_velocity" />
+ <field field_ref="utau" name="tauuo" long_name="surface_downward_x_stress" />
+ </field_group>
+
+ <field_group id="groupV" >
+ <field field_ref="voce" name="vo" long_name="sea_water_y_velocity" />
+ <field field_ref="ssv" name="vos" long_name="sea_surface_y_velocity" />
+ <field field_ref="vtau" name="tauvo" long_name="surface_downward_y_stress" />
+ </field_group>
+
+ <field_group id="groupW" >
+ <field field_ref="woce" name="wo" long_name="ocean vertical velocity" />
+ </field_group>
+
+ <!-- TMB diagnostic output -->
+ <field_group id="1h_grid_T_tmb" grid_ref="grid_T_2D" operation="instant">
+ <field id="top_temp" name="votemper_top" unit="degC" />
+ <field id="mid_temp" name="votemper_mid" unit="degC" />
+ <field id="bot_temp" name="votemper_bot" unit="degC" />
+ <field id="top_sal" name="vosaline_top" unit="psu" />
+ <field id="mid_sal" name="vosaline_mid" unit="psu" />
+ <field id="bot_sal" name="vosaline_bot" unit="psu" />
+ <field id="sshnmasked" name="sossheig" unit="m" />
+ </field_group>
+
+ <field_group id="1h_grid_U_tmb" grid_ref="grid_U_2D" operation="instant">
+ <field id="top_u" name="vozocrtx_top" unit="m/s" />
+ <field id="mid_u" name="vozocrtx_mid" unit="m/s" />
+ <field id="bot_u" name="vozocrtx_bot" unit="m/s" />
+ <field id="baro_u" name="vobtcrtx" unit="m/s" />
+ </field_group>
+
+ <field_group id="1h_grid_V_tmb" grid_ref="grid_V_2D" operation="instant">
+ <field id="top_v" name="vomecrty_top" unit="m/s" />
+ <field id="mid_v" name="vomecrty_mid" unit="m/s" />
+ <field id="bot_v" name="vomecrty_bot" unit="m/s" />
+ <field id="baro_v" name="vobtcrty" unit="m/s" />
+ </field_group>
+
+ <!-- 25h diagnostic output -->
+ <field_group id="25h_grid_T" grid_ref="grid_T_3D" 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" unit="m" />
+ </field_group>
+
+ <field_group id="25h_grid_U" grid_ref="grid_U_3D" operation="instant" >
+ <field id="vozocrtx25h" name="i current 25h mean" unit="m/s" />
+ </field_group>
+
+ <field_group id="25h_grid_V" grid_ref="grid_V_3D" operation="instant">
+ <field id="vomecrty25h" name="j current 25h mean" unit="m/s" />
+ </field_group>
+
+ <field_group id="25h_grid_W" grid_ref="grid_W_3D" operation="instant">
+ <field id="vomecrtz25h" name="k current 25h mean" unit="m/s" />
+ <field id="avt25h" name="vertical diffusivity25h mean" unit="m2/s" />
+ <field id="avm25h" name="vertical viscosity 25h mean" unit="m2/s" />
+ <field id="tke25h" name="turbulent kinetic energy 25h mean" />
+ <field id="mxln25h" name="mixing length 25h mean" unit="m" />
+ </field_group>
+
+ <!--
+============================================================================================================
+ -->
+ <!-- output variables for my configuration (example) -->
+
+ <field_group id="myvarOCE" >
+ <!-- grid T -->
+ <field field_ref="e3t" name="e3t" long_name="vertical scale factor" />
+ <field field_ref="sst" name="tos" long_name="sea_surface_temperature" />
+ <field field_ref="sss" name="sos" long_name="sea_surface_salinity" />
+ <field field_ref="ssh" name="zos" long_name="sea_surface_height_above_geoid" />
+
+ <!-- grid U -->
+ <field field_ref="e3u" name="e3u" long_name="vertical scale factor" />
+ <field field_ref="ssu" name="uos" long_name="sea_surface_x_velocity" />
+
+ <!-- grid V -->
+ <field field_ref="e3v" name="e3v" long_name="vertical scale factor" />
+ <field field_ref="ssv" name="vos" long_name="sea_surface_y_velocity" />
+ </field_group>
+
+ </field_definition>
diff --git a/EXP_Apr19/field_def_nemo-pisces.xml b/EXP_Apr19/field_def_nemo-pisces.xml
new file mode 100755
index 0000000000000000000000000000000000000000..3642fecd45165685616e1350b326fec4a0e93d3f
--- /dev/null
+++ b/EXP_Apr19/field_def_nemo-pisces.xml
@@ -0,0 +1,298 @@
+<?xml version="1.0"?>
+ <!-- $id$ -->
+
+ <!--
+============================================================================================================
+= definition of all existing variables =
+= DO NOT CHANGE =
+============================================================================================================
+ -->
+ <field_definition level="1" prec="4" operation="average" enabled=".TRUE." default_value="1.e20" > <!-- time step automaticaly defined -->
+
+ <!--
+============================================================================================================
+ Biogeochemistry model variables
+============================================================================================================
+ -->
+
+ <!-- ptrc on T grid -->
+
+ <field_group id="ptrc_T" grid_ref="grid_T_3D">
+ <!-- PISCES standard : variables available with ln_p4z -->
+ <field id="DIC" long_name="Dissolved inorganic Concentration" unit="mmol/m3" />
+ <field id="DIC_e3t" long_name="DIC * e3t" unit="mmol/m2" > DIC * e3t </field >
+ <field id="Alkalini" long_name="Total Alkalinity Concentration" unit="mmol/m3" />
+ <field id="Alkalini_e3t" long_name="Alkalini * e3t" unit="mmol/m2" > Alkalini * e3t </field >
+ <field id="O2" long_name="Oxygen Concentration" unit="mmol/m3" />
+ <field id="O2_e3t" long_name="O2 * e3t" unit="mmol/m2" > O2 * e3t </field >
+ <field id="CaCO3" long_name="Calcite Concentration" unit="mmol/m3" />
+ <field id="CaCO3_e3t" long_name="CaCO3 * e3t" unit="mmol/m2" > CaCO3 * e3t </field >
+ <field id="PO4" long_name="Phosphate Concentration" unit="mmol/m3" />
+ <field id="PO4_e3t" long_name="PO4 * e3t" unit="mmol/m2" > PO4 * e3t </field >
+ <field id="POC" long_name="Small organic carbon Concentration" unit="mmol/m3" />
+ <field id="POC_e3t" long_name="POC * e3t" unit="mmol/m2" > POC * e3t </field >
+ <field id="Si" long_name="Silicate Concentration" unit="mmol/m3" />
+ <field id="Si_e3t" long_name="Si * e3t" unit="mmol/m2" > Si * e3t </field >
+ <field id="PHY" long_name="(Nano)Phytoplankton Concentration" unit="mmol/m3" />
+ <field id="PHY_e3t" long_name="PHY * e3t" unit="mmol/m2" > PHY * e3t </field >
+ <field id="ZOO" long_name="(Micro)Zooplankton Concentration" unit="mmol/m3" />
+ <field id="ZOO_e3t" long_name="ZOO2 * e3t" unit="mmol/m2" > ZOO * e3t </field >
+ <field id="DOC" long_name="Dissolved organic Concentration" unit="mmol/m3" />
+ <field id="DOC_e3t" long_name="DOC * e3t" unit="mmol/m2" > DOC * e3t </field >
+ <field id="PHY2" long_name="Diatoms Concentration" unit="mmol/m3" />
+ <field id="PHY2_e3t" long_name="PHY2 * e3t" unit="mmol/m2" > PHY2 * e3t </field >
+ <field id="ZOO2" long_name="Mesozooplankton Concentration" unit="mmol/m3" />
+ <field id="ZOO2_e3t" long_name="ZOO2 * e3t" unit="mmol/m2" > ZOO2 * e3t </field >
+ <field id="DSi" long_name="Diatoms Silicate Concentration" unit="mmol/m3" />
+ <field id="DSi_e3t" long_name="Dsi * e3t" unit="mmol/m2" > DSi * e3t </field >
+ <field id="Fer" long_name="Dissolved Iron Concentration" unit="mmol/m3" />
+ <field id="Fer_e3t" long_name="Fer * e3t" unit="mmol/m2" > Fer * e3t </field >
+ <field id="BFe" long_name="Big iron particles Concentration" unit="mmol/m3" />
+ <field id="BFe_e3t" long_name="BFe * e3t" unit="mmol/m2" > BFe * e3t </field >
+ <field id="GOC" long_name="Big organic carbon Concentration" unit="mmol/m3" />
+ <field id="GOC_e3t" long_name="GOC * e3t" unit="mmol/m2" > GOC * e3t </field >
+ <field id="SFe" long_name="Small iron particles Concentration" unit="mmol/m3" />
+ <field id="SFe_e3t" long_name="SFe * e3t" unit="mmol/m2" > SFe * e3t </field >
+ <field id="DFe" long_name="Diatoms iron Concentration" unit="mmol/m3" />
+ <field id="DFe_e3t" long_name="DFe * e3t" unit="mmol/m2" > DFe * e3t </field >
+ <field id="GSi" long_name="Sinking biogenic Silicate Concentration" unit="mmol/m3" />
+ <field id="GSi_e3t" long_name="GSi * e3t" unit="mmol/m2" > GSi * e3t </field >
+ <field id="NFe" long_name="Nano iron Concentration" unit="mmol/m3" />
+ <field id="NFe_e3t" long_name="NFe * e3t" unit="mmol/m2" > NFe * e3t </field >
+ <field id="NCHL" long_name="Nano chlorophyl Concentration" unit="mg/m3" />
+ <field id="NCHL_e3t" long_name="NCHL * e3t" unit="mmol/m2" > NCHL * e3t </field >
+ <field id="DCHL" long_name="Diatoms chlorophyl Concentration" unit="mg/m3" />
+ <field id="DCHL_e3t" long_name="DCHL * e3t" unit="mmol/m2" > DCHL * e3t </field >
+ <field id="NO3" long_name="Nitrate Concentration" unit="mmol/m3" />
+ <field id="NO3_e3t" long_name="NO3 * e3t" unit="mmol/m2" > NO3 * e3t </field >
+ <field id="NH4" long_name="Ammonium Concentration" unit="mmol/m3" />
+ <field id="NH4_e3t" long_name="NH4 * e3t" unit="mmol/m2" > NH4 * e3t </field >
+
+ <!-- PISCES quota : variables available with ln_p5z -->
+
+ <field id="DON" long_name="Dissolved organic N Concentration" unit="mmol/m3" />
+ <field id="DON_e3t" long_name="DON * e3t" unit="mmol/m2" > DON * e3t </field >
+ <field id="DOP" long_name="Dissolved organic P Concentration" unit="mmol/m3" />
+ <field id="DOP_e3t" long_name="DOP * e3t" unit="mmol/m2" > DOP * e3t </field >
+ <field id="PON" long_name="Small PON Concentration" unit="mmol/m3" />
+ <field id="PON_e3t" long_name="PON * e3t" unit="mmol/m2" > PON * e3t </field >
+ <field id="POP" long_name="Small POP Concentration" unit="mmol/m3" />
+ <field id="POP_e3t" long_name="POP * e3t" unit="mmol/m2" > POP * e3t </field >
+ <field id="GON" long_name="Big PON Concentration" unit="mmol/m3" />
+ <field id="GON_e3t" long_name="GON * e3t" unit="mmol/m2" > GON * e3t </field >
+ <field id="GOP" long_name="Big POP Concentration" unit="mmol/m3" />
+ <field id="GOP_e3t" long_name="GOP * e3t" unit="mmol/m2" > GOP * e3t </field >
+ <field id="PHYN" long_name="Nanophytoplankton N biomass" unit="mmol/m3" />
+ <field id="PHYN_e3t" long_name="PHYN * e3t" unit="mmol/m2" > PHYN * e3t </field >
+ <field id="PHYP" long_name="Nanophytoplankton P biomass" unit="mmol/m3" />
+ <field id="PHYP_e3t" long_name="PHYP * e3t" unit="mmol/m2" > PHYP * e3t </field >
+ <field id="DIAN" long_name="Diatoms N biomass" unit="mmol/m3" />
+ <field id="DIAN_e3t" long_name="DIAN * e3t" unit="mmol/m2" > DIAN * e3t </field >
+ <field id="DIAP" long_name="Diatoms P biomass" unit="mmol/m3" />
+ <field id="DIAP_e3t" long_name="DIAP * e3t" unit="mmol/m2" > DIAP * e3t </field >
+ <field id="PIC" long_name="Picophytoplankton C biomass" unit="mmol/m3" />
+ <field id="PIC_e3t" long_name="PIC * e3t" unit="mmol/m2" > PIC * e3t </field >
+ <field id="PICN" long_name="Picophytoplankton N biomass" unit="mmol/m3" />
+ <field id="PICN_e3t" long_name="PICN * e3t" unit="mmol/m2" > PICN * e3t </field >
+ <field id="PICP" long_name="Picophytoplankton P biomass" unit="mmol/m3" />
+ <field id="PICP_e3t" long_name="PICP * e3t" unit="mmol/m2" > PICP * e3t </field >
+ <field id="PFe" long_name="Picophytoplankton Fe biomass" unit="mmol/m3" />
+ <field id="PFe_e3t" long_name="PFe * e3t" unit="mmol/m2" > PFe * e3t </field >
+ <field id="PCHL" long_name="Picophytoplankton Chl biomass" unit="mg/m3" />
+ <field id="PCHL_e3t" long_name="PCHL * e3t" unit="mmol/m2" > PCHL * e3t </field >
+
+ <!-- PISCES with ligand parametisation : variables available namelist paramter ln_ligand -->
+ <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 >
+ <field id="LFe" long_name="Lithogenic iron concentration" unit="mmol/m3" />
+ <field id="LFe_e3t" long_name="LFe * e3t" unit="mmol/m2" > LFe * 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 >
+
+ <!-- CFC11 : variables available with ln_cfc11 -->
+ <field id="CFC11" long_name="Chlorofluoro carbon11 Concentration" unit="umol/m3" />
+ <field id="CFC11_e3t" long_name="CFC11 * e3t" unit="umol/m2" > CFC11 * e3t </field >
+
+ <!-- CFC12 : variables available with ln_cfc12 -->
+ <field id="CFC12" long_name="Chlorofluoro carbon12 Concentration" unit="umol/m3" />
+ <field id="CFC12_e3t" long_name="CFC12 * e3t" unit="umol/m2" > CFC12 * e3t </field >
+
+ <!-- SF6 : variables available with ln_sf6 -->
+ <field id="SF6" long_name="Sulfur hexafluoride Concentration" unit="umol/m3" />
+ <field id="SF6_e3t" long_name="SF6 * e3t" unit="umol/m2" > SF6 * e3t </field >
+
+ <!-- C14 : variables available with ln_c14 -->
+ <field id="RC14" long_name="Radiocarbon ratio" unit="-" />
+ <field id="RC14_e3t" long_name="RC14 * e3t" unit="m" > RC14 * e3t </field >
+
+ <!-- AGE : variables available with ln_age -->
+ <field id="Age" long_name="Sea water age since surface contact" unit="yr" />
+ <field id="Age_e3t" long_name="Age * e3t" unit="yr * m" > Age * e3t </field >
+
+ </field_group>
+
+ <!-- PISCES additional diagnostics on T grid -->
+
+ <field_group id="diad_T" grid_ref="grid_T_2D">
+ <field id="PH" long_name="PH" unit="1" grid_ref="grid_T_3D" />
+ <field id="CO3" long_name="Bicarbonates" unit="mol/m3" grid_ref="grid_T_3D" />
+ <field id="CO3sat" long_name="CO3 saturation" unit="mol/m3" grid_ref="grid_T_3D" />
+ <field id="PAR" long_name="Photosynthetically Available Radiation" unit="W/m2" grid_ref="grid_T_3D" />
+ <field id="PARDM" long_name="Daily mean PAR" unit="W/m2" grid_ref="grid_T_3D" />
+ <field id="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="Fe2" long_name="Iron II concentration" unit="nmol/m3" grid_ref="grid_T_3D" />
+ <field id="Fe3" long_name="Iron III concentration" unit="nmol/m3" grid_ref="grid_T_3D" />
+ <field id="FeL1" long_name="Complexed Iron concentration with L1" unit="nmol/m3" grid_ref="grid_T_3D" />
+ <field id="FeL2" long_name="Complexed Iron concentration with L2" unit="nmol/m3" grid_ref="grid_T_3D" />
+ <field id="FeP" long_name="Precipitated Iron III" unit="nmol/m3" grid_ref="grid_T_3D" />
+ <field id="TL1" long_name="Total L1 concentration" unit="nmol/m3" grid_ref="grid_T_3D" />
+ <field id="TL2" long_name="Total L2 concentration" unit="nmol/m3" grid_ref="grid_T_3D" />
+ <field id="pdust" long_name="dust concentration" unit="g/m3" />
+ <field id="Totlig" long_name="Total ligand concentation" unit="nmol/m3" grid_ref="grid_T_3D" />
+ <field id="Biron" long_name="Bioavailable iron" unit="nmol/m3" grid_ref="grid_T_3D" />
+ <field id="Sdenit" long_name="Nitrate reduction in the sediments" unit="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="Dpo2" long_name="Delta O2" unit="uatm" />
+ <field id="Heup" long_name="Euphotic layer depth" unit="m" />
+ <field id="Irondep" long_name="Iron deposition from dust" unit="mol/m2/s" />
+ <field id="Ironsed" long_name="Iron deposition from sediment" unit="mol/m2/s" grid_ref="grid_T_3D" />
+
+ <!-- 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" />
+ <field id="INTPPPHYN" long_name="Vertically integrated primary production by nanophy" unit="mol/m2/s" />
+ <field id="INTPPPHYD" long_name="Vertically integrated primary production by diatom" unit="mol/m2/s" />
+ <field id="INTPP" long_name="Vertically integrated primary production by phyto" unit="mol/m2/s" />
+ <field id="INTPNEW" long_name="Vertically integrated new primary production" unit="mol/m2/s" />
+ <field id="INTPBFE" long_name="Vertically integrated of biogenic iron production" unit="mol/m2/s" />
+ <field id="INTPBSI" long_name="Vertically integrated of biogenic Si production" unit="mol/m2/s" />
+ <field id="INTPCAL" long_name="Vertically integrated of calcite production" unit="mol/m2/s" />
+
+ <!-- PISCES light : variables available with key_pisces_reduced -->
+ <field id="FNO3PHY" long_name="FNO3PHY" unit="" grid_ref="grid_T_3D" />
+ <field id="FNH4PHY" long_name="FNH4PHY" unit="" grid_ref="grid_T_3D" />
+ <field id="FNH4NO3" long_name="FNH4NO3" unit="" grid_ref="grid_T_3D" />
+ <field id="TNO3PHY" long_name="TNO3PHY" unit="" />
+ <field id="TNH4PHY" long_name="TNH4PHY" unit="" />
+ <field id="TPHYDOM" long_name="TPHYDOM" unit="" />
+ <field id="TPHYNH4" long_name="TPHYNH4" unit="" />
+ <field id="TPHYZOO" long_name="TPHYZOO" unit="" />
+ <field id="TPHYDET" long_name="TPHYDET" unit="" />
+ <field id="TDETZOO" long_name="TDETZOO" unit="" />
+ <field id="TZOODET" long_name="TZOODET" unit="" />
+ <field id="TZOOBOD" long_name="TZOOBOD" unit="" />
+ <field id="TZOONH4" long_name="TZOONH4" unit="" />
+ <field id="TZOODOM" long_name="TZOODOM" unit="" />
+ <field id="TNH4NO3" long_name="TNH4NO3" unit="" />
+ <field id="TDOMNH4" long_name="TDOMNH4" unit="" />
+ <field id="TDETNH4" long_name="TDETNH4" unit="" />
+ <field id="TPHYTOT" long_name="TPHYTOT" unit="" />
+ <field id="TZOOTOT" long_name="TZOOTOT" unit="" />
+ <field id="SEDPOC" long_name="SEDPOC" unit="" />
+ <field id="TDETSED" long_name="TDETSED" unit="" />
+
+ <!-- CFC11 : variables available with ln_cfc11 -->
+ <field id="qtr_CFC11" long_name="Air-sea flux of CFC-11" unit="mol/m2/s" />
+ <field id="qint_CFC11" long_name="Cumulative air-sea flux of CFC-11" unit="mol/m2" />
+
+ <!-- CFC12 : variables available with ln_cfc12 -->
+ <field id="qtr_CFC12" long_name="Air-sea flux of CFC12" unit="mol/m2/s" />
+ <field id="qint_CFC12" long_name="Cumulative air-sea flux of CFC12" unit="mol/m2" />
+
+ <!-- SF6 : variables available with ln_sf6 -->
+ <field id="qtr_SF6" long_name="Air-sea flux of SF6" unit="mol/m2/s" />
+ <field id="qint_SF6" long_name="Cumulative air-sea flux of SF6" unit="mol/m2" />
+
+ <!-- C14 : variables available with ln_c14 -->
+ <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="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" />
+ </field_group>
+
+ <field_group id="tracer_scalar" grid_ref="grid_T_2D" >
+ <!-- PISCES scalar -->
+ <field id="pno3tot" long_name="Global mean nitrate concentration" unit="mol/m3" />
+ <field id="ppo4tot" long_name="global mean phosphorus concentration" unit="mol/m3" />
+ <field id="psiltot" long_name="Global mean silicate concentration" unit="mol/m3" />
+ <field id="palktot" long_name="Global mean alkalinity concentration" unit="mol/m3" />
+ <field id="pfertot" long_name="Global mean iron concentration" unit="mol/m3" />
+ <field id="tcflx" long_name="Total Flux of Carbon out of the ocean" unit="mol/s" />
+ <field id="tcflxcum" long_name="Cumulative total Flux of Carbon out of the ocean" unit="mol/s" />
+ <field id="tcexp" long_name="Total Carbon export at 100m" unit="mol/s" />
+ <field id="tintpp" long_name="Global total integrated primary production" unit="mol/s" />
+ <field id="tnfix" long_name="Global total nitrogen fixation" unit="mol/s" />
+ <field id="tdenit" long_name="Total denitrification" unit="mol/s" />
+ <!-- C14 scalar -->
+ <field id="AtmCO2" long_name="Global atmospheric CO2" unit="ppm" />
+ <field id="AtmC14" long_name="Global atmospheric DeltaC14" unit="permil" />
+ <field id="K_C14" long_name="Global 14C/C exchange velocity" unit="m/yr" />
+ <field id="K_CO2" long_name="Global CO2 piston velocity" unit="cm/h" />
+ <field id="C14Inv" long_name="global Radiocarbon ocean inventory" unit="10^26 atoms" />
+ </field_group>
+
+ </field_definition>
diff --git a/EXP_Apr19/file_def_nemo-opa.xml b/EXP_Apr19/file_def_nemo-opa.xml
new file mode 100755
index 0000000000000000000000000000000000000000..0c764355f476753dab24fcc1735e284fc5564ac8
--- /dev/null
+++ b/EXP_Apr19/file_def_nemo-opa.xml
@@ -0,0 +1,55 @@
+<?xml version="1.0"?>
+
+ <!--
+============================================================================================================
+= output files definition =
+= Define your own files =
+= put the variables you want... =
+============================================================================================================
+ -->
+
+ <file_definition type="multiple_file" name="@expname@_@freq@_@startdate@_@enddate@" sync_freq="10d" min_digits="4">
+
+ <file_group id="15mi" output_freq="15mi" output_level="10" enabled=".TRUE." > <!-- 1 time step files -->
+ <file id="file1" name_suffix="_grid_T" description="ocean T grid variables" >
+ <field field_ref="ssh" name="zos" long_name="sea_surface_height_above_geoid" operation="instant" />
+ </file>
+
+ <file id="file2" name_suffix="_grid_U" description="ocean U grid variables" >
+ <field field_ref="ssu" name="uos" long_name="sea_water_surface_x_velocity" operation="instant" />
+ <field field_ref="uwnd" name="uwnd" long_name="u_component_of_wind" operation="instant" />
+ </file>
+
+ <file id="file3" name_suffix="_grid_V" description="ocean V grid variables" >
+ <field field_ref="ssv" name="vos" long_name="sea_water_surface_y_velocity" operation="instant" />
+ <field field_ref="vwnd" name="vwnd" long_name="v_component_of_wind" operation="instant" />
+ </file>
+
+ </file_group>
+
+
+
+ <file_group id="1h" output_freq="1h" output_level="10" enabled=".TRUE." />
+ <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="25h_mean" output_freq="1d" output_level="10" enabled=".TRUE."/>
+ <file_group id="1d" output_freq="1d" output_level="10" enabled=".TRUE."/> <!-- 1d files -->
+ <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="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_definition>
+
+
diff --git a/EXP_Apr19/file_def_nemo.xml b/EXP_Apr19/file_def_nemo.xml
new file mode 100755
index 0000000000000000000000000000000000000000..cf06a038a66766810b2f6d44155c0793525ed55f
--- /dev/null
+++ b/EXP_Apr19/file_def_nemo.xml
@@ -0,0 +1,151 @@
+<?xml version="1.0"?>
+
+ <!--
+============================================================================================================
+= output files definition =
+= Define your own files =
+= put the variables you want... =
+============================================================================================================
+ -->
+
+ <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 id="file19" name_suffix="_SSH" description="ocean T grid variables" >
+ <field field_ref="ssh" name="zos" />
+ <field field_ref="sst" name="tos" />
+ </file>
+ </file_group>
+
+ <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="1d" output_freq="1d" output_level="10" enabled=".TRUE."/> <!-- 1d files -->
+ <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 id="file20" name_suffix="_Tides" description="tidal harmonics" >
+ <field field_ref="M2x" name="M2x" long_name="M2 Elevation harmonic real part" />
+ <field field_ref="M2y" name="M2y" long_name="M2 Elevation harmonic imaginary part" />
+ <field field_ref="M2x_u" name="M2x_u" long_name="M2 current barotrope along i-axis harmonic real part " />
+ <field field_ref="M2y_u" name="M2y_u" long_name="M2 current barotrope along i-axis harmonic imaginary part " />
+ <field field_ref="M2x_v" name="M2x_v" long_name="M2 current barotrope along j-axis harmonic real part " />
+ <field field_ref="M2y_v" name="M2y_v" long_name="M2 current barotrope along j-axis harmonic imaginary part " />
+ </file>
+<!--
+ <file id="file11" name_suffix="_Tides_T" >
+ <field field_ref="M2x" name="M2_x_elev" />
+ <field field_ref="M2y" name="M2_y_elev" />
+ <field field_ref="S2x" name="S2_x_elev" />
+ <field field_ref="S2y" name="S2_y_elev" />
+ <field field_ref="N2x" name="N2_x_elev" />
+ <field field_ref="N2y" name="N2_y_elev" />
+ <field field_ref="K2x" name="K2_x_elev" />
+ <field field_ref="K2y" name="K2_y_elev" />
+ <field field_ref="K1x" name="K1_x_elev" />
+ <field field_ref="K1y" name="K1_y_elev" />
+ <field field_ref="M4x" name="M4_x_elev" />
+ <field field_ref="M4y" name="M4_y_elev" />
+ <field field_ref="Q1x" name="Q1_x_elev" />
+ <field field_ref="Q1y" name="Q1_y_elev" />
+ <field field_ref="O1x" name="O1_x_elev" />
+ <field field_ref="O1y" name="O1_y_elev" />
+ <field field_ref="P1x" name="P1_x_elev" />
+ <field field_ref="P1y" name="P1_y_elev" />
+ <field field_ref="S1x" name="S1_x_elev" />
+ <field field_ref="S1y" name="S1_y_elev" />
+ <field field_ref="2N2x" name="a2N2_x_elev" />
+ <field field_ref="2N2y" name="a2N2_y_elev" />
+ <field field_ref="MU2x" name="MU2_x_elev" />
+ <field field_ref="MU2y" name="MU2_y_elev" />
+ <field field_ref="NU2x" name="NU2_x_elev" />
+ <field field_ref="NU2y" name="NU2_y_elev" />
+ <field field_ref="L2x" name="L2_x_elev" />
+ <field field_ref="L2y" name="L2_y_elev" />
+ <field field_ref="T2x" name="T2_x_elev" />
+ <field field_ref="T2y" name="T2_y_elev" />
+ </file>
+
+ <file id="file12" name_suffix="_Tides_U" >
+ <field field_ref="M2x_u" name="M2_x_u" />
+ <field field_ref="M2y_u" name="M2_y_u" />
+ <field field_ref="S2x_u" name="S2_x_u" />
+ <field field_ref="S2y_u" name="S2_y_u" />
+ <field field_ref="N2x_u" name="N2_x_u" />
+ <field field_ref="N2y_u" name="N2_y_u" />
+ <field field_ref="K2x_u" name="K2_x_u" />
+ <field field_ref="K2y_u" name="K2_y_u" />
+ <field field_ref="K1x_u" name="K1_x_u" />
+ <field field_ref="K1y_u" name="K1_y_u" />
+ <field field_ref="M4x_u" name="M4_x_u" />
+ <field field_ref="M4y_u" name="M4_y_u" />
+ <field field_ref="Q1x_u" name="Q1_x_u" />
+ <field field_ref="Q1y_u" name="Q1_y_u" />
+ <field field_ref="O1x_u" name="O1_x_u" />
+ <field field_ref="O1y_u" name="O1_y_u" />
+ <field field_ref="P1x_u" name="P1_x_u" />
+ <field field_ref="P1y_u" name="P1_y_u" />
+ <field field_ref="S1x_u" name="S1_x_u" />
+ <field field_ref="S1y_u" name="S1_y_u" />
+ <field field_ref="2N2x_u" name="a2N2_x_u" />
+ <field field_ref="2N2y_u" name="a2N2_y_u" />
+ <field field_ref="MU2x_u" name="MU2_x_u" />
+ <field field_ref="MU2y_u" name="MU2_y_u" />
+ <field field_ref="NU2x_u" name="NU2_x_u" />
+ <field field_ref="NU2y_u" name="NU2_y_u" />
+ <field field_ref="L2x_u" name="L2_x_u" />
+ <field field_ref="L2y_u" name="L2_y_u" />
+ <field field_ref="T2x_u" name="T2_x_u" />
+ <field field_ref="T2y_u" name="T2_y_u" />
+ </file>
+
+ <file id="file13" name_suffix="_Tides_V" >
+ <field field_ref="M2x_v" name="M2_x_v" />
+ <field field_ref="M2y_v" name="M2_y_v" />
+ <field field_ref="S2x_v" name="S2_x_v" />
+ <field field_ref="S2y_v" name="S2_y_v" />
+ <field field_ref="N2x_v" name="N2_x_v" />
+ <field field_ref="N2y_v" name="N2_y_v" />
+ <field field_ref="K2x_v" name="K2_x_v" />
+ <field field_ref="K2y_v" name="K2_y_v" />
+ <field field_ref="K1x_v" name="K1_x_v" />
+ <field field_ref="K1y_v" name="K1_y_v" />
+ <field field_ref="M4x_v" name="M4_x_v" />
+ <field field_ref="M4y_v" name="M4_y_v" />
+ <field field_ref="Q1x_v" name="Q1_x_v" />
+ <field field_ref="Q1y_v" name="Q1_y_v" />
+ <field field_ref="O1x_v" name="O1_x_v" />
+ <field field_ref="O1y_v" name="O1_y_v" />
+ <field field_ref="P1x_v" name="P1_x_v" />
+ <field field_ref="P1y_v" name="P1_y_v" />
+ <field field_ref="S1x_v" name="S1_x_v" />
+ <field field_ref="S1y_v" name="S1_y_v" />
+ <field field_ref="2N2x_v" name="a2N2_x_v" />
+ <field field_ref="2N2y_v" name="a2N2_y_v" />
+ <field field_ref="MU2x_v" name="MU2_x_v" />
+ <field field_ref="MU2y_v" name="MU2_y_v" />
+ <field field_ref="NU2x_v" name="NU2_x_v" />
+ <field field_ref="NU2y_v" name="NU2_y_v" />
+ <field field_ref="L2x_v" name="L2_x_v" />
+ <field field_ref="L2y_v" name="L2_y_v" />
+ <field field_ref="T2x_v" name="T2_x_v" />
+ <field field_ref="T2y_v" name="T2_y_v" />
+ </file>
+-->
+ </file_group>
+
+ <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_definition>
+
diff --git a/EXP_Apr19/iodef.xml b/EXP_Apr19/iodef.xml
new file mode 100755
index 0000000000000000000000000000000000000000..9a1f846b72da76dd0491bf9947325406d84f4798
--- /dev/null
+++ b/EXP_Apr19/iodef.xml
@@ -0,0 +1,27 @@
+<?xml version="1.0"?>
+<simulation>
+
+<!-- ============================================================================================ -->
+<!-- XIOS context -->
+<!-- ============================================================================================ -->
+
+ <context id="xios">
+
+ <variable_definition>
+
+ <variable id="info_level" type="int">0</variable>
+ <variable id="using_server" type="bool">true</variable>
+ <variable id="using_oasis" type="bool">false</variable>
+ <variable id="oasis_codes_id" type="string" >oceanx</variable>
+
+ </variable_definition>
+
+ </context>
+
+<!-- ============================================================================================ -->
+<!-- NEMO CONTEXT add and suppress the components you need -->
+<!-- ============================================================================================ -->
+
+ <context id="nemo" src="./context_nemo.xml"/> <!-- NEMO -->
+
+</simulation>
diff --git a/EXP_Apr19/namelist_cfg b/EXP_Apr19/namelist_cfg
new file mode 100755
index 0000000000000000000000000000000000000000..8ba50b2b8e4d85e1da6b79dbbc51e415bdf36686
--- /dev/null
+++ b/EXP_Apr19/namelist_cfg
@@ -0,0 +1,407 @@
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!! NEMO/OPA : AMM12 configuration namelist used to overwrite defaults values defined in SHARED/namelist_ref
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!-----------------------------------------------------------------------
+&namrun ! parameters of the run
+!-----------------------------------------------------------------------
+ cn_exp = "BoBEAS" ! experience name
+ nn_it000 = 1 ! 57601 ! first time step ! CHNAGE THIS TO END OF LAST RUN + 1
+ nn_itend = 100 ! 218880 crashed in Apr16 ! 138240 ! last time step (std 1 day = 144) (for dt = 6 min, so 240*dt = 1 day) ! CHANGE THIS TO THIS + 56 DAYS
+ nn_date0 = 20190401 !20110717 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
+ nn_time0 = 0 ! initial time of day in hhmm
+ nn_leapy = 1 ! Leap year calendar (1) or not (0)
+ ln_rstart = .false. ! start from rest (F) or from a restart file (T)
+ nn_euler = 1 ! = 0 : start with forward time step if ln_rstart=T
+ nn_rstctl = 2 ! restart control ==> activated only if ln_rstart=T
+ ! ! = 0 nn_date0 read in namelist ; nn_it000 : read in namelist
+ ! ! = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart
+ ! ! = 2 nn_date0 read in restart ; nn_it000 : check consistancy between namelist and restart
+ cn_ocerst_in = "INDIAN_00138240_restart" ! suffix of ocean restart name (input) ! CHANGE THIS TO LAST RESTART FILE NUMBER
+ cn_ocerst_indir = "./restarts" ! directory from which to read input ocean restarts
+ cn_ocerst_out = "restart" ! suffix of ocean restart name (output)
+ cn_ocerst_outdir= "./restarts" ! directory in which to write output ocean restarts
+! nn_euler = 1 ! = 0 : start with forward time step if ln_rstart=.true.
+! cn_ocerst_in = "ammsurge_restart_oce" ! suffix of ocean restart name (input)
+! cn_ocerst_out = "restart_oce_out" ! suffix of ocean restart name (input)
+ nn_istate = 0 ! output the initial state (1) or not (0)
+ nn_stock = 100 ! 40320 =28days ! frequency of creation of a restart file (modulo referenced to 1) ! 28 DAYS
+ nn_write = 10 ! 86400 ! frequency of write in the output file (modulo referenced to nit000)
+/
+!-----------------------------------------------------------------------
+&namcfg ! parameters of the configuration
+!-----------------------------------------------------------------------
+ ln_read_cfg = .true. ! (=T) read the domain configuration file
+ ! (=F) user defined configuration ==>>> see usrdef(_...) modules
+ cn_domcfg = "domain_cfg.nc" ! domain configuration filename
+/
+!-----------------------------------------------------------------------
+!&namzgr ! vertical coordinate
+!-----------------------------------------------------------------------
+! ln_zps = .false. ! z-coordinate - partial steps (T/F)
+! ln_sco = .true. ! s- or hybrid z-s-coordinate (T/F)
+!/
+!-----------------------------------------------------------------------
+!&namzgr_sco ! s-coordinate or hybrid z-s-coordinate
+!-----------------------------------------------------------------------
+! rn_sbot_min = 6.0 ! minimum depth of s-bottom surface (>0) (m)
+! rn_hc = 0.0 ! critical depth for transition to stretched coordinates
+! rn_rmax = 1.0 ! maximum cut-off r-value allowed (0<r_max<1)
+! rn_theta = 0.0 ! surface control parameter (0<=theta<=20)
+!/
+!-----------------------------------------------------------------------
+&namdom ! space and time domain (bathymetry, mesh, timestep)
+!-----------------------------------------------------------------------
+! ln_2d = .false. ! (=T) run in 2D barotropic mode (no tracer processes or vertical diffusion)
+ rn_rdt = 60. ! time step for the dynamics (and tracer if nn_acc=0)
+/
+
+!-----------------------------------------------------------------------
+&namtsd ! data : Temperature & Salinity
+!-----------------------------------------------------------------------
+! ! 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_tem = 'initcd_votemper.nc', -1 ,'votemper', .false. , .true. , 'yearly' , '' , '' , ''
+ sn_sal = 'initcd_vosaline.nc', -1 ,'vosaline', .false. , .true. , 'yearly' , '' , '' , ''
+ sn_dep = 'initcd_depth.nc' , -12 ,'gdept_4D', .false. , .true. , 'yearly' , '' , '' , ''
+ sn_msk = 'sosie_initcd_mask.nc', -12 ,'mask', .false. , .true. , 'yearly' , '' , '' , ''
+
+ !
+ cn_dir = './ICS/' ! root directory for the location of the runoff files
+ ln_tsd_init = .true. ! Initialisation of ocean T & S with T & S input data (T) or not (F)
+ ln_tsd_interp = .true. ! Interpolation of T & S in the verticalinput data (T) or not (F)
+ ln_tsd_tradmp = .false. ! damping of ocean T & S toward T & S input data (T) or not (F)
+/
+!-----------------------------------------------------------------------
+&namsbc ! Surface Boundary Condition (surface module)
+!-----------------------------------------------------------------------
+ nn_fsbc = 1 ! frequency of surface boundary condition computation
+ ! (also = the frequency of sea-ice model call)
+ ln_usr = .false.
+ ln_blk = .true.
+ ln_apr_dyn = .false. ! Patm gradient added in ocean & ice Eqs. (T => fill namsbc_apr )
+ nn_ice = 0 ! =0 no ice boundary condition ,
+ ln_rnf = .true. ! Runoffs (T => fill namsbc_rnf)
+ ln_ssr = .false. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr)
+ ln_traqsr = .false. ! Light penetration in the ocean (T => fill namtra_qsr)
+ nn_fwb = 0 ! FreshWater Budget: =0 unchecked
+/
+
+
+!-----------------------------------------------------------------------
+&namsbc_usr ! namsbc_surge surge model fluxes
+!-----------------------------------------------------------------------
+ ln_use_sbc = .false. ! (T) to turn on surge fluxes (wind and pressure only)
+ ! (F) for no fluxes (ie tide only case)
+
+!
+! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation !
+! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing !
+! sn_wndi = 'CARIB_u10_weights_bicubic_atmos.nc' , 1 ,'x_wind', .true. , .false. , 'daily' ,'weights_era5_INDIAN_bicubic.nc' , ''
+! sn_wndj = 'CARIB_v10_weights_bicubic_atmos.nc' , 1 ,'y_wind', .true. , .false. , 'daily' ,'weights_era5_INDIAN_bicubic.nc' , ''
+ cn_dir = './FORCING/' ! root directory for the location of the bulk files
+ rn_vfac = 1. ! multiplicative factor for ocean/ice velocity
+ ! in the calculation of the wind stress (0.=absolute winds or 1.=relative winds)
+ rn_charn_const = 0.0275
+/
+
+!-----------------------------------------------------------------------
+&namtra_qsr ! penetrative solar radiation
+!-----------------------------------------------------------------------
+/
+
+!-----------------------------------------------------------------------
+&namsbc_apr ! Atmospheric pressure used as ocean forcing or in bulk
+!-----------------------------------------------------------------------
+! ! 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_apr= 'ERA5_MSL', 1 , 'air_pressure_at_sea_level' , .true. , .false., 'daily' , 'weights_era5_INDIAN_bicubic.nc' , '' , 'ERA5_LSM'
+ cn_dir = './FORCING/'! root directory for the location of the bulk files
+ rn_pref = 101200. ! reference atmospheric pressure [N/m2]/
+ ln_ref_apr = .false. ! ref. pressure: global mean Patm (T) or a constant (F)
+ ln_apr_obc = .true. ! inverse barometer added to OBC ssh data
+/
+!-----------------------------------------------------------------------
+&namberg ! iceberg parameters (default: No iceberg)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namlbc ! lateral momentum boundary condition
+!-----------------------------------------------------------------------
+ rn_shlat = 0 ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat
+ ! free slip ! partial slip ! no slip ! strong slip
+/
+
+!-----------------------------------------------------------------------
+&nam_tide ! tide parameters
+!-----------------------------------------------------------------------
+ ln_tide = .true.
+ rdttideramp = 0
+ clname(1) = 'M2' ! name of constituent
+ clname(2) = 'K2'
+ clname(3) = 'S2'
+ clname(4) = 'N2'
+ clname(5) = 'Q1' ! name of constituent
+ clname(6) = 'O1'
+ clname(7) = 'P1'
+ clname(8) = 'K1'
+! clname(5) = '2N2'
+! clname(6) = 'MU2'
+! clname(7) = 'N2'
+! clname(8) = 'NU2'
+! clname(9) = 'M2'
+! clname(10) = 'L2'
+! clname(11) = 'T2'
+! clname(12) = 'S2'
+! clname(13) = 'K2'
+! clname(14) = 'M4'
+! clname(15) = '2MK6'
+! clname(16) = '2MS6'
+! clname(17) = '2SM2'
+! clname(18) = '3M2S2'
+! clname(19) = 'Lam2'
+! clname(20) = 'M6'
+! clname(21) = 'MK3'
+! clname(22) = 'MN4'
+! clname(23) = 'MNS2'
+! clname(24) = 'MO3'
+! clname(25) = 'MS4'
+! clname(26) = 'MSN2'
+/
+!-----------------------------------------------------------------------
+&nambdy ! unstructured open boundaries
+!-----------------------------------------------------------------------
+ ln_bdy = .true.
+ nb_bdy = 1 ! number of open boundary sets
+ cn_coords_file = 'coordinates.bdy.nc' ! bdy coordinates files
+ cn_dyn2d = 'flather' ! was flather
+ nn_dyn2d_dta = 0 ! = 0, bdy data are equal to the initial state
+ ! = 1, bdy data are read in 'bdydata .nc' files
+ ! = 2, use tidal harmonic forcing data from files
+ ! = 3, use external data AND tidal harmonic forcing
+
+ cn_dyn3d = 'frs' !
+ nn_dyn3d_dta = 0 ! = 0, bdy data are equal to the initial state
+ ! = 1, bdy data are read in 'bdydata .nc' files
+
+
+ cn_tra = 'frs' ! 'frs'
+ nn_tra_dta = 0 ! = 0, bdy data are equal to the initial state
+ ! = 1, bdy data are read in 'bdydata .nc' files
+ nn_rimwidth = 1 ! width of the relaxation zone
+ ln_tra_dmp = .false.
+/
+!-----------------------------------------------------------------------
+&nambdy_dta ! open boundaries - external data
+!-----------------------------------------------------------------------
+! ! 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 !
+ bn_ssh = 'INDIAN_bdy_SSH' , 24 , 'sossheig', .true. , .false. , 'monthly' , '' , '' , ''
+ bn_u2d = 'INDIAN_bdy_U' , 24 , 'vobtcrtx', .true. , .false. , 'monthly' , '' , '' , ''
+ bn_v2d = 'INDIAN_bdy_V' , 24 , 'vobtcrty', .true. , .false. , 'monthly' , '' , '' , ''
+ bn_u3d = 'INDIAN_bdy_U' , 24 , 'vozocrtx', .true. , .false. , 'monthly' , '' , '' , ''
+ bn_v3d = 'INDIAN_bdy_V' , 24 , 'vomecrty', .true. , .false. , 'monthly' , '' , '' , ''
+ bn_tem = 'INDIAN_bdy_TEMP', 24 , 'votemper', .true. , .false. , 'monthly' , '' , '' , ''
+ bn_sal = 'INDIAN_bdy_SAL' , 24 , 'vosaline', .true. , .false. , 'monthly' , '' , '' , ''
+/
+!-----------------------------------------------------------------------
+&nambdy_tide ! tidal forcing at open boundaries
+!-----------------------------------------------------------------------
+ filtide = 'TIDES/INDIAN_bdytide_rotT_' ! file name root of tidal forcing files
+ ln_bdytide_2ddta = .false.
+ ln_bdytide_conj = .false. !
+/
+!-----------------------------------------------------------------------
+&nambfr ! bottom friction
+!-----------------------------------------------------------------------
+ nn_bfr = 2 ! type of bottom friction : = 0 : free slip, = 1 : linear friction
+ ! = 2 : nonlinear friction
+ rn_bfri2 = 2.4e-3 ! bottom drag coefficient (non linear case)
+ rn_bfeb2 = 0.0e0 ! bottom turbulent kinetic energy background (m2/s2)
+ ln_loglayer = .false. ! loglayer bottom friction (only effect when nn_bfr = 2)
+ rn_bfrz0 = 0.003 ! bottom roughness (only effect when ln_loglayer = .true.)
+/
+!-----------------------------------------------------------------------
+&nambbc ! bottom temperature boundary condition
+!-----------------------------------------------------------------------
+ ln_trabbc = .false. ! Apply a geothermal heating at the ocean bottom
+/
+!-----------------------------------------------------------------------
+&nambbl ! bottom boundary layer scheme
+!-----------------------------------------------------------------------
+ nn_bbl_ldf = 0 ! diffusive bbl (=1) or not (=0)
+/
+!-----------------------------------------------------------------------
+&nameos ! ocean physical parameters
+!-----------------------------------------------------------------------
+ ln_teos10 = .true. ! = Use TEOS-10 equation of state
+/
+!-----------------------------------------------------------------------
+&namtra_adv ! advection scheme for tracer
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param)
+!-----------------------------------------------------------------------
+/
+!----------------------------------------------------------------------------------
+&namtra_ldf ! lateral diffusion scheme for tracers
+!----------------------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namtra_ldfeiv ! eddy induced velocity param. (default: NO)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namtra_dmp ! tracer: T & S newtonian damping
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namdyn_adv ! formulation of the momentum advection
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namdyn_vor ! option of physics/algorithm (not control by CPP keys)
+!-----------------------------------------------------------------------
+ ln_dynvor_een = .true. ! energy & enstrophy scheme
+/
+!-----------------------------------------------------------------------
+&namdyn_hpg ! Hydrostatic pressure gradient option
+!-----------------------------------------------------------------------
+ ln_hpg_zps = .false. ! z-coordinate - partial steps (interpolation)
+ ln_hpg_sco = .true. ! s-coordinate (Standard Jacobian scheme)
+/
+!-----------------------------------------------------------------------
+&namdyn_spg ! surface pressure gradient (CPP key only)
+!-----------------------------------------------------------------------
+ ln_dynspg_ts = .true. ! split-explicit free surface
+ ln_bt_auto = .true. ! Set nn_baro automatically to be just below
+ ! a user defined maximum courant number (rn_bt_cmax)
+/
+!-----------------------------------------------------------------------
+&namdyn_ldf ! lateral diffusion on momentum
+!-----------------------------------------------------------------------
+ ! ! Type of the operator :
+ ln_dynldf_blp = .true. ! bilaplacian operator
+ ln_dynldf_lap = .false. ! bilaplacian operator
+ ! ! Direction of action :
+ ln_dynldf_lev = .true. ! iso-level
+ ! Coefficient
+ rn_ahm_0 = 60.0 ! horizontal laplacian eddy viscosity [m2/s]
+ rn_bhm_0 = -1.0e+9 ! horizontal bilaplacian eddy viscosity [m4/s]
+/
+!-----------------------------------------------------------------------
+&namzdf ! vertical physics
+!-----------------------------------------------------------------------
+ rn_avm0 = 0.1e-6 ! vertical eddy viscosity [m2/s] (background Kz if not "key_zdfcst")
+ rn_avt0 = 0.1e-6 ! vertical eddy diffusivity [m2/s] (background Kz if not "key_zdfcst")
+ ln_zdfevd = .false. ! enhanced vertical diffusion (evd) (T) or not (F)
+ nn_evdm = 1 ! evd apply on tracer (=0) or on tracer and momentum (=1)
+/
+!-----------------------------------------------------------------------
+&namzdf_ric ! richardson number dependent vertical diffusion ("key_zdfric" )
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion ("key_zdftke")
+!-----------------------------------------------------------------------
+/
+!------------------------------------------------------------------------
+&namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionally:
+!------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")
+/
+!-----------------------------------------------------------------------
+&namzdf_gls ! GLS vertical diffusion ("key_zdfgls")
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namzdf_ddm ! double diffusive mixing parameterization ("key_zdfddm")
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namzdf_tmx ! tidal mixing parameterization ("key_zdftmx")
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namsol ! elliptic solver / island / free surface
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&nammpp ! Massively Parallel Processing ("key_mpp_mpi)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namctl ! Control prints & Benchmark
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namnc4 ! netcdf4 chunking and compression settings ("key_netcdf4")
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namtrd ! diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra")
+! ! or mixed-layer trends or barotropic vorticity ("key_trdmld" or "key_trdvor")
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namflo ! float parameters ("key_float")
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namptr ! Poleward Transport Diagnostic
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namhsb ! Heat and salt budgets
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namdiu ! Cool skin and warm layer models (default F)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&nam_diaharm ! Harmonic analysis of tidal constituents ('key_diaharm')
+!-----------------------------------------------------------------------
+ nit000_han = 1 !43201 ! First time step used for harmonic analysis
+ nitend_han = 43680 !86400 ! Last time step used for harmonic analysis
+ nstep_han = 5 ! Time step frequency for harmonic analysis
+ tname(1) = 'M2' ! Name of tidal constituents
+ tname(2) = 'K2'
+ tname(3) = 'S2'
+/
+!-----------------------------------------------------------------------
+&namdct ! transports through sections
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&nam_diatmb ! Top Middle Bottom Output (default F)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&nam_dia25h ! 25h Mean Output (default F)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namobs ! observation usage switch ('key_diaobs')
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&nam_asminc ! assimilation increments ('key_asminc')
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namsbc_wave ! External fields from wave model
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namdyn_nept ! Neptune effect (simplified: lateral and vertical diffusions removed)
+!-----------------------------------------------------------------------
+/
+!-----------------------------------------------------------------------
+&namwad ! Wetting and Drying namelist
+!-----------------------------------------------------------------------
+ ln_wd = .false. !: key to turn on/off wetting/drying (T: on, F: off)
+ rn_wdmin1=0.1 !: minimum water depth on dried cells
+ rn_wdmin2 = 0.01 !: tolerrance of minimum water depth on dried cells
+ rn_wdld = 20.0 !: land elevation below which wetting/drying will be considered
+ nn_wdit = 10 !: maximum number of iteration for W/D limiter
+/
diff --git a/EXP_Apr19/namelist_ref b/EXP_Apr19/namelist_ref
new file mode 100755
index 0000000000000000000000000000000000000000..d301714c2823a2b470b4cc49773e571e16b7a8fe
--- /dev/null
+++ b/EXP_Apr19/namelist_ref
@@ -0,0 +1,1200 @@
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!! namelist_ref
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!! NEMO/OPA : 1 - run manager (namrun)
+!! namelists 2 - Domain (namcfg, namzgr, namdom, namtsd, namcrs, namc1d, namc1d_uvd)
+!! 3 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_sas)
+!! namsbc_cpl, namtra_qsr, namsbc_rnf,
+!! namsbc_apr, namsbc_ssr, namsbc_alb, namsbc_wave)
+!! 4 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide)
+!! 5 - bottom boundary (nambfr, nambbc, nambbl)
+!! 6 - Tracer (nameos, namtra_adv, namtra_ldf, namtra_ldfeiv, namtra_dmp)
+!! 7 - dynamics (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf)
+!! 8 - Verical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_ddm, namzdf_tmx, namzdf_tmx_new)
+!! 9 - diagnostics (namnc4, namtrd, namspr, namflo, namhsb, namsto)
+!! 10 - miscellaneous (nammpp, namctl)
+!! 11 - Obs & Assim (namobs, nam_asminc)
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+
+!!======================================================================
+!! *** Run management namelists ***
+!!======================================================================
+!! namrun parameters of the run
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namrun ! parameters of the run
+!-----------------------------------------------------------------------
+ nn_no = 0 ! job number (no more used...)
+ cn_exp = "ORCA2" ! experience name
+ nn_it000 = 1 ! first time step
+ nn_itend = 5475 ! last time step (std 5475)
+ nn_date0 = 010101 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
+ nn_time0 = 0 ! initial time of day in hhmm
+ nn_leapy = 0 ! Leap year calendar (1) or not (0)
+ ln_rstart = .false. ! start from rest (F) or from a restart file (T)
+ nn_euler = 1 ! = 0 : start with forward time step if ln_rstart=T
+ nn_rstctl = 0 ! restart control ==> activated only if ln_rstart=T
+ ! ! = 0 nn_date0 read in namelist ; nn_it000 : read in namelist
+ ! ! = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart
+ ! ! = 2 nn_date0 read in restart ; nn_it000 : check consistancy between namelist and restart
+ cn_ocerst_in = "restart" ! suffix of ocean restart name (input)
+ cn_ocerst_indir = "." ! directory from which to read input ocean restarts
+ cn_ocerst_out = "restart" ! suffix of ocean restart name (output)
+ cn_ocerst_outdir= "." ! directory in which to write output ocean restarts
+ ln_iscpl = .false. ! cavity evolution forcing or coupling to ice sheet model
+ nn_istate = 0 ! output the initial state (1) or not (0)
+ ln_rst_list = .false. ! output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F)
+ nn_stock = 5475 ! frequency of creation of a restart file (modulo referenced to 1)
+ nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written
+ nn_write = 5475 ! frequency of write in the output file (modulo referenced to nn_it000)
+ ln_mskland = .false. ! mask land points in NetCDF outputs (costly: + ~15%)
+ ln_cfmeta = .false. ! output additional data to netCDF files required for compliance with the CF metadata standard
+ ln_clobber = .true. ! clobber (overwrite) an existing file
+ nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (works only with iom_nf90 routines)
+/
+!
+!!======================================================================
+!! *** Domain namelists ***
+!!======================================================================
+!! namcfg parameters of the configuration
+!! namdom space and time domain (bathymetry, mesh, timestep)
+!! namwad Wetting and drying (default F)
+!! namtsd data: temperature & salinity
+!! namcrs coarsened grid (for outputs and/or TOP) ("key_crs")
+!! namc1d 1D configuration options ("key_c1d")
+!! namc1d_dyndmp 1D newtonian damping applied on currents ("key_c1d")
+!! namc1d_uvd 1D data (currents) ("key_c1d")
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namcfg ! parameters of the configuration
+!-----------------------------------------------------------------------
+ ln_read_cfg = .false. ! (=T) read the domain configuration file
+ ! ! (=F) user defined configuration ==>>> see usrdef(_...) modules
+ cn_domcfg = "domain_cfg" ! domain configuration filename
+ !
+ ln_write_cfg= .false. ! (=T) create the domain configuration file
+ cn_domcfg_out = "domain_cfg_out" ! newly created domain configuration filename
+ !
+ ln_use_jattr = .false. ! use (T) the file attribute: open_ocean_jstart, if present
+ ! ! in netcdf input files, as the start j-row for reading
+/
+!-----------------------------------------------------------------------
+&namdom ! space and time domain (bathymetry, mesh, timestep)
+!-----------------------------------------------------------------------
+ ln_linssh = .false. ! =T linear free surface ==>> model level are fixed in time
+ nn_closea = 0 ! remove (=0) or keep (=1) closed seas and lakes (ORCA)
+ !
+ nn_msh = 0 ! create (>0) a mesh file or not (=0)
+ rn_isfhmin = 1.00 ! treshold (m) to discriminate grounding ice to floating ice
+ !
+ rn_rdt = 5760. ! time step for the dynamics (and tracer if nn_acc=0)
+ rn_atfp = 0.1 ! asselin time filter parameter
+ !
+ ln_crs = .false. ! Logical switch for coarsening module
+/
+!-----------------------------------------------------------------------
+&namtsd ! data : Temperature & Salinity
+!-----------------------------------------------------------------------
+! ! 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_tem = './ICS/initcd_votemper.nc', -1 ,'votemper', .false. , .true. , 'yearly' , '' , '' , ''
+ sn_sal = './ICS/initcd_vosaline.nc', -1 ,'vosaline', .false. , .true. , 'yearly' , '' , '' , ''
+
+
+ !
+ cn_dir = './' ! root directory for the location of the runoff files
+ ln_tsd_init = .true. ! Initialisation of ocean T & S with T & S input data (T) or not (F)
+ ln_tsd_tradmp = .true. ! damping of ocean T & S toward T & S input data (T) or not (F)
+/
+!-----------------------------------------------------------------------
+&namwad ! Wetting and drying (default F)
+!-----------------------------------------------------------------------
+ ln_wd = .false. ! T/F activation of wetting and drying
+ rn_wdmin1 = 0.1 ! Minimum wet depth on dried cells
+ rn_wdmin2 = 0.01 ! Tolerance of min wet depth on dried cells
+ rn_wdld = 20.0 ! Land elevation below which wetting/drying is allowed
+ nn_wdit = 10 ! Max iterations for W/D limiter
+/
+!-----------------------------------------------------------------------
+&namcrs ! coarsened grid (for outputs and/or TOP) ("key_crs")
+!-----------------------------------------------------------------------
+ nn_factx = 3 ! Reduction factor of x-direction
+ nn_facty = 3 ! Reduction factor of y-direction
+ nn_binref = 0 ! Bin centering preference: NORTH or EQUAT
+ ! 0, coarse grid is binned with preferential treatment of the north fold
+ ! 1, coarse grid is binned with centering at the equator
+ ! Symmetry with nn_facty being odd-numbered. Asymmetry with even-numbered nn_facty.
+ nn_msh_crs = 1 ! create (=1) a mesh file or not (=0)
+ nn_crs_kz = 0 ! 0, MEAN of volume boxes
+ ! 1, MAX of boxes
+ ! 2, MIN of boxes
+ ln_crs_wn = .true. ! wn coarsened (T) or computed using horizontal divergence ( F )
+/
+!-----------------------------------------------------------------------
+&namc1d ! 1D configuration options ("key_c1d")
+!-----------------------------------------------------------------------
+ rn_lat1d = 50 ! Column latitude (default at PAPA station)
+ rn_lon1d = -145 ! Column longitude (default at PAPA station)
+ ln_c1d_locpt= .true. ! Localization of 1D config in a grid (T) or independant point (F)
+/
+!-----------------------------------------------------------------------
+&namc1d_dyndmp ! U & V newtonian damping ("key_c1d")
+!-----------------------------------------------------------------------
+ ln_dyndmp = .false. ! add a damping term (T) or not (F)
+/
+!-----------------------------------------------------------------------
+&namc1d_uvd ! data: U & V currents ("key_c1d")
+!-----------------------------------------------------------------------
+! ! 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_ucur = 'ucurrent' , -1 ,'u_current', .false. , .true. , 'monthly' , '' , 'Ume' , ''
+ sn_vcur = 'vcurrent' , -1 ,'v_current', .false. , .true. , 'monthly' , '' , 'Vme' , ''
+!
+ cn_dir = './' ! root directory for the location of the files
+ ln_uvd_init = .false. ! Initialisation of ocean U & V with U & V input data (T) or not (F)
+ ln_uvd_dyndmp = .false. ! damping of ocean U & V toward U & V input data (T) or not (F)
+/
+
+!!======================================================================
+!! *** Surface Boundary Condition namelists ***
+!!======================================================================
+!! namsbc surface boundary condition
+!! namsbc_flx flux formulation (ln_flx =T)
+!! namsbc_blk Bulk formulae formulation (ln_blk =T)
+!! namsbc_cpl CouPLed formulation ("key_oasis3" )
+!! namsbc_sas Stand-Alone Surface module
+!! namtra_qsr penetrative solar radiation (ln_traqsr =T)
+!! namsbc_rnf river runoffs (ln_rnf =T)
+!! namsbc_isf ice shelf melting/freezing (nn_isf >0)
+!! namsbc_iscpl coupling option between land ice model and ocean
+!! namsbc_apr Atmospheric Pressure (ln_apr_dyn =T)
+!! namsbc_ssr sea surface restoring term (for T and/or S) (ln_ssr =T)
+!! namsbc_alb albedo parameters
+!! namsbc_wave external fields from wave model (ln_wave =T)
+!! namberg iceberg floats (ln_icebergs=T)
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namsbc ! Surface Boundary Condition (surface module)
+!-----------------------------------------------------------------------
+ nn_fsbc = 5 ! frequency of surface boundary condition computation
+ ! (also = the frequency of sea-ice & iceberg model call)
+ ! Type of air-sea fluxes
+ ln_usr = .false. ! user defined formulation (T => check usrdef_sbc)
+ ln_flx = .false. ! flux formulation (T => fill namsbc_flx )
+ ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk )
+ ! Type of coupling (Ocean/Ice/Atmosphere) :
+ ln_cpl = .false. ! atmosphere coupled formulation ( requires key_oasis3 )
+ ln_mixcpl = .false. ! forced-coupled mixed formulation ( requires key_oasis3 )
+ nn_components = 0 ! configuration of the opa-sas OASIS coupling
+ ! =0 no opa-sas OASIS coupling: default single executable configuration
+ ! =1 opa-sas OASIS coupling: multi executable configuration, OPA component
+ ! =2 opa-sas OASIS coupling: multi executable configuration, SAS component
+ nn_limflx = -1 ! LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used)
+ ! =-1 Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled
+ ! = 0 Average per-category fluxes (forced and coupled mode)
+ ! = 1 Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled
+ ! = 2 Redistribute a single flux over categories (coupled mode only)
+ ! Sea-ice :
+ nn_ice = 3 ! =0 no ice boundary condition ,
+ ! =1 use observed ice-cover ,
+ ! =2 to 4 : ice-model used (LIM2, LIM3 or CICE) ("key_lim3", "key_lim2", or "key_cice")
+ nn_ice_embd = 1 ! =0 levitating ice (no mass exchange, concentration/dilution effect)
+ ! =1 levitating ice with mass and salt exchange but no presure effect
+ ! =2 embedded sea-ice (full salt and mass exchanges and pressure)
+ ! Misc. options of sbc :
+ ln_traqsr = .false. ! Light penetration in the ocean (T => fill namtra_qsr)
+ ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave
+ ln_rnf = .true. ! runoffs (T => fill namsbc_rnf)
+ ln_ssr = .true. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr)
+ nn_fwb = 2 ! FreshWater Budget: =0 unchecked
+ ! =1 global mean of e-p-r set to zero at each time step
+ ! =2 annual global mean of e-p-r set to zero
+ ln_apr_dyn = .false. ! Patm gradient added in ocean & ice Eqs. (T => fill namsbc_apr )
+ ln_isf = .false. ! ice shelf (T => fill namsbc_isf)
+ ln_wave = .false. ! Activate coupling with wave (T => fill namsbc_wave)
+ ln_cdgw = .false. ! Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave)
+ ln_sdw = .false. ! Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave)
+ ln_tauoc = .false. ! Activate ocean stress modified by external wave induced stress (T => ln_wave=.true. & fill namsbc_wave)
+ ln_stcor = .false. ! Activate Stokes Coriolis term (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave)
+ nn_lsm = 1 ! =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) ,
+ ! =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field)
+/
+!-----------------------------------------------------------------------
+&namsbc_flx ! surface boundary condition : flux formulation
+!-----------------------------------------------------------------------
+! ! 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_utau = 'utau' , 24 , 'utau' , .false. , .false., 'yearly' , '' , '' , ''
+ sn_vtau = 'vtau' , 24 , 'vtau' , .false. , .false., 'yearly' , '' , '' , ''
+ sn_qtot = 'qtot' , 24 , 'qtot' , .false. , .false., 'yearly' , '' , '' , ''
+ sn_qsr = 'qsr' , 24 , 'qsr' , .false. , .false., 'yearly' , '' , '' , ''
+ sn_emp = 'emp' , 24 , 'emp' , .false. , .false., 'yearly' , '' , '' , ''
+
+ cn_dir = './fluxes/' ! root directory for the location of the flux files
+/
+!-----------------------------------------------------------------------
+&namsbc_blk ! namsbc_blk generic Bulk formula (ln_blk = T)
+!-----------------------------------------------------------------------
+! ! 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_humi= 'ERA5_SPH', 1, 'SPH', .true., .false., 'yearly', 'weights_era5_INDIAN_bicubic.nc', '', 'ERA5_LSM'
+ sn_prec= 'ERA5_TP', 1, 'TP', .true., .false., 'yearly', 'weights_era5_INDIAN_bicubic.nc', '', 'ERA5_LSM'
+ sn_qlw = 'ERA5_STRD', 1, 'STRD', .true., .false., 'yearly', 'weights_era5_INDIAN_bicubic.nc', '', 'ERA5_LSM'
+ sn_qsr= 'ERA5_SSRD', 1, 'SSRD', .true., .false., 'yearly', 'weights_era5_INDIAN_bicubic.nc', '', 'ERA5_LSM'
+ sn_snow= 'ERA5_SF', 1, 'SF', .true., .false., 'yearly', 'weights_era5_INDIAN_bicubic.nc', '', 'ERA5_LSM'
+ sn_tair= 'ERA5_T2M', 1, 'T2M', .true., .false., 'yearly', 'weights_era5_INDIAN_bicubic.nc', '', 'ERA5_LSM'
+ sn_wndi= 'ERA5_U10', 1, 'U10', .true., .false., 'yearly', 'weights_era5_INDIAN_bicubic.nc', 'Uwnd', 'ERA5_LSM'
+ sn_wndj= 'ERA5_V10', 1, 'V10', .true., .false., 'yearly', 'weights_era5_INDIAN_bicubic.nc', 'Vwnd', 'ERA5_LSM'
+ sn_slp= 'SPH_ERA5_SP', 1, 'SP', .true., .false., 'yearly', 'weights_era5_INDIAN_bicubic.nc', '', 'ERA5_LSM'
+! sn_tdif= 'SPH_ERA5_D2M', 1, 'D2M', .true., .false., 'yearly', 'weights_era5_INDIAN_bicubic.nc', '', 'ERA5_LSM'
+ ! ! bulk algorithm :
+ ln_NCAR = .true. ! "NCAR" algorithm (Large and Yeager 2008)
+ ln_COARE_3p0= .false. ! "COARE 3.0" algorithm (Fairall et al. 2003)
+ ln_COARE_3p5= .false. ! "COARE 3.5" algorithm (Edson et al. 2013)
+ ln_ECMWF = .false. ! "ECMWF" algorithm (IFS cycle 31)
+ !
+ cn_dir = './fluxes/' ! root directory for the location of the bulk files
+ ln_taudif = .false. ! HF tau contribution: use "mean of stress module - module of the mean stress" data
+ rn_zqt = 2. ! Air temperature and humidity reference height (m)
+ rn_zu = 10. ! Wind vector reference height (m)
+ rn_pfac = 1. ! multiplicative factor for precipitation (total & snow)
+ rn_efac = 1. ! multiplicative factor for evaporation (0. or 1.)
+ rn_vfac = 1. ! multiplicative factor for ocean/ice velocity
+ ! in the calculation of the wind stress (0.=absolute winds or 1.=relative winds)
+ ln_Cd_L12 = .false. ! Modify the drag ice-atm and oce-atm depending on ice concentration
+ ! This parameterization is from Lupkes et al. (JGR 2012)
+/
+!-----------------------------------------------------------------------
+&namsbc_cpl ! coupled ocean/atmosphere model ("key_oasis3")
+!-----------------------------------------------------------------------
+! ! description ! multiple ! vector ! vector ! vector !
+! ! ! categories ! reference ! orientation ! grids !
+! send
+ sn_snd_temp = 'weighted oce and ice' , 'no' , '' , '' , ''
+ sn_snd_alb = 'weighted ice' , 'no' , '' , '' , ''
+ sn_snd_thick = 'none' , 'no' , '' , '' , ''
+ sn_snd_crt = 'none' , 'no' , 'spherical' , 'eastward-northward' , 'T'
+ sn_snd_co2 = 'coupled' , 'no' , '' , '' , ''
+ sn_snd_crtw = 'none' , 'no' , '' , '' , 'U,V'
+ sn_snd_ifrac = 'none' , 'no' , '' , '' , ''
+ sn_snd_wlev = 'coupled' , 'no' , '' , '' , ''
+! receive
+ sn_rcv_w10m = 'none' , 'no' , '' , '' , ''
+ sn_rcv_taumod = 'coupled' , 'no' , '' , '' , ''
+ sn_rcv_tau = 'oce only' , 'no' , 'cartesian' , 'eastward-northward', 'U,V'
+ sn_rcv_dqnsdt = 'coupled' , 'no' , '' , '' , ''
+ sn_rcv_qsr = 'oce and ice' , 'no' , '' , '' , ''
+ sn_rcv_qns = 'oce and ice' , 'no' , '' , '' , ''
+ sn_rcv_emp = 'conservative' , 'no' , '' , '' , ''
+ sn_rcv_rnf = 'coupled' , 'no' , '' , '' , ''
+ sn_rcv_cal = 'coupled' , 'no' , '' , '' , ''
+ sn_rcv_co2 = 'coupled' , 'no' , '' , '' , ''
+ sn_rcv_hsig = 'none' , 'no' , '' , '' , ''
+ sn_rcv_iceflx = 'none' , 'no' , '' , '' , ''
+ sn_rcv_mslp = 'none' , 'no' , '' , '' , ''
+ sn_rcv_phioc = 'none' , 'no' , '' , '' , ''
+ sn_rcv_sdrfx = 'none' , 'no' , '' , '' , ''
+ sn_rcv_sdrfy = 'none' , 'no' , '' , '' , ''
+ sn_rcv_wper = 'none' , 'no' , '' , '' , ''
+ sn_rcv_wnum = 'none' , 'no' , '' , '' , ''
+ sn_rcv_wstrf = 'none' , 'no' , '' , '' , ''
+ sn_rcv_wdrag = 'none' , 'no' , '' , '' , ''
+!
+ nn_cplmodel = 1 ! Maximum number of models to/from which NEMO is potentialy sending/receiving data
+ ln_usecplmask = .false. ! use a coupling mask file to merge data received from several models
+ ! ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel)
+/
+!-----------------------------------------------------------------------
+&namsbc_sas ! Stand Alone Surface boundary condition
+!-----------------------------------------------------------------------
+! ! 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 !
+ l_sasread = .TRUE. ! Read fields in a file if .TRUE. , or initialize to 0. in sbcssm.F90 if .FALSE.
+ sn_usp = 'sas_grid_U', 120 , 'vozocrtx', .true. , .true. , 'yearly' , '' , '' , ''
+ sn_vsp = 'sas_grid_V', 120 , 'vomecrty', .true. , .true. , 'yearly' , '' , '' , ''
+ sn_tem = 'sas_grid_T', 120 , 'sosstsst', .true. , .true. , 'yearly' , '' , '' , ''
+ sn_sal = 'sas_grid_T', 120 , 'sosaline', .true. , .true. , 'yearly' , '' , '' , ''
+ sn_ssh = 'sas_grid_T', 120 , 'sossheig', .true. , .true. , 'yearly' , '' , '' , ''
+ sn_e3t = 'sas_grid_T', 120 , 'e3t_m' , .true. , .true. , 'yearly' , '' , '' , ''
+ sn_frq = 'sas_grid_T', 120 , 'frq_m' , .true. , .true. , 'yearly' , '' , '' , ''
+
+ ln_3d_uve = .true. ! specify whether we are supplying a 3D u,v and e3 field
+ ln_read_frq = .false. ! specify whether we must read frq or not
+ cn_dir = './' ! root directory for the location of the bulk files are
+/
+!-----------------------------------------------------------------------
+&namtra_qsr ! penetrative solar radiation (ln_traqsr=T)
+!-----------------------------------------------------------------------
+! ! 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_chl ='chlorophyll', -1 , 'CHLA' , .true. , .true. , 'yearly' , '' , '' , ''
+
+ cn_dir = './' ! root directory for the location of the runoff files
+ ln_qsr_rgb = .true. ! RGB (Red-Green-Blue) light penetration
+ ln_qsr_2bd = .false. ! 2 bands light penetration
+ ln_qsr_bio = .false. ! bio-model light penetration
+ nn_chldta = 1 ! RGB : Chl data (=1) or cst value (=0)
+ rn_abs = 0.58 ! RGB & 2 bands: fraction of light (rn_si1)
+ rn_si0 = 0.35 ! RGB & 2 bands: shortess depth of extinction
+ rn_si1 = 23.0 ! 2 bands: longest depth of extinction
+ ln_qsr_ice = .false. ! light penetration for ice-model LIM3
+/
+!-----------------------------------------------------------------------
+&namsbc_rnf ! runoffs namelist surface boundary condition (ln_rnf=T)
+!-----------------------------------------------------------------------
+! ! 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_rnf = 'river_test' , -1 , 'rorunoff', .true. , .true. , 'yearly' , '' , '' , ''
+ sn_cnf = 'river_test' , 0 , 'socoefr0', .false. , .true. , 'yearly' , '' , '' , ''
+ sn_s_rnf = 'runoffs' , 24 , 'rosaline', .true. , .true. , 'yearly' , '' , '' , ''
+ sn_t_rnf = 'runoffs' , 24 , 'rotemper', .true. , .true. , 'yearly' , '' , '' , ''
+ sn_dep_rnf = 'runoffs' , 0 , 'rodepth' , .false. , .true. , 'yearly' , '' , '' , ''
+
+ cn_dir = './' ! root directory for the location of the runoff files
+ ln_rnf_mouth= .false. ! specific treatment at rivers mouths
+ rn_hrnf = 15.e0 ! depth over which enhanced vertical mixing is used (ln_rnf_mouth=T)
+ rn_avt_rnf = 1.e-3 ! value of the additional vertical mixing coef. [m2/s] (ln_rnf_mouth=T)
+ rn_rfact = 1.e0 ! multiplicative factor for runoff
+ ln_rnf_depth= .false. ! read in depth information for runoff
+ ln_rnf_tem = .false. ! read in temperature information for runoff
+ ln_rnf_sal = .false. ! read in salinity information for runoff
+ ln_rnf_depth_ini = .false. ! compute depth at initialisation from runoff file
+ rn_rnf_max = 5.735e-4 ! max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true )
+ rn_dep_max = 150. ! depth over which runoffs is spread ( ln_rnf_depth_ini = .true )
+ nn_rnf_depth_file = 0 ! create (=1) a runoff depth file or not (=0)
+/
+!-----------------------------------------------------------------------
+&namsbc_isf ! Top boundary layer (ISF) (nn_isf >0)
+!-----------------------------------------------------------------------
+! ! 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 !
+! nn_isf == 4
+ sn_fwfisf = 'rnfisf' , -12 ,'sowflisf', .false. , .true. , 'yearly' , '' , '' , ''
+! nn_isf == 3
+ sn_rnfisf = 'rnfisf' , -12 ,'sofwfisf', .false. , .true. , 'yearly' , '' , '' , ''
+! nn_isf == 2 and 3
+ sn_depmax_isf='rnfisf' , -12 ,'sozisfmax', .false. , .true. , 'yearly' , '' , '' , ''
+ sn_depmin_isf='rnfisf' , -12 ,'sozisfmin', .false. , .true. , 'yearly' , '' , '' , ''
+! nn_isf == 2
+ sn_Leff_isf = 'rnfisf' , -12 ,'Leff' , .false. , .true. , 'yearly' , '' , '' , ''
+!
+! for all case
+ nn_isf = 1 ! ice shelf melting/freezing
+ ! 1 = presence of ISF 2 = bg03 parametrisation
+ ! 3 = rnf file for isf 4 = ISF fwf specified
+ ! option 1 and 4 need ln_isfcav = .true. (domzgr)
+! only for nn_isf = 1 or 2
+ rn_gammat0 = 1.e-4 ! gammat coefficient used in blk formula
+ rn_gammas0 = 1.e-4 ! gammas coefficient used in blk formula
+! only for nn_isf = 1 or 4
+ rn_hisf_tbl = 30. ! thickness of the top boundary layer (Losh et al. 2008)
+ ! ! 0 => thickness of the tbl = thickness of the first wet cell
+! only for nn_isf = 1
+ nn_isfblk = 1 ! 1 ISOMIP like: 2 equations formulation (Hunter et al., 2006)
+ ! ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015)
+ nn_gammablk = 1 ! 0 = cst Gammat (= gammat/s)
+ ! ! 1 = velocity dependend Gamma (u* * gammat/s) (Jenkins et al. 2010)
+ ! ! 2 = velocity and stability dependent Gamma (Holland et al. 1999)
+/
+!-----------------------------------------------------------------------
+&namsbc_iscpl ! land ice / ocean coupling option
+!-----------------------------------------------------------------------
+ nn_drown = 10 ! number of iteration of the extrapolation loop (fill the new wet cells)
+ ln_hsb = .false. ! activate conservation module (conservation exact after a time of rn_fiscpl)
+ nn_fiscpl = 43800 ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency)
+/
+!-----------------------------------------------------------------------
+&namsbc_apr ! Atmospheric pressure used as ocean forcing (ln_apr_dyn =T)
+!-----------------------------------------------------------------------
+! ! 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_apr = 'patm' , -1 ,'somslpre', .true. , .true. , 'yearly' , '' , '' , ''
+
+ cn_dir = './' ! root directory for the location of the bulk files
+ rn_pref = 101000. ! reference atmospheric pressure [N/m2]/
+ ln_ref_apr = .false. ! ref. pressure: global mean Patm (T) or a constant (F)
+ ln_apr_obc = .false. ! inverse barometer added to OBC ssh data
+/
+!-----------------------------------------------------------------------
+&namsbc_ssr ! surface boundary condition : sea surface restoring (ln_ssr=T)
+!-----------------------------------------------------------------------
+! ! 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_sst = 'sst_data', 24 , 'sst' , .false. , .false., 'yearly' , '' , '' , ''
+ sn_sss = 'sss_data', -1 , 'sss' , .true. , .true. , 'yearly' , '' , '' , ''
+
+ cn_dir = './' ! root directory for the location of the runoff files
+ nn_sstr = 0 ! add a retroaction term in the surface heat flux (=1) or not (=0)
+ nn_sssr = 2 ! add a damping term in the surface freshwater flux (=2)
+ ! or to SSS only (=1) or no damping term (=0)
+ rn_dqdt = -40. ! magnitude of the retroaction on temperature [W/m2/K]
+ rn_deds = -166.67 ! magnitude of the damping on salinity [mm/day]
+ ln_sssr_bnd = .true. ! flag to bound erp term (associated with nn_sssr=2)
+ rn_sssr_bnd = 4.e0 ! ABS(Max/Min) value of the damping erp term [mm/day]
+/
+!-----------------------------------------------------------------------
+&namsbc_alb ! albedo parameters
+!-----------------------------------------------------------------------
+ nn_ice_alb = 1 ! parameterization of ice/snow albedo
+ ! 0: Shine & Henderson-Sellers (JGR 1985), giving clear-sky albedo
+ ! 1: "home made" based on Brandt et al. (JClim 2005) and Grenfell & Perovich (JGR 2004),
+ ! giving cloud-sky albedo
+ rn_alb_sdry = 0.85 ! dry snow albedo : 0.80 (nn_ice_alb = 0); 0.85 (nn_ice_alb = 1); obs 0.85-0.87 (cloud-sky)
+ rn_alb_smlt = 0.75 ! melting snow albedo : 0.65 ( '' ) ; 0.75 ( '' ) ; obs 0.72-0.82 ( '' )
+ rn_alb_idry = 0.60 ! dry ice albedo : 0.72 ( '' ) ; 0.60 ( '' ) ; obs 0.54-0.65 ( '' )
+ rn_alb_imlt = 0.50 ! bare puddled ice albedo : 0.53 ( '' ) ; 0.50 ( '' ) ; obs 0.49-0.58 ( '' )
+/
+!-----------------------------------------------------------------------
+&namsbc_wave ! External fields from wave model (ln_wave=T)
+!-----------------------------------------------------------------------
+! ! 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_cdg = 'sdw_wave' , 1 , 'drag_coeff' , .true. , .false. , 'daily' , '' , '' , ''
+ sn_usd = 'sdw_wave' , 1 , 'u_sd2d' , .true. , .false. , 'daily' , '' , '' , ''
+ sn_vsd = 'sdw_wave' , 1 , 'v_sd2d' , .true. , .false. , 'daily' , '' , '' , ''
+ sn_hsw = 'sdw_wave' , 1 , 'hs' , .true. , .false. , 'daily' , '' , '' , ''
+ sn_wmp = 'sdw_wave' , 1 , 'wmp' , .true. , .false. , 'daily' , '' , '' , ''
+ sn_wnum = 'sdw_wave' , 1 , 'wave_num' , .true. , .false. , 'daily' , '' , '' , ''
+ sn_tauoc = 'sdw_wave' , 1 , 'wave_stress', .true. , .false. , 'daily' , '' , '' , ''
+!
+ cn_dir = './' ! root directory for the location of drag coefficient files
+/
+!-----------------------------------------------------------------------
+&namberg ! iceberg parameters (default: No iceberg)
+!-----------------------------------------------------------------------
+ ln_icebergs = .false. ! iceberg floats or not
+ ln_bergdia = .true. ! Calculate budgets
+ nn_verbose_level = 1 ! Turn on more verbose output if level > 0
+ nn_verbose_write = 15 ! Timesteps between verbose messages
+ nn_sample_rate = 1 ! Timesteps between sampling for trajectory storage
+ ! Initial mass required for an iceberg of each class
+ rn_initial_mass = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11
+ ! Proportion of calving mass to apportion to each class
+ rn_distribution = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02
+ ! Ratio between effective and real iceberg mass (non-dim)
+ ! i.e. number of icebergs represented at a point
+ rn_mass_scaling = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1
+ ! thickness of newly calved bergs (m)
+ rn_initial_thickness = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250.
+ rn_rho_bergs = 850. ! Density of icebergs
+ rn_LoW_ratio = 1.5 ! Initial ratio L/W for newly calved icebergs
+ ln_operator_splitting = .true. ! Use first order operator splitting for thermodynamics
+ rn_bits_erosion_fraction = 0. ! Fraction of erosion melt flux to divert to bergy bits
+ rn_sicn_shift = 0. ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1)
+ ln_passive_mode = .false. ! iceberg - ocean decoupling
+ nn_test_icebergs = 10 ! Create test icebergs of this class (-1 = no)
+ ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2)
+ rn_test_box = 108.0, 116.0, -66.0, -58.0
+ rn_speed_limit = 0. ! CFL speed limit for a berg
+
+! ! 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_icb = 'calving', -1 , 'calvingmask', .true. , .true. , 'yearly' , '' , '' , ''
+
+ cn_dir = './'
+/
+
+!!======================================================================
+!! *** Lateral boundary condition ***
+!!======================================================================
+!! namlbc lateral momentum boundary condition
+!! namagrif agrif nested grid ( read by child model only ) ("key_agrif")
+!! nam_tide Tidal forcing
+!! nambdy Unstructured open boundaries
+!! nambdy_dta Unstructured open boundaries - external data
+!! nambdy_tide tidal forcing at open boundaries
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namlbc ! lateral momentum boundary condition
+!-----------------------------------------------------------------------
+ ! ! free slip ! partial slip ! no slip ! strong slip
+ rn_shlat = 2. ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat
+ ln_vorlat = .false. ! consistency of vorticity boundary condition with analytical Eqs.
+/
+!-----------------------------------------------------------------------
+&namagrif ! AGRIF zoom ("key_agrif")
+!-----------------------------------------------------------------------
+ nn_cln_update = 3 ! baroclinic update frequency
+ ln_spc_dyn = .true. ! use 0 as special value for dynamics
+ rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [m2/s]
+ rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [m2/s]
+ ln_chk_bathy = .FALSE. !
+/
+!-----------------------------------------------------------------------
+&nam_tide ! tide parameters
+!-----------------------------------------------------------------------
+ ln_tide = .false.
+ ln_tide_pot = .true. ! use tidal potential forcing
+ ln_tide_ramp= .false. !
+ rdttideramp = 0. !
+ clname(1) = 'DUMMY' ! name of constituent - all tidal components must be set in namelist_cfg
+/
+!-----------------------------------------------------------------------
+&nambdy ! unstructured open boundaries
+!-----------------------------------------------------------------------
+ ln_bdy = .false. ! Use unstructured open boundaries
+ nb_bdy = 0 ! number of open boundary sets
+ ln_coords_file = .true. ! =T : read bdy coordinates from file
+ cn_coords_file = 'coordinates.bdy.nc' ! bdy coordinates files
+ ln_mask_file = .false. ! =T : read mask from file
+ cn_mask_file = '' ! name of mask file (if ln_mask_file=.TRUE.)
+ cn_dyn2d = 'none' !
+ nn_dyn2d_dta = 0 ! = 0, bdy data are equal to the initial state
+ ! = 1, bdy data are read in 'bdydata .nc' files
+ ! = 2, use tidal harmonic forcing data from files
+ ! = 3, use external data AND tidal harmonic forcing
+ cn_dyn3d = 'none' !
+ nn_dyn3d_dta = 0 ! = 0, bdy data are equal to the initial state
+ ! = 1, bdy data are read in 'bdydata .nc' files
+ cn_tra = 'none' !
+ nn_tra_dta = 0 ! = 0, bdy data are equal to the initial state
+ ! = 1, bdy data are read in 'bdydata .nc' files
+ cn_ice_lim = 'none' !
+ nn_ice_lim_dta = 0 ! = 0, bdy data are equal to the initial state
+ ! = 1, bdy data are read in 'bdydata .nc' files
+ rn_ice_tem = 270. ! lim3 only: arbitrary temperature of incoming sea ice
+ rn_ice_sal = 10. ! lim3 only: -- salinity --
+ rn_ice_age = 30. ! lim3 only: -- age --
+
+ ln_tra_dmp =.false. ! open boudaries conditions for tracers
+ ln_dyn3d_dmp =.false. ! open boundary condition for baroclinic velocities
+ rn_time_dmp = 1. ! Damping time scale in days
+ rn_time_dmp_out = 1. ! Outflow damping time scale
+ nn_rimwidth = 10 ! width of the relaxation zone
+ ln_vol = .false. ! total volume correction (see nn_volctl parameter)
+ nn_volctl = 1 ! = 0, the total water flux across open boundaries is zero
+ nb_jpk_bdy = -1 ! number of levels in the bdy data (set < 0 if consistent with planned run)
+/
+!-----------------------------------------------------------------------
+&nambdy_dta ! open boundaries - external data
+!-----------------------------------------------------------------------
+! ! 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 !
+ bn_ssh = 'INDIAN_bdyT_u2d', 1 , 'sossheig', .true. , .false. , 'monthly' , '' , '' , ''
+ bn_u2d = 'INDIAN_bdyU_u2d', 1 , 'vobtcrtx', .true. , .false. , 'monthly' , '' , '' , ''
+ bn_v2d = 'INDIAN_bdyV_u2d', 1 , 'vobtcrty', .true. , .false. , 'monthly' , '' , '' , ''
+ bn_u3d = 'INDIAN_bdyU_u3d', 1 , 'vozocrtx', .true. , .false. , 'monthly' , '' , '' , ''
+ bn_v3d = 'INDIAN_bdyV_u3d', 1 , 'vomecrty', .true. , .false. , 'monthly' , '' , '' , ''
+ bn_tem = 'INDIAN_bdyT_tra', 1 , 'votemper', .true. , .false. , 'monthly' , '' , '' , ''
+ bn_sal = 'INDIAN_bdyT_tra', 1 , 'vosaline', .true. , .false. , 'monthly' , '' , '' , ''
+! for lim2
+! bn_frld = 'amm12_bdyT_ice', 24 , 'ileadfra', .true. , .false. , 'daily' , '' , '' , ''
+! bn_hicif = 'amm12_bdyT_ice', 24 , 'iicethic', .true. , .false. , 'daily' , '' , '' , ''
+! bn_hsnif = 'amm12_bdyT_ice', 24 , 'isnowthi', .true. , .false. , 'daily' , '' , '' , ''
+! for lim3
+! bn_a_i = 'amm12_bdyT_ice', 24 , 'ileadfra', .true. , .false. , 'daily' , '' , '' , ''
+! bn_ht_i = 'amm12_bdyT_ice', 24 , 'iicethic', .true. , .false. , 'daily' , '' , '' , ''
+! bn_ht_s = 'amm12_bdyT_ice', 24 , 'isnowthi', .true. , .false. , 'daily' , '' , '' , ''
+
+ cn_dir = 'OBC/' ! root directory for the location of the bulk files
+ ln_full_vel = .false. !
+/
+!-----------------------------------------------------------------------
+&nambdy_tide ! tidal forcing at open boundaries
+!-----------------------------------------------------------------------
+ filtide = 'bdydta/amm12_bdytide_' ! file name root of tidal forcing files
+ ln_bdytide_2ddta = .false. !
+ ln_bdytide_conj = .false. !
+/
+
+!!======================================================================
+!! *** Bottom boundary condition ***
+!!======================================================================
+!! nambfr bottom friction
+!! nambbc bottom temperature boundary condition
+!! nambbl bottom boundary layer scheme ("key_trabbl")
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&nambfr ! bottom friction (default: linear)
+!-----------------------------------------------------------------------
+ nn_bfr = 1 ! type of bottom friction : = 0 : free slip, = 1 : linear friction
+ ! = 2 : nonlinear friction
+ rn_bfri1 = 4.e-4 ! bottom drag coefficient (linear case)
+ rn_bfri2 = 1.e-3 ! bottom drag coefficient (non linear case). Minimum coeft if ln_loglayer=T
+ rn_bfri2_max= 1.e-1 ! max. bottom drag coefficient (non linear case and ln_loglayer=T)
+ rn_bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m2/s2)
+ rn_bfrz0 = 3.e-3 ! bottom roughness [m] if ln_loglayer=T
+ ln_bfr2d = .false. ! horizontal variation of the bottom friction coef (read a 2D mask file )
+ rn_bfrien = 50. ! local multiplying factor of bfr (ln_bfr2d=T)
+ rn_tfri1 = 4.e-4 ! top drag coefficient (linear case)
+ rn_tfri2 = 2.5e-3 ! top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T
+ rn_tfri2_max= 1.e-1 ! max. top drag coefficient (non linear case and ln_loglayer=T)
+ rn_tfeb2 = 0.0 ! top turbulent kinetic energy background (m2/s2)
+ rn_tfrz0 = 3.e-3 ! top roughness [m] if ln_loglayer=T
+ ln_tfr2d = .false. ! horizontal variation of the top friction coef (read a 2D mask file )
+ rn_tfrien = 50. ! local multiplying factor of tfr (ln_tfr2d=T)
+
+ ln_bfrimp = .true. ! implicit bottom friction (requires ln_zdfexp = .false. if true)
+ ln_loglayer = .false. ! logarithmic formulation (non linear case)
+/
+!-----------------------------------------------------------------------
+&nambbc ! bottom temperature boundary condition (default: NO)
+!-----------------------------------------------------------------------
+! ! 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_qgh ='geothermal_heating.nc', -12. , 'heatflow', .false. , .true. , 'yearly' , '' , '' , ''
+ !
+ ln_trabbc = .false. ! Apply a geothermal heating at the ocean bottom
+ nn_geoflx = 2 ! geothermal heat flux: = 0 no flux
+ ! = 1 constant flux
+ ! = 2 variable flux (read in geothermal_heating.nc in mW/m2)
+ rn_geoflx_cst = 86.4e-3 ! Constant value of geothermal heat flux [W/m2]
+ cn_dir = './' ! root directory for the location of the runoff files
+/
+!-----------------------------------------------------------------------
+&nambbl ! bottom boundary layer scheme ("key_trabbl")
+!-----------------------------------------------------------------------
+ nn_bbl_ldf = 1 ! diffusive bbl (=1) or not (=0)
+ nn_bbl_adv = 0 ! advective bbl (=1/2) or not (=0)
+ rn_ahtbbl = 1000. ! lateral mixing coefficient in the bbl [m2/s]
+ rn_gambbl = 10. ! advective bbl coefficient [s]
+/
+
+!!======================================================================
+!! Tracer (T & S ) namelists
+!!======================================================================
+!! nameos equation of state
+!! namtra_adv advection scheme
+!! namtra_adv_mle mixed layer eddy param. (Fox-Kemper param.)
+!! namtra_ldf lateral diffusion scheme
+!! namtra_ldfeiv eddy induced velocity param.
+!! namtra_dmp T & S newtonian damping
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&nameos ! ocean physical parameters
+!-----------------------------------------------------------------------
+ ln_teos10 = .false. ! = Use TEOS-10 equation of state
+ ln_eos80 = .false. ! = Use EOS80 equation of state
+ ln_seos = .false. ! = Use simplified equation of state (S-EOS)
+ !
+ ! ! S-EOS coefficients (ln_seos=T):
+ ! ! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS
+ rn_a0 = 1.6550e-1 ! thermal expension coefficient
+ rn_b0 = 7.6554e-1 ! saline expension coefficient
+ rn_lambda1 = 5.9520e-2 ! cabbeling coeff in T^2 (=0 for linear eos)
+ rn_lambda2 = 7.4914e-4 ! cabbeling coeff in S^2 (=0 for linear eos)
+ rn_mu1 = 1.4970e-4 ! thermobaric coeff. in T (=0 for linear eos)
+ rn_mu2 = 1.1090e-5 ! thermobaric coeff. in S (=0 for linear eos)
+ rn_nu = 2.4341e-3 ! cabbeling coeff in T*S (=0 for linear eos)
+/
+!-----------------------------------------------------------------------
+&namtra_adv ! advection scheme for tracer (default: NO advection)
+!-----------------------------------------------------------------------
+ ln_traadv_cen = .false. ! 2nd order centered scheme
+ nn_cen_h = 4 ! =2/4, horizontal 2nd order CEN / 4th order CEN
+ nn_cen_v = 4 ! =2/4, vertical 2nd order CEN / 4th order COMPACT
+ ln_traadv_fct = .false. ! FCT scheme
+ nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order
+ nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order
+ nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping
+ ! ! (number of sub-timestep = nn_fct_zts)
+ ln_traadv_mus = .false. ! MUSCL scheme
+ ln_mus_ups = .false. ! use upstream scheme near river mouths
+ ln_traadv_ubs = .false. ! UBS scheme
+ nn_ubs_v = 2 ! =2 , vertical 2nd order FCT / COMPACT 4th order
+ ln_traadv_qck = .true. ! QUICKEST scheme
+/
+!-----------------------------------------------------------------------
+&namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param) (default: NO)
+!-----------------------------------------------------------------------
+ ln_mle = .false. ! (T) use the Mixed Layer Eddy (MLE) parameterisation
+ rn_ce = 0.06 ! magnitude of the MLE (typical value: 0.06 to 0.08)
+ nn_mle = 1 ! MLE type: =0 standard Fox-Kemper ; =1 new formulation
+ rn_lf = 5.e+3 ! typical scale of mixed layer front (meters) (case rn_mle=0)
+ rn_time = 172800. ! time scale for mixing momentum across the mixed layer (seconds) (case rn_mle=0)
+ rn_lat = 20. ! reference latitude (degrees) of MLE coef. (case rn_mle=1)
+ nn_mld_uv = 0 ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max)
+ nn_conv = 0 ! =1 no MLE in case of convection ; =0 always MLE
+ rn_rho_c_mle= 0.01 ! delta rho criterion used to calculate MLD for FK
+/
+!-----------------------------------------------------------------------
+&namtra_ldf ! lateral diffusion scheme for tracers (default: NO diffusion)
+!-----------------------------------------------------------------------
+ ! ! Operator type:
+ ! ! no diffusion: set ln_traldf_lap=..._blp=F
+ ln_traldf_lap = .false. ! laplacian operator
+ ln_traldf_blp = .false. ! bilaplacian operator
+ !
+ ! ! Direction of action:
+ ln_traldf_lev = .false. ! iso-level
+ ln_traldf_hor = .false. ! horizontal (geopotential)
+ ln_traldf_iso = .false. ! iso-neutral (standard operator)
+ ln_traldf_triad = .false. ! iso-neutral (triad operator)
+ !
+ ! ! iso-neutral options:
+ ln_traldf_msc = .false. ! Method of Stabilizing Correction (both operators)
+ rn_slpmax = 0.01 ! slope limit (both operators)
+ ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only)
+ rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only)
+ ln_botmix_triad = .false. ! lateral mixing on bottom (triad only)
+ !
+ ! ! Coefficients:
+ nn_aht_ijk_t = 0 ! space/time variation of eddy coef
+ ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file
+ ! ! = 0 constant
+ ! ! = 10 F(k) =ldf_c1d
+ ! ! = 20 F(i,j) =ldf_c2d
+ ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation
+ ! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d
+ ! ! = 31 F(i,j,k,t)=F(local velocity and grid-spacing)
+ rn_aht_0 = 2000. ! lateral eddy diffusivity (lap. operator) [m2/s]
+ rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s]
+/
+!-----------------------------------------------------------------------
+&namtra_ldfeiv ! eddy induced velocity param. (default: NO)
+!-----------------------------------------------------------------------
+ ln_ldfeiv =.false. ! use eddy induced velocity parameterization
+ ln_ldfeiv_dia =.false. ! diagnose eiv stream function and velocities
+ rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s]
+ nn_aei_ijk_t = 21 ! space/time variation of the eiv coeficient
+ ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file
+ ! ! = 0 constant
+ ! ! = 10 F(k) =ldf_c1d
+ ! ! = 20 F(i,j) =ldf_c2d
+ ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation
+ ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d
+/
+!-----------------------------------------------------------------------
+&namtra_dmp ! tracer: T & S newtonian damping (default: NO)
+!-----------------------------------------------------------------------
+ ln_tradmp = .false. ! add a damping termn (T) or not (F)
+ nn_zdmp = 0 ! vertical shape =0 damping throughout the water column
+ ! =1 no damping in the mixing layer (kz criteria)
+ ! =2 no damping in the mixed layer (rho crieria)
+ cn_resto ='resto.nc' ! Name of file containing restoration coeff. field (use dmp_tools to create this)
+/
+
+!!======================================================================
+!! *** Dynamics namelists ***
+!!======================================================================
+!! namdyn_adv formulation of the momentum advection
+!! namdyn_vor advection scheme
+!! namdyn_hpg hydrostatic pressure gradient
+!! namdyn_spg surface pressure gradient
+!! namdyn_ldf lateral diffusion scheme
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namdyn_adv ! formulation of the momentum advection (default: vector form)
+!-----------------------------------------------------------------------
+ ln_dynadv_vec = .true. ! vector form (T) or flux form (F)
+ nn_dynkeg = 0 ! scheme for grad(KE): =0 C2 ; =1 Hollingsworth correction
+ ln_dynadv_cen2= .false. ! flux form - 2nd order centered scheme
+ ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme
+ ln_dynzad_zts = .false. ! Use (T) sub timestepping for vertical momentum advection
+/
+!-----------------------------------------------------------------------
+&nam_vvl ! vertical coordinate options (default: zstar)
+!-----------------------------------------------------------------------
+ ln_vvl_zstar = .true. ! zstar vertical coordinate
+ ln_vvl_ztilde = .false. ! ztilde vertical coordinate: only high frequency variations
+ ln_vvl_layer = .false. ! full layer vertical coordinate
+ ln_vvl_ztilde_as_zstar = .false. ! ztilde vertical coordinate emulating zstar
+ ln_vvl_zstar_at_eqtor = .false. ! ztilde near the equator
+ rn_ahe3 = 0.0e0 ! thickness diffusion coefficient
+ rn_rst_e3t = 30.e0 ! ztilde to zstar restoration timescale [days]
+ rn_lf_cutoff = 5.0e0 ! cutoff frequency for low-pass filter [days]
+ rn_zdef_max = 0.9e0 ! maximum fractional e3t deformation
+ ln_vvl_dbg = .true. ! debug prints (T/F)
+/
+!-----------------------------------------------------------------------
+&namdyn_vor ! Vorticity / Coriolis scheme (default: NO)
+!-----------------------------------------------------------------------
+ ln_dynvor_ene = .false. ! enstrophy conserving scheme
+ ln_dynvor_ens = .false. ! energy conserving scheme
+ ln_dynvor_mix = .false. ! mixed scheme
+ ln_dynvor_een = .false. ! energy & enstrophy scheme
+ nn_een_e3f = 1 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1)
+ ln_dynvor_msk = .false. ! vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) ! PLEASE DO NOT ACTIVATE
+/
+!-----------------------------------------------------------------------
+&namdyn_hpg ! Hydrostatic pressure gradient option (default: zps)
+!-----------------------------------------------------------------------
+ ln_hpg_zco = .false. ! z-coordinate - full steps
+ ln_hpg_zps = .false. ! z-coordinate - partial steps (interpolation)
+ ln_hpg_sco = .false. ! s-coordinate (standard jacobian formulation)
+ ln_hpg_isf = .false. ! s-coordinate (sco ) adapted to isf
+ ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial)
+ ln_hpg_prj = .false. ! s-coordinate (Pressure Jacobian scheme)
+/
+!-----------------------------------------------------------------------
+&namdyn_spg ! surface pressure gradient (default: NO)
+!-----------------------------------------------------------------------
+ ln_dynspg_exp = .false. ! explicit free surface
+ ln_dynspg_ts = .false. ! split-explicit free surface
+ ln_bt_fw = .true. ! Forward integration of barotropic Eqs.
+ ln_bt_av = .true. ! Time filtering of barotropic variables
+ nn_bt_flt = 1 ! Time filter choice = 0 None
+ ! ! = 1 Boxcar over nn_baro sub-steps
+ ! ! = 2 Boxcar over 2*nn_baro " "
+ ln_bt_auto = .true. ! Number of sub-step defined from:
+ rn_bt_cmax = 0.8 ! =T : the Maximum Courant Number allowed
+ nn_baro = 30 ! =F : the number of sub-step in rn_rdt seconds
+/
+!-----------------------------------------------------------------------
+&namdyn_ldf ! lateral diffusion on momentum (default: NO)
+!-----------------------------------------------------------------------
+ ! ! Type of the operator :
+ ! ! no diffusion: set ln_dynldf_lap=..._blp=F
+ ln_dynldf_lap = .false. ! laplacian operator
+ ln_dynldf_blp = .false. ! bilaplacian operator
+ ! ! Direction of action :
+ ln_dynldf_lev = .false. ! iso-level
+ ln_dynldf_hor = .false. ! horizontal (geopotential)
+ ln_dynldf_iso = .false. ! iso-neutral
+ ! ! Coefficient
+ nn_ahm_ijk_t = 0 ! space/time variation of eddy coef
+ ! ! =-30 read in eddy_viscosity_3D.nc file
+ ! ! =-20 read in eddy_viscosity_2D.nc file
+ ! ! = 0 constant
+ ! ! = 10 F(k)=c1d
+ ! ! = 20 F(i,j)=F(grid spacing)=c2d
+ ! ! = 30 F(i,j,k)=c2d*c1d
+ ! ! = 31 F(i,j,k)=F(grid spacing and local velocity)
+ ! ! = 32 F(i,j,k)=F(local gridscale and deformation rate)
+ ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km)
+ rn_ahm_0 = 40000. ! horizontal laplacian eddy viscosity [m2/s]
+ rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s]
+ rn_bhm_0 = 1.e+12 ! horizontal bilaplacian eddy viscosity [m4/s]
+ ! ! Smagorinsky settings (nn_ahm_ijk_t = 32) :
+ rn_csmc = 3.5 ! Smagorinsky constant of proportionality
+ rn_minfac = 1.0 ! multiplier of theorectical lower limit
+ rn_maxfac = 1.0 ! multiplier of theorectical upper limit
+/
+
+!!======================================================================
+!! Tracers & Dynamics vertical physics namelists
+!!======================================================================
+!! namzdf vertical physics
+!! namzdf_ric richardson number dependent vertical mixing ("key_zdfric")
+!! namzdf_tke TKE dependent vertical mixing ("key_zdftke")
+!! namzdf_gls GLS vertical mixing ("key_zdfgls")
+!! namzdf_ddm double diffusive mixing parameterization ("key_zdfddm")
+!! namzdf_tmx tidal mixing parameterization ("key_zdftmx")
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namzdf ! vertical physics
+!-----------------------------------------------------------------------
+ rn_avm0 = 1.2e-4 ! vertical eddy viscosity [m2/s] (background Kz if not "key_zdfcst")
+ rn_avt0 = 1.2e-5 ! vertical eddy diffusivity [m2/s] (background Kz if not "key_zdfcst")
+ nn_avb = 0 ! profile for background avt & avm (=1) or not (=0)
+ nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0)
+ ln_zdfevd = .true. ! enhanced vertical diffusion (evd) (T) or not (F)
+ nn_evdm = 0 ! evd apply on tracer (=0) or on tracer and momentum (=1)
+ rn_avevd = 100. ! evd mixing coefficient [m2/s]
+ ln_zdfnpc = .false. ! Non-Penetrative Convective algorithm (T) or not (F)
+ nn_npc = 1 ! frequency of application of npc
+ nn_npcp = 365 ! npc control print frequency
+ ln_zdfexp = .false. ! time-stepping: split-explicit (T) or implicit (F) time stepping
+ nn_zdfexp = 3 ! number of sub-timestep for ln_zdfexp=T
+ ln_zdfqiao = .false. ! Enhanced wave vertical mixing Qiao (2010) (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave)
+/
+!-----------------------------------------------------------------------
+&namzdf_ric ! richardson number dependent vertical diffusion ("key_zdfric" )
+!-----------------------------------------------------------------------
+ rn_avmri = 100.e-4 ! maximum value of the vertical viscosity
+ rn_alp = 5. ! coefficient of the parameterization
+ nn_ric = 2 ! coefficient of the parameterization
+ rn_ekmfc = 0.7 ! Factor in the Ekman depth Equation
+ rn_mldmin = 1.0 ! minimum allowable mixed-layer depth estimate (m)
+ rn_mldmax = 1000.0 ! maximum allowable mixed-layer depth estimate (m)
+ rn_wtmix = 10.0 ! vertical eddy viscosity coeff [m2/s] in the mixed-layer
+ rn_wvmix = 10.0 ! vertical eddy diffusion coeff [m2/s] in the mixed-layer
+ ln_mldw = .true. ! Flag to use or not the mixed layer depth param.
+/
+!-----------------------------------------------------------------------
+&namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion ("key_zdftke")
+!-----------------------------------------------------------------------
+ rn_ediff = 0.1 ! coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) )
+ rn_ediss = 0.7 ! coef. of the Kolmogoroff dissipation
+ rn_ebb = 67.83 ! coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T)
+ rn_emin = 1.e-6 ! minimum value of tke [m2/s2]
+ rn_emin0 = 1.e-4 ! surface minimum value of tke [m2/s2]
+ rn_bshear = 1.e-20 ! background shear (>0) currently a numerical threshold (do not change it)
+ nn_mxl = 2 ! mixing length: = 0 bounded by the distance to surface and bottom
+ ! = 1 bounded by the local vertical scale factor
+ ! = 2 first vertical derivative of mixing length bounded by 1
+ ! = 3 as =2 with distinct disspipative an mixing length scale
+ nn_pdl = 1 ! Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm)
+ ln_mxl0 = .true. ! surface mixing length scale = F(wind stress) (T) or not (F)
+ rn_mxl0 = 0.04 ! surface buoyancy lenght scale minimum value
+ ln_lc = .true. ! Langmuir cell parameterisation (Axell 2002)
+ rn_lc = 0.15 ! coef. associated to Langmuir cells
+ nn_etau = 1 ! penetration of tke below the mixed layer (ML) due to near intertial waves
+ ! = 0 no penetration
+ ! = 1 add a tke source below the ML
+ ! = 2 add a tke source just at the base of the ML
+ ! = 3 as = 1 applied on HF part of the stress (ln_cpl=T)
+ rn_efr = 0.05 ! fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2)
+ nn_htau = 1 ! type of exponential decrease of tke penetration below the ML
+ ! = 0 constant 10 m length scale
+ ! = 1 0.5m at the equator to 30m poleward of 40 degrees
+/
+!-----------------------------------------------------------------------
+&namzdf_gls ! GLS vertical diffusion ("key_zdfgls")
+!-----------------------------------------------------------------------
+ rn_emin = 1.e-7 ! minimum value of e [m2/s2]
+ rn_epsmin = 1.e-12 ! minimum value of eps [m2/s3]
+ ln_length_lim = .true. ! limit on the dissipation rate under stable stratification (Galperin et al., 1988)
+ rn_clim_galp = 0.267 ! galperin limit
+ ln_sigpsi = .true. ! Activate or not Burchard 2001 mods on psi schmidt number in the wb case
+ rn_crban = 100. ! Craig and Banner 1994 constant for wb tke flux
+ rn_charn = 70000. ! Charnock constant for wb induced roughness length
+ rn_hsro = 0.02 ! Minimum surface roughness
+ rn_frac_hs = 1.3 ! Fraction of wave height as roughness (if nn_z0_met=2)
+ nn_z0_met = 2 ! Method for surface roughness computation (0/1/2/3)
+ ! ! =3 requires ln_wave=T
+ nn_bc_surf = 1 ! surface condition (0/1=Dir/Neum)
+ nn_bc_bot = 1 ! bottom condition (0/1=Dir/Neum)
+ nn_stab_func = 2 ! stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB)
+ nn_clos = 1 ! predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen)
+/
+!-----------------------------------------------------------------------
+&namzdf_ddm ! double diffusive mixing parameterization ("key_zdfddm")
+!-----------------------------------------------------------------------
+ rn_avts = 1.e-4 ! maximum avs (vertical mixing on salinity)
+ rn_hsbfr = 1.6 ! heat/salt buoyancy flux ratio
+/
+!-----------------------------------------------------------------------
+&namzdf_tmx ! tidal mixing parameterization ("key_zdftmx")
+!-----------------------------------------------------------------------
+ rn_htmx = 500. ! vertical decay scale for turbulence (meters)
+ rn_n2min = 1.e-8 ! threshold of the Brunt-Vaisala frequency (s-1)
+ rn_tfe = 0.333 ! tidal dissipation efficiency
+ rn_me = 0.2 ! mixing efficiency
+ ln_tmx_itf = .true. ! ITF specific parameterisation
+ rn_tfe_itf = 1. ! ITF tidal dissipation efficiency
+/
+!-----------------------------------------------------------------------
+&namzdf_tmx_new ! internal wave-driven mixing parameterization ("key_zdftmx_new" & "key_zdfddm")
+!-----------------------------------------------------------------------
+ nn_zpyc = 1 ! pycnocline-intensified dissipation scales as N (=1) or N^2 (=2)
+ ln_mevar = .true. ! variable (T) or constant (F) mixing efficiency
+ ln_tsdiff = .true. ! account for differential T/S mixing (T) or not (F)
+/
+!!======================================================================
+!! *** Miscellaneous namelists ***
+!!======================================================================
+!! nammpp Massively Parallel Processing ("key_mpp_mpi)
+!! namctl Control prints
+!! namsto Stochastic parametrization of EOS
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&nammpp ! Massively Parallel Processing ("key_mpp_mpi)
+!-----------------------------------------------------------------------
+ cn_mpi_send = 'I' ! mpi send/recieve type ='S', 'B', or 'I' for standard send,
+ ! buffer blocking send or immediate non-blocking sends, resp.
+ nn_buffer = 0 ! size in bytes of exported buffer ('B' case), 0 no exportation
+ ln_nnogather= .false. ! activate code to avoid mpi_allgather use at the northfold
+ jpni = 0 ! jpni number of processors following i (set automatically if < 1)
+ jpnj = 0 ! jpnj number of processors following j (set automatically if < 1)
+ jpnij = 0 ! jpnij number of local domains (set automatically if < 1)
+/
+!-----------------------------------------------------------------------
+&namctl ! Control prints
+!-----------------------------------------------------------------------
+ ln_ctl = .false. ! trends control print (expensive!)
+ nn_print = 0 ! level of print (0 no extra print)
+ nn_ictls = 0 ! start i indice of control sum (use to compare mono versus
+ nn_ictle = 0 ! end i indice of control sum multi processor runs
+ nn_jctls = 0 ! start j indice of control over a subdomain)
+ nn_jctle = 0 ! end j indice of control
+ nn_isplt = 1 ! number of processors in i-direction
+ nn_jsplt = 1 ! number of processors in j-direction
+ nn_timing = 0 ! timing by routine activated (=1) creates timing.output file, or not (=0)
+ nn_diacfl = 0 ! Write out CFL diagnostics (=1) in cfl_diagnostics.ascii, or not (=0)
+/
+!-----------------------------------------------------------------------
+&namsto ! Stochastic parametrization of EOS (default: NO)
+!-----------------------------------------------------------------------
+ ln_sto_eos = .false. ! stochastic equation of state
+ nn_sto_eos = 1 ! number of independent random walks
+ rn_eos_stdxy= 1.4 ! random walk horz. standard deviation (in grid points)
+ rn_eos_stdz = 0.7 ! random walk vert. standard deviation (in grid points)
+ rn_eos_tcor = 1440. ! random walk time correlation (in timesteps)
+ nn_eos_ord = 1 ! order of autoregressive processes
+ nn_eos_flt = 0 ! passes of Laplacian filter
+ rn_eos_lim = 2.0 ! limitation factor (default = 3.0)
+ ln_rststo = .false. ! start from mean parameter (F) or from restart file (T)
+ ln_rstseed = .true. ! read seed of RNG from restart file
+ cn_storst_in = "restart_sto" ! suffix of stochastic parameter restart file (input)
+ cn_storst_out = "restart_sto" ! suffix of stochastic parameter restart file (output)
+/
+
+!!======================================================================
+!! *** Diagnostics namelists ***
+!!======================================================================
+!! namtrd dynamics and/or tracer trends (default F)
+!! namptr Poleward Transport Diagnostics (default F)
+!! namhsb Heat and salt budgets (default F)
+!! namdiu Cool skin and warm layer models (default F)
+!! namdiu Cool skin and warm layer models (default F)
+!! namflo float parameters ("key_float")
+!! nam_diaharm Harmonic analysis of tidal constituents ("key_diaharm")
+!! namdct transports through some sections ("key_diadct")
+!! nam_diatmb Top Middle Bottom Output (default F)
+!! nam_dia25h 25h Mean Output (default F)
+!! namnc4 netcdf4 chunking and compression settings ("key_netcdf4")
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namtrd ! trend diagnostics (default F)
+!-----------------------------------------------------------------------
+ ln_glo_trd = .false. ! (T) global domain averaged diag for T, T^2, KE, and PE
+ ln_dyn_trd = .false. ! (T) 3D momentum trend output
+ ln_dyn_mxl = .false. ! (T) 2D momentum trends averaged over the mixed layer (not coded yet)
+ ln_vor_trd = .false. ! (T) 2D barotropic vorticity trends (not coded yet)
+ ln_KE_trd = .false. ! (T) 3D Kinetic Energy trends
+ ln_PE_trd = .false. ! (T) 3D Potential Energy trends
+ ln_tra_trd = .false. ! (T) 3D tracer trend output
+ ln_tra_mxl = .false. ! (T) 2D tracer trends averaged over the mixed layer (not coded yet)
+ nn_trd = 365 ! print frequency (ln_glo_trd=T) (unit=time step)
+/
+!!gm nn_ctls = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk)
+!!gm rn_ucf = 1. ! unit conversion factor (=1 -> /seconds ; =86400. -> /day)
+!!gm cn_trdrst_in = "restart_mld" ! suffix of ocean restart name (input)
+!!gm cn_trdrst_out = "restart_mld" ! suffix of ocean restart name (output)
+!!gm ln_trdmld_restart = .false. ! restart for ML diagnostics
+!!gm ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S
+!!gm
+!-----------------------------------------------------------------------
+&namptr ! Poleward Transport Diagnostic (default F)
+!-----------------------------------------------------------------------
+ ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F)
+ ln_subbas = .false. ! Atlantic/Pacific/Indian basins computation (T) or not
+/
+!-----------------------------------------------------------------------
+&namhsb ! Heat and salt budgets (default F)
+!-----------------------------------------------------------------------
+ ln_diahsb = .false. ! check the heat and salt budgets (T) or not (F)
+/
+!-----------------------------------------------------------------------
+&namdiu ! Cool skin and warm layer models (default F)
+!-----------------------------------------------------------------------
+ ln_diurnal = .false. !
+ ln_diurnal_only = .false. !
+/
+!-----------------------------------------------------------------------
+&namflo ! float parameters ("key_float")
+!-----------------------------------------------------------------------
+ jpnfl = 1 ! total number of floats during the run
+ jpnnewflo = 0 ! number of floats for the restart
+ ln_rstflo = .false. ! float restart (T) or not (F)
+ nn_writefl = 75 ! frequency of writing in float output file
+ nn_stockfl = 5475 ! frequency of creation of the float restart file
+ ln_argo = .false. ! Argo type floats (stay at the surface each 10 days)
+ ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T)
+ ! ! or computed with Blanke' scheme (F)
+ ln_ariane = .true. ! Input with Ariane tool convention(T)
+ ln_flo_ascii= .true. ! Output with Ariane tool netcdf convention(F) or ascii file (T)
+/
+!-----------------------------------------------------------------------
+&nam_diaharm ! Harmonic analysis of tidal constituents ("key_diaharm")
+!-----------------------------------------------------------------------
+ nit000_han = 1 ! First time step used for harmonic analysis
+ nitend_han = 75 ! Last time step used for harmonic analysis
+ nstep_han = 15 ! Time step frequency for harmonic analysis
+ tname(1) = 'M2' ! Name of tidal constituents
+ tname(2) = 'K1'
+/
+!-----------------------------------------------------------------------
+&namdct ! transports through some sections ("key_diadct")
+!-----------------------------------------------------------------------
+ nn_dct = 15 ! time step frequency for transports computing
+ nn_dctwri = 15 ! time step frequency for transports writing
+ nn_secdebug= 112 ! 0 : no section to debug
+ ! ! -1 : debug all section
+ ! ! 0 < n : debug section number n
+/
+!-----------------------------------------------------------------------
+&nam_diatmb ! Top Middle Bottom Output (default F)
+!-----------------------------------------------------------------------
+ ln_diatmb = .false. ! Choose Top Middle and Bottom output or not
+/
+!-----------------------------------------------------------------------
+&nam_dia25h ! 25h Mean Output (default F)
+!-----------------------------------------------------------------------
+ ln_dia25h = .false. ! Choose 25h mean output or not
+/
+!-----------------------------------------------------------------------
+&namnc4 ! netcdf4 chunking and compression settings ("key_netcdf4")
+!-----------------------------------------------------------------------
+ nn_nchunks_i= 4 ! number of chunks in i-dimension
+ nn_nchunks_j= 4 ! number of chunks in j-dimension
+ nn_nchunks_k= 31 ! number of chunks in k-dimension
+ ! ! setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which
+ ! ! is optimal for postprocessing which works exclusively with horizontal slabs
+ ln_nc4zip = .true. ! (T) use netcdf4 chunking and compression
+ ! ! (F) ignore chunking information and produce netcdf3-compatible files
+/
+
+!!======================================================================
+!! *** Observation & Assimilation ***
+!!======================================================================
+!! namobs observation and model comparison
+!! nam_asminc assimilation increments ('key_asminc')
+!!======================================================================
+!
+!-----------------------------------------------------------------------
+&namobs ! observation usage switch
+!-----------------------------------------------------------------------
+ ln_diaobs = .false. ! Logical switch for the observation operator
+ ln_t3d = .false. ! Logical switch for T profile observations
+ ln_s3d = .false. ! Logical switch for S profile observations
+ ln_sla = .false. ! Logical switch for SLA observations
+ ln_sst = .false. ! Logical switch for SST observations
+ ln_sic = .false. ! Logical switch for Sea Ice observations
+ ln_vel3d = .false. ! Logical switch for velocity observations
+ ln_altbias = .false. ! Logical switch for altimeter bias correction
+ ln_nea = .false. ! Logical switch for rejection of observations near land
+ ln_grid_global = .true. ! Logical switch for global distribution of observations
+ ln_grid_search_lookup = .false. ! Logical switch for obs grid search w/lookup table
+ ln_ignmis = .true. ! Logical switch for ignoring missing files
+ ln_s_at_t = .false. ! Logical switch for computing model S at T obs if not there
+ ln_sstnight = .false. ! Logical switch for calculating night-time average for SST obs
+! All of the *files* variables below are arrays. Use namelist_cfg to add more files
+ cn_profbfiles = 'profiles_01.nc' ! Profile feedback input observation file names
+ cn_slafbfiles = 'sla_01.nc' ! SLA feedback input observation file names
+ cn_sstfbfiles = 'sst_01.nc' ! SST feedback input observation file names
+ cn_sicfbfiles = 'sic_01.nc' ! SIC feedback input observation file names
+ cn_velfbfiles = 'vel_01.nc' ! Velocity feedback input observation file names
+ cn_altbiasfile = 'altbias.nc' ! Altimeter bias input file name
+ cn_gridsearchfile='gridsearch.nc' ! Grid search file name
+ rn_gridsearchres = 0.5 ! Grid search resolution
+ rn_dobsini = 00010101.000000 ! Initial date in window YYYYMMDD.HHMMSS
+ rn_dobsend = 00010102.000000 ! Final date in window YYYYMMDD.HHMMSS
+ nn_1dint = 0 ! Type of vertical interpolation method
+ nn_2dint = 0 ! Type of horizontal interpolation method
+ nn_msshc = 0 ! MSSH correction scheme
+ rn_mdtcorr = 1.61 ! MDT correction
+ rn_mdtcutoff = 65.0 ! MDT cutoff for computed correction
+ nn_profdavtypes = -1 ! Profile daily average types - array
+ ln_sstbias = .false. !
+ cn_sstbias_files = 'sstbias.nc' !
+/
+!-----------------------------------------------------------------------
+&nam_asminc ! assimilation increments ('key_asminc')
+!-----------------------------------------------------------------------
+ ln_bkgwri = .false. ! Logical switch for writing out background state
+ ln_trainc = .false. ! Logical switch for applying tracer increments
+ ln_dyninc = .false. ! Logical switch for applying velocity increments
+ ln_sshinc = .false. ! Logical switch for applying SSH increments
+ ln_asmdin = .false. ! Logical switch for Direct Initialization (DI)
+ ln_asmiau = .false. ! Logical switch for Incremental Analysis Updating (IAU)
+ nitbkg = 0 ! Timestep of background in [0,nitend-nit000-1]
+ nitdin = 0 ! Timestep of background for DI in [0,nitend-nit000-1]
+ nitiaustr = 1 ! Timestep of start of IAU interval in [0,nitend-nit000-1]
+ nitiaufin = 15 ! Timestep of end of IAU interval in [0,nitend-nit000-1]
+ niaufn = 0 ! Type of IAU weighting function
+ ln_salfix = .false. ! Logical switch for ensuring that the sa > salfixmin
+ salfixmin = -9999 ! Minimum salinity after applying the increments
+ nn_divdmp = 0 ! Number of iterations of divergence damping operator
+/
diff --git a/EXP_Apr19/run_script.pbs b/EXP_Apr19/run_script.pbs
new file mode 100755
index 0000000000000000000000000000000000000000..8d6cc0c20a75afd53cd9950c1085e65b2460efc5
--- /dev/null
+++ b/EXP_Apr19/run_script.pbs
@@ -0,0 +1,41 @@
+#!/bin/bash
+# ---------------------------
+#===============================================================
+# CLUSTER BITS
+#===============================================================
+#PBS -N BoBEAS
+#PBS -l select=21
+#PBS -l walltime=00:20:00
+#PBS -A n01-ACCORD
+#PBS -j oe
+#PBS -r n
+# mail alert at (b)eginning, (e)nd and (a)bortion of execution
+#PBS -m bea
+#PBS -M jelt@noc.ac.uk
+
+module swap PrgEnv-cray PrgEnv-intel
+module load cray-netcdf-hdf5parallel
+module load cray-hdf5-parallel
+
+export PBS_O_WORKDIR=$(readlink -f $PBS_O_WORKDIR)
+# Change to the direcotry that the job was submitted from
+cd $PBS_O_WORKDIR
+
+
+# Set the number of threads to 1
+# This prevents any system libraries from automatically
+# using threading.
+export OMP_NUM_THREADS=1
+# Change to the directory that the job was submitted from
+ulimit -s unlimited
+ulimit -c unlimited
+
+export NEMOproc=480 #550
+export XIOSproc=1
+
+#===============================================================
+# LAUNCH JOB
+#===============================================================
+echo `date` : Launch Job
+aprun -b -n 5 -N 5 ./xios_server.exe : -n $NEMOproc -N 24 ./opa
+exit