- Timestamp:
- 2017-07-04T17:46:48+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/wrk_OMP_test_for_Silvia/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfric.F90
r8056 r8279 164 164 ! !== avm and avt = F(Richardson number) ==! 165 165 DO jk = 2, jpkm1 166 DO jj = k_Jstr, k_Jend167 DO ji = k_Jstr, k_Iend! coefficient = F(richardson number) (avm-weighted Ri)166 DO jj = tnldj, tnlej 167 DO ji = tnldi, tnlei ! coefficient = F(richardson number) (avm-weighted Ri) 168 168 zcfRi = 1._wp / ( 1._wp + rn_alp * MAX( 0._wp , p_avm(ji,jj,jk) * rn2(ji,jj,jk) / ( p_sh2(ji,jj,jk) + 1.e-20 ) ) ) 169 169 zav = rn_avmri * zcfRi**nn_ric … … 180 180 IF( ln_mldw ) THEN !== set a minimum value in the Ekman layer ==! 181 181 ! 182 DO jj = k_Jstr, k_Jend!* Ekman depth183 DO ji = k_Jstr, k_Iend182 DO jj = tnldj, tnlej !* Ekman depth 183 DO ji = tnldi, tnlei 184 184 zustar = SQRT( taum(ji,jj) * r1_rau0 ) 185 185 zed = rn_ekmfc * zustar / ( ABS( ff_t(ji,jj) ) + rsmall ) ! Ekman depth … … 188 188 END DO 189 189 DO jk = 2, jpkm1 !* minimum mixing coeff. within the Ekman layer 190 DO jj = k_Jstr, k_Jend191 DO ji = k_Jstr, k_Iend190 DO jj = tnldj, tnlej 191 DO ji = tnldi, tnlei 192 192 IF( pdept(ji,jj,jk) < zh_ekm(ji,jj) ) THEN 193 193 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.