/[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.3 - (hide 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 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

  ViewVC Help
Powered by ViewVC 1.1.22