high_res_cube/matlab-topology-generator/driver.m

% Create exch2 communication map and schedule for a cube sphere grid with 
% constant tile size tnx x tny.

% Use red-green-blue shorthand for cube index space specification
% In this notation cube faces are laid out as shown below
%
%             
%                         f5(nr,ng)  f6(nb,ng)
%              f3(ng,nb)  f4(nr,nb)
%   f1(nb,nr)  f2(nb,nr)
% 
nr=64;  nb=32; ng=128;
%nr=576;  nb=576; ng=576;
%nr=672;  nb=672; ng=672;
nr=32; nb=32; ng=32;
%nr=510; nb=510; ng=510;
%nr=30; nb=30; ng=30;

% Choose tile subgrid sizes for each face.
% nr,nb,ng must be integer multiples of tnx and tny.
tnx=85;tny=85;
%tnx=10;tny=10;
%tnx=16;tny=16;
tnx=16;tny=32;
%tnx=192;tny=64;

% Make list of domains. Assume MITgcm standard cube layout, three color path labeling and 
% global indexing convention.
clear domain ndomains domain_nx domain_ny
[ndomains,domain,domain_nx,domain_ny] = exch2_setup_cs6_domains(nr,nb,ng);

% Now create basic tile definitions for each domain with their offsets within the domain
% tn[xy]        :: tile extents in x and y
% tbase[xy]     :: offset of tile local coords from domain coords
%               :: tdom[xy] = tlocal[xy] + tbase[xy]
% t[xy]globallo :: global composite domain coordinate associated with
%               :: tlocal[xy]=(1,1)
% mydomain      :: domain number the tile belongs to
% tileid        :: identifier number for tile
% tx            :: tile x coordinate within domain tiling
% ty            :: tile y coordinate within domain tiling
clear tile
[tile,ntiles,ierr,domain]= ...
 exch2_setup_cs6_tiles(tnx,tny,domain,ndomains);

% Set neighbor domains for each tile
[tile] = exch2_setup_cs6_get_neighbor_domains(tile, domain, ntiles);

% Let try and figure out what points I send my edges to. We do this by a search procedure 
% rather than a functional relationship. The search procedure visits each edge of each tile 
% in turn.  For internal edges (edges that don't cross a domain boundary) the index range at 
% +/-1 in the normal direction to the edge is searched for. This identifies all the tiles that % border this tile.
[tile] = exch2_setup_cs6_get_internal_neighbor_tiles(       tile, domain, ntiles);
[tile] = exch2_setup_cs6_get_internal_neighbor_index_ranges(tile, domain, ntiles);

[tile] = exch2_setup_cs6_get_external_neighbor_tiles(       tile, domain, ntiles);
% Draw a picture of the domain and its tiles in the standard cube layout
exch2_setup_cs6_plot(domain, tile, tnx, tny);

% Write attributes for a tile
exch2_setup_cs6_print( domain, tile, 1)
1	% Create exch2 communication map and schedule for a cube sphere grid with
2	% constant tile size tnx x tny.
3
4	% Use red-green-blue shorthand for cube index space specification
5	% In this notation cube faces are laid out as shown below
6	%
7	%
8	% f5(nr,ng) f6(nb,ng)
9	% f3(ng,nb) f4(nr,nb)
10	% f1(nb,nr) f2(nb,nr)
11	%
12	nr=64; nb=32; ng=128;
13	%nr=576; nb=576; ng=576;
14	%nr=672; nb=672; ng=672;
15	nr=32; nb=32; ng=32;
16	%nr=510; nb=510; ng=510;
17	%nr=30; nb=30; ng=30;
18
19	% Choose tile subgrid sizes for each face.
20	% nr,nb,ng must be integer multiples of tnx and tny.
21	tnx=85;tny=85;
22	%tnx=10;tny=10;
23	%tnx=16;tny=16;
24	tnx=16;tny=32;
25	%tnx=192;tny=64;
26
27	% Make list of domains. Assume MITgcm standard cube layout, three color path labeling and
28	% global indexing convention.
29	clear domain ndomains domain_nx domain_ny
30	[ndomains,domain,domain_nx,domain_ny] = exch2_setup_cs6_domains(nr,nb,ng);
31
32	% Now create basic tile definitions for each domain with their offsets within the domain
33	% tn[xy] :: tile extents in x and y
34	% tbase[xy] :: offset of tile local coords from domain coords
35	% :: tdom[xy] = tlocal[xy] + tbase[xy]
36	% t[xy]globallo :: global composite domain coordinate associated with
37	% :: tlocal[xy]=(1,1)
38	% mydomain :: domain number the tile belongs to
39	% tileid :: identifier number for tile
40	% tx :: tile x coordinate within domain tiling
41	% ty :: tile y coordinate within domain tiling
42	clear tile
43	[tile,ntiles,ierr,domain]= ...
44	exch2_setup_cs6_tiles(tnx,tny,domain,ndomains);
45
46	% Set neighbor domains for each tile
47	[tile] = exch2_setup_cs6_get_neighbor_domains(tile, domain, ntiles);
48
49	% Let try and figure out what points I send my edges to. We do this by a search procedure
50	% rather than a functional relationship. The search procedure visits each edge of each tile
51	% in turn. For internal edges (edges that don't cross a domain boundary) the index range at
52	% +/-1 in the normal direction to the edge is searched for. This identifies all the tiles that % border this tile.
53	[tile] = exch2_setup_cs6_get_internal_neighbor_tiles( tile, domain, ntiles);
54	[tile] = exch2_setup_cs6_get_internal_neighbor_index_ranges(tile, domain, ntiles);
55
56	[tile] = exch2_setup_cs6_get_external_neighbor_tiles( tile, domain, ntiles);
57	% Draw a picture of the domain and its tiles in the standard cube layout
58	exch2_setup_cs6_plot(domain, tile, tnx, tny);
59
60	% Write attributes for a tile
61	exch2_setup_cs6_print( domain, tile, 1)