Changeset 4019 for LMDZ6

Timestamp:

Nov 25, 2021, 2:59:52 PM (3 years ago)

Author:

evignon

Message:

Correction (esperee) du probleme des nan dans le calcul de l_mix dans yamada4

File:

• LMDZ6/trunk/libf/phylmd/yamada4.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/yamada4.F90

@@ -986 +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
@@ -1101 +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
@@ -1167 +1168 @@
         ENDDO
      ENDDO
-   ENDIF                                                                        ! pbl_lmixmin_alpha
+   ENDIF                                                                           ! pbl_lmixmin_alpha
 
 !==================================================================================
@@ -1183 +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