Changeset 3355 for LMDZ6/branches/DYNAMICO-conv/libf/phylmd/cv3_routines.F90

Timestamp:

Jun 29, 2018, 11:12:44 AM (6 years ago)

Author:

Laurent Fairhead

Message:

Commiting variables initialisation as found by Yann
LF

File:

• LMDZ6/branches/DYNAMICO-conv/libf/phylmd/cv3_routines.F90 (modified) (1 diff)

LMDZ6/branches/DYNAMICO-conv/libf/phylmd/cv3_routines.F90

@@ -3747 +3747 @@
     END DO ! cld
 
+!ym BIG Warning : it seems that the k loop is missing !!!
+!ym Strong advice to check this
+!ym add a k loop temporary 
+
 ! (particular case: no detraining level is found)                              ! cld
-    DO il = 1, ncum                                                            ! cld
-      IF (i<=inb(il) .AND. nent(il,i)==0 .AND. iflag(il)<=1) THEN              ! cld
-        qcond(il, i) = qcond(il, i) + (1.-ep(il,i))*clw(il, i)                 ! cld
-        qtment(il, i) = qent(il,k,i) + qtment(il,i)                          ! cld
-        nqcond(il, i) = nqcond(il, i) + 1.                                     ! cld
-      END IF                                                                   ! cld
-    END DO                                                                     ! cld
+    DO k = i + 1, nl
+      DO il = 1, ncum        !ym k loop added                                                      ! cld
+        IF (i<=inb(il) .AND. nent(il,i)==0 .AND. iflag(il)<=1) THEN              ! cld
+          qcond(il, i) = qcond(il, i) + (1.-ep(il,i))*clw(il, i)                 ! cld
+!jyg<   Bug correction 20180620
+!      PROBLEM: Should not qent(il,i,i) be taken into account even if nent(il,i)/=0?
+!!        qtment(il, i) = qent(il,k,i) + qtment(il,i)
+! cld
+          qtment(il, i) = qent(il,i,i) + qtment(il,i)
+! cld
+!>jyg
+
+          nqcond(il, i) = nqcond(il, i) + 1.                                     ! cld
+        END IF                                                                   ! cld
+      END DO
+    ENDDO                                                                     ! cld
 
     DO il = 1, ncum                                                            ! cld