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

Timestamp:

2005-09-30T12:20:54+02:00 (19 years ago)

Author:

opalod

Message:

nemo_v1_update_017:RB: added a new solver (nsolv=4) corresponding to solsor with extra outer halo.

File:

• trunk/NEMO/OPA_SRC/SOL/sol_oce.F90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/SOL/sol_oce.F90

@@ -24 +24 @@
    !! ----------------------------------
    INTEGER , PUBLIC ::      & !!: namsol   elliptic solver / island / free surface
-      nsolv    =    1 ,     &  !: = 1/2/3 type of elliptic solver
+      nsolv    =    1 ,     &  !: = 1/2/3/4 type of elliptic solver
       nsol_arp =    0 ,     &  !: = 0/1 absolute/relative precision convergence test
       nmin     =  300 ,     &  !: minimum of iterations for the SOR solver
@@ -55 +55 @@
       rr               !: coefficient  =(rn,rn)
 
-   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,4) ::   &  !:
+   REAL(wp), PUBLIC, DIMENSION(1-jpr2di:jpi+jpr2di,1-jpr2dj:jpj+jpr2dj,4) ::   &  !:
       gcp              !: barotropic matrix extra-diagonal elements
 
-   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   &  !:
+   REAL(wp), PUBLIC, DIMENSION(1-jpr2di:jpi+jpr2di,1-jpr2dj:jpj+jpr2dj) ::   &  !:
       gcx, gcxb,    &  !: now, before solution of the elliptic equation
       gcdprc,       &  !: inverse diagonal preconditioning matrix