/[MITgcm]/MITgcm_contrib/high_res_cube/README_ice
ViewVC logotype

Annotation of /MITgcm_contrib/high_res_cube/README_ice

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


Revision 1.41 - (hide annotations) (download)
Fri Jul 9 23:54:52 2004 UTC (19 years, 11 months ago) by dimitri
Branch: MAIN
Changes since 1.40: +55 -30 lines
updated for checkpoint54b_post

1 dimitri 1.1 Instructions for setting up a cube-sphere integration with sea-ice.
2 dimitri 1.41 Date last tested: July 9, 2004.
3     Code and input used for integrating cube9.
4 dimitri 1.8
5    
6     ##########################################################################
7     ##########################################################################
8 dimitri 1.14 # getting the code from anonymous CVS server
9 dimitri 1.1
10 dimitri 1.38 cvs co MITgcm_contrib/high_res_cube/README_ice
11     cvs co MITgcm_contrib/high_res_cube/code-mods
12     cvs co MITgcm_contrib/high_res_cube/input
13     cvs co MITgcm_contrib/high_res_cube/results
14     cvs co MITgcm_contrib/edge_problem
15 dimitri 1.41 cvs co -r checkpoint54b_post MITgcm_code
16     cvs co -r checkpoint54b_post MITgcm/verification/global_ocean.cs32x15
17     cvs co -r checkpoint54b_post MITgcm/verification/testreport
18 dimitri 1.14
19    
20     ##########################################################################
21     ##########################################################################
22     # 32*32*6*15, 1-cpu cube sphere
23 dimitri 1.3
24 dimitri 1.1 cd MITgcm/verification/global_ocean.cs32x15/code
25 dimitri 1.8 \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
26 dimitri 1.33 \cp ../../../utils/exch2/code-mods/s12t_16x32/* .
27 dimitri 1.37 cd ../build
28     \rm *
29     ../../../tools/genmake2 -mods=../code
30     make depend
31     make
32     cd ..
33     mkdir run
34     cd run
35     \rm *
36     \cp ../input/* .
37     \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
38     ../build/mitgcmuv >& output.txt
39 dimitri 1.3
40 dimitri 1.37 comparison output is in:
41     ../../../../MITgcm_contrib/high_res_cube/results/output.txt
42    
43     to use matlab to look at the output
44     matlab
45    
46     for fld={'Eta','VICE','UICE','HEFF','AREA','Qsw','Qnet','FV','FU','EmPmR'}
47     tmp=permute(readbin([fld{1} '.0000000020.data'],[32,6,32]),[1 3 2]);
48     cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
49     end
50    
51     heff=permute(readbin(['HEFF.0000000020.data'],[32,6,32]),[1 3 2]);
52     mask=0*heff; mask(find(heff>.001))=1;
53     uice=mask.*permute(readbin(['UICE.0000000020.data'],[32,6,32]),[1 3 2]);
54     vice=mask.*permute(readbin(['VICE.0000000020.data'],[32,6,32]),[1 3 2]);
55     clf
56     subplot(221), mypcolor(uice(:,:,3)'), colorbar, title('uice, tile 3')
57     subplot(222), mypcolor(uice(:,:,6)'), colorbar, title('uice, tile 6')
58     subplot(223), mypcolor(vice(:,:,3)'), colorbar, title('vice, tile 3')
59     subplot(224), mypcolor(vice(:,:,6)'), colorbar, title('vice, tile 6')
60    
61    
62     ##########################################################################
63     ##########################################################################
64     # 32*32*6*15, 17-tile, 1-cpu cube sphere experiment
65    
66     cd MITgcm/verification/global_ocean.cs32x15/code
67     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
68     \cp ../code_alt/code.176t_8x4/* .
69 dimitri 1.25 cd ../build
70 dimitri 1.37 \rm *
71 dimitri 1.10 ../../../tools/genmake2 -mods=../code
72 dimitri 1.1 make depend
73     make
74     cd ..
75     mkdir run
76     cd run
77 dimitri 1.37 \rm *
78 dimitri 1.8 \cp ../input/* .
79     \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
80 dimitri 1.3 ../build/mitgcmuv >& output.txt
81    
82     comparison output is in:
83     ../../../../MITgcm_contrib/high_res_cube/results/output.txt
84    
85     to use matlab to look at the output
86     matlab
87 dimitri 1.17
88 dimitri 1.5 for fld={'Eta','VICE','UICE','HEFF','AREA','Qsw','Qnet','FV','FU','EmPmR'}
89 dimitri 1.13 tmp=permute(readbin([fld{1} '.0000000020.data'],[32,6,32]),[1 3 2]);
90 dimitri 1.4 cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
91     end
92 dimitri 1.8
93 dimitri 1.17 heff=permute(readbin(['HEFF.0000000020.data'],[32,6,32]),[1 3 2]);
94     mask=0*heff; mask(find(heff>.001))=1;
95     uice=mask.*permute(readbin(['UICE.0000000020.data'],[32,6,32]),[1 3 2]);
96     vice=mask.*permute(readbin(['VICE.0000000020.data'],[32,6,32]),[1 3 2]);
97     clf
98     subplot(221), mypcolor(uice(:,:,3)'), colorbar, title('uice, tile 3')
99     subplot(222), mypcolor(uice(:,:,6)'), colorbar, title('uice, tile 6')
100     subplot(223), mypcolor(vice(:,:,3)'), colorbar, title('vice, tile 3')
101     subplot(224), mypcolor(vice(:,:,6)'), colorbar, title('vice, tile 6')
102    
103 dimitri 1.8
104     ##########################################################################
105     ##########################################################################
106     # 32*32*6*15, 2-cpu cube sphere on the altices with efc
107    
108     cd MITgcm/verification/global_ocean.cs32x15/code
109     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
110 dimitri 1.33 \cp ../../../utils/exch2/code-mods/s12t_16x32/* .
111 dimitri 1.8 \cp SIZE.h_mpi SIZE.h
112    
113 dimitri 1.25 cd ../build
114 dimitri 1.14 \rm *
115 dimitri 1.41 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi_altix
116 dimitri 1.8 make depend
117 dimitri 1.41 make -j
118 dimitri 1.8
119     cd ..
120     mkdir run
121     cd run
122 dimitri 1.14 \rm *
123 dimitri 1.8 \cp ../input/* .
124     \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
125     mpirun -np 2 ../build/mitgcmuv
126    
127    
128     ##########################################################################
129     ##########################################################################
130 dimitri 1.35 # 510*510*6*50, 216-tile, 54-cpu cube sphere on the altices with efc
131 dimitri 1.14 # requires input files in sudirectories run_template and ncep_rgau
132     # presently residing under orion:/tmp1/dmenem/cube
133 dimitri 1.8
134     cd MITgcm/verification/global_ocean.cs32x15/code
135     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
136     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* .
137     \cp SIZE.h_54 SIZE.h
138 dimitri 1.11 \cp packages.conf_hr packages.conf
139 dimitri 1.23
140 dimitri 1.25 cd ../build
141 dimitri 1.23 \rm *
142 dimitri 1.41 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi_altix
143 dimitri 1.23 make depend
144 dimitri 1.41 make -j
145 dimitri 1.23
146     cd ..
147     mkdir run
148     cd run
149     \rm *
150     \cp ../input/* .
151     \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
152     qsub job_orion
153    
154    
155     ##########################################################################
156     ##########################################################################
157 dimitri 1.35 # 510*510*6*50, 54-tile, 54-cpu cube sphere on the altices with efc
158 dimitri 1.36 # starting with pickup file from end of 1992
159 dimitri 1.35 # requires input files in sudirectories run_template and ncep_rgau
160    
161     cd MITgcm/verification/global_ocean.cs32x15/code
162     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
163     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s54t_170x170/* .
164     \cp packages.conf_hr packages.conf
165    
166     cd ../build
167     \rm *
168 dimitri 1.40 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi_altix
169 dimitri 1.35 make depend
170 dimitri 1.41 make -j
171 dimitri 1.35
172     cd ..
173     mkdir run
174     cd run
175     \rm *
176     \cp ../input/* .
177     \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
178 dimitri 1.41 qsub job_altix92_54t
179 dimitri 1.35
180    
181     ##########################################################################
182     ##########################################################################
183 dimitri 1.41 # 510*510*6*50, 216-cpu cube sphere on altix with efc
184 dimitri 1.23 # requires input files in sudirectories run_template and ncep_rgau
185 dimitri 1.41 # these can be obained from lou:/u/menemenl/cube
186 dimitri 1.23
187     cd MITgcm/verification/global_ocean.cs32x15/code
188     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
189 dimitri 1.41 \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* .
190     \cp SIZE.h_216 SIZE.h
191 dimitri 1.23 \cp packages.conf_hr packages.conf
192 dimitri 1.8
193 dimitri 1.25 cd ../build
194 dimitri 1.14 \rm *
195 dimitri 1.41 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi_altix
196 dimitri 1.8 make depend
197 dimitri 1.41 make -j
198 dimitri 1.8
199     cd ..
200     mkdir run
201     cd run
202 dimitri 1.14 \rm *
203 dimitri 1.8 \cp ../input/* .
204     \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
205 dimitri 1.41 qsub job_altix92
206    
207    
208    
209    
210    
211    
212    
213    
214    
215    
216    
217    
218    
219    
220    
221    
222    
223    
224    
225    
226    
227    
228    
229    
230    
231 dimitri 1.35
232    
233     ##########################################################################
234     ##########################################################################
235 dimitri 1.41 # 510*510*6*50, 205-tile, 41-cpu cube sphere on the altices with efc
236 dimitri 1.35 # requires input files in sudirectories run_template and ncep_rgau
237     # presently residing under orion:/tmp1/dmenem/cube
238    
239     cd MITgcm/verification/global_ocean.cs32x15/code
240     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
241 dimitri 1.41 \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s205t_85x85/* .
242     \cp SIZE.h_41 SIZE.h
243 dimitri 1.35 \cp packages.conf_hr packages.conf
244    
245     cd ../build
246     \rm *
247 dimitri 1.41 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi_altix
248 dimitri 1.35 make depend
249 dimitri 1.41 make -j
250 dimitri 1.35
251     cd ..
252     mkdir run
253     cd run
254     \rm *
255     \cp ../input/* .
256     \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
257 dimitri 1.41 qsub job_orion41
258 dimitri 1.9
259    
260     ##########################################################################
261     ##########################################################################
262 dimitri 1.41 # 510*510*6*50, 1500-tile, 50-cpu cube sphere on the altices with efc
263 dimitri 1.17 # requires input files in sudirectories run_template and ncep_rgau
264 dimitri 1.41 # presently residing under orion:/tmp1/dmenem/cube
265 dimitri 1.9
266     cd MITgcm/verification/global_ocean.cs32x15/code
267     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
268 dimitri 1.41 \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s1500t_17x51/* .
269     \cp SIZE.h_50 SIZE.h
270 dimitri 1.12 \cp packages.conf_hr packages.conf
271 dimitri 1.9
272 dimitri 1.25 cd ../build
273 dimitri 1.15 \rm *
274 dimitri 1.13 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi_altix
275 dimitri 1.9 make depend
276 dimitri 1.41 make -j
277 dimitri 1.9
278     cd ..
279     mkdir run
280     cd run
281 dimitri 1.15 \rm *
282     \cp ../input/* .
283     \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
284 dimitri 1.41 qsub job_orion
285 dimitri 1.17
286    
287     ##########################################################################
288     ##########################################################################
289     # matlab script for generating pChkptFreq
290    
291     nTimeSteps=26352;
292 dimitri 1.34 for niter0=216:26352:(26352*12)
293 dimitri 1.17 pChkptFreq=(niter0+nTimeSteps)*1200;
294     day=pChkptFreq/24/60/60;
295     disp(['niter0=' int2str(niter0) ',day=' int2str(day) ', ' ...
296     datestr(datenum(1992,1,1)+day) ', pChkptFreq=' int2str(pChkptFreq)])
297     end
298    
299     niter0 = 216 , day = 369 , 04-Jan-1993, pChkptFreq = 31881600
300     niter0 = 26568 , day = 735 , 05-Jan-1994, pChkptFreq = 63504000
301     niter0 = 52920 , day = 1101, 06-Jan-1995, pChkptFreq = 95126400
302     niter0 = 79272 , day = 1467, 07-Jan-1996, pChkptFreq = 126748800
303     niter0 = 105624, day = 1833, 07-Jan-1997, pChkptFreq = 158371200
304     niter0 = 131976, day = 2199, 08-Jan-1998, pChkptFreq = 189993600
305     niter0 = 158328, day = 2565, 09-Jan-1999, pChkptFreq = 221616000
306     niter0 = 184680, day = 2931, 10-Jan-2000, pChkptFreq = 253238400
307     niter0 = 211032, day = 3297, 10-Jan-2001, pChkptFreq = 284860800
308     niter0 = 237384, day = 3663, 11-Jan-2002, pChkptFreq = 316483200
309     niter0 = 263736, day = 4029, 12-Jan-2003, pChkptFreq = 348105600
310 dimitri 1.34 niter0 = 290088, day = 4395, 13-Jan-2004, pChkptFreq = 379728000

  ViewVC Help
Powered by ViewVC 1.1.22