Changeset 2580


Ignore:
Timestamp:
2011-02-08T09:01:02+01:00 (10 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:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist_ice_lim2

    r2528 r2580  
    4848   c_rhg       =  20.0     !  2nd bulk-rhelogy parameter 
    4949   etamn       =   0.0e+07 !  minimun value for viscosity 
    50    creepl      =   2.0e-08 !  creep limit 
     50   creepl      =   1.0e-08 !  creep limit 
    5151   ecc         =   2.0     !  eccentricity of the elliptical yield curve 
    5252   ahi0        = 350.e0    !  horizontal eddy diffusivity coefficient for sea-ice [m2/s] 
    53    nevp        =   360     !  number of EVP subcycling iterations 
    54    telast      =   3600    !  timescale for EVP elastic waves 
     53   nevp        =   120     !  number of EVP subcycling iterations 
     54   telast      =   9600    !  timescale for EVP elastic waves 
    5555   alphaevp    =   1.0     !  coefficient for the solution of EVP int. stresses 
    5656/ 
  • trunk/NEMOGCM/NEMO/LIM_SRC_3/limrhg.F90

    r2528 r2580  
    190190            zpresh(ji,jj) = tms(ji,jj) *  strength(ji,jj) * 0.5_wp 
    191191#endif 
     192#if defined key_lim2 
     193            zpresh(ji,jj) = tms(ji,jj) *  pstar * vt_i(ji,jj) * EXP( -c_rhg * (1. - at_i(ji,jj) ) ) 
     194#endif 
    192195            ! tmi = 1 where there is ice or on land 
    193196            tmi(ji,jj)    = 1._wp - ( 1._wp - MAX( 0._wp , SIGN ( 1._wp , vt_i(ji,jj) - epsd ) ) ) * tms(ji,jj) 
Note: See TracChangeset for help on using the changeset viewer.