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	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