Changeset 4171 for LMDZ6/branches/LMDZ-ECRAD/libf/phylmd/yamada4.F90

Timestamp:

Jun 17, 2022, 4:24:49 PM (2 years ago)

Author:

lguez

Message:

Sync latest trunk changes to branch LMDZ-ECRAD.

Location:

LMDZ6/branches/LMDZ-ECRAD

Files:

• . (modified) (1 prop)
• libf/phylmd/yamada4.F90 (modified) (7 diffs)

LMDZ6/branches/LMDZ-ECRAD/libf/phylmd/yamada4.F90

@@ -6 +6 @@
   USE dimphy
   USE ioipsl_getin_p_mod, ONLY : getin_p
-  USE phys_local_var_mod, only: tke_dissip
+  USE phys_local_var_mod, only: tke_dissip,wprime
   
   IMPLICIT NONE
@@ -726 +726 @@
       lyam(1:ngrid, 2:klev)*5.17*kn(1:ngrid, 2:klev)*n2(1:ngrid, 2:klev)/ &
       sqrt(q2(1:ngrid,2:klev))
-
+ 
     t2yam(1:ngrid, 2:klev) = 9.1*kn(1:ngrid, 2:klev)* &
       dtetadz(1:ngrid, 2:klev)**2/sqrt(q2(1:ngrid,2:klev))* &
@@ -750 +750 @@
 
 !============================================================================
-! Diagnostique de la dissipation
+! Diagnostique de la dissipation et vitesse verticale
 !============================================================================
 
 ! Diagnostics
  tke_dissip(1:ngrid,:,nsrf)=0.
-! DO k=2,klev
-!    DO ig=1,ngrid
-!       jg=ni(ig)
-!       tke_dissip(jg,k,nsrf)=dissip(ig,k)
-!    ENDDO
-! ENDDO
+ wprime(1:ngrid,:,nsrf)=0.
+ DO k=2,klev
+    DO ig=1,ngrid
+       jg=ni(ig)
+       wprime(jg,k,nsrf)=sqrt(MAX(1./3*q2(ig,k),0.))
+       tke_dissip(jg,k,nsrf)=dissip(ig,k)
+    ENDDO
+ ENDDO
  
 !=============================================================================
@@ -984 +986 @@
  INTEGER            ngrid              ! Nombre de points concern??s sur l'horizontal
  INTEGER            iflag_pbl          ! Choix du sch??ma de turbulence
- REAL            pbl_lmixmin_alpha  ! on active ou non le calcul de la longueur de melange minimum
+ REAL               pbl_lmixmin_alpha  ! on active ou non le calcul de la longueur de melange minimum en fonction du relief
  REAL               lmixmin            ! Minimum absolu de la longueur de m??lange
  REAL               zlay(klon, klev)   ! altitude du centre de la couche
@@ -1099 +1101 @@
   END IF
 
-!=================================================================================
-!  CALCUL d'une longueur de melange en fonctions de la topographie sous maille: l2
-! si plb_lmixmin_alpha=TRUE et si on se trouve sur de la terre ( pas actif sur les 
+!===========================================================================================
+!  CALCUL d'une longueur de melange minimum en fonctions de la topographie sous maille: l2
+! si pbl_lmixmin_alpha=TRUE et si on se trouve sur de la terre ( pas actif sur les 
 ! glacier, la glace de mer et les oc??ans)
-!=================================================================================
+!===========================================================================================
 
    l2(1:ngrid,:)=0.0
    l_mixmin(1:ngrid,:,nsrf)=0.
    l_mix(1:ngrid,:,nsrf)=0.
+   hlim(1:ngrid)=0.
 
    IF (nsrf .EQ. 1) THEN
@@ -1165 +1168 @@
         ENDDO
      ENDDO
-   ENDIF                                                                        ! pbl_lmixmin_alpha
+   ENDIF                                                                           ! pbl_lmixmin_alpha
 
 !==================================================================================
@@ -1181 +1184 @@
 ! Diagnostics
 
- DO k=2,klev
+ DO k=1,klev+1
     DO ig=1,ngrid
        jg=ni(ig)
        l_mix(jg,k,nsrf)=lmix(ig,k)
-       l_mixmin(jg,k,nsrf)=l2(ig,k)
+       l_mixmin(jg,k,nsrf)=MAX(l2(ig,k),lmixmin)
     ENDDO
  ENDDO
- DO ig=1,ngrid
-    jg=ni(ig)
-    l_mix(jg,1,nsrf)=hlim(ig)
- ENDDO