Changeset 11867
- Timestamp:
- 2019-11-06T13:08:50+01:00 (5 years ago)
- Location:
- NEMO/branches/2019/dev_r11842_SI3-10_EAP
- Files:
-
- 1 added
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11842_SI3-10_EAP/cfgs/SHARED/namelist_ice_ref
r11731 r11867 97 97 rn_relast = 0.333 ! ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast 98 98 ! advised value: 1/3 (rn_nevp=120) or 1/9 (rn_nevp=300) 99 ln_rhg_EAP = .false. ! EVP rheology 99 100 / 100 101 !------------------------------------------------------------------------------ -
NEMO/branches/2019/dev_r11842_SI3-10_EAP/src/ICE/icedyn_rhg.F90
r11536 r11867 17 17 USE ice ! sea-ice: variables 18 18 USE icedyn_rhg_evp ! sea-ice: EVP rheology 19 USE icedyn_rhg_eap ! sea-ice: EAP rheology 19 20 USE icectl ! sea-ice: control prints 20 21 ! … … 33 34 ! ! associated indices: 34 35 INTEGER, PARAMETER :: np_rhgEVP = 1 ! EVP rheology 35 !!INTEGER, PARAMETER :: np_rhgEAP = 2 ! EAP rheology36 INTEGER, PARAMETER :: np_rhgEAP = 2 ! EAP rheology 36 37 37 38 ! ** namelist (namrhg) ** 38 39 LOGICAL :: ln_rhg_EVP ! EVP rheology 40 LOGICAL :: ln_rhg_EAP ! EVP rheology 39 41 ! 40 42 !! * Substitutions … … 81 83 CALL ice_dyn_rhg_evp( kt, stress1_i, stress2_i, stress12_i, shear_i, divu_i, delta_i ) 82 84 ! 85 ! !----------------------------! 86 CASE( np_rhgEAP ) ! Elasto-Anisotropic-Plastic ! 87 ! !----------------------------! 88 CALL ice_dyn_rhg_eap( kt, stress1_i, stress2_i, stress12_i, shear_i, divu_i, delta_i ) 83 89 END SELECT 84 90 ! 85 IF( lrst_ice ) THEN !* write EVP fields in the restart file 86 IF( ln_rhg_EVP ) CALL rhg_evp_rst( 'WRITE', kt ) 91 IF( lrst_ice ) THEN 92 IF( ln_rhg_EVP ) CALL rhg_evp_rst( 'WRITE', kt ) !* write EVP fields in the restart file 93 IF( ln_rhg_EAP ) CALL rhg_eap_rst( 'WRITE', kt ) !* write EAP fields in the restart file 87 94 ENDIF 88 95 ! … … 110 117 INTEGER :: ios, ioptio ! Local integer output status for namelist read 111 118 !! 112 NAMELIST/namdyn_rhg/ ln_rhg_EVP, ln_aEVP, rn_creepl, rn_ecc , nn_nevp, rn_relast 119 NAMELIST/namdyn_rhg/ ln_rhg_EVP, ln_aEVP, rn_creepl, rn_ecc , nn_nevp, rn_relast, ln_rhg_EAP 113 120 !!------------------------------------------------------------------- 114 121 ! … … 132 139 WRITE(numout,*) ' number of iterations for subcycling nn_nevp = ', nn_nevp 133 140 WRITE(numout,*) ' ratio of elastic timescale over ice time step rn_relast = ', rn_relast 141 WRITE(numout,*) ' rheology EAP (icedyn_rhg_eap) ln_rhg_EAP = ', ln_rhg_EAP 134 142 ENDIF 135 143 ! … … 137 145 ioptio = 0 138 146 IF( ln_rhg_EVP ) THEN ; ioptio = ioptio + 1 ; nice_rhg = np_rhgEVP ; ENDIF 139 !!IF( ln_rhg_EAP ) THEN ; ioptio = ioptio + 1 ; nice_rhg = np_rhgEAP ; ENDIF147 IF( ln_rhg_EAP ) THEN ; ioptio = ioptio + 1 ; nice_rhg = np_rhgEAP ; ENDIF 140 148 IF( ioptio /= 1 ) CALL ctl_stop( 'ice_dyn_rhg_init: choose one and only one ice rheology' ) 141 149 ! 142 150 IF( ln_rhg_EVP ) CALL rhg_evp_rst( 'READ' ) !* read or initialize all required files 151 IF( ln_rhg_EAP ) CALL rhg_eap_rst( 'READ' ) !* read or initialize all required files 143 152 ! 144 153 END SUBROUTINE ice_dyn_rhg_init
Note: See TracChangeset
for help on using the changeset viewer.