Changeset 4714

Timestamp:

Oct 5, 2023, 8:04:29 AM (15 months ago)

Author:

evignon

Message:

ajout nouvelle option du calcul de la longueur de melange en fonction du cisaillement de vent
dans atke

File:

• LMDZ6/trunk/libf/phylmd/lmdz_atke_exchange_coeff.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_atke_exchange_coeff.F90

@@ -20 +20 @@
 !
 ! Main assumptions of the model :
-! (1) dry atmosphere
-! (2) horizontal homogeneity (Dx=Dy=0.)
+! (1) horizontal homogeneity (Dx=Dy=0.)
 !=======================================================================
 
@@ -205 +204 @@
    ENDDO
 
+ELSE IF (iflag_atke_lmix .EQ. 3 ) THEN
+! add effect of wind shear on lstrat following grisogono 2010, qjrms
+DO ilay=2,nlay
+      DO igrid=1,ngrid
+          l_exchange(igrid,ilay) = kappa*l0*z_interf(igrid,ilay) / (kappa*z_interf(igrid,ilay) + l0)
+          IF (N2(igrid,ilay) .GT. 0. .AND. shear2(igrid,ilay) .GT. 0.) THEN
+             lstrat=clmix*sqrt(tke(igrid,ilay))/sqrt(shear2(igrid,ilay))*(1.0+Ri(igrid,ilay)/(2.*Prandtl(igrid,ilay)))
+             lstrat=max(lstrat,lmin)
+             !Inverse interpolation, Van de Wiel et al. 2010   
+             l_exchange(igrid,ilay)=(1./(l_exchange(igrid,ilay))+1./(lstrat))**(-1.0)
+          ENDIF
+      ENDDO
+   ENDDO
+