Changeset 352 for trunk/NEMO/OPA_SRC/mppini.F90

Timestamp:

2005-12-12T15:18:10+01:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_033 : CT : Switch to IOIPSL-3-0 new library

File:

• trunk/NEMO/OPA_SRC/mppini.F90 (modified) (5 diffs)

trunk/NEMO/OPA_SRC/mppini.F90

@@ -16 +16 @@
    USE sol_oce         ! ocean elliptic solver
    USE lib_mpp         ! distribued memory computing library
+   USE ioipsl
 
    IMPLICIT NONE
@@ -51 +52 @@
       !!   9.0  !  04-01  (G. Madec, J.M. Molines)  F90 : free form, north fold jpni >1
       !!----------------------------------------------------------------------
+
       ! No mpp computation
-
       nimpp  = 1
       njmpp  = 1
@@ -64 +65 @@
       nbondi = 2
       nbondj = 2
+      nidom  = FLIO_DOM_NONE
+
       IF(lwp) THEN
          WRITE(numout,*)
@@ -80 +83 @@
           nstop = nstop + 1
       ENDIF
+
 
    END SUBROUTINE mpp_init
@@ -499 +503 @@
       !!   9.0  !  04-03  (G. Madec)  MPP-IOIPSL 
       !!----------------------------------------------------------------------
-      USE ioipsl
-      INTEGER, DIMENSION(4) ::   &
-         iglo, iloc, iabsf, iabsl, ihals, ihale   ! ???
-      !!----------------------------------------------------------------------
+      !! Local declarations
+
+      INTEGER, DIMENSION(2) ::   &
+         iglo, iloc, iabsf, iabsl, ihals, ihale, idid
+      INTEGER ::   &
+         inbondi, inbondj
+      !!----------------------------------------------------------------------
+
+      ! The domain is splitted only horizontally along i- or/and j- direction
+      ! So we need at the most only 1D arrays with 2 elements
       iglo(1) = jpiglo
       iglo(2) = jpjglo
-      iglo(3) = jpk
-      iglo(4) = 1
       iloc(1) = nlci
       iloc(2) = nlcj
-      iloc(3) = jpk
-      iloc(4) = 1
       iabsf(1) = nimppt(narea)
       iabsf(2) = njmppt(narea)
-      iabsf(3) = 1
-      iabsf(4) = 1
       iabsl(:) = iabsf(:) + iloc(:) - 1
       ihals(1) = jpreci
       ihals(2) = jprecj
-      ihals(3) = 0
-      ihals(4) = 0
       ihale(1) = jpreci
       ihale(2) = jprecj
-      ihale(3) = 0
-      ihale(4) = 0
-      IF( nbondi == -1 .OR. nbondi == 2 )   ihals(1) = 0
-      IF( nbondi ==  1 .OR. nbondi == 2 )   ihale(1) = 0
-      IF( nbondj == -1 .OR. nbondj == 2 )   ihals(2) = 0
-      IF( nbondj ==  1 .OR. nbondj == 2 )   ihale(2) = 0
+      idid(1) = 1
+      idid(2) = 2
+
+      inbondj = -1                                   ! general case
+      IF( narea   >  jpni          )   inbondj = 0      ! first row of processor
+      IF( narea   >  (jpnj-1)*jpni )   inbondj = 1      ! last  row of processor
+      IF( jpnj    == 1             )   inbondj = 2      ! one processor only in j-direction
+      IF( inbondj == -1 .OR. inbondj == 2 )   ihals(2) = 0
+      IF( inbondj ==  1 .OR. inbondj == 2 )   ihale(2) = 0
+         
+      inbondi = 0                                    ! 
+      IF( MOD( narea, jpni ) == 1 )   inbondi = -1      !
+      IF( MOD( narea, jpni ) == 0 )   inbondi =  1      !
+      IF( jpni               == 1 )   inbondi =  2      ! one processor only in i-direction
+      IF( inbondi == -1 .OR. inbondi == 2 )   ihals(1) = 0
+      IF( inbondi ==  1 .OR. inbondi == 2 )   ihale(1) = 0
+
+
       IF(lwp) THEN
-          WRITE(numout,*) 'mpp_init_ioipsl :   iloc  = ', iloc (1), iloc (2), iloc (3), iloc (4)
-          WRITE(numout,*) '~~~~~~~~~~~~~~~     iabsf = ', iabsf(1), iabsf(2), iabsf(3), iabsf(4)
-          WRITE(numout,*) '                    ihals = ', ihals(1), ihals(2), ihals(3), ihals(4)
-          WRITE(numout,*) '                    ihale = ', ihale(1), ihale(2), ihale(3), ihale(4)
-      ENDIF
-
-      CALL ioipsl_inimpp( jpnij, nproc, iglo, iloc, iabsf, iabsl, ihals, ihale )
+          WRITE(numout,*) 'mpp_init_ioipsl :   iloc  = ', iloc (1), iloc (2)
+          WRITE(numout,*) '~~~~~~~~~~~~~~~     iabsf = ', iabsf(1), iabsf(2)
+          WRITE(numout,*) '                    ihals = ', ihals(1), ihals(2)
+          WRITE(numout,*) '                    ihale = ', ihale(1), ihale(2)
+      ENDIF
+
+      CALL flio_dom_set ( jpnij, nproc, idid, iglo, iloc, iabsf, iabsl, ihals, ihale, 'BOX', nidom)
 
    END SUBROUTINE mpp_init_ioipsl