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/15)
  if mod(Nx,i)==0, px=[px i]; end
end
for i=2:(Ny/15)
  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
         384          48           8          60         272
         384          96           4          30         544
         384         192           2          15        1088
         408           3         136         960          16
         408           6          68         480          32
         408          12          34         240          64
         408          24          17         120         128
         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

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% determine possible combinations of tile sizes
clear all, close all, Nx=5760; Ny=4352; 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         192         272
         510          30          17         192         256
         512          32          16         180         272
         544          32          17         180         256
         576          36          16         160         272
         612          36          17         160         256
         640          40          16         144         272
         680          40          17         144         256
         720          45          16         128         272
         765          45          17         128         256
         768          24          32         240         136
         816          24          34         240         128
         960          30          32         192         136
        1020          30          34         192         128
        1024          32          32         180         136
        1088          32          34         180         128
        1152          36          32         160         136
        1224          36          34         160         128
        1280          40          32         144         136
        1360          40          34         144         128
        1440          45          32         128         136
        1530          45          34         128         128
        1536          48          32         120         136
        1632          48          34         120         128
        1920          60          32          96         136
        2040          60          34          96         128
        2048          64          32          90         136
        2176          64          34          90         128
        2304          72          32          80         136
        2448          72          34          80         128
        2560          80          32          72         136
        2720          80          34          72         128
        2880          90          32          64         136
        3060          90          34          64         128
        3072          48          64         120          68
        3264          48          68         120          64
        3840          60          64          96          68
        4080          60          68          96          64
        4096          64          64          90          68
1	dimitri	1.1	% determine possible combinations of tile sizes
2			clear all, close all, Nx=2880; Ny=2176; px=[]; py=[];
3	dimitri	1.4	for i=2:(Nx/15)
4	dimitri	1.1	if mod(Nx,i)==0, px=[px i]; end
5			end
6	dimitri	1.4	for i=2:(Ny/15)
7	dimitri	1.1	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	dimitri	1.4	384 48 8 60 272
28			384 96 4 30 544
29			384 192 2 15 1088
30			408 3 136 960 16
31			408 6 68 480 32
32			408 12 34 240 64
33			408 24 17 120 128
34	dimitri	1.1	480 30 16 96 136
35			510 30 17 96 128
36			512 32 16 90 136
37			544 32 17 90 128
38			576 36 16 80 136
39			612 36 17 80 128
40			640 40 16 72 136
41			680 40 17 72 128
42			720 45 16 64 136
43			765 45 17 64 128
44			768 24 32 120 68
45			816 24 34 120 64
46			960 30 32 96 68
47			1020 30 34 96 64
48			1024 32 32 90 68
49	dimitri	1.2	1088 32 34 90 64
50			1152 36 32 80 68
51			1224 36 34 80 64
52			1280 40 32 72 68
53			1360 40 34 72 64
54			1440 45 32 64 68
55			1530 45 34 64 64
56			1536 48 32 60 68
57			1632 48 34 60 64
58			1920 60 32 48 68
59			2040 60 34 48 64
60	dimitri	1.3
61			%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
62			%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
63
64			% determine possible combinations of tile sizes
65			clear all, close all, Nx=5760; Ny=4352; px=[]; py=[];
66			for i=2:(Nx/30)
67			if mod(Nx,i)==0, px=[px i]; end
68			end
69			for i=2:(Ny/30)
70			if mod(Ny,i)==0, py=[py i]; end
71			end
72			n=0;
73			for i=1:length(px)
74			for j=1:length(py)
75			n=n+1;
76			pnx(n)=px(i);
77			snx(n)=Nx/pnx(n);
78			pny(n)=py(j);
79			sny(n)=Ny/pny(n);
80			pnxpny(n)=pnx(n)*pny(n);
81			end
82			end
83			[pnxpny i]=sort(pnxpny);
84			pnx=pnx(i); pny=pny(i);
85			snx=snx(i); sny=sny(i);
86			[pnxpny' pnx' pny' snx' sny']
87
88			npx*npy npx npy snx sny
89			480 30 16 192 272
90			510 30 17 192 256
91			512 32 16 180 272
92			544 32 17 180 256
93			576 36 16 160 272
94			612 36 17 160 256
95			640 40 16 144 272
96			680 40 17 144 256
97			720 45 16 128 272
98			765 45 17 128 256
99			768 24 32 240 136
100			816 24 34 240 128
101			960 30 32 192 136
102			1020 30 34 192 128
103			1024 32 32 180 136
104			1088 32 34 180 128
105			1152 36 32 160 136
106			1224 36 34 160 128
107			1280 40 32 144 136
108			1360 40 34 144 128
109			1440 45 32 128 136
110			1530 45 34 128 128
111			1536 48 32 120 136
112			1632 48 34 120 128
113			1920 60 32 96 136
114			2040 60 34 96 128
115			2048 64 32 90 136
116			2176 64 34 90 128
117			2304 72 32 80 136
118			2448 72 34 80 128
119			2560 80 32 72 136
120			2720 80 34 72 128
121			2880 90 32 64 136
122			3060 90 34 64 128
123			3072 48 64 120 68
124			3264 48 68 120 64
125			3840 60 64 96 68
126			4080 60 68 96 64
127			4096 64 64 90 68