Changeset 4632 for LMDZ6/trunk/libf/phylmd/atke_exchange_coeff_mod.F90
- Timestamp:
- Jul 18, 2023, 3:27:56 PM (15 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/atke_exchange_coeff_mod.F90
r4631 r4632 176 176 ENDDO 177 177 178 ELSE IF (iflag_atke_lmix .EQ. 2 ) THEN 179 ! add effect of wind shear on lstrat following grisogono and belusic 2008, qjrms, eq 2 180 DO ilay=2,nlay 181 DO igrid=1,ngrid 182 l_exchange(igrid,ilay) = kappa*l0*z_interf(igrid,ilay) / (kappa*z_interf(igrid,ilay) + l0) 183 IF (N2(igrid,ilay) .GT. 0.) THEN 184 lstrat=clmix*sqrt(tke(igrid,ilay))/(2.*max(sqrt(shear2(igrid,ilay)),1E-10)*(1.+sqrt(Ri(igrid,ilay))/2.)) 185 IF (lstrat .LT. l_exchange(igrid,ilay)) THEN 186 l_exchange(igrid,ilay)=max(lstrat,lmin) 187 ENDIF 188 ENDIF 189 ENDDO 190 ENDDO 191 192 193 178 194 ELSE 179 ! default: neglect effect of local stratification 195 ! default: neglect effect of local stratification and shear 180 196 181 197 DO ilay=2,nlay+1 … … 217 233 218 234 ! full implicit scheme resolved with a second order polynomial equation 219 ! + exact resolution for very stable cases 235 ! + exact resolution for very stable cases (iflag_atke_lmix=1) 220 236 DO ilay=2,nlay 221 237 DO igrid=1,ngrid
Note: See TracChangeset
for help on using the changeset viewer.