2 c---------------------------------------------------------------------
3 c---------------------------------------------------------------------
7 c---------------------------------------------------------------------
8 c The following include file is generated automatically by the
9 c "setparams" utility. It defines
10 c maxcells: the square root of the maximum number of processors
11 c problem_size: 12, 64, 102, 162 (for class T, A, B, C)
12 c dt_default: default time step for this problem size if no
14 c niter_default: default number of iterations for this problem size
15 c---------------------------------------------------------------------
19 integer ncells, grid_points(3)
20 common /global/ ncells, grid_points
22 double precision tx1, tx2, tx3, ty1, ty2, ty3, tz1, tz2, tz3,
23 > dx1, dx2, dx3, dx4, dx5, dy1, dy2, dy3, dy4,
24 > dy5, dz1, dz2, dz3, dz4, dz5, dssp, dt,
25 > ce(5,13), dxmax, dymax, dzmax, xxcon1, xxcon2,
26 > xxcon3, xxcon4, xxcon5, dx1tx1, dx2tx1, dx3tx1,
27 > dx4tx1, dx5tx1, yycon1, yycon2, yycon3, yycon4,
28 > yycon5, dy1ty1, dy2ty1, dy3ty1, dy4ty1, dy5ty1,
29 > zzcon1, zzcon2, zzcon3, zzcon4, zzcon5, dz1tz1,
30 > dz2tz1, dz3tz1, dz4tz1, dz5tz1, dnxm1, dnym1,
31 > dnzm1, c1c2, c1c5, c3c4, c1345, conz1, c1, c2,
32 > c3, c4, c5, c4dssp, c5dssp, dtdssp, dttx1, bt,
33 > dttx2, dtty1, dtty2, dttz1, dttz2, c2dttx1,
34 > c2dtty1, c2dttz1, comz1, comz4, comz5, comz6,
35 > c3c4tx3, c3c4ty3, c3c4tz3, c2iv, con43, con16
37 common /constants/ tx1, tx2, tx3, ty1, ty2, ty3, tz1, tz2, tz3,
38 > dx1, dx2, dx3, dx4, dx5, dy1, dy2, dy3, dy4,
39 > dy5, dz1, dz2, dz3, dz4, dz5, dssp, dt,
40 > ce, dxmax, dymax, dzmax, xxcon1, xxcon2,
41 > xxcon3, xxcon4, xxcon5, dx1tx1, dx2tx1, dx3tx1,
42 > dx4tx1, dx5tx1, yycon1, yycon2, yycon3, yycon4,
43 > yycon5, dy1ty1, dy2ty1, dy3ty1, dy4ty1, dy5ty1,
44 > zzcon1, zzcon2, zzcon3, zzcon4, zzcon5, dz1tz1,
45 > dz2tz1, dz3tz1, dz4tz1, dz5tz1, dnxm1, dnym1,
46 > dnzm1, c1c2, c1c5, c3c4, c1345, conz1, c1, c2,
47 > c3, c4, c5, c4dssp, c5dssp, dtdssp, dttx1, bt,
48 > dttx2, dtty1, dtty2, dttz1, dttz2, c2dttx1,
49 > c2dtty1, c2dttz1, comz1, comz4, comz5, comz6,
50 > c3c4tx3, c3c4ty3, c3c4tz3, c2iv, con43, con16
52 integer EAST, WEST, NORTH, SOUTH,
55 parameter (EAST=2000, WEST=3000, NORTH=4000, SOUTH=5000,
56 > BOTTOM=6000, TOP=7000)
58 integer cell_coord (3,maxcells), cell_low (3,maxcells),
59 > cell_high (3,maxcells), cell_size(3,maxcells),
60 > predecessor(3), slice (3,maxcells),
61 > grid_size (3), successor(3),
62 > start (3,maxcells), end (3,maxcells)
63 common /partition/ cell_coord, cell_low, cell_high, cell_size,
64 > grid_size, successor, predecessor, slice,
67 integer IMAX, JMAX, KMAX, MAX_CELL_DIM, BUF_SIZE, IMAXP, JMAXP
69 parameter (MAX_CELL_DIM = (problem_size/maxcells)+1)
71 parameter (IMAX=MAX_CELL_DIM,JMAX=MAX_CELL_DIM,KMAX=MAX_CELL_DIM)
72 parameter (IMAXP=IMAX/2*2+1,JMAXP=JMAX/2*2+1)
74 c---------------------------------------------------------------------
75 c +1 at end to avoid zero length arrays for 1 node
76 c---------------------------------------------------------------------
77 parameter (BUF_SIZE=MAX_CELL_DIM*MAX_CELL_DIM*(maxcells-1)*60*2+1)
80 > u (-2:IMAXP+1,-2:JMAXP+1,-2:KMAX+1, 5,maxcells),
81 > us (-1:IMAX, -1:JMAX, -1:KMAX, maxcells),
82 > vs (-1:IMAX, -1:JMAX, -1:KMAX, maxcells),
83 > ws (-1:IMAX, -1:JMAX, -1:KMAX, maxcells),
84 > qs (-1:IMAX, -1:JMAX, -1:KMAX, maxcells),
85 > ainv (-1:IMAX, -1:JMAX, -1:KMAX, maxcells),
86 > rho_i (-1:IMAX, -1:JMAX, -1:KMAX, maxcells),
87 > speed (-1:IMAX, -1:JMAX, -1:KMAX, maxcells),
88 > square (-1:IMAX, -1:JMAX, -1:KMAX, maxcells),
89 > rhs ( 0:IMAXP-1, 0:JMAXP-1, 0:KMAX-1, 5,maxcells),
90 > forcing ( 0:IMAXP-1, 0:JMAXP-1, 0:KMAX-1, 5,maxcells),
91 > lhs ( 0:IMAXP-1, 0:JMAXP-1, 0:KMAX-1,15,maxcells),
92 > in_buffer(BUF_SIZE), out_buffer(BUF_SIZE)
93 common /fields/ u, us, vs, ws, qs, ainv, rho_i, speed, square,
94 > rhs, forcing, lhs, in_buffer, out_buffer
96 double precision cv(-2:MAX_CELL_DIM+1), rhon(-2:MAX_CELL_DIM+1),
97 > rhos(-2:MAX_CELL_DIM+1), rhoq(-2:MAX_CELL_DIM+1),
98 > cuf(-2:MAX_CELL_DIM+1), q(-2:MAX_CELL_DIM+1),
99 > ue(-2:MAX_CELL_DIM+1,5), buf(-2:MAX_CELL_DIM+1,5)
100 common /work_1d/ cv, rhon, rhos, rhoq, cuf, q, ue, buf
102 integer west_size, east_size, bottom_size, top_size,
103 > north_size, south_size, start_send_west,
104 > start_send_east, start_send_south, start_send_north,
105 > start_send_bottom, start_send_top, start_recv_west,
106 > start_recv_east, start_recv_south, start_recv_north,
107 > start_recv_bottom, start_recv_top
108 common /box/ west_size, east_size, bottom_size,
109 > top_size, north_size, south_size,
110 > start_send_west, start_send_east, start_send_south,
111 > start_send_north, start_send_bottom, start_send_top,
112 > start_recv_west, start_recv_east, start_recv_south,
113 > start_recv_north, start_recv_bottom, start_recv_top