Changeset 5122


Ignore:
Timestamp:
2015-03-04T15:34:45+01:00 (6 years ago)
Author:
vancop
Message:

fixes to verify SETTE tests. Evt ok, except SAS (no solver) and ORCA_AGRIF_LIM (compiler crash)

Location:
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/CONFIG/SHARED/namelist_ice_lim2_ref

    r5054 r5122  
    4848   c_rhg       =  20.0     !  2nd bulk-rhelogy parameter 
    4949   etamn       =   0.0e+07 !  minimun value for viscosity 
    50    creepl      =   1.0e-08 !  creep limit 
    51    ecc         =   2.0     !  eccentricity of the elliptical yield curve 
     50   rn_creepl   =   1.0e-08 !  creep limit 
     51   rn_ecc      =   2.0     !  eccentricity of the elliptical yield curve 
    5252   ahi0        = 350.e0    !  horizontal eddy diffusivity coefficient for sea-ice [m2/s] 
    53    nevp        =   120     !  number of EVP subcycling iterations 
     53   nn_nevp     =   120     !  number of EVP subcycling iterations 
    5454   telast      =   9600    !  timescale for EVP elastic waves 
    5555   alphaevp    =   1.0     !  coefficient for the solution of EVP int. stresses 
  • branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_2/ice_2.F90

    r5058 r5122  
    3535   INTEGER , PUBLIC ::   nbiter      !: number of sub-time steps for relaxation 
    3636   INTEGER , PUBLIC ::   nbitdr      !: maximum number of iterations for relaxation 
    37    INTEGER , PUBLIC ::   nevp        !: number of EVP subcycling iterations 
     37   INTEGER , PUBLIC ::   nn_nevp     !: number of EVP subcycling iterations 
    3838   INTEGER , PUBLIC ::   telast      !: timescale for EVP elastic waves 
    3939   REAL(wp), PUBLIC ::   epsd        !: tolerance parameter for dynamic 
     
    4747   REAL(wp), PUBLIC ::   c_rhg       !: second bulk-rhelogy parameter 
    4848   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 
     49   REAL(wp), PUBLIC ::   rn_creepl   !: creep limit 
     50   REAL(wp), PUBLIC ::   rn_ecc      !: eccentricity of the elliptical yield curve 
    5151   REAL(wp), PUBLIC ::   ahi0        !: sea-ice hor. eddy diffusivity coeff. (m2/s) 
    5252   REAL(wp), PUBLIC ::   alphaevp    !: coefficient for the solution of EVP int. stresses 
    5353 
    54    REAL(wp), PUBLIC ::   usecc2                !:  = 1.0 / ( ecc * ecc ) 
     54   REAL(wp), PUBLIC ::   usecc2                !:  = 1.0 / ( rn_ecc * rn_ecc ) 
    5555   REAL(wp), PUBLIC ::   rhoco                 !: = rau0 * cw 
    5656   REAL(wp), PUBLIC ::   sangvg, cangvg        !: sin and cos of the turning angle for ocean stress 
  • branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_2/limdyn_2.F90

    r5053 r5122  
    227227      NAMELIST/namicedyn/ epsd, alpha,     & 
    228228         &                dm, nbiter, nbitdr, om, resl, cw, angvg, pstar,   & 
    229          &                c_rhg, etamn, creepl, ecc, ahi0,                  & 
    230          &                nevp, telast, alphaevp 
     229         &                c_rhg, etamn, rn_creepl, rn_ecc, ahi0,                  & 
     230         &                nn_nevp, telast, alphaevp 
    231231      !!------------------------------------------------------------------- 
    232232                     
     
    256256         WRITE(numout,*) '       second bulk-rhelogy parameter                    c_rhg  = ', c_rhg 
    257257         WRITE(numout,*) '       minimun value for viscosity                      etamn  = ', etamn 
    258          WRITE(numout,*) '       creep limit                                      creepl = ', creepl 
    259          WRITE(numout,*) '       eccentricity of the elliptical yield curve       ecc    = ', ecc 
     258         WRITE(numout,*) '       creep limit                                      rn_creepl = ', rn_creepl 
     259         WRITE(numout,*) '       eccentricity of the elliptical yield curve       rn_ecc = ', rn_ecc 
    260260         WRITE(numout,*) '       horizontal diffusivity coeff. for sea-ice        ahi0   = ', ahi0 
    261          WRITE(numout,*) '       number of iterations for subcycling nevp   = ', nevp 
     261         WRITE(numout,*) '       number of iterations for subcycling              nn_nevp= ', nn_nevp 
    262262         WRITE(numout,*) '       timescale for elastic waves telast = ', telast 
    263263         WRITE(numout,*) '       coefficient for the solution of int. stresses alphaevp = ', alphaevp 
     
    271271 
    272272      !  Initialization 
    273       usecc2 = 1.0 / ( ecc * ecc ) 
     273      usecc2 = 1.0 / ( rn_ecc * rn_ecc ) 
    274274      rhoco  = rau0 * cw 
    275275      angvg  = angvg * rad      ! convert angvg from degree to radian 
  • branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_2/limrhg_2.F90

    r3680 r5122  
    266266 
    267267               !  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) 
    269269 
    270270               !-  Computation of viscosities. 
  • branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limrhg.F90

    r5070 r5122  
    102102      !!                 and charge ellipse. 
    103103      !!                 The user should make sure that the parameters 
    104       !!                 nevp, elastic time scale and creepl maintain stress state 
     104      !!                 nn_nevp, elastic time scale and rn_creepl maintain stress state 
    105105      !!                 on the charge ellipse for plastic flow 
    106106      !!                 e.g. in the Canadian Archipelago 
  • branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90

    r5078 r5122  
    655655      WRITE(*,*) 'sbc_ice_lim: You should not have seen this print! error?', kt, kblk 
    656656   END SUBROUTINE sbc_ice_lim 
     657   SUBROUTINE sbc_lim_init                 ! Dummy routine 
     658   END SUBROUTINE sbc_lim_init 
    657659#endif 
    658660 
Note: See TracChangeset for help on using the changeset viewer.