- Timestamp:
- Oct 5, 2023, 8:04:29 AM (15 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/lmdz_atke_exchange_coeff.F90
r4688 r4714 20 20 ! 21 21 ! Main assumptions of the model : 22 ! (1) dry atmosphere 23 ! (2) horizontal homogeneity (Dx=Dy=0.) 22 ! (1) horizontal homogeneity (Dx=Dy=0.) 24 23 !======================================================================= 25 24 … … 205 204 ENDDO 206 205 206 ELSE IF (iflag_atke_lmix .EQ. 3 ) THEN 207 ! add effect of wind shear on lstrat following grisogono 2010, qjrms 208 DO ilay=2,nlay 209 DO igrid=1,ngrid 210 l_exchange(igrid,ilay) = kappa*l0*z_interf(igrid,ilay) / (kappa*z_interf(igrid,ilay) + l0) 211 IF (N2(igrid,ilay) .GT. 0. .AND. shear2(igrid,ilay) .GT. 0.) THEN 212 lstrat=clmix*sqrt(tke(igrid,ilay))/sqrt(shear2(igrid,ilay))*(1.0+Ri(igrid,ilay)/(2.*Prandtl(igrid,ilay))) 213 lstrat=max(lstrat,lmin) 214 !Inverse interpolation, Van de Wiel et al. 2010 215 l_exchange(igrid,ilay)=(1./(l_exchange(igrid,ilay))+1./(lstrat))**(-1.0) 216 ENDIF 217 ENDDO 218 ENDDO 219 207 220 208 221
Note: See TracChangeset
for help on using the changeset viewer.