Changeset 5469

Timestamp:

Jan 10, 2025, 1:19:15 PM (9 hours ago)

Author:

jyg

Message:

Bug fix in subroutine cv3_tracer in cv3_routines.f90.
This is an old bug (>10y).

File:

• LMDZ6/trunk/libf/phylmd/cv3_routines.f90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/cv3_routines.f90

@@ -5022 +5022 @@
         d1a(i, j) = d1a(i, j) + ment(i, k, j)*ep(i, k)*(1.-sigij(i,k,j))
         IF (k<=j) THEN
-          dam(i, j) = dam(i, j) + ment(i, k, j)*epm(i, k, j)*(1.-ep(i,k))*(1.-sigij(i,k,j))
+!!jyg          dam(i, j) = dam(i, j) + ment(i, k, j)*epm(i, k, j)*(1.-ep(i,k))*(1.-sigij(i,k,j))
+          dam(i, j) = dam(i, j) + ment(i, k, j)*epm(i, j, k)*(1.-ep(i,k))*(1.-sigij(i,k,j))
           phi2(i, j, k) = phi(i, j, k)*epm(i, j, k)
         END IF