[MITgcm-devel] ini_cg3d

[Jean-Michel Campin]

Hi,

I have a question about  the preconditioner in "ini_cg3d.F":
I don't see clearly why
aU is set to aV3d(I,J,K-1,bi,bj) the 2nd time
lines 291,295:
           IF ( K .NE. 1 ) THEN
            aU = aV3d(I,J,K-1,bi,bj)
           ELSE
            aU = 0
           ENDIF
whereas it is set to aV3d(I,J,K,bi,bj) the 1rst time (which seems
more logical to me)
lines 234,238:
           IF ( K .NE. 1 ) THEN
            aU = aV3d(I,J,K,bi,bj)
           ELSE
            aU = 0
           ENDIF

It might not make a big difference anyway, since the 2nd time
is just to reset the preconditioner to unity if aC is zero.

Jean-Michel