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

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

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