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

Contents of /MITgcm_contrib/high_res_cube/README_ice

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


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

1 Instructions for setting up a cube-sphere integration with sea-ice.
2 Date last tested: July 9, 2004.
3 Code and input used for integrating cube9.
4
5
6 ##########################################################################
7 ##########################################################################
8 # getting the code from anonymous CVS server
9
10 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 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
19
20 ##########################################################################
21 ##########################################################################
22 # 32*32*6*15, 1-cpu cube sphere
23
24 cd MITgcm/verification/global_ocean.cs32x15/code
25 \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
26 \cp ../../../utils/exch2/code-mods/s12t_16x32/* .
27 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
40 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 cd ../build
70 \rm *
71 ../../../tools/genmake2 -mods=../code
72 make depend
73 make
74 cd ..
75 mkdir run
76 cd run
77 \rm *
78 \cp ../input/* .
79 \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
80 ../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
88 for fld={'Eta','VICE','UICE','HEFF','AREA','Qsw','Qnet','FV','FU','EmPmR'}
89 tmp=permute(readbin([fld{1} '.0000000020.data'],[32,6,32]),[1 3 2]);
90 cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
91 end
92
93 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
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 \cp ../../../utils/exch2/code-mods/s12t_16x32/* .
111 \cp SIZE.h_mpi SIZE.h
112
113 cd ../build
114 \rm *
115 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi_altix
116 make depend
117 make -j
118
119 cd ..
120 mkdir run
121 cd run
122 \rm *
123 \cp ../input/* .
124 \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
125 mpirun -np 2 ../build/mitgcmuv
126
127
128 ##########################################################################
129 ##########################################################################
130 # 510*510*6*50, 216-tile, 54-cpu cube sphere on the altices with efc
131 # requires input files in sudirectories run_template and ncep_rgau
132 # presently residing under orion:/tmp1/dmenem/cube
133
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 \cp packages.conf_hr packages.conf
139
140 cd ../build
141 \rm *
142 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi_altix
143 make depend
144 make -j
145
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 # 510*510*6*50, 54-tile, 54-cpu cube sphere on the altices with efc
158 # starting with pickup file from end of 1992
159 # 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 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi_altix
169 make depend
170 make -j
171
172 cd ..
173 mkdir run
174 cd run
175 \rm *
176 \cp ../input/* .
177 \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
178 qsub job_altix92_54t
179
180
181 ##########################################################################
182 ##########################################################################
183 # 510*510*6*50, 216-cpu cube sphere on altix with efc
184 # requires input files in sudirectories run_template and ncep_rgau
185 # these can be obained from lou:/u/menemenl/cube
186
187 cd MITgcm/verification/global_ocean.cs32x15/code
188 \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
189 \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* .
190 \cp SIZE.h_216 SIZE.h
191 \cp packages.conf_hr packages.conf
192
193 cd ../build
194 \rm *
195 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi_altix
196 make depend
197 make -j
198
199 cd ..
200 mkdir run
201 cd run
202 \rm *
203 \cp ../input/* .
204 \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
205 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
232
233 ##########################################################################
234 ##########################################################################
235 # 510*510*6*50, 205-tile, 41-cpu cube sphere on the altices with efc
236 # 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 \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s205t_85x85/* .
242 \cp SIZE.h_41 SIZE.h
243 \cp packages.conf_hr packages.conf
244
245 cd ../build
246 \rm *
247 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi_altix
248 make depend
249 make -j
250
251 cd ..
252 mkdir run
253 cd run
254 \rm *
255 \cp ../input/* .
256 \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
257 qsub job_orion41
258
259
260 ##########################################################################
261 ##########################################################################
262 # 510*510*6*50, 1500-tile, 50-cpu cube sphere on the altices with efc
263 # requires input files in sudirectories run_template and ncep_rgau
264 # presently residing under orion:/tmp1/dmenem/cube
265
266 cd MITgcm/verification/global_ocean.cs32x15/code
267 \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
268 \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s1500t_17x51/* .
269 \cp SIZE.h_50 SIZE.h
270 \cp packages.conf_hr packages.conf
271
272 cd ../build
273 \rm *
274 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi_altix
275 make depend
276 make -j
277
278 cd ..
279 mkdir run
280 cd run
281 \rm *
282 \cp ../input/* .
283 \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
284 qsub job_orion
285
286
287 ##########################################################################
288 ##########################################################################
289 # matlab script for generating pChkptFreq
290
291 nTimeSteps=26352;
292 for niter0=216:26352:(26352*12)
293 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 niter0 = 290088, day = 4395, 13-Jan-2004, pChkptFreq = 379728000

  ViewVC Help
Powered by ViewVC 1.1.22