Changeset 2580

Timestamp:

2011-02-08T09:01:02+01:00 (13 years ago)

Author:

cbricaud

Message:

Corrections for EVP rheology in LIM2: add missing lines in limrhg.F90 and change parameters for ORCA2 in namelist_ice_lim2

Location:

trunk/NEMOGCM

Files:

• CONFIG/ORCA2_LIM/EXP00/namelist_ice_lim2 (modified) (1 diff)
• NEMO/LIM_SRC_3/limrhg.F90 (modified) (1 diff)

trunk/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist_ice_lim2

@@ -48 +48 @@
    c_rhg       =  20.0     !  2nd bulk-rhelogy parameter
    etamn       =   0.0e+07 !  minimun value for viscosity
-   creepl      =   2.0e-08 !  creep limit
+   creepl      =   1.0e-08 !  creep limit
    ecc         =   2.0     !  eccentricity of the elliptical yield curve
    ahi0        = 350.e0    !  horizontal eddy diffusivity coefficient for sea-ice [m2/s]
-   nevp        =   360     !  number of EVP subcycling iterations
-   telast      =   3600    !  timescale for EVP elastic waves
+   nevp        =   120     !  number of EVP subcycling iterations
+   telast      =   9600    !  timescale for EVP elastic waves
    alphaevp    =   1.0     !  coefficient for the solution of EVP int. stresses
 /

trunk/NEMOGCM/NEMO/LIM_SRC_3/limrhg.F90

@@ -190 +190 @@
             zpresh(ji,jj) = tms(ji,jj) *  strength(ji,jj) * 0.5_wp
 #endif
+#if defined key_lim2
+            zpresh(ji,jj) = tms(ji,jj) *  pstar * vt_i(ji,jj) * EXP( -c_rhg * (1. - at_i(ji,jj) ) )
+#endif
             ! tmi = 1 where there is ice or on land
             tmi(ji,jj)    = 1._wp - ( 1._wp - MAX( 0._wp , SIGN ( 1._wp , vt_i(ji,jj) - epsd ) ) ) * tms(ji,jj)