- Timestamp:
- 2021-03-09T09:37:00+01:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2021/dev_r14393_HPC-03_Mele_Comm_Cleanup/src/OCE/ZDF/zdfric.F90
r14072 r14601 156 156 ! 157 157 ! !== avm and avt = F(Richardson number) ==! 158 DO_3D( 1, 0, 1, 0, 2, jpkm1 ) ! coefficient = F(richardson number) (avm-weighted Ri) 158 ! [comm_cleanup] ! DO_3D( 1, 0, 1, 0, 2, jpkm1 ) ! coefficient = F(richardson number) (avm-weighted Ri) 159 DO_3D( nn_hls, nn_hls-1, nn_hls, nn_hls-1, 2, jpkm1 ) ! coefficient = F(richardson number) (avm-weighted Ri) 159 160 zcfRi = 1._wp / ( 1._wp + rn_alp * MAX( 0._wp , avm(ji,jj,jk) * rn2(ji,jj,jk) / ( p_sh2(ji,jj,jk) + 1.e-20 ) ) ) 160 161 zav = rn_avmri * zcfRi**nn_ric … … 169 170 IF( ln_mldw ) THEN !== set a minimum value in the Ekman layer ==! 170 171 ! 171 DO_2D( 0, 0, 0, 0 ) !* Ekman depth 172 ! [comm_cleanup] ! DO_2D( 0, 0, 0, 0 ) !* Ekman depth 173 DO_2D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1 ) 172 174 zustar = SQRT( taum(ji,jj) * r1_rho0 ) 173 175 zhek = rn_ekmfc * zustar / ( ABS( ff_t(ji,jj) ) + rsmall ) ! Ekman depth 174 176 zh_ekm(ji,jj) = MAX( rn_mldmin , MIN( zhek , rn_mldmax ) ) ! set allowed range 175 177 END_2D 176 DO_3D( 0, 0, 0, 0, 2, jpkm1 ) !* minimum mixing coeff. within the Ekman layer 178 ! [comm_cleanup] ! DO_3D( 0, 0, 0, 0, 2, jpkm1 ) !* minimum mixing coeff. within the Ekman layer 179 DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 2, jpkm1 ) !* minimum mixing coeff. within the Ekman layer 177 180 IF( gdept(ji,jj,jk,Kmm) < zh_ekm(ji,jj) ) THEN 178 181 p_avm(ji,jj,jk) = MAX( p_avm(ji,jj,jk), rn_wvmix ) * wmask(ji,jj,jk)
Note: See TracChangeset
for help on using the changeset viewer.