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

Diff of /MITgcm_contrib/quarter_degree_global/size.m

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

revision 1.2 by dimitri, Wed Dec 15 16:45:03 2004 UTC revision 1.4 by dimitri, Sun Oct 15 17:05:05 2006 UTC
# Line 1  Line 1 
1  % determine possible combinations of tile sizes  % determine possible combinations of tile sizes
2  clear all, close all, Nx=2880; Ny=2176; px=[]; py=[];  clear all, close all, Nx=2880; Ny=2176; px=[]; py=[];
3  for i=2:(Nx/30)  for i=2:(Nx/15)
4    if mod(Nx,i)==0, px=[px i]; end    if mod(Nx,i)==0, px=[px i]; end
5  end  end
6  for i=2:(Ny/30)  for i=2:(Ny/15)
7    if mod(Ny,i)==0, py=[py i]; end    if mod(Ny,i)==0, py=[py i]; end
8  end  end
9  n=0;  n=0;
# Line 24  snx=snx(i); sny=sny(i); Line 24  snx=snx(i); sny=sny(i);
24    
25    
26       npx*npy         npx         npy         snx         sny       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           480          30          16          96         136
35           510          30          17          96         128           510          30          17          96         128
36           512          32          16          90         136           512          32          16          90         136
# Line 50  snx=snx(i); sny=sny(i); Line 57  snx=snx(i); sny=sny(i);
57          1632          48          34          60          64          1632          48          34          60          64
58          1920          60          32          48          68          1920          60          32          48          68
59          2040          60          34          48          64          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
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22