Changeset 13601 for NEMO/trunk/cfgs

Timestamp:

2020-10-14T17:59:34+02:00 (4 years ago)

Author:

clem

Message:

trunk: rewrite heat budget of sea ice to make it perfectly conservative by construction. Also, activating ln_icediachk now gives an ascii file (icedrift_diagnostics.ascii) containing mass, salt and heat global conservation issues (if any). In addition, 2D drift files can be outputed (icedrift_mass…) and field_def is changed accordingly. Note that advection schemes are not yet commited since there is still a restartability issue that I do not understand.

File:

• NEMO/trunk/cfgs/SHARED/field_def_nemo-ice.xml (modified) (3 diffs)

NEMO/trunk/cfgs/SHARED/field_def_nemo-ice.xml

@@ -78 +78 @@
           <field id="isig2"        long_name="2nd principal stress component for EVP rhg"                                                                        unit=""     />
           <field id="isig3"        long_name="convergence measure for EVP rheology (must be around 1)"                                                           unit=""     />
+          <!-- clem: sig1 sig2 sig 3 will be replaced by sig1 and sig2 in a following commit -->         
+          <field id="sig1_pnorm"   long_name="P-normalized 1st principal stress component"                                                                       unit=""     />
+          <field id="sig2_pnorm"   long_name="P-normalized 2nd principal stress component"                                                                       unit=""     />
           <field id="normstr"      long_name="Average normal stress in sea ice"                        standard_name="average_normal_stress"                     unit="N/m"  />
           <field id="sheastr"      long_name="Maximum shear stress in sea ice"                         standard_name="maximum_shear_stress"                      unit="N/m"  />
@@ -165 +168 @@
 
           <!-- diags -->
-          <field id="hfxdhc"       long_name="Heat content variation in snow and ice (neg = ice cooling)"   unit="W/m2" />
+          <field id="hfxdhc"        long_name="Heat content variation in snow and ice (neg = ice cooling)"   unit="W/m2" />
+          <!-- available if ln_icediachk=T -->
+          <field id="icedrift_mass" long_name="Ice mass drift (conservation check)"   unit="kg/m2/s" />
+          <field id="icedrift_salt" long_name="Ice salt drift (conservation check)"   unit="kg/m2/s" />
+          <field id="icedrift_heat" long_name="Ice heat drift (conservation check)"   unit="W/m2"    />
 
      <!-- sbcssm variables -->
@@ -459 +466 @@
           <field field_ref="vfxthin"          name="vfxthin" />
    
-     <!-- diag error for negative ice volume after advection -->
-     <field field_ref="iceneg_pres"      name="sineg_pres" />
-     <field field_ref="iceneg_volu"      name="sineg_volu" />
-     <field field_ref="iceneg_hfx"       name="sineg_hfx"  />
         </field_group>