Changeset 5123 for trunk/NEMOGCM/NEMO/LIM_SRC_2
- Timestamp:
- 2015-03-04T17:06:03+01:00 (9 years ago)
- Location:
- trunk/NEMOGCM/NEMO/LIM_SRC_2
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/LIM_SRC_2/ice_2.F90
r5056 r5123 35 35 INTEGER , PUBLIC :: nbiter !: number of sub-time steps for relaxation 36 36 INTEGER , PUBLIC :: nbitdr !: maximum number of iterations for relaxation 37 INTEGER , PUBLIC :: n evp!: number of EVP subcycling iterations37 INTEGER , PUBLIC :: nn_nevp !: number of EVP subcycling iterations 38 38 INTEGER , PUBLIC :: telast !: timescale for EVP elastic waves 39 39 REAL(wp), PUBLIC :: epsd !: tolerance parameter for dynamic … … 47 47 REAL(wp), PUBLIC :: c_rhg !: second bulk-rhelogy parameter 48 48 REAL(wp), PUBLIC :: etamn !: minimun value for viscosity 49 REAL(wp), PUBLIC :: creepl!: creep limit50 REAL(wp), PUBLIC :: ecc!: eccentricity of the elliptical yield curve49 REAL(wp), PUBLIC :: rn_creepl !: creep limit 50 REAL(wp), PUBLIC :: rn_ecc !: eccentricity of the elliptical yield curve 51 51 REAL(wp), PUBLIC :: ahi0 !: sea-ice hor. eddy diffusivity coeff. (m2/s) 52 52 REAL(wp), PUBLIC :: alphaevp !: coefficient for the solution of EVP int. stresses 53 REAL(wp), PUBLIC :: hminrhg = 0.001_wp !: clem : ice volume (a*h in m) below which ice velocity is set to ocean velocity 54 55 REAL(wp), PUBLIC :: usecc2 !: = 1.0 / ( ecc * ecc ) 53 54 REAL(wp), PUBLIC :: usecc2 !: = 1.0 / ( rn_ecc * rn_ecc ) 56 55 REAL(wp), PUBLIC :: rhoco !: = rau0 * cw 57 56 REAL(wp), PUBLIC :: sangvg, cangvg !: sin and cos of the turning angle for ocean stress -
trunk/NEMOGCM/NEMO/LIM_SRC_2/limdyn_2.F90
r4624 r5123 227 227 NAMELIST/namicedyn/ epsd, alpha, & 228 228 & dm, nbiter, nbitdr, om, resl, cw, angvg, pstar, & 229 & c_rhg, etamn, creepl,ecc, ahi0, &230 & n evp, telast, alphaevp, hminrhg229 & c_rhg, etamn, rn_creepl, rn_ecc, ahi0, & 230 & nn_nevp, telast, alphaevp 231 231 !!------------------------------------------------------------------- 232 232 … … 256 256 WRITE(numout,*) ' second bulk-rhelogy parameter c_rhg = ', c_rhg 257 257 WRITE(numout,*) ' minimun value for viscosity etamn = ', etamn 258 WRITE(numout,*) ' creep limit creepl = ',creepl259 WRITE(numout,*) ' eccentricity of the elliptical yield curve ecc = ',ecc258 WRITE(numout,*) ' creep limit rn_creepl = ', rn_creepl 259 WRITE(numout,*) ' eccentricity of the elliptical yield curve rn_ecc = ', rn_ecc 260 260 WRITE(numout,*) ' horizontal diffusivity coeff. for sea-ice ahi0 = ', ahi0 261 WRITE(numout,*) ' number of iterations for subcycling nevp = ',nevp261 WRITE(numout,*) ' number of iterations for subcycling nn_nevp= ', nn_nevp 262 262 WRITE(numout,*) ' timescale for elastic waves telast = ', telast 263 263 WRITE(numout,*) ' coefficient for the solution of int. stresses alphaevp = ', alphaevp 264 WRITE(numout,*) ' min ice thickness for rheology calculations hminrhg = ', hminrhg265 264 ENDIF 266 265 ! … … 272 271 273 272 ! Initialization 274 usecc2 = 1.0 / ( ecc *ecc )273 usecc2 = 1.0 / ( rn_ecc * rn_ecc ) 275 274 rhoco = rau0 * cw 276 275 angvg = angvg * rad ! convert angvg from degree to radian -
trunk/NEMOGCM/NEMO/LIM_SRC_2/limrhg_2.F90
r3680 r5123 266 266 267 267 ! Creep limit depends on the size of the grid. 268 zdelta = MAX( SQRT( ztrace2 + ( ztrace2 - 4._wp * zdeter ) * usecc2 ), creepl)268 zdelta = MAX( SQRT( ztrace2 + ( ztrace2 - 4._wp * zdeter ) * usecc2 ), rn_creepl) 269 269 270 270 !- Computation of viscosities.
Note: See TracChangeset
for help on using the changeset viewer.