Changeset 10312 for NEMO/branches/2018/dev_r9947_SI3_advection/cfgs

Timestamp:

2018-11-15T10:46:57+01:00 (5 years ago)

Author:

clem

Message:

add the parameterization of landfast ice from Lemieux2015 and 2016 with the addition of isotropic tensile strength

Location:

NEMO/branches/2018/dev_r9947_SI3_advection/cfgs/SHARED

Files:

• field_def_nemo-ice.xml (modified) (1 diff)
• namelist_ice_ref (modified) (1 diff)

NEMO/branches/2018/dev_r9947_SI3_advection/cfgs/SHARED/field_def_nemo-ice.xml

@@ -29 +29 @@
           <field id="icemask15"    long_name="Ice mask (0 if ice conc. lower than 15%, 1 otherwise)"   standard_name="sea_ice_mask15"                            unit="" />
      <field id="icepres"      long_name="Fraction of time steps with sea ice"                     standard_name="sea_ice_time_fraction"                     unit="" />
-     
+          <field id="fasticepres"  long_name="Fraction of time steps with landfast ice"                standard_name="fast_ice_time_fraction"                    unit="" />
+ 
      <!-- general fields -->
           <field id="icemass"      long_name="Sea-ice mass per area"                                   standard_name="sea_ice_amount"                            unit="kg/m2"/>

NEMO/branches/2018/dev_r9947_SI3_advection/cfgs/SHARED/namelist_ice_ref

@@ -52 +52 @@
    ln_dynRHGADV     = .false.         !  dyn.: no ridge/raft & no corrections  (rheology + advection)
    ln_dynADV        = .false.         !  dyn.: only advection w prescribed vel.(rn_uvice + advection)
-      rn_uice       =   0.00001       !        prescribed ice u-velocity
-      rn_vice       =   0.            !        prescribed ice v-velocity
-   rn_ishlat        =   2.            !  free slip (0) ; partial slip (0-2) ; no slip (2) ; strong slip (>2)
-   ln_landfast      = .false.         !  landfast ice parameterization (T or F)                           
-      rn_gamma      =   0.15          !     fraction of ocean depth that ice must reach to initiate landfast
-                                      !        recommended range: [0.1 ; 0.25]
-      rn_icebfr     =  10.            !     maximum bottom stress per unit area of contact [N/m2]                 
-                                      !        a very large value ensures ice velocity=0 even with a small contact area
-                                      !        recommended range: ?? (should be greater than atm-ice stress => >0.1 N/m2)
-      rn_lfrelax    =   1.e-5         !     relaxation time scale to reach static friction [s-1]
+      rn_uice       =   0.5           !        prescribed ice u-velocity
+      rn_vice       =   0.5           !        prescribed ice v-velocity
+   rn_ishlat        =   2.            !  lbc : free slip (0) ; partial slip (0-2) ; no slip (2) ; strong slip (>2)
+   ln_landfast_L16  = .false.         !  landfast: parameterization from Lemieux 2016
+   ln_landfast_home = .false.         !  landfast: parameterization from "home made"
+      rn_depfra     =   0.125         !        fraction of ocean depth that ice must reach to initiate landfast
+                                      !          recommended range: [0.1 ; 0.25] - L16=0.125 - home=0.15
+      rn_icebfr     =  15.            !        ln_landfast_L16:  maximum bottom stress per unit volume [N/m3]
+                                      !        ln_landfast_home: maximum bottom stress per unit area of contact [N/m2]
+                                      !          recommended range: ?? L16=15 - home=10
+      rn_lfrelax    =   1.e-5         !        relaxation time scale to reach static friction [s-1]
+      rn_tensile    =   0.2           !        isotropic tensile strength
 /
 !------------------------------------------------------------------------------