Changeset 13536 for NEMO/branches/2020/SI3-03_VP_rheology/cfgs

Timestamp:

2020-09-29T11:22:04+02:00 (4 years ago)

Author:

vancop

Message:

VP rheology, first code version finalised

Location:

NEMO/branches/2020/SI3-03_VP_rheology/cfgs/SHARED

Files:

• field_def_nemo-ice.xml (modified) (2 diffs)
• namelist_ice_ref (modified) (2 diffs)

NEMO/branches/2020/SI3-03_VP_rheology/cfgs/SHARED/field_def_nemo-ice.xml

@@ -75 +75 @@
           <field id="utau_bi"      long_name="X-component of ocean bottom stress on sea ice -landfast" standard_name="ocean_bottom_upward_x_stress"              unit="N/m2" />
           <field id="vtau_bi"      long_name="Y-component of ocean bottom stress on sea ice -landfast" standard_name="ocean_bottom_upward_y_stress"              unit="N/m2" />
-          <field id="isig1"        long_name="1st principal stress component for EVP rhg"                                                                        unit=""     />
-          <field id="isig2"        long_name="2nd principal stress component for EVP rhg"                                                                        unit=""     />
-          <field id="isig3"        long_name="convergence measure for EVP rheology (must be around 1)"                                                           unit=""     />
+          <field id="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"  />
@@ -398 +397 @@
      <field field_ref="normstr"          name="normstr" />
      <field field_ref="sheastr"          name="sheastr" />
-     <field field_ref="isig1"            name="isig1"   />
-     <field field_ref="isig2"            name="isig2"   />
-     <field field_ref="isig3"            name="isig3"   />
+     <field field_ref="sig1_pnorm"       name="sig1_pnorm"   />
+     <field field_ref="sig2_pnorm"       name="sig2_pnorm"   />
      
      <!-- heat fluxes -->

NEMO/branches/2020/SI3-03_VP_rheology/cfgs/SHARED/namelist_ice_ref

@@ -91 +91 @@
 &namdyn_rhg     !   Ice rheology
 !------------------------------------------------------------------------------
-   ln_rhg_EVP       = .true.          !  EVP rheology
+   ln_rhg_EVP       = .false.         !  EVP rheology
       ln_aEVP       = .true.          !     adaptive rheology (Kimmritz et al. 2016 & 2017)
       rn_creepl     =   2.0e-9        !     creep limit [1/s]
@@ -99 +99 @@
                                       !        advised value: 1/3 (rn_nevp=120) or 1/9 (rn_nevp=300)
    ln_rhg_chkcvg    = .false.         !  check convergence of rheology (outputs: file ice_cvg.nc & variable uice_cvg)
+   ln_rhg_VP        = .true.          !  VP rheology
+   nn_nout_vp       = 10              !     number of outer iterations
+   nn_ninn_vp       = 1500            !     number of inner iterations
+   ln_zebra_vp      = .false.         !     activate zebra solver
+   rn_relaxu_vp     = 0.95            !     relaxation factor for u-velocity
+   rn_relaxv_vp     = 0.95            !     relaxation factor for v-velocity
+   rn_uerr_max_vp   = 0.80            !     maximum error on velocity
+   rn_uerr_min_vp   = 1.0e-04         !     velocity to decide convergence
+   nn_cvgchk_vp     = 5               !     iteration step for convergence check
 /
 !------------------------------------------------------------------------------