/[MITgcm]/MITgcm_contrib/quarter_degree_global/size.m
ViewVC logotype

Contents of /MITgcm_contrib/quarter_degree_global/size.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.3 - (show annotations) (download)
Thu Mar 10 01:45:01 2005 UTC (20 years, 4 months ago) by dimitri
Branch: MAIN
Changes since 1.2: +68 -0 lines
quarter_degree_global/size.m

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

  ViewVC Help
Powered by ViewVC 1.1.22