New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 4147 for branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/LIM_SRC_2/ice_2.F90 – NEMO

Ignore:
Timestamp:
2013-11-04T12:51:55+01:00 (10 years ago)
Author:
cetlod
Message:

merge in dev_LOCEAN_2013, the 1st development branch dev_r3853_CNRS9_Confsetting, from its starting point ( r3853 ) on the trunk: see ticket #1169

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/LIM_SRC_2/ice_2.F90

    r3625 r4147  
    2222   REAL(wp), PUBLIC ::   rdt_ice      !: ice time step 
    2323 
    24    !                                                                     !!* namicerun read in iceini  * 
    25    CHARACTER(len=32)     , PUBLIC ::   cn_icerst_in  = "restart_ice_in"   !: suffix of ice restart name (input) 
    26    CHARACTER(len=32)     , PUBLIC ::   cn_icerst_out = "restart_ice"      !: suffix of ice restart name (output) 
    27    LOGICAL               , PUBLIC ::   ln_limdyn     = .TRUE.             !: flag for ice dynamics (T) or not (F) 
    28    LOGICAL               , PUBLIC ::   ln_limdmp     = .FALSE.            !: Ice damping 
    29    LOGICAL               , PUBLIC ::   ln_nicep      = .TRUE.             !: flag grid points output (T) or not (F) 
    30    REAL(wp)              , PUBLIC ::   hsndif        = 0._wp              !: snow temp. computation (0) or not (9999) 
    31    REAL(wp)              , PUBLIC ::   hicdif        = 0._wp              !: ice  temp. computation (0) or not (9999) 
    32    REAL(wp), DIMENSION(2), PUBLIC ::   acrit = (/ 1.e-6_wp , 1.e-6_wp /)  !: minimum lead fraction in the 2 hemisphere 
    33     
    34    !                                          !!* ice-dynamic namelist (namicedyn) * 
    35    INTEGER , PUBLIC ::   nbiter =    1         !: number of sub-time steps for relaxation 
    36    INTEGER , PUBLIC ::   nbitdr =  250         !: maximum number of iterations for relaxation 
    37    INTEGER , PUBLIC ::   nevp   =  360         !: number of EVP subcycling iterations 
    38    INTEGER , PUBLIC ::   telast = 3600         !: timescale for EVP elastic waves 
    39    REAL(wp), PUBLIC ::   epsd   = 1.0e-20_wp   !: tolerance parameter for dynamic 
    40    REAL(wp), PUBLIC ::   alpha  = 0.5_wp       !: coefficient for semi-implicit coriolis 
    41    REAL(wp), PUBLIC ::   dm     = 0.6e+03_wp   !: diffusion constant for dynamics 
    42    REAL(wp), PUBLIC ::   om     = 0.5_wp       !: relaxation constant 
    43    REAL(wp), PUBLIC ::   resl   = 5.0e-05_wp   !: maximum value for the residual of relaxation 
    44    REAL(wp), PUBLIC ::   cw     = 5.0e-03_wp   !: drag coefficient for oceanic stress 
    45    REAL(wp), PUBLIC ::   angvg  = 0._wp        !: turning angle for oceanic stress 
    46    REAL(wp), PUBLIC ::   pstar  = 1.0e+04_wp   !: first bulk-rheology parameter 
    47    REAL(wp), PUBLIC ::   c_rhg  = 20._wp       !: second bulk-rhelogy parameter 
    48    REAL(wp), PUBLIC ::   etamn  = 0._wp        !: minimun value for viscosity 
    49    REAL(wp), PUBLIC ::   creepl = 2.e-08_wp    !: creep limit 
    50    REAL(wp), PUBLIC ::   ecc    = 2._wp        !: eccentricity of the elliptical yield curve 
    51    REAL(wp), PUBLIC ::   ahi0   = 350._wp      !: sea-ice hor. eddy diffusivity coeff. (m2/s) 
    52    REAL(wp), PUBLIC ::   alphaevp = 1._wp      !: coefficient for the solution of EVP int. stresses 
     24   !                                                   !!* namicerun read in iceini  * 
     25   CHARACTER(len=32)     , PUBLIC ::   cn_icerst_in     !: suffix of ice restart name (input) 
     26   CHARACTER(len=32)     , PUBLIC ::   cn_icerst_out    !: suffix of ice restart name (output) 
     27   LOGICAL               , PUBLIC ::   ln_limdyn        !: flag for ice dynamics (T) or not (F) 
     28   LOGICAL               , PUBLIC ::   ln_limdmp        !: Ice damping 
     29   REAL(wp)              , PUBLIC ::   hsndif           !: snow temp. computation (0) or not (9999) 
     30   REAL(wp)              , PUBLIC ::   hicdif           !: ice temp. computation (0) or not (9999) 
     31   REAL(wp), DIMENSION(2), PUBLIC ::   acrit            !: minimum lead fraction in the 2 hemisphere 
     32   ! 
     33   LOGICAL               , PUBLIC ::   ln_nicep      = .TRUE.     !: flag grid points output (T) or not (F) 
     34   !                                !!* ice-dynamic namelist (namicedyn) * 
     35   INTEGER , PUBLIC ::   nbiter      !: number of sub-time steps for relaxation 
     36   INTEGER , PUBLIC ::   nbitdr      !: maximum number of iterations for relaxation 
     37   INTEGER , PUBLIC ::   nevp        !: number of EVP subcycling iterations 
     38   INTEGER , PUBLIC ::   telast      !: timescale for EVP elastic waves 
     39   REAL(wp), PUBLIC ::   epsd        !: tolerance parameter for dynamic 
     40   REAL(wp), PUBLIC ::   alpha       !: coefficient for semi-implicit coriolis 
     41   REAL(wp), PUBLIC ::   dm          !: diffusion constant for dynamics 
     42   REAL(wp), PUBLIC ::   om          !: relaxation constant 
     43   REAL(wp), PUBLIC ::   resl        !: maximum value for the residual of relaxation 
     44   REAL(wp), PUBLIC ::   cw          !: drag coefficient for oceanic stress 
     45   REAL(wp), PUBLIC ::   angvg       !: turning angle for oceanic stress 
     46   REAL(wp), PUBLIC ::   pstar       !: first bulk-rheology parameter 
     47   REAL(wp), PUBLIC ::   c_rhg       !: second bulk-rhelogy parameter 
     48   REAL(wp), PUBLIC ::   etamn       !: minimun value for viscosity 
     49   REAL(wp), PUBLIC ::   creepl      !: creep limit 
     50   REAL(wp), PUBLIC ::   ecc         !: eccentricity of the elliptical yield curve 
     51   REAL(wp), PUBLIC ::   ahi0        !: sea-ice hor. eddy diffusivity coeff. (m2/s) 
     52   REAL(wp), PUBLIC ::   alphaevp    !: coefficient for the solution of EVP int. stresses 
    5353 
    5454   REAL(wp), PUBLIC ::   usecc2                !:  = 1.0 / ( ecc * ecc ) 
Note: See TracChangeset for help on using the changeset viewer.