1 c---------------------------------------------------------------------
2 c---------------------------------------------------------------------
6 c---------------------------------------------------------------------
7 c---------------------------------------------------------------------
11 c---------------------------------------------------------------------
12 c The following include file is generated automatically by the
13 c "setparams" utility. It defines
14 c maxcells: the square root of the maximum number of processors
15 c problem_size: 12, 64, 102, 162 (for class T, A, B, C)
16 c dt_default: default time step for this problem size if no
18 c niter_default: default number of iterations for this problem size
19 c---------------------------------------------------------------------
23 integer aa, bb, cc, BLOCK_SIZE
24 parameter (aa=1, bb=2, cc=3, BLOCK_SIZE=5)
26 integer ncells, grid_points(3)
27 double precision elapsed_time
28 common /global/ elapsed_time, ncells, grid_points
30 double precision tx1, tx2, tx3, ty1, ty2, ty3, tz1, tz2, tz3,
31 > dx1, dx2, dx3, dx4, dx5, dy1, dy2, dy3, dy4,
32 > dy5, dz1, dz2, dz3, dz4, dz5, dssp, dt,
33 > ce(5,13), dxmax, dymax, dzmax, xxcon1, xxcon2,
34 > xxcon3, xxcon4, xxcon5, dx1tx1, dx2tx1, dx3tx1,
35 > dx4tx1, dx5tx1, yycon1, yycon2, yycon3, yycon4,
36 > yycon5, dy1ty1, dy2ty1, dy3ty1, dy4ty1, dy5ty1,
37 > zzcon1, zzcon2, zzcon3, zzcon4, zzcon5, dz1tz1,
38 > dz2tz1, dz3tz1, dz4tz1, dz5tz1, dnxm1, dnym1,
39 > dnzm1, c1c2, c1c5, c3c4, c1345, conz1, c1, c2,
40 > c3, c4, c5, c4dssp, c5dssp, dtdssp, dttx1, bt,
41 > dttx2, dtty1, dtty2, dttz1, dttz2, c2dttx1,
42 > c2dtty1, c2dttz1, comz1, comz4, comz5, comz6,
43 > c3c4tx3, c3c4ty3, c3c4tz3, c2iv, con43, con16
45 common /constants/ tx1, tx2, tx3, ty1, ty2, ty3, tz1, tz2, tz3,
46 > dx1, dx2, dx3, dx4, dx5, dy1, dy2, dy3, dy4,
47 > dy5, dz1, dz2, dz3, dz4, dz5, dssp, dt,
48 > ce, dxmax, dymax, dzmax, xxcon1, xxcon2,
49 > xxcon3, xxcon4, xxcon5, dx1tx1, dx2tx1, dx3tx1,
50 > dx4tx1, dx5tx1, yycon1, yycon2, yycon3, yycon4,
51 > yycon5, dy1ty1, dy2ty1, dy3ty1, dy4ty1, dy5ty1,
52 > zzcon1, zzcon2, zzcon3, zzcon4, zzcon5, dz1tz1,
53 > dz2tz1, dz3tz1, dz4tz1, dz5tz1, dnxm1, dnym1,
54 > dnzm1, c1c2, c1c5, c3c4, c1345, conz1, c1, c2,
55 > c3, c4, c5, c4dssp, c5dssp, dtdssp, dttx1, bt,
56 > dttx2, dtty1, dtty2, dttz1, dttz2, c2dttx1,
57 > c2dtty1, c2dttz1, comz1, comz4, comz5, comz6,
58 > c3c4tx3, c3c4ty3, c3c4tz3, c2iv, con43, con16
60 integer EAST, WEST, NORTH, SOUTH,
63 parameter (EAST=2000, WEST=3000, NORTH=4000, SOUTH=5000,
64 > BOTTOM=6000, TOP=7000)
66 integer cell_coord (3,maxcells), cell_low (3,maxcells),
67 > cell_high (3,maxcells), cell_size(3,maxcells),
68 > predecessor(3), slice (3,maxcells),
69 > grid_size (3), successor(3) ,
70 > start (3,maxcells), end (3,maxcells)
71 common /partition/ cell_coord, cell_low, cell_high, cell_size,
72 > grid_size, successor, predecessor, slice,
75 integer IMAX, JMAX, KMAX, MAX_CELL_DIM, BUF_SIZE
77 parameter (MAX_CELL_DIM = (problem_size/maxcells)+1)
79 parameter (IMAX=MAX_CELL_DIM,JMAX=MAX_CELL_DIM,KMAX=MAX_CELL_DIM)
81 parameter (BUF_SIZE=MAX_CELL_DIM*MAX_CELL_DIM*(maxcells-1)*60+1)
84 > us ( -1:IMAX, -1:JMAX, -1:KMAX, maxcells),
85 > vs ( -1:IMAX, -1:JMAX, -1:KMAX, maxcells),
86 > ws ( -1:IMAX, -1:JMAX, -1:KMAX, maxcells),
87 > qs ( -1:IMAX, -1:JMAX, -1:KMAX, maxcells),
88 > rho_i ( -1:IMAX, -1:JMAX, -1:KMAX, maxcells),
89 > square ( -1:IMAX, -1:JMAX, -1:KMAX, maxcells),
90 > forcing (5, 0:IMAX-1, 0:JMAX-1, 0:KMAX-1, maxcells),
91 > u (5, -2:IMAX+1,-2:JMAX+1,-2:KMAX+1, maxcells),
92 > rhs (5, -1:IMAX-1,-1:JMAX-1,-1:KMAX-1, maxcells),
93 > lhsc (5,5,-1:IMAX-1,-1:JMAX-1,-1:KMAX-1, maxcells),
94 > backsub_info (5, 0:MAX_CELL_DIM, 0:MAX_CELL_DIM, maxcells),
95 > in_buffer(BUF_SIZE), out_buffer(BUF_SIZE)
96 common /fields/ u, us, vs, ws, qs, rho_i, square,
97 > rhs, forcing, lhsc, in_buffer, out_buffer,
100 double precision cv(-2:MAX_CELL_DIM+1), rhon(-2:MAX_CELL_DIM+1),
101 > rhos(-2:MAX_CELL_DIM+1), rhoq(-2:MAX_CELL_DIM+1),
102 > cuf(-2:MAX_CELL_DIM+1), q(-2:MAX_CELL_DIM+1),
103 > ue(-2:MAX_CELL_DIM+1,5), buf(-2:MAX_CELL_DIM+1,5)
104 common /work_1d/ cv, rhon, rhos, rhoq, cuf, q, ue, buf
106 integer west_size, east_size, bottom_size, top_size,
107 > north_size, south_size, start_send_west,
108 > start_send_east, start_send_south, start_send_north,
109 > start_send_bottom, start_send_top, start_recv_west,
110 > start_recv_east, start_recv_south, start_recv_north,
111 > start_recv_bottom, start_recv_top
112 common /box/ west_size, east_size, bottom_size,
113 > top_size, north_size, south_size,
114 > start_send_west, start_send_east, start_send_south,
115 > start_send_north, start_send_bottom, start_send_top,
116 > start_recv_west, start_recv_east, start_recv_south,
117 > start_recv_north, start_recv_bottom, start_recv_top
119 double precision tmp_block(5,5), b_inverse(5,5), tmp_vec(5)
120 common /work_solve/ tmp_block, b_inverse, tmp_vec
123 c These are used by btio
125 integer collbuf_nodes, collbuf_size, iosize, eltext,
126 $ combined_btype, fp, idump, record_length, element,
127 $ combined_ftype, idump_sub, rd_interval
128 common /btio/ collbuf_nodes, collbuf_size, iosize, eltext,
129 $ combined_btype, fp, idump, record_length,
130 $ idump_sub, rd_interval
131 double precision sum(niter_default), xce_sub(5)
132 common /btio/ sum, xce_sub
134 common /btio/ iseek, element, combined_ftype