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 4659 for branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3 – NEMO

Ignore:
Timestamp:
2014-06-04T19:18:57+02:00 (10 years ago)
Author:
clem
Message:

cleaning namelist parameters in LIM3 for NEMO3.4

Location:
branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/ice.F90

    r4649 r4659  
    166166 
    167167   !                                          !!** ice-dynamic namelist (namicedyn) ** 
    168    INTEGER , PUBLIC ::   nbiter = 1            !: number of sub-time steps for relaxation 
    169    INTEGER , PUBLIC ::   nbitdr = 250          !: maximum number of iterations for relaxation 
    170168   INTEGER , PUBLIC ::   nevp   = 400          !: number of iterations for subcycling 
    171169   INTEGER , PUBLIC ::   nlay_i = 5            !: number of layers in the ice 
    … …  
    173171   !                                          !!** ice-dynamic namelist (namicedyn) ** 
    174172   REAL(wp), PUBLIC ::   epsd   = 1.0e-20_wp   !: tolerance parameter for dynamic 
    175    REAL(wp), PUBLIC ::   alpha  = 0.5_wp       !: coefficient for semi-implicit coriolis 
    176    REAL(wp), PUBLIC ::   dm     = 0.6e+03_wp   !: diffusion constant for dynamics 
    177173   REAL(wp), PUBLIC ::   om     = 0.5_wp       !: relaxation constant 
    178    REAL(wp), PUBLIC ::   resl   = 5.0e-05_wp   !: maximum value for the residual of relaxation 
    179174   REAL(wp), PUBLIC ::   cw     = 5.0e-03_wp   !: drag coefficient for oceanic stress 
    180175   REAL(wp), PUBLIC ::   angvg  = 0._wp        !: turning angle for oceanic stress 
    181176   REAL(wp), PUBLIC ::   pstar  = 1.0e+04_wp   !: determines ice strength (N/M), Hibler JPO79 
    182177   REAL(wp), PUBLIC ::   c_rhg  = 20._wp       !: determines changes in ice strength 
    183    REAL(wp), PUBLIC ::   etamn  = 0.0e+07_wp   !: minimun value for viscosity : has to be 0 
    184178   REAL(wp), PUBLIC ::   creepl = 2.0e-08_wp   !: creep limit : has to be under 1.0e-9 
    185179   REAL(wp), PUBLIC ::   ecc    = 2._wp        !: eccentricity of the elliptical yield curve 
    … …  
    220214   REAL(wp), PUBLIC ::   Craft     = 5.0_wp        !: coefficient for smoothness of the hyperbolic tangent in rafting 
    221215   REAL(wp), PUBLIC ::   ridge_por = 0.0_wp        !: initial porosity of ridges (0.3 regular value) 
    222    REAL(wp), PUBLIC ::   sal_max_ridge = 15.0_wp   !: maximum ridged ice salinity (ppt) 
    223    REAL(wp), PUBLIC ::   betas    = 1.0_wp         !: coef. for partitioning of snowfall between leads and sea ice 
    224    REAL(wp), PUBLIC ::   kappa_i  = 1.0_wp         !: coef. for the extinction of radiation Grenfell et al. (2006) [1/m] 
     216   REAL(wp), PUBLIC ::   betas     = 1.0_wp        !: coef. for partitioning of snowfall between leads and sea ice 
     217   REAL(wp), PUBLIC ::   kappa_i   = 1.0_wp        !: coef. for the extinction of radiation Grenfell et al. (2006) [1/m] 
    225218   REAL(wp), PUBLIC ::   nconv_i_thd = 50_wp       !: maximal number of iterations for heat diffusion 
    226219   REAL(wp), PUBLIC ::   maxer_i_thd = 1.0e-4_wp   !: maximal tolerated error (C) for heat diffusion 
    … …  
    228221   !                                              !!** ice-mechanical redistribution namelist (namiceitdme) 
    229222   INTEGER , PUBLIC ::   ridge_scheme_swi = 0      !: scheme used for ice ridging 
    230    INTEGER , PUBLIC ::   raftswi          = 1      !: rafting of ice or not                         
     223   INTEGER , PUBLIC ::   raft_swi         = 1      !: rafting of ice or not                         
    231224   INTEGER , PUBLIC ::   partfun_swi      = 1      !: participation function: =0 Thorndike et al. (1975), =1 Lipscomb et al. (2007) 
    232    INTEGER , PUBLIC ::   transfun_swi     = 0      !: transfer function: =0 Hibler 1980, =1 Lipscomb et al. 2007 
    233225   INTEGER , PUBLIC ::   brinstren_swi    = 0      !: use brine volume to diminish ice strength 
    234226 

  • branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limdyn.F90

    r4649 r4659  
    240240      !! ** input   :   Namelist namicedyn 
    241241      !!------------------------------------------------------------------- 
    242       NAMELIST/namicedyn/ epsd, alpha,     & 
    243          &                dm, nbiter, nbitdr, om, resl, cw, angvg, pstar,   & 
    244          &                c_rhg, etamn, creepl, ecc, ahi0, & 
     242      NAMELIST/namicedyn/ epsd, om, cw, angvg, pstar,   & 
     243         &                c_rhg, creepl, ecc, ahi0,     & 
    245244         &                nevp, telast, alphaevp, hminrhg 
    246245      !!------------------------------------------------------------------- 
    … …  
    254253         WRITE(numout,*) '~~~~~~~~~~~~' 
    255254         WRITE(numout,*) '   tolerance parameter                              epsd   = ', epsd 
    256          WRITE(numout,*) '   coefficient for semi-implicit coriolis           alpha  = ', alpha 
    257          WRITE(numout,*) '   diffusion constant for dynamics                  dm     = ', dm 
    258          WRITE(numout,*) '   number of sub-time steps for relaxation          nbiter = ', nbiter 
    259          WRITE(numout,*) '   maximum number of iterations for relaxation      nbitdr = ', nbitdr 
    260255         WRITE(numout,*) '   relaxation constant                              om     = ', om 
    261          WRITE(numout,*) '   maximum value for the residual of relaxation     resl   = ', resl 
    262256         WRITE(numout,*) '   drag coefficient for oceanic stress              cw     = ', cw 
    263257         WRITE(numout,*) '   turning angle for oceanic stress                 angvg  = ', angvg 
    264258         WRITE(numout,*) '   first bulk-rheology parameter                    pstar  = ', pstar 
    265259         WRITE(numout,*) '   second bulk-rhelogy parameter                    c_rhg  = ', c_rhg 
    266          WRITE(numout,*) '   minimun value for viscosity                      etamn  = ', etamn 
    267260         WRITE(numout,*) '   creep limit                                      creepl = ', creepl 
    268261         WRITE(numout,*) '   eccentricity of the elliptical yield curve       ecc    = ', ecc 

  • branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90

    r4649 r4659  
    618618      !!---------------------------------------------------------------------! 
    619619      INTEGER ::   ji,jj, jl    ! dummy loop indices 
    620       INTEGER ::   krdg_index   !  
    621620      REAL(wp) ::   Gstari, astari, hi, hrmean, zdummy   ! local scalar 
    622621      REAL(wp), POINTER, DIMENSION(:,:)   ::   zworka    ! temporary array used here 
    … …  
    694693      !----------------------------------------------------------------- 
    695694 
    696       krdg_index = 1 
    697  
    698       IF( krdg_index == 0 ) THEN       !--- Linear formulation (Thorndike et al., 1975) 
     695      IF( partfun_swi == 0 ) THEN       !--- Linear formulation (Thorndike et al., 1975) 
    699696         DO jl = 0, ice_cat_bounds(1,2)       ! only undeformed ice participates 
    700697            DO jj = 1, jpj  
    … …  
    724721         END DO 
    725722         ! 
    726       ENDIF ! krdg_index 
    727  
    728       IF( raftswi == 1 ) THEN      ! Ridging and rafting ice participation functions 
     723      ENDIF ! partfun_swi 
     724 
     725      IF( raft_swi == 1 ) THEN      ! Ridging and rafting ice participation functions 
    729726         ! 
    730727         DO jl = 1, jpl 
    … …  
    742739         END DO ! jl 
    743740 
    744       ELSE  ! raftswi = 0 
     741      ELSE  ! raft_swi = 0 
    745742         ! 
    746743         DO jl = 1, jpl 
    … …  
    750747      ENDIF 
    751748 
    752       IF ( raftswi == 1 ) THEN 
     749      IF ( raft_swi == 1 ) THEN 
    753750 
    754751         IF( MAXVAL(aridge + araft - athorn(:,:,1:jpl)) .GT. epsi10 ) THEN 
    … …  
    13481345      !! ** input   :   Namelist namiceitdme 
    13491346      !!------------------------------------------------------------------- 
    1350       NAMELIST/namiceitdme/ ridge_scheme_swi, Cs, Cf, fsnowrdg, fsnowrft,&  
    1351          Gstar, astar,                                & 
    1352          Hstar, raftswi, hparmeter, Craft, ridge_por, & 
    1353          sal_max_ridge,  partfun_swi, transfun_swi,   & 
    1354          brinstren_swi 
     1347      NAMELIST/namiceitdme/ ridge_scheme_swi, Cs, Cf, fsnowrdg, fsnowrft,              &  
     1348        &                   Gstar, astar, Hstar, raft_swi, hparmeter, Craft, ridge_por, & 
     1349        &                   partfun_swi, brinstren_swi 
    13551350      !!------------------------------------------------------------------- 
    13561351      ! 
    … …  
    13701365         WRITE(numout,*)'   Equivalent to G* for an exponential part function       astar           ', astar 
    13711366         WRITE(numout,*)'   Quantity playing a role in max ridged ice thickness     Hstar           ', Hstar 
    1372          WRITE(numout,*)'   Rafting of ice sheets or not                            raftswi         ', raftswi 
     1367         WRITE(numout,*)'   Rafting of ice sheets or not                            raft_swi        ', raft_swi 
    13731368         WRITE(numout,*)'   Parmeter thickness (threshold between ridge-raft)       hparmeter       ', hparmeter 
    13741369         WRITE(numout,*)'   Rafting hyperbolic tangent coefficient                  Craft           ', Craft   
    13751370         WRITE(numout,*)'   Initial porosity of ridges                              ridge_por       ', ridge_por 
    1376          WRITE(numout,*)'   Maximum salinity of ridging ice                         sal_max_ridge   ', sal_max_ridge 
    13771371         WRITE(numout,*)'   Switch for part. function (0) linear (1) exponential    partfun_swi     ', partfun_swi 
    1378          WRITE(numout,*)'   Switch for tran. function (0) linear (1) exponential    transfun_swi    ', transfun_swi 
    13791372         WRITE(numout,*)'   Switch for including brine volume in ice strength comp. brinstren_swi   ', brinstren_swi 
    13801373      ENDIF 

  • branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90

    r4649 r4659  
    610610      !!------------------------------------------------------------------- 
    611611      NAMELIST/namicethd/ hmelt , hiccrit, fraz_swi, maxfrazb, vfrazb, Cfrazb,   & 
    612          &                hicmin, hiclim,                                        & 
    613          &                sbeta  , parlat, hakspl, hibspl, exld,                 & 
    614          &                hakdif, hnzst  , thth  , parsub, alphs, betas,         &  
     612         &                hiclim, hnzst, parsub, betas,                          &  
    615613         &                kappa_i, nconv_i_thd, maxer_i_thd, thcon_i_swi 
    616614      !!------------------------------------------------------------------- 
    … …  
    629627         WRITE(numout,*)'   Namelist of ice parameters for ice thermodynamic computation ' 
    630628         WRITE(numout,*)'      maximum melting at the bottom                           hmelt        = ', hmelt 
    631          WRITE(numout,*)'      ice thick. for lateral accretion in NH (SH)             hiccrit(1/2) = ', hiccrit 
     629         WRITE(numout,*)'      ice thick. for lateral accretion                        hiccrit      = ', hiccrit 
    632630         WRITE(numout,*)'      Frazil ice thickness as a function of wind or not       fraz_swi     = ', fraz_swi 
    633631         WRITE(numout,*)'      Maximum proportion of frazil ice collecting at bottom   maxfrazb     = ', maxfrazb 
    634632         WRITE(numout,*)'      Thresold relative drift speed for collection of frazil  vfrazb       = ', vfrazb 
    635633         WRITE(numout,*)'      Squeezing coefficient for collection of frazil          Cfrazb       = ', Cfrazb 
    636          WRITE(numout,*)'      ice thick. corr. to max. energy stored in brine pocket  hicmin       = ', hicmin   
    637634         WRITE(numout,*)'      minimum ice thickness                                   hiclim       = ', hiclim  
    638635         WRITE(numout,*)'      numerical carac. of the scheme for diffusion in ice ' 
    639          WRITE(numout,*)'      Cranck-Nicholson (=0.5), implicit (=1), explicit (=0)   sbeta        = ', sbeta 
    640          WRITE(numout,*)'      percentage of energy used for lateral ablation          parlat       = ', parlat 
    641          WRITE(numout,*)'      slope of distr. for Hakkinen-Mellor lateral melting     hakspl       = ', hakspl   
    642          WRITE(numout,*)'      slope of distribution for Hibler lateral melting        hibspl       = ', hibspl 
    643          WRITE(numout,*)'      exponent for leads-closure rate                         exld         = ', exld 
    644          WRITE(numout,*)'      coefficient for diffusions of ice and snow              hakdif       = ', hakdif 
    645          WRITE(numout,*)'      threshold thick. for comp. of eq. thermal conductivity  zhth         = ', thth  
    646636         WRITE(numout,*)'      thickness of the surf. layer in temp. computation       hnzst        = ', hnzst 
    647637         WRITE(numout,*)'      switch for snow sublimation  (=1) or not (=0)           parsub       = ', parsub   
    648          WRITE(numout,*)'      coefficient for snow density when snow ice formation    alphs        = ', alphs 
    649638         WRITE(numout,*)'      coefficient for ice-lead partition of snowfall          betas        = ', betas 
    650639         WRITE(numout,*)'      extinction radiation parameter in sea ice (1.0)         kappa_i      = ', kappa_i 
    … …  
    655644      ENDIF 
    656645      ! 
    657       rcdsn = hakdif * rcdsn  
    658       rcdic = hakdif * rcdic 
    659       ! 
    660646   END SUBROUTINE lim_thd_init 
    661647 

  • branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limthd_lac.F90

    r4649 r4659  
    157157 
    158158      ! Default new ice thickness  
    159       hicol(:,:) = hiccrit(1) 
    160  
    161       IF( fraz_swi == 1._wp ) THEN 
     159      hicol(:,:) = hiccrit 
     160 
     161      IF( fraz_swi == 1 ) THEN 
    162162 
    163163         !-------------------- 
    … …  
    322322         !---------------------- 
    323323         DO ji = 1, nbpac 
    324             zh_newice(ji) = hiccrit(1) 
    325          END DO 
    326          IF( fraz_swi == 1.0 ) zh_newice(:) = hicol_b(:) 
     324            zh_newice(ji) = hiccrit 
     325         END DO 
     326         IF( fraz_swi == 1 ) zh_newice(:) = hicol_b(:) 
    327327 
    328328         !---------------------- 

  • branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/thd_ice.F90

    r4649 r4659  
    2020   !                                         !!! ** ice-thermo namelist (namicethd) ** 
    2121   REAL(wp), PUBLIC ::   hmelt   = -0.15     !: maximum melting at the bottom; active only for one category 
    22    REAL(wp), PUBLIC ::   hicmin  = 0.2       !: (REMOVE) 
    2322   REAL(wp), PUBLIC ::   hiclim  = 0.05      !: minimum ice thickness 
    24    REAL(wp), PUBLIC ::   sbeta   = 1.0       !: numerical scheme for diffusion in ice  (REMOVE) 
    25    REAL(wp), PUBLIC ::   parlat  = 0.0       !: (REMOVE) 
    26    REAL(wp), PUBLIC ::   hakspl  = 0.5       !: (REMOVE) 
    27    REAL(wp), PUBLIC ::   hibspl  = 0.5       !: (REMOVE) 
    28    REAL(wp), PUBLIC ::   exld    = 2.0       !: (REMOVE) 
    29    REAL(wp), PUBLIC ::   hakdif  = 1.0       !: (REMOVE) 
    30    REAL(wp), PUBLIC ::   thth    = 0.2       !: (REMOVE) 
    3123   REAL(wp), PUBLIC ::   hnzst   = 0.1       !: thick. of the surf. layer in temp. comp. 
    3224   REAL(wp), PUBLIC ::   parsub  = 1.0       !: switch for snow sublimation or not 
    33    REAL(wp), PUBLIC ::   alphs   = 1.0       !: coef. for snow density when snow-ice formation 
    34    REAL(wp), PUBLIC ::   fraz_swi= 1.0       !: use of frazil ice collection in function of wind (1.0) or not (0.0) 
    3525   REAL(wp), PUBLIC ::   maxfrazb= 0.7       !: maximum portion of frazil ice collecting at the ice bottom 
    3626   REAL(wp), PUBLIC ::   vfrazb  = 0.41667   !: threshold drift speed for collection of bottom frazil ice 
    3727   REAL(wp), PUBLIC ::   Cfrazb  = 5.0       !: squeezing coefficient for collection of bottom frazil ice 
     28   REAL(wp), PUBLIC ::   hiccrit = 0.3       !: ice th. for lateral accretion in the NH (SH) (m) 
    3829 
    39    REAL(wp), PUBLIC, DIMENSION(2) ::   hiccrit = (/0.3,0.3/)  !: ice th. for lateral accretion in the NH (SH) (m) 
     30   INTEGER , PUBLIC ::   fraz_swi= 1         !: use of frazil ice collection in function of wind (1) or not (0) 
    4031 
    4132   !!----------------------------- 
Note: See TracChangeset for help on using the changeset viewer.