Changeset 12114


Ignore:
Timestamp:
2019-12-07T15:06:58+01:00 (6 months ago)
Author:
cetlod
Message:

Minor updates

Location:
NEMO/branches/2019/dev_r11643_ENHANCE-11_CEthe_Shaconemo_diags
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11643_ENHANCE-11_CEthe_Shaconemo_diags/cfgs/SHARED/namelist_ice_ref

    r11998 r12114  
    6060      rn_icebfr     =  15.            !        maximum bottom stress per unit volume [N/m3] 
    6161      rn_lfrelax    =   1.e-5         !        relaxation time scale to reach static friction [s-1] 
    62       rn_tensile    =   0.05          !        isotropic tensile strength [0-0.5??] 
     62      rn_tensile    =   0.2           !        isotropic tensile strength [0-0.5??] 
    6363/ 
    6464!------------------------------------------------------------------------------ 
     
    9191!------------------------------------------------------------------------------ 
    9292   ln_rhg_EVP       = .true.          !  EVP rheology 
    93       ln_aEVP       = .true.          !     adaptive rheology (Kimmritz et al. 2016 & 2017) 
     93      ln_aEVP       = .false.         !     adaptive rheology (Kimmritz et al. 2016 & 2017) 
    9494      rn_creepl     =   2.0e-9        !     creep limit [1/s] 
    9595      rn_ecc        =   2.0           !     eccentricity of the elliptical yield curve           
  • NEMO/branches/2019/dev_r11643_ENHANCE-11_CEthe_Shaconemo_diags/cfgs/SHARED/namelist_ref

    r11998 r12114  
    392392      ln_sssr_bnd =  .true.   !  flag to bound erp term (associated with nn_sssr=2) 
    393393      rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day] 
     394      nn_sssr_ice =   1       ! control of sea surface restoring under sea-ice 
     395                              ! 0 = no restoration under ice : * (1-icefrac) 
     396                              ! 1 = restoration everywhere  
     397                              ! >1 = enhanced restoration under ice : 1+(nn_icedmp-1)*icefrac 
    394398 
    395399   cn_dir      = './'      !  root directory for the SST/SSS data location 
  • NEMO/branches/2019/dev_r11643_ENHANCE-11_CEthe_Shaconemo_diags/src/OCE/SBC/sbcssr.F90

    r11995 r12114  
    4343   LOGICAL         ::   ln_sssr_bnd     ! flag to bound erp term  
    4444   REAL(wp)        ::   rn_sssr_bnd     ! ABS(Max./Min.) value of erp term [mm/day] 
    45    INTEGER         ::   nn_icedmp       ! Control of restoring under ice 
     45   INTEGER         ::   nn_sssr_ice     ! Control of restoring under ice 
    4646 
    4747   REAL(wp) , ALLOCATABLE, DIMENSION(:) ::   buffer   ! Temporary buffer for exchange 
     
    102102            ENDIF 
    103103            ! 
    104             IF( nn_sssr /= 0 .AND. nn_icedmp /= 1 ) THEN 
    105               ! use fraction of ice ( fr_i ) to adjust relaxation under ice if nn_icedmp .ne. 1 
    106               ! n.b. coefice is initialised and fixed to 1._wp if nn_icedmp = 1 
     104            IF( nn_sssr /= 0 .AND. nn_sssr_ice /= 1 ) THEN 
     105              ! use fraction of ice ( fr_i ) to adjust relaxation under ice if nn_sssr_ice .ne. 1 
     106              ! n.b. coefice is initialised and fixed to 1._wp if nn_sssr_ice = 1 
    107107               DO jj = 1, jpj 
    108108                  DO ji = 1, jpi 
    109                      SELECT CASE ( nn_icedmp ) 
     109                     SELECT CASE ( nn_sssr_ice ) 
    110110                       CASE ( 0 )    ;  coefice(ji,jj) = 1._wp - fr_i(ji,jj)              ! no/reduced damping under ice 
    111                        CASE  DEFAULT ;  coefice(ji,jj) = 1._wp +(nn_icedmp-1)*fr_i(ji,jj) ! reinforced damping (x nn_icedmp) under ice ) 
     111                       CASE  DEFAULT ;  coefice(ji,jj) = 1._wp + ( nn_sssr_ice - 1 ) * fr_i(ji,jj) ! reinforced damping (x nn_sssr_ice) under ice ) 
    112112                     END SELECT 
    113113                  END DO 
     
    170170      TYPE(FLD_N) ::   sn_sst, sn_sss        ! informations about the fields to be read 
    171171      NAMELIST/namsbc_ssr/ cn_dir, nn_sstr, nn_sssr, rn_dqdt, rn_deds, sn_sst, & 
    172               & sn_sss, ln_sssr_bnd, rn_sssr_bnd, nn_icedmp 
     172              & sn_sss, ln_sssr_bnd, rn_sssr_bnd, nn_sssr_ice 
    173173      INTEGER     ::  ios 
    174174      !!---------------------------------------------------------------------- 
     
    198198         WRITE(numout,*) '         flag to bound erp term                 ln_sssr_bnd = ', ln_sssr_bnd 
    199199         WRITE(numout,*) '         ABS(Max./Min.) erp threshold           rn_sssr_bnd = ', rn_sssr_bnd, ' mm/day' 
    200          WRITE(numout,*) '      Cntrl of surface restoration under ice nn_icedmp      = ', nn_icedmp 
     200         WRITE(numout,*) '      Cntrl of surface restoration under ice nn_sssr_ice    = ', nn_sssr_ice 
    201201         WRITE(numout,*) '          ( 0 = no restoration under ice)' 
    202202         WRITE(numout,*) '          ( 1 = restoration everywhere  )' 
     
    232232      ENDIF 
    233233      ! 
    234       coefice(:,:) = 1._wp         !  Initialise coefice to 1._wp ; will not need to be changed if nn_icedmp=1 
     234      coefice(:,:) = 1._wp         !  Initialise coefice to 1._wp ; will not need to be changed if nn_sssr_ice=1 
    235235      !                            !* Initialize qrp and erp if no restoring  
    236236      IF( nn_sstr /= 1                   )   qrp(:,:) = 0._wp 
Note: See TracChangeset for help on using the changeset viewer.