Changeset 4863 for trunk/NEMOGCM/NEMO/LIM_SRC_3/limdyn.F90

Timestamp:

2014-11-18T15:45:57+01:00 (9 years ago)

Author:

clem

Message:

LIM3: change elastic timescale to ratio of elastic timescale over ice time step

File:

• trunk/NEMOGCM/NEMO/LIM_SRC_3/limdyn.F90 (modified) (3 diffs)

trunk/NEMOGCM/NEMO/LIM_SRC_3/limdyn.F90

@@ -243 +243 @@
       NAMELIST/namicedyn/ epsd, om, cw, angvg, pstar,   &
          &                c_rhg, creepl, ecc, ahi0,     &
-         &                nevp, telast, alphaevp, hminrhg
+         &                nevp, relast, alphaevp, hminrhg
       !!-------------------------------------------------------------------
 
@@ -269 +269 @@
          WRITE(numout,*) '   horizontal diffusivity coeff. for sea-ice        ahi0   = ', ahi0
          WRITE(numout,*) '   number of iterations for subcycling              nevp   = ', nevp
-         WRITE(numout,*) '   timescale for elastic waves                      telast = ', telast
+         WRITE(numout,*) '   ratio of elastic timescale over ice time step    relast = ', relast
          WRITE(numout,*) '   coefficient for the solution of int. stresses  alphaevp = ', alphaevp
          WRITE(numout,*) '   min ice thickness for rheology calculations     hminrhg = ', hminrhg
@@ -287 +287 @@
       pstarh = pstar * 0.5_wp
 
+      ! elastic damping
+      telast = relast * rdt_ice
+
       !  Diffusion coefficients.
       ahiu(:,:) = ahi0 * umask(:,:,1)