Changeset 2502

Timestamp:

May 4, 2016, 11:51:01 AM (8 years ago)

Author:

jyg

Message:

Bug fix in cv3p1_closure.F90 and
cv3p2_closure.F90: initialization loops extended
up to nd.

Location:

LMDZ5/trunk/libf/phylmd

Files:

• cv3p1_closure.F90 (modified) (3 diffs)
• cv3p2_closure.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/cv3p1_closure.F90

@@ -113 +113 @@
   IF (prt_level>=20) PRINT *, 'cv3p1_param nloc ncum nd icb inb nl', nloc, &
     ncum, nd, icb(nloc), inb(nloc), nl
-  DO k = 1, nl
+  DO k = 1, nd          !jyg: initialization up to nd
     DO il = 1, ncum
       m(il, k) = 0.0
@@ -120 +120 @@
 
 !CR: initializations for erosion of adiabatic ascent 
-  DO k = 1,nl
+  DO k = 1,nd           !jyg: initialization up to nd
     DO il = 1, ncum
         mad(il,k)=0.
@@ -628 +628 @@
 
 !Verification sum(me)=sum(m)
-  DO k = 1,nl+1
+  DO k = 1,nd                         !jyg: initialization up to nd
     DO il = 1, ncum
        md(il,k)=0.

LMDZ5/trunk/libf/phylmd/cv3p2_closure.F90

@@ -126 +126 @@
   DO k = 1, nl
     DO il = 1, ncum
-      m(il, k) = 0.0
       rhodp(il,k) = 0.007*p(il, k)*(ph(il,k)-ph(il,k+1))/tv(il, k)
     END DO
   END DO
 
-!CR: initializations for erosion of adiabatic ascent 
-  DO k = 1,nl
+!CR+jyg: initializations (up to nd) for erosion of adiabatic ascent and of m and wlim
+  DO k = 1,nd
     DO il = 1, ncum
         mad(il,k)=0.
@@ -138 +137 @@
         betalim(il,k)=1.
         wlim(il,k)=0.
+        m(il, k) = 0.0
     ENDDO
   ENDDO
@@ -711 +711 @@
 
 !Verification sum(me)=sum(m)
-  DO k = 1,nl+1
+  DO k = 1,nd
     DO il = 1, ncum
        md(il,k)=0.