Changeset 2605 for LMDZ5/trunk/libf

Timestamp:

Jul 28, 2016, 11:46:18 AM (8 years ago)

Author:

jyg

Message:

Adding some initialization in cv3_routines.F90.
Should help guarantee 1+1=2.

File:

• LMDZ5/trunk/libf/phylmd/cv3_routines.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/cv3_routines.F90

@@ -238 +238 @@
 ! gz = phi at the full levels (same as p).
 
-  DO i = 1, len
-    gz(i, 1) = 0.0
-  END DO
+!!  DO i = 1, len                    !jyg
+!!    gz(i, 1) = 0.0                 !jyg
+!!  END DO                           !jyg
+    gz(:,:) = 0.                     !jyg: initialization of the whole array
 ! ori      do 140 k=2,nlp
   DO k = 2, nl ! convect3
@@ -409 +410 @@
     IF (ok_new_feed) THEN
       IF (iter==niter) THEN
-        DO k = minorig, nl
+        DO i = 1,len                         !jyg
+          pfeedmin(i) = ph(i,minorig+1)      !jyg
+        ENDDO                                !jyg
+        DO k = minorig+1, nl                 !jyg
+!!        DO k = minorig, nl                 !jyg
           DO i = 1, len
             IF (ph(i,k)>=plclfeed(i)) pfeedmin(i) = ph(i, k)