Changeset 120 for trunk/NEMO/OPA_SRC/SOL/sol_oce.F90

Timestamp:

2004-06-29T17:39:08+02:00 (20 years ago)

Author:

opalod

Message:

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

File:

• trunk/NEMO/OPA_SRC/SOL/sol_oce.F90 (modified) (1 diff)

trunk/NEMO/OPA_SRC/SOL/sol_oce.F90

@@ -23 +23 @@
    INTEGER , PUBLIC ::      & !!: namsol   elliptic solver / island / free surface
       nsolv    =    1 ,     &  !: = 1/2/3 type of elliptic solver
-      nsol_arp =    1 ,     &  !: = 1/0   absolute/relative precision convergence test
-      nmin     =  300 ,     &  !: minimum of iterations for the solver
-      nmax     =  800 ,     &  !: maximum of iterations for the solver
-      nmod     =   10 ,     &  !: frequency of test for the solver
+      nsol_arp =    0 ,     &  !: = 0/1 absolute/relative precision convergence test
+      nmin     =  300 ,     &  !: minimum of iterations for the SOR solver
+      nmax     =  800 ,     &  !: maximum of iterations for the SOR solver
+      nmod     =   10 ,     &  !: frequency of test for the SOR solver
       nmisl    = 4000          !: maximum pcg iterations for island
      
    REAL(wp), PUBLIC ::      & !!: namsol   elliptic solver / island / free surface
-      eps    = 1.e-6_wp ,   &  !: absolute precision of the solver
-      sor    = 1.92_wp  ,   &  !: optimal coefficient for sor solver
-      epsisl = 1.e-10_wp,   &  !: absolute precision on stream function solver
-      rnu    = 1.0_wp          !: strength of the additional force used in free surface
+      eps    =  1.e-6_wp ,  &  !: absolute precision of the solver
+      resmax = 1.e-14_wp ,  &  !: absolute precision for the SOR solver
+      sor    =   1.92_wp ,  &  !: optimal coefficient for the SOR solver
+      epsisl = 1.e-10_wp ,  &  !: absolute precision on stream function solver
+      rnu    =    1.0_wp       !: strength of the additional force used in free surface
 
    CHARACTER(len=1), PUBLIC ::   &  !:
       c_solver_pt = 'T'        !: nature of grid-points T (S) for free surface case
-      !                        !                       F (G) for rigid-lid case
+      !                        !                        F (G) for rigid-lid case
 
    INTEGER , PUBLIC ::   &  !: