MITgcm_contrib/quarter_degree_global/size.m

% determine possible combinations of tile sizes
clear all, close all, Nx=2880; Ny=2176; px=[]; py=[];
for i=2:(Nx/30)
  if mod(Nx,i)==0, px=[px i]; end
end
for i=2:(Ny/30)
  if mod(Ny,i)==0, py=[py i]; end
end
n=0;
for i=1:length(px)
  for j=1:length(py)
    n=n+1;
    pnx(n)=px(i);
    snx(n)=Nx/pnx(n);
    pny(n)=py(j);
    sny(n)=Ny/pny(n);
    pnxpny(n)=pnx(n)*pny(n);
  end
end
[pnxpny i]=sort(pnxpny);
pnx=pnx(i); pny=pny(i);
snx=snx(i); sny=sny(i);
[pnxpny' pnx' pny' snx' sny']


     npx*npy         npx         npy         snx         sny
         480          30          16          96         136
         510          30          17          96         128
         512          32          16          90         136
         544          32          17          90         128
         576          36          16          80         136
         612          36          17          80         128
         640          40          16          72         136
         680          40          17          72         128
         720          45          16          64         136
         765          45          17          64         128
         768          24          32         120          68
         816          24          34         120          64
         960          30          32          96          68
        1020          30          34          96          64
        1024          32          32          90          68
        1088          32          34          90          64
        1152          36          32          80          68
        1224          36          34          80          64
        1280          40          32          72          68
        1360          40          34          72          64
        1440          45          32          64          68
        1530          45          34          64          64
        1536          48          32          60          68
        1632          48          34          60          64
        1920          60          32          48          68
        2040          60          34          48          64
1	% determine possible combinations of tile sizes
2	clear all, close all, Nx=2880; Ny=2176; px=[]; py=[];
3	for i=2:(Nx/30)
4	if mod(Nx,i)==0, px=[px i]; end
5	end
6	for i=2:(Ny/30)
7	if mod(Ny,i)==0, py=[py i]; end
8	end
9	n=0;
10	for i=1:length(px)
11	for j=1:length(py)
12	n=n+1;
13	pnx(n)=px(i);
14	snx(n)=Nx/pnx(n);
15	pny(n)=py(j);
16	sny(n)=Ny/pny(n);
17	pnxpny(n)=pnx(n)*pny(n);
18	end
19	end
20	[pnxpny i]=sort(pnxpny);
21	pnx=pnx(i); pny=pny(i);
22	snx=snx(i); sny=sny(i);
23	[pnxpny' pnx' pny' snx' sny']
24
25
26	npx*npy npx npy snx sny
27	480 30 16 96 136
28	510 30 17 96 128
29	512 32 16 90 136
30	544 32 17 90 128
31	576 36 16 80 136
32	612 36 17 80 128
33	640 40 16 72 136
34	680 40 17 72 128
35	720 45 16 64 136
36	765 45 17 64 128
37	768 24 32 120 68
38	816 24 34 120 64
39	960 30 32 96 68
40	1020 30 34 96 64
41	1024 32 32 90 68
42	1088 32 34 90 64
43	1152 36 32 80 68
44	1224 36 34 80 64
45	1280 40 32 72 68
46	1360 40 34 72 64
47	1440 45 32 64 68
48	1530 45 34 64 64
49	1536 48 32 60 68
50	1632 48 34 60 64
51	1920 60 32 48 68
52	2040 60 34 48 64