- Timestamp:
- 2017-11-24T17:56:51+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/NEMO/LIM_SRC_3/icedyn_rhg.F90
r8586 r8813 78 78 CASE( np_rhgEVP ) ! Elasto-Viscous-Plastic ! 79 79 ! !------------------------! 80 CALL ice_dyn_rhg_evp( kt, stress1_i, stress2_i, stress12_i, u_ice, v_ice,shear_i, divu_i, delta_i )80 CALL ice_dyn_rhg_evp( kt, stress1_i, stress2_i, stress12_i, shear_i, divu_i, delta_i ) 81 81 ! 82 82 END SELECT … … 107 107 INTEGER :: ios, ioptio ! Local integer output status for namelist read 108 108 !! 109 NAMELIST/namdyn_rhg/ ln_rhg_EVP, rn_creepl, rn_ecc , nn_nevp, rn_relast109 NAMELIST/namdyn_rhg/ ln_rhg_EVP, ln_aEVP, rn_creepl, rn_ecc , nn_nevp, rn_relast 110 110 !!------------------------------------------------------------------- 111 111 ! … … 125 125 WRITE(numout,*) ' Namelist namdyn_rhg:' 126 126 WRITE(numout,*) ' rheology EVP (icedyn_rhg_evp) ln_rhg_EVP = ', ln_rhg_EVP 127 WRITE(numout,*) ' use adaptive EVP (aEVP) ln_aEVP = ', ln_aEVP 127 128 WRITE(numout,*) ' creep limit rn_creepl = ', rn_creepl 128 129 WRITE(numout,*) ' eccentricity of the elliptical yield curve rn_ecc = ', rn_ecc … … 137 138 IF( ioptio /= 1 ) CALL ctl_stop( 'ice_dyn_rhg_init: choose one and only one ice rheology' ) 138 139 ! 139 IF( ln_rhg_EVP ) CALL rhg_evp_rst( 'READ' ) !* read or initialize all required files140 IF( ln_rhg_EVP ) CALL rhg_evp_rst( 'READ' ) !* read or initialize all required files 140 141 ! 141 142 END SUBROUTINE ice_dyn_rhg_init
Note: See TracChangeset
for help on using the changeset viewer.