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 13346 – NEMO

Changeset 13346


Ignore:
Timestamp:
2020-07-27T14:52:36+02:00 (4 years ago)
Author:
clem
Message:

change namelist parameter ln_rhg_chkcvg into nn_rhg_chkcvg to ensure more flexibility

Location:
NEMO/releases/r4.0/r4.0-HEAD
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • NEMO/releases/r4.0/r4.0-HEAD/cfgs/SHARED/namelist_ice_ref

    r13343 r13346  
    9898      rn_relast     =   0.333         !     ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast  
    9999                                      !        advised value: 1/3 (nn_nevp=100) or 1/9 (nn_nevp=300) 
    100    ln_rhg_chkcvg    = .false.         !  check convergence of rheology (outputs: file ice_cvg.nc & variable uice_cvg) 
     100   nn_rhg_chkcvg    =   0             !  check convergence of rheology 
     101                                      !     = 0  no check 
     102                                      !     = 1  check at the main time step (output xml: uice_cvg) 
     103                                      !     = 2  check at both main and rheology time steps (additional output: ice_cvg.nc) 
     104                                      !          this option 2 asks a lot of communications between cpu 
    101105/ 
    102106!------------------------------------------------------------------------------ 
  • NEMO/releases/r4.0/r4.0-HEAD/src/ICE/ice.F90

    r13284 r13346  
    155155   INTEGER , PUBLIC ::   nn_nevp          !: number of iterations for subcycling 
    156156   REAL(wp), PUBLIC ::   rn_relast        !: ratio => telast/rdt_ice (1/3 or 1/9 depending on nb of subcycling nevp)  
    157    LOGICAL , PUBLIC ::   ln_rhg_chkcvg    !: check ice rheology convergence  
     157   INTEGER , PUBLIC ::   nn_rhg_chkcvg    !: check ice rheology convergence  
    158158   ! 
    159159   !                                     !!** ice-advection namelist (namdyn_adv) ** 
  • NEMO/releases/r4.0/r4.0-HEAD/src/ICE/icedyn_rhg.F90

    r13284 r13346  
    110110      INTEGER ::   ios, ioptio   ! Local integer output status for namelist read 
    111111      !! 
    112       NAMELIST/namdyn_rhg/  ln_rhg_EVP, ln_aEVP, rn_creepl, rn_ecc , nn_nevp, rn_relast, ln_rhg_chkcvg 
     112      NAMELIST/namdyn_rhg/  ln_rhg_EVP, ln_aEVP, rn_creepl, rn_ecc , nn_nevp, rn_relast, nn_rhg_chkcvg 
    113113      !!------------------------------------------------------------------- 
    114114      ! 
     
    132132         WRITE(numout,*) '         number of iterations for subcycling               nn_nevp       = ', nn_nevp 
    133133         WRITE(numout,*) '         ratio of elastic timescale over ice time step     rn_relast     = ', rn_relast 
    134          WRITE(numout,*) '      check convergence of rheology                        ln_rhg_chkcvg = ', ln_rhg_chkcvg 
     134         WRITE(numout,*) '      check convergence of rheology                        nn_rhg_chkcvg = ', nn_rhg_chkcvg 
     135         IF    ( nn_rhg_chkcvg == 0 ) THEN   ;   WRITE(numout,*) '         no check' 
     136         ELSEIF( nn_rhg_chkcvg == 1 ) THEN   ;   WRITE(numout,*) '         check cvg at the main time step' 
     137         ELSEIF( nn_rhg_chkcvg == 2 ) THEN   ;   WRITE(numout,*) '         check cvg at both main and rheology time steps' 
     138         ENDIF 
    135139      ENDIF 
    136140      ! 
  • NEMO/releases/r4.0/r4.0-HEAD/src/ICE/icedyn_rhg_evp.F90

    r13284 r13346  
    362362         ! 
    363363         ! convergence test 
    364          IF( ln_rhg_chkcvg ) THEN 
     364         IF( nn_rhg_chkcvg == 1 .OR. nn_rhg_chkcvg == 2 ) THEN 
    365365            DO jj = 1, jpj 
    366366               DO ji = 1, jpi 
     
    715715 
    716716         ! convergence test 
    717          IF( ln_rhg_chkcvg )   CALL rhg_cvg( kt, jter, nn_nevp, u_ice, v_ice, zu_ice, zv_ice ) 
     717         IF( nn_rhg_chkcvg == 2 )   CALL rhg_cvg( kt, jter, nn_nevp, u_ice, v_ice, zu_ice, zv_ice ) 
    718718         ! 
    719719         !                                                ! ==================== ! 
     
    890890      ! 
    891891      ! --- convergence tests --- ! 
    892       IF( ln_rhg_chkcvg ) THEN 
     892      IF( nn_rhg_chkcvg == 1 .OR. nn_rhg_chkcvg == 2 ) THEN 
    893893         IF( iom_use('uice_cvg') ) THEN 
    894894            IF( ln_aEVP ) THEN   ! output: beta * ( u(t=nn_nevp) - u(t=nn_nevp-1) ) 
Note: See TracChangeset for help on using the changeset viewer.