Changeset 3577 for LMDZ6

Timestamp:

Oct 7, 2019, 3:30:19 PM (5 years ago)

Author:

fhourdin

Message:

Corrections d'erreurs identifiees en debug ancienne physique sur Jean-Zay
Vive le poihl interactif.

Location:

LMDZ6/trunk/libf/phylmd

Files:

• coef_diff_turb_mod.F90 (modified) (1 diff)
• cv30_routines.F90 (modified) (6 diffs)
• physiq_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/coef_diff_turb_mod.F90

@@ -182 +182 @@
        ycoefh(1:knon,2:klev)=ykmn(1:knon,2:klev)
                 
+    ELSE
+       ! No TKE for Standard Physics
+       yq2=0.
     ENDIF !(iflag_pbl.ge.3)
 

LMDZ6/trunk/libf/phylmd/cv30_routines.F90

@@ -3050 +3050 @@
   ! variables pour tracer dans precip de l'AA et des mel
   ! local variables:
-  INTEGER i, j, k
+  INTEGER i, j, k, nam1
   REAL epm(nloc, na, na)
 
+  nam1=na-1 ! Introduced because ep is not defined for j=na
   ! variables d'Emanuel : du second indice au troisieme
   ! --->    tab(i,k,j) -> de l origine k a l arrivee j
@@ -3082 +3083 @@
   ! fraction deau condensee dans les melanges convertie en precip : epm
   ! et eau condensée précipitée dans masse d'air saturé : l_m*dM_m/dzdz.dzdz
-  DO j = 1, na
+  DO j = 1, nam1
     DO k = 1, j - 1
       DO i = 1, ncum
@@ -3095 +3096 @@
   END DO
 
-  DO j = 1, na
-    DO k = 1, na
+  DO j = 1, nam1
+    DO k = 1, nam1
       DO i = 1, ncum
         IF (k>=icb(i) .AND. k<=inb(i)) THEN
@@ -3106 +3107 @@
   END DO
 
-  DO j = 1, na
+  DO j = 1, nam1
     DO k = 1, j - 1
       DO i = 1, ncum
@@ -3117 +3118 @@
 
   ! matrices pour calculer la tendance des concentrations dans cvltr.F90
-  DO j = 1, na
-    DO k = 1, na
+  DO j = 1, nam1
+    DO k = 1, nam1
       DO i = 1, ncum
         da(i, j) = da(i, j) + (1.-sij(i,k,j))*ment(i, k, j)
@@ -3127 +3128 @@
   END DO
 
-  DO j = 1, na
+  DO j = 1, nam1
     DO k = 1, j - 1
       DO i = 1, ncum

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -1701 +1701 @@
 
 !jyg<
-       IF (klon_glo==1) THEN
+       IF (iflag_pbl<=1) THEN
+          ! No TKE for Standard Physics
+          pbl_tke(:,:,:)=0.
+
+       ELSE IF (klon_glo==1) THEN
           pbl_tke(:,:,is_ave) = 0.
           DO nsrf=1,nbsrf