/[MITgcm]/MITgcm/pkg/exch2/exch2_get_scal_bounds.F
ViewVC logotype

Annotation of /MITgcm/pkg/exch2/exch2_get_scal_bounds.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.1 - (hide annotations) (download)
Sat May 30 21:22:13 2009 UTC (15 years ago) by jmc
Branch: MAIN
CVS Tags: checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint64, checkpoint65, checkpoint62, checkpoint63, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65o, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint62c, checkpoint62b, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61q, checkpoint61z, checkpoint61x, checkpoint61y, HEAD 
simplify argument list of S/R exch2_get_uv_bounds and exch2_get_scal_bounds
 (which replaces exch2_get_recv_bounds & exch2_get_send_bounds) by including
 Topology header file (and common blocks);
implement EXCH_IGNORE_CORNERS in S/R exch2_get_scal_bounds.
1 jmc 1.1 C $Header: /u/gcmpack/MITgcm/pkg/exch2/exch2_get_send_bounds.F,v 1.3 2007/10/09 00:03:33 jmc Exp $
2     C $Name: $
3    
4     #include "CPP_EEOPTIONS.h"
5    
6     CBOP 0
7     C !ROUTINE: EXCH2_GET_SCAL_BOUNDS
8    
9     C !INTERFACE:
10     SUBROUTINE EXCH2_GET_SCAL_BOUNDS(
11     I fCode, eWdth, updateCorners,
12     I tgTile, tgNb,
13     O tIlo, tIhi, tJlo, tJhi,
14     O tiStride, tjStride,
15     I myThid )
16    
17     C !DESCRIPTION:
18     C Return the index range & step of the part of the array (overlap-region)
19     C which is going to be updated by the exchange with 1 neighbour.
20     C Scalar field version (as opposed to 2-components vector field).
21     C Note: presently only coded for grid-centered (mass point) scalar.
22     C needs to extend it to grid-corner (vorticity point) scalar.
23    
24     C !USES:
25     IMPLICIT NONE
26    
27     #include "SIZE.h"
28     #include "W2_EXCH2_SIZE.h"
29     #include "W2_EXCH2_TOPOLOGY.h"
30    
31     C !INPUT/OUTPUT PARAMETERS:
32     C fCode :: field code (position on staggered grid)
33     C eWdth :: width of data region to exchange
34     C updateCorners :: flag, do update corner halo region if true
35     C tgTile :: target tile
36     C tgNb :: target Neighbour entry
37     C tIlo, tIhi :: index range in I that will be filled in target array
38     C tJlo, tJhi :: index range in J that will be filled in target array
39     C tiStride :: index step in I that will be filled in target array
40     C tjStride :: index step in J that will be filled in target array
41     C myThid :: my Thread Id. number
42    
43     CHARACTER*2 fCode
44     INTEGER eWdth
45     LOGICAL updateCorners
46     INTEGER tgTile, tgNb
47     INTEGER tIlo, tIhi, tJlo, tJhi
48     INTEGER tiStride, tjStride
49     INTEGER myThid
50    
51     C !LOCAL VARIABLES:
52     CEOP
53    
54     C--- Initialise index range from Topology values:
55     tIlo = exch2_iLo(tgNb,tgTile)
56     tIhi = exch2_iHi(tgNb,tgTile)
57     tJlo = exch2_jLo(tgNb,tgTile)
58     tJhi = exch2_jHi(tgNb,tgTile)
59    
60     C--- Expand index range according to exchange-Width "eWdth"
61     IF ( tIlo .EQ. tIhi .AND. tIlo .EQ. 0 ) THEN
62     C Filling a west edge overlap
63     tIlo = 1-eWdth
64     tiStride = 1
65     IF ( tJlo .LE. tJhi ) THEN
66     tjStride = 1
67     ELSE
68     tjStride =-1
69     ENDIF
70     IF ( updateCorners ) THEN
71     tJlo = tJlo-tjStride*(eWdth-1)
72     tJhi = tJhi+tjStride*(eWdth-1)
73     ELSE
74     tJlo = tJlo+tjStride
75     tJhi = tJhi-tjStride
76     ENDIF
77     ENDIF
78     IF ( tIlo .EQ. tIhi .AND. tIlo .GT. 1 ) THEN
79     C Filling an east edge overlap
80     tIhi = tIhi+eWdth-1
81     tiStride = 1
82     IF ( tJlo .LE. tJhi ) THEN
83     tjStride = 1
84     ELSE
85     tjStride =-1
86     ENDIF
87     IF ( updateCorners ) THEN
88     tJlo = tJlo-tjStride*(eWdth-1)
89     tJhi = tJhi+tjStride*(eWdth-1)
90     ELSE
91     tJlo = tJlo+tjStride
92     tJhi = tJhi-tjStride
93     ENDIF
94     ENDIF
95     IF ( tJlo .EQ. tJhi .AND. tJlo .EQ. 0 ) THEN
96     C Filling a south edge overlap
97     tJlo = 1-eWdth
98     tjStride = 1
99     IF ( tIlo .LE. tIhi ) THEN
100     tiStride = 1
101     ELSE
102     tiStride =-1
103     ENDIF
104     IF ( updateCorners ) THEN
105     tIlo = tIlo-tiStride*(eWdth-1)
106     tIhi = tIhi+tiStride*(eWdth-1)
107     ELSE
108     tIlo = tIlo+tiStride
109     tIhi = tIhi-tiStride
110     ENDIF
111     ENDIF
112     IF ( tJlo .EQ. tJhi .AND. tJlo .GT. 1 ) THEN
113     C Filling a north edge overlap
114     tJhi = tJhi+eWdth-1
115     tjStride = 1
116     IF ( tIlo .LE. tIhi ) THEN
117     tiStride = 1
118     ELSE
119     tiStride =-1
120     ENDIF
121     IF ( updateCorners ) THEN
122     tIlo = tIlo-tiStride*(eWdth-1)
123     tIhi = tIhi+tiStride*(eWdth-1)
124     ELSE
125     tIlo = tIlo+tiStride
126     tIhi = tIhi-tiStride
127     ENDIF
128     ENDIF
129    
130     RETURN
131     END
  ViewVC Help
Powered by ViewVC 1.1.22