Changeset 111 for trunk/NEMO/OPA_SRC/SOL/solmat.F90

Timestamp:

2004-06-28T15:22:55+02:00 (20 years ago)

Author:

opalod

Message:

CT : UPDATE070 : Optimisation of the Red-Black SOR algorithm convergence

File:

• trunk/NEMO/OPA_SRC/SOL/solmat.F90 (modified) (8 diffs)

trunk/NEMO/OPA_SRC/SOL/solmat.F90

@@ -222 +222 @@
          DO jj = 1, jpj
             DO ji = 1, jpi
-               IF( bmask(ji,jj) /= 0. )   gcdprc(ji,jj) = 1.e0 / gcdmat(ji,jj)
+               IF( bmask(ji,jj) /= 0.e0 )   gcdprc(ji,jj) = 1.e0 / gcdmat(ji,jj)
             END DO
          END DO
@@ -240 +240 @@
          DO jj = 1, jpj
             DO ji = 1, jpi
-               IF( bmask(ji,jj) /= 0. ) THEN
-                  gcdprc(ji,jj) = 1.
+               IF( bmask(ji,jj) /= 0.e0 ) THEN
+                  gcdprc(ji,jj) = 1.e0
                ELSE
-                  gcdprc(ji,jj) = 0.
+                  gcdprc(ji,jj) = 0.e0
                ENDIF
             END DO
@@ -531 +531 @@
             DO jj = 1, jpj
                IF( iimask(1,jj) /= 0 ) THEN
-                  gcp(1,jj,2) = 0.
+                  gcp(1,jj,2) = 0.e0
                   gcp(1,jj,1) = zdemi * gcp(1,jj,1)
                   gcp(1,jj,4) = zdemi * gcp(1,jj,4)
@@ -546 +546 @@
             DO jj = 1, jpj
                IF( iimask(iiend,jj) /= 0 ) THEN
-                  gcp(iiend,jj,3) = 0.
+                  gcp(iiend,jj,3) = 0.e0
                   gcp(iiend,jj,1) = zdemi * gcp(iiend,jj,1)
                   gcp(iiend,jj,4) = zdemi * gcp(iiend,jj,4)
@@ -564 +564 @@
          DO ji = 1, jpi
             IF( iimask(ji,1) /= 0 ) THEN
-               gcp(ji,1,1) = 0.
+               gcp(ji,1,1) = 0.e0
                gcp(ji,1,2) = zdemi * gcp(ji,1,2)
                gcp(ji,1,3) = zdemi * gcp(ji,1,3)
@@ -580 +580 @@
          DO ji = 1, jpi
             IF( iimask(ji,ijend) /= 0 ) THEN
-               gcp(ji,ijend,4) = 0.
+               gcp(ji,ijend,4) = 0.e0
                gcp(ji,ijend,2) = zdemi * gcp(ji,ijend,2) 
                gcp(ji,ijend,3) = zdemi * gcp(ji,ijend,3)
@@ -682 +682 @@
       CALL feti_vsub(noeuds,wfeti(may),wfeti(maz),wfeti(maz))
 
-      zres2 = 0.
+      zres2 = 0.e0
       DO jl = 1, noeuds
          zres2 = zres2 + wfeti(may+jl-1) * wfeti(may+jl-1)
@@ -688 +688 @@
       CALL mpp_sum(zres2,1,zres)
 
-      res2 = 0.
+      res2 = 0.e0
       DO jl = 1, noeuds
          res2 = res2 + wfeti(maz+jl-1) * wfeti(maz+jl-1)