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

Annotation of /MITgcm_contrib/quarter_degree_global/size.m

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


Revision 1.4 - (hide annotations) (download)
Sun Oct 15 17:05:05 2006 UTC (18 years, 9 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 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

  ViewVC Help
Powered by ViewVC 1.1.22