Changeset 4055

Timestamp:

2013-10-09T16:39:46+02:00 (11 years ago)

Author:

clem

Message:

correction for LIM2, see ticket:#1154

Location:

branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_2

Files:

• ice_2.F90 (modified) (1 diff)
• limdyn_2.F90 (modified) (2 diffs)

branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_2/ice_2.F90

@@ -51 +51 @@
    REAL(wp), PUBLIC ::   ahi0   = 350._wp      !: sea-ice hor. eddy diffusivity coeff. (m2/s)
    REAL(wp), PUBLIC ::   alphaevp = 1._wp      !: coefficient for the solution of EVP int. stresses
+   REAL(wp), PUBLIC ::   hminrhg = 0.05_wp     !: clem : ice thickness (in m) below which ice velocity is set to ocean velocity
 
    REAL(wp), PUBLIC ::   usecc2                !:  = 1.0 / ( ecc * ecc )

branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_2/limdyn_2.F90

@@ -227 +227 @@
          &                dm, nbiter, nbitdr, om, resl, cw, angvg, pstar,   &
          &                c_rhg, etamn, creepl, ecc, ahi0,                  &
-         &                nevp, telast,alphaevp
+         &                nevp, telast, alphaevp, hminrhg
       !!-------------------------------------------------------------------
 
@@ -255 +255 @@
          WRITE(numout,*) '       timescale for elastic waves telast = ', telast
          WRITE(numout,*) '       coefficient for the solution of int. stresses alphaevp = ', alphaevp
+         WRITE(numout,*) '       min ice thickness for rheology calculations     hminrhg = ', hminrhg
       ENDIF
       !