Changeset 15259 for NEMO/branches/UKMO/NEMO_4.0.4_EAP_rheology/src/ICE
- Timestamp:
- 2021-09-14T16:41:21+02:00 (3 years ago)
- Location:
- NEMO/branches/UKMO/NEMO_4.0.4_EAP_rheology/src/ICE
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0.4_EAP_rheology/src/ICE/icedyn_rdgrft.F90
r14998 r15259 187 187 ! closing_net = rate at which open water area is removed + ice area removed by ridging 188 188 ! - ice area added in new ridges 189 closing_net(ji) = rn_csrdg * 0.5_wp * ( zdelt(ji) - ABS( zdivu(ji) ) ) - MIN( zdivu(ji), 0._wp )190 189 IF( ln_rhg_EVP ) closing_net(ji) = rn_csrdg * 0.5_wp * ( zdelt(ji) - ABS( zdivu(ji) ) ) - MIN( zdivu(ji), 0._wp ) 191 190 IF( ln_rhg_EAP ) closing_net(ji) = zconv(ji) … … 773 772 ! !--------------------------------------------------! 774 773 strength(:,:) = rn_pstar * SUM( v_i(:,:,:), dim=3 ) * EXP( -rn_crhg * ( 1._wp - SUM( a_i(:,:,:), dim=3 ) ) ) 775 ! ismooth = 1 ! original code 776 ismooth = 0 ! try for EAP stability 774 ismooth = 1 ! original code 775 ! ismooth = 0 ! try for EAP stability 776 ! 777 777 ! !--------------------------------------------------! 778 778 ELSE ! Zero strength ! -
NEMO/branches/UKMO/NEMO_4.0.4_EAP_rheology/src/ICE/icedyn_rhg_eap.F90
r14998 r15259 799 799 800 800 ! convergence test 801 IF( nn_rhg_chkcvg == 2 ) CALL rhg_cvg ( kt, jter, nn_nevp, u_ice, v_ice, zu_ice, zv_ice )801 IF( nn_rhg_chkcvg == 2 ) CALL rhg_cvg_eap( kt, jter, nn_nevp, u_ice, v_ice, zu_ice, zv_ice ) 802 802 ! 803 803 ! ! ==================== ! … … 1044 1044 1045 1045 1046 SUBROUTINE rhg_cvg ( kt, kiter, kitermax, pu, pv, pub, pvb )1046 SUBROUTINE rhg_cvg_eap( kt, kiter, kitermax, pu, pv, pub, pvb ) 1047 1047 !!---------------------------------------------------------------------- 1048 !! *** ROUTINE rhg_cvg ***1048 !! *** ROUTINE rhg_cvg_eap *** 1049 1049 !! 1050 1050 !! ** Purpose : check convergence of oce rheology … … 1072 1072 IF( lwp ) THEN 1073 1073 WRITE(numout,*) 1074 WRITE(numout,*) 'rhg_cvg : ice rheology convergence control'1074 WRITE(numout,*) 'rhg_cvg_eap : ice rheology convergence control' 1075 1075 WRITE(numout,*) '~~~~~~~' 1076 1076 ENDIF … … 1111 1111 ENDIF 1112 1112 1113 END SUBROUTINE rhg_cvg 1113 END SUBROUTINE rhg_cvg_eap 1114 1114 1115 1115 -
NEMO/branches/UKMO/NEMO_4.0.4_EAP_rheology/src/ICE/iceupdate.F90
r14075 r15259 146 146 ! 147 147 ! the non-solar is simply derived from the solar flux 148 qns(ji,jj) = qt_oce_ai(ji,jj) - zqsr148 qns(ji,jj) = qt_oce_ai(ji,jj) - qsr(ji,jj) 149 149 150 150 ! Mass flux at the atm. surface
Note: See TracChangeset
for help on using the changeset viewer.