Back to home page

MITgcm

 
 

    


File indexing completed on 2018-03-02 18:39:44 UTC

view on githubraw file Latest commit add29e06 on 2018-01-31 20:35:05 UTC
6d0e054297 Jean*0001 CBOP
                0002 C     !ROUTINE: W2_EXCH2_PARAMS.h
                0003 C     !INTERFACE:
                0004 C     #include W2_EXCH2_PARAMS.h
                0005 
                0006 C     !DESCRIPTION:
                0007 C     *==========================================================*
                0008 C     | W2_EXCH2_PARAMS.h
017b6b2289 Jean*0009 C     | o Header file defining WRAPPER2-EXCH2 topology parameters
6d0e054297 Jean*0010 C     *==========================================================*
                0011 CEOP
                0012 
017b6b2289 Jean*0013 C--   COMMON /W2_EXCH2_PARM_I/ Integer valued parameters used by W2-EXCH2
                0014 C     preDefTopol :: pre-defined Topology selector:
                0015 C                 :: = 0 : topology defined from processing "data.exch2";
                0016 C                 :: = 1 : simple, single facet topology;
                0017 C                 :: = 2 : customized topology (w2_set_myown_facets)
                0018 C                 :: = 3 : 6-face Cube (3 face-dims: nRed, nGreen, nBlue).
                0019 C     nFacets     :: Number of facets (or domains) in this topology
                0020 C     facet_dims  :: facet pair of dimensions (n1x,n1y, n2x,n2y, ...)
d4ac4b712d Jean*0021 C     nTiles      :: Number of tiles in this topology <- moved to W2_EXCH2_TOPOLOGY.h
017b6b2289 Jean*0022 C     nBlankTiles :: Number of "Blank-Tiles"
                0023 C     blankList   :: List of "Blank-Tiles" (non active)
                0024 C--
                0025 C     W2_mapIO    :: select option for global-IO mapping:
                0026 C             =-1 :: old format, put facets 1 after the other in the X dir.;
                0027 C                 :: this is not necessary "compact".
                0028 C             = 1 :: compact format, 1 facet after the other (mostly in Y dir.)
                0029 C                 :: but needs to fold some facet (domain) if too large
                0030 C             = 0 :: compact format (= 1 long line), one facet after the other.
                0031 C--
                0032 C     W2_oUnit    :: output fortran Unit for printing Std messages
                0033 C     W2_printMsg :: select option for information messages printing
                0034 C             < 0 :: open and print to "w2_tile_topology" log file
d6ea3164dc Jean*0035 C             = 0 :: print the minimum, to StdOut
017b6b2289 Jean*0036 C             = 1 :: no duplicated print from different processes (only Master)
                0037 C             = 2 :: all processes do print (used to check).
                0038       INTEGER preDefTopol
                0039       INTEGER nFacets
                0040       INTEGER facet_dims(2*W2_maxNbFacets)
                0041       INTEGER nBlankTiles
                0042       INTEGER blankList(W2_maxNbTiles)
                0043       INTEGER W2_mapIO
                0044       INTEGER W2_oUnit, W2_printMsg
                0045       COMMON /W2_EXCH2_PARM_I/
                0046      &        preDefTopol,
                0047      &        nFacets, facet_dims,
d4ac4b712d Jean*0048      &        nBlankTiles, blankList,
017b6b2289 Jean*0049      &        W2_mapIO,
                0050      &        W2_oUnit, W2_printMsg
046fd16d1c Andr*0051 
6fc6b32488 Jean*0052 C--   COMMON /W2_EXCH2_PARM_L/ Logical valued parameters used by W2-EXCH2
                0053 C     W2_useE2ioLayOut :: =T: use Exch2 global-IO Layout; =F: use model default
                0054       LOGICAL W2_useE2ioLayOut
                0055       COMMON /W2_EXCH2_PARM_L/
                0056      &        W2_useE2ioLayOut
                0057 
017b6b2289 Jean*0058 C--   COMMON /W2_EXCH2_BUILD_I/ Integer variables used to build topology
                0059 C     facet_owns  :: Range of tiles this facet "owns"
                0060 C     facet_pij   \  ::
                0061 C     facet_oi     } :: indices correspondence matrix (facet_pij) & offsets:
                0062 C     facet_oj    /  ::
                0063 C-with:  suffix "so" for indices of source facet j ;
                0064 C        suffix "tg" for indices of target facet jj= INT(facet_link(i,j))
                0065 C      pij(:,i,j) : matrix which gives so indices when applied to tg indices
                0066 C        iso = pij(1)*itg + pij(2)*jtg + oi
                0067 C        jso = pij(3)*itg + pij(4)*jtg + oj
                0068 C-----
                0069       INTEGER facet_owns(2,W2_maxNbFacets)
                0070       INTEGER facet_pij(4,4,W2_maxNbFacets)
                0071       INTEGER facet_oi(4,W2_maxNbFacets)
                0072       INTEGER facet_oj(4,W2_maxNbFacets)
                0073       COMMON /W2_EXCH2_BUILD_I/
                0074      &        facet_owns,
                0075      &        facet_pij, facet_oi, facet_oj
046fd16d1c Andr*0076 
017b6b2289 Jean*0077 C--   COMMON /W2_EXCH2_PARM_R/ Real*4 valued parameters used by W2-EXCH2
                0078 C--   topology defined from processing "data.exch2" (preDefTopol=0):
                0079 C     facet_link  :: Face-Edge connectivity map:
                0080 C       facet_link(i,j)=XX.1 : face(j)-edge(i) (i=1,2,3,4 <==> N,S,E,W)
                0081 C       is connected to Northern edge of face "XX" ; similarly,
                0082 C       = XX.2 : to Southern.E, XX.3 = Eastern.E, XX.4 = Western.E of face "XX".
                0083       Real*4  facet_link( 4, W2_maxNbFacets )
                0084       COMMON /W2_EXCH2_PARM_R/ facet_link
046fd16d1c Andr*0085 
017b6b2289 Jean*0086 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|