9 |
% f3(ng,nb) f4(nr,nb) |
% f3(ng,nb) f4(nr,nb) |
10 |
% f1(nb,nr) f2(nb,nr) |
% f1(nb,nr) f2(nb,nr) |
11 |
% |
% |
12 |
nr=64; nb=32; ng=128; |
%nr=64; nb=32; ng=128; |
13 |
%nr=576; nb=576; ng=576; |
%nr=576; nb=576; ng=576; |
14 |
%nr=672; nb=672; ng=672; |
%nr=672; nb=672; ng=672; |
15 |
nr=32; nb=32; ng=32; |
nr=32; nb=32; ng=32; |
18 |
|
|
19 |
% Choose tile subgrid sizes for each face. |
% Choose tile subgrid sizes for each face. |
20 |
% nr,nb,ng must be integer multiples of tnx and tny. |
% nr,nb,ng must be integer multiples of tnx and tny. |
21 |
tnx=85;tny=85; |
%tnx=85;tny=85; |
22 |
%tnx=10;tny=10; |
%tnx=10;tny=10; |
23 |
%tnx=16;tny=16; |
%tnx=16;tny=16; |
24 |
tnx=16;tny=32; |
tnx=16;tny=32; |
25 |
|
%tnx=32;tny=32; |
26 |
|
%tnx=32;tny=8; |
27 |
%tnx=192;tny=64; |
%tnx=192;tny=64; |
28 |
|
tnx=8;tny=4; |
29 |
|
|
30 |
% Make list of domains. Assume MITgcm standard cube layout, three color path labeling and |
% Make list of domains. Assume MITgcm standard cube layout, three color path labeling and |
31 |
% global indexing convention. |
% global indexing convention. |
57 |
[tile] = exch2_setup_cs6_get_internal_neighbor_index_ranges(tile, domain, ntiles); |
[tile] = exch2_setup_cs6_get_internal_neighbor_index_ranges(tile, domain, ntiles); |
58 |
|
|
59 |
[tile] = exch2_setup_cs6_get_external_neighbor_tiles( tile, domain, ntiles); |
[tile] = exch2_setup_cs6_get_external_neighbor_tiles( tile, domain, ntiles); |
60 |
% Draw a picture of the domain and its tiles in the standard cube layout |
|
61 |
exch2_setup_cs6_plot(domain, tile, tnx, tny); |
% Draw a picture of the full domain and its tiles in the standard cube layout |
62 |
|
% exch2_setup_cs6_plot(domain, tile, tnx, tny); |
63 |
|
|
64 |
|
% Squeeze the blank tiles out |
65 |
|
[tile, domain] = exch2_setup_squeeze_blanks(domain, tile, tnx, tny); |
66 |
|
|
67 |
% Write attributes for a tile |
% Write attributes for a tile |
68 |
exch2_setup_cs6_print( domain, tile, 1) |
exch2_setup_cs6_print( domain, tile, tnx, tny, 1) |