Changeset 13346
- Timestamp:
- 2020-07-27T14:52:36+02:00 (4 years ago)
- 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 98 98 rn_relast = 0.333 ! ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast 99 99 ! 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 101 105 / 102 106 !------------------------------------------------------------------------------ -
NEMO/releases/r4.0/r4.0-HEAD/src/ICE/ice.F90
r13284 r13346 155 155 INTEGER , PUBLIC :: nn_nevp !: number of iterations for subcycling 156 156 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 convergence157 INTEGER , PUBLIC :: nn_rhg_chkcvg !: check ice rheology convergence 158 158 ! 159 159 ! !!** ice-advection namelist (namdyn_adv) ** -
NEMO/releases/r4.0/r4.0-HEAD/src/ICE/icedyn_rhg.F90
r13284 r13346 110 110 INTEGER :: ios, ioptio ! Local integer output status for namelist read 111 111 !! 112 NAMELIST/namdyn_rhg/ ln_rhg_EVP, ln_aEVP, rn_creepl, rn_ecc , nn_nevp, rn_relast, ln_rhg_chkcvg112 NAMELIST/namdyn_rhg/ ln_rhg_EVP, ln_aEVP, rn_creepl, rn_ecc , nn_nevp, rn_relast, nn_rhg_chkcvg 113 113 !!------------------------------------------------------------------- 114 114 ! … … 132 132 WRITE(numout,*) ' number of iterations for subcycling nn_nevp = ', nn_nevp 133 133 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 135 139 ENDIF 136 140 ! -
NEMO/releases/r4.0/r4.0-HEAD/src/ICE/icedyn_rhg_evp.F90
r13284 r13346 362 362 ! 363 363 ! convergence test 364 IF( ln_rhg_chkcvg) THEN364 IF( nn_rhg_chkcvg == 1 .OR. nn_rhg_chkcvg == 2 ) THEN 365 365 DO jj = 1, jpj 366 366 DO ji = 1, jpi … … 715 715 716 716 ! 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 ) 718 718 ! 719 719 ! ! ==================== ! … … 890 890 ! 891 891 ! --- convergence tests --- ! 892 IF( ln_rhg_chkcvg) THEN892 IF( nn_rhg_chkcvg == 1 .OR. nn_rhg_chkcvg == 2 ) THEN 893 893 IF( iom_use('uice_cvg') ) THEN 894 894 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.