Changeset 911

Timestamp:

06/17/19 00:31:27 (5 years ago)

Author:

dubos

Message:

devel : fix bug in tridiagonal solver for vertical diffusion in Venus test case

File:

• codes/icosagcm/devel/src/initial/etat0_venus.f90 (modified) (2 diffs)

codes/icosagcm/devel/src/initial/etat0_venus.f90

@@ -168 +168 @@
     DO ij=1,ngrid
        X_ij = 1./B(ij,1)
-       C(ij,2) = -A(ij,2) * X_ij
+       C(ij,1)  = -A(ij,2) * X_ij
        xu(ij,1) = Ru(ij,1) * X_ij
        xv(ij,1) = Rv(ij,1) * X_ij
@@ -175 +175 @@
        !DIR$ SIMD
        DO ij=1,ngrid
-          X_ij = 1./( B(ij,l) + A(ij,l)*C(ij,l) )
-          C(ij,l+1) = -A(ij,l+1) * X_ij ! zero for l=llm
+          X_ij = 1./( B(ij,l) + A(ij,l)*C(ij,l-1) )
+          C(ij,l)  = -A(ij,l+1) * X_ij ! zero for l=llm
           xu(ij,l) = (Ru(ij,l)+A(ij,l)*xu(ij,l-1)) * X_ij
           xv(ij,l) = (Rv(ij,l)+A(ij,l)*xv(ij,l-1)) * X_ij