/[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.4 - (show annotations) (download)
Sun Oct 15 17:05:05 2006 UTC (17 years, 6 months ago) by dimitri
Branch: MAIN
CVS Tags: HEAD
Changes since 1.3: +9 -2 lines
adding some old stuff related to eight/BULK integration

1 % determine possible combinations of tile sizes
2 clear all, close all, Nx=2880; Ny=2176; px=[]; py=[];
3 for i=2:(Nx/15)
4 if mod(Nx,i)==0, px=[px i]; end
5 end
6 for i=2:(Ny/15)
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 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 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 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
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

  ViewVC Help
Powered by ViewVC 1.1.22