1 |
C $Header: /u/gcmpack/MITgcm_contrib/high_res_cube/cube92_latest/code/SIZE.h,v 1.2 2019/07/31 19:47:11 dimitri Exp $ |
2 |
CBOP |
3 |
C !ROUTINE: SIZE.h |
4 |
C !INTERFACE: |
5 |
C include SIZE.h |
6 |
C !DESCRIPTION: \bv |
7 |
C *==========================================================* |
8 |
C | SIZE.h Declare size of underlying computational grid. |
9 |
C *==========================================================* |
10 |
C | The design here supports a three-dimensional model grid |
11 |
C | with indices I,J and K. The three-dimensional domain |
12 |
C | is comprised of nPx*nSx blocks (or tiles) of size sNx |
13 |
C | along the first (left-most index) axis, nPy*nSy blocks |
14 |
C | of size sNy along the second axis and one block of size |
15 |
C | Nr along the vertical (third) axis. |
16 |
C | Blocks/tiles have overlap regions of size OLx and OLy |
17 |
C | along the dimensions that are subdivided. |
18 |
C *==========================================================* |
19 |
C \ev |
20 |
C |
21 |
C Voodoo numbers controlling data layout: |
22 |
C sNx :: Number of X points in tile. |
23 |
C sNy :: Number of Y points in tile. |
24 |
C OLx :: Tile overlap extent in X. |
25 |
C OLy :: Tile overlap extent in Y. |
26 |
C nSx :: Number of tiles per process in X. |
27 |
C nSy :: Number of tiles per process in Y. |
28 |
C nPx :: Number of processes to use in X. |
29 |
C nPy :: Number of processes to use in Y. |
30 |
C Nx :: Number of points in X for the full domain. |
31 |
C Ny :: Number of points in Y for the full domain. |
32 |
C Nr :: Number of points in vertical direction. |
33 |
CEOP |
34 |
INTEGER sNx |
35 |
INTEGER sNy |
36 |
INTEGER OLx |
37 |
INTEGER OLy |
38 |
INTEGER nSx |
39 |
INTEGER nSy |
40 |
INTEGER nPx |
41 |
INTEGER nPy |
42 |
INTEGER Nx |
43 |
INTEGER Ny |
44 |
INTEGER Nr |
45 |
PARAMETER ( |
46 |
& sNx = 34, |
47 |
& sNy = 102, |
48 |
& OLx = 8, |
49 |
& OLy = 8, |
50 |
& nSx = 1, |
51 |
& nSy = 1, |
52 |
& nPx = 450, |
53 |
& nPy = 1, |
54 |
& Nx = sNx*nSx*nPx, |
55 |
& Ny = sNy*nSy*nPy, |
56 |
& Nr = 50) |
57 |
|
58 |
C MAX_OLX :: Set to the maximum overlap region size of any array |
59 |
C MAX_OLY that will be exchanged. Controls the sizing of exch |
60 |
C routine buffers. |
61 |
INTEGER MAX_OLX |
62 |
INTEGER MAX_OLY |
63 |
PARAMETER ( MAX_OLX = OLx, |
64 |
& MAX_OLY = OLy ) |
65 |
|