1 |
Instructions for setting up a quarter degree integration |
2 |
forced 1. by NCEP fields and 2. by SIO fields. |
3 |
|
4 |
########################################################################## |
5 |
# getting the code from anonymous CVS server |
6 |
|
7 |
cvs co MITgcm_contrib/quarter_degree_global |
8 |
cvs co MITgcm_code |
9 |
cvs co MITgcm/verification/global_with_exf |
10 |
cvs co MITgcm/verification/testreport |
11 |
|
12 |
|
13 |
############################################################################# |
14 |
# 1-CPU verification/global_with_exf including pkg/diagnostics |
15 |
|
16 |
cd MITgcm/verification/global_with_exf/build |
17 |
\rm * |
18 |
\cp ../code/* . |
19 |
\cp ../../../../MITgcm_contrib/quarter_degree_global/code/packages.conf . |
20 |
\cp ../../../../MITgcm_contrib/quarter_degree_global/code/ECCO_CPPOPTIONS.h . |
21 |
\cp ../../../../MITgcm_contrib/quarter_degree_global/code/CPP_OPTIONS.h . |
22 |
\cp ../../../../MITgcm_contrib/quarter_degree_global/code/DIAGNOSTICS_SIZE.h . |
23 |
\cp ../../../../MITgcm_contrib/quarter_degree_global/code/SIZE.h_test1 . |
24 |
\mv SIZE.h_test1 SIZE.h |
25 |
../../../tools/genmake2 |
26 |
make depend |
27 |
make -j |
28 |
cd .. |
29 |
mkdir run |
30 |
cd run |
31 |
\rm * |
32 |
\cp ../build/mitgcmuv . |
33 |
\ln -sf ../input/* . |
34 |
\rm data |
35 |
\cp ../../../../MITgcm_contrib/quarter_degree_global/input/data_test data |
36 |
\cp ../../../../MITgcm_contrib/quarter_degree_global/input/data.pkg . |
37 |
\cp ../../../../MITgcm_contrib/quarter_degree_global/input/data.kpp . |
38 |
\cp ../../../../MITgcm_contrib/quarter_degree_global/input/data.diagnostics . |
39 |
mitgcmuv >& output.txt |
40 |
|
41 |
|
42 |
there is comparison output in |
43 |
../../../../MITgcm_contrib/quarter_degree_global/results/output.txt |
44 |
|
45 |
|
46 |
############################################################################# |
47 |
# mpi test verification/global_with_exf including pkg/diagnostics on columbia |
48 |
|
49 |
cd MITgcm/verification/global_with_exf/build |
50 |
\rm * |
51 |
\cp ../code/* . |
52 |
\cp ../../../../MITgcm_contrib/quarter_degree_global/code/packages.conf . |
53 |
\cp ../../../../MITgcm_contrib/quarter_degree_global/code/ECCO_CPPOPTIONS.h . |
54 |
\cp ../../../../MITgcm_contrib/quarter_degree_global/code/CPP_OPTIONS.h . |
55 |
\cp ../../../../MITgcm_contrib/quarter_degree_global/code/SIZE.h_test . |
56 |
\mv SIZE.h_test SIZE.h |
57 |
../../../tools/genmake2 -of ../../../tools/build_options/linux_ia64_ifort+mpi_altix_nas |
58 |
make depend |
59 |
make -j |
60 |
cd .. |
61 |
mkdir run |
62 |
cd run |
63 |
\rm * |
64 |
\cp ../build/mitgcmuv . |
65 |
\ln -sf ../input/* . |
66 |
\rm data |
67 |
\cp ../../../../MITgcm_contrib/quarter_degree_global/input/data_test data |
68 |
mpirun -np 2 mitgcmuv |
69 |
|
70 |
there is comparison output in |
71 |
../../../../MITgcm_contrib/quarter_degree_global/results/STDOUT.0000 |
72 |
|
73 |
|
74 |
########################################################################## |
75 |
# 1440x1088 204-cpu (120x12)x(64x17) 1/4-deg integration on altix with efc |
76 |
# requires input files in subdirectory iter69 to be at same level as MITgcm |
77 |
# and MITgcm_contrib. A copy of iter69 resides on columbia10:/nobackup10 |
78 |
# The integration length is 4018 days, 1992-2002. |
79 |
# First three days need to be integrated using deltat=200. |
80 |
|
81 |
cd MITgcm |
82 |
mkdir bin exe |
83 |
cd bin |
84 |
\rm * |
85 |
\cp ../../MITgcm_contrib/quarter_degree_global/code/* . |
86 |
../tools/genmake2 -of ../tools/build_options/linux_ia64_ifort+mpi_altix_nas |
87 |
make depend |
88 |
make -j |
89 |
cd ../exe |
90 |
\rm * |
91 |
\cp ../bin/mitgcmuv . |
92 |
\cp ../../MITgcm_contrib/quarter_degree_global/input/* . |
93 |
\ln -sf ../../iter69/* . |
94 |
qsub job_altix |
95 |
|
96 |
|
97 |
########################################################################## |
98 |
# 2880x2176 480-cpu (96x30)x(136x16) 1/4-deg integration on altix with efc |
99 |
# requires input files in sudirectory iter69 to be at same level as MITgcm |
100 |
# and MITgcm_contrib. A copy of iter69 with 1/8-deg files resides on |
101 |
# columbia18:/nobackup18/menemenl/eight. |
102 |
# The integration length is 3 days. |
103 |
# First three days need to be integrated using deltat=100. |
104 |
|
105 |
cd MITgcm |
106 |
mkdir bin exe |
107 |
cd bin |
108 |
\rm * |
109 |
\cp ../../MITgcm_contrib/quarter_degree_global/code/* . |
110 |
\cp SIZE.h_8_480 SIZE.h |
111 |
../tools/genmake2 -of ../tools/build_options/linux_ia64_efc+mpi_altix |
112 |
make depend |
113 |
make -j |
114 |
cd ../exe |
115 |
\rm * |
116 |
\cp ../bin/mitgcmuv . |
117 |
\cp ../../MITgcm_contrib/quarter_degree_global/input/* . |
118 |
\ln -sf ../../iter69/* . |
119 |
\cp data.exf_clim_8 data.exf_clim |
120 |
\cp data_SIO_8 data |
121 |
qsub job_altix |
122 |
|
123 |
|
124 |
########################################################################## |
125 |
# 2880x2176 960-cpu (96x30)x(68x32) 1/4-deg integration on altix with efc |
126 |
# requires input files in sudirectory iter69 to be at same level as MITgcm |
127 |
# and MITgcm_contrib. A copy of iter69 with 1/8-deg files resides on |
128 |
# columbia18:/nobackup18/menemenl/eight. |
129 |
# The integration length is 3 days. |
130 |
# First three days need to be integrated using deltat=100. |
131 |
|
132 |
module purge |
133 |
module load modules scsl.1.5.0.0 intel-comp.8.0.070 mpt.1.11-85 |
134 |
setenv MPI_DSM_DISTRIBUTE |
135 |
cd MITgcm |
136 |
mkdir bin exe |
137 |
cd bin |
138 |
\rm * |
139 |
\cp ../../MITgcm_contrib/quarter_degree_global/code/* . |
140 |
\cp SIZE.h_8_960 SIZE.h |
141 |
../tools/genmake2 -of ../tools/build_options/linux_ia64_ifort+mpi_altix_nas |
142 |
make depend |
143 |
make -j 4 |
144 |
|
145 |
cd ../exe |
146 |
\cp ../bin/mitgcmuv . |
147 |
\cp ../../MITgcm_contrib/quarter_degree_global/input/* . |
148 |
\ln -sf ../../iter69/* . |
149 |
\cp data.exf_clim_8 data.exf_clim |
150 |
\cp data_SIO_8 data |
151 |
modpickup_8 |
152 |
|
153 |
/opt/sgi/mpt/1.11-85/bin/mpirun.real -v -a a2k \ |
154 |
-d /nobackup/menemenl/eighth/MITgcm/exe \ |
155 |
c19 480 runon 11-490 mitgcmuv960 : \ |
156 |
c18 480 runon 11-490 mitgcmuv960 \ |
157 |
>&! output.run & |
158 |
|
159 |
|
160 |
|
161 |
/opt/sgi/mpt/1.11-85/bin/mpirun.real -v -a a2k \ |
162 |
-d /nobackup/menemenl/eighth/MITgcm/exe \ |
163 |
c17 480 runon 11-490 mitgcmuv1920 : \ |
164 |
c18 480 runon 11-490 mitgcmuv1920 : \ |
165 |
c19 480 runon 11-490 mitgcmuv1920 : \ |
166 |
c20 480 runon 11-490 mitgcmuv1920 \ |
167 |
>&! output.run & |
168 |
|
169 |
|
170 |
|
171 |
|
172 |
|
173 |
|
174 |
fails for 1/8 |
175 |
> viscA4GridMin=3.e-4, |
176 |
> viscA4GridMax=4.e-3, |
177 |
> viscC4Leith=0.05, |
178 |
> deltat=300 |
179 |
fails: |
180 |
< viscA4GridMin=6.e-4, |
181 |
< viscA4GridMax=8.e-3, |
182 |
< viscC4Leith=0.1, |
183 |
fails: |
184 |
< viscA4GridMin=6.e-4, |
185 |
< viscA4GridMax=8.e-3, |
186 |
< viscC4Leith=0.2, |
187 |
fails: |
188 |
< viscA4GridMin=1.e-3, |
189 |
< viscA4GridMax=1.e-2, |
190 |
< viscC4Leith=0.2, |
191 |
fails: |
192 |
< viscA4GridMin=1.e-3, |
193 |
< viscA4GridMax=1.e-2, |
194 |
< viscC4Leith=0.3, |
195 |
fails: |
196 |
< viscA4GridMin=1.2e-3, |
197 |
< viscA4GridMax=1.6e-2, |
198 |
< viscC4Leith=0.2, |
199 |
fails: |
200 |
< viscA4GridMin=1.2e-3, |
201 |
< viscA4GridMax=1.6e-2, |
202 |
< viscC4Leith=0.4, |
203 |
fails |
204 |
< viscA4GridMin=1.2e-3, |
205 |
< viscA4GridMax=1.2e-2, |
206 |
< viscC4Leith=1.2, |
207 |
fails: |
208 |
< viscA4GridMin=2e-3, |
209 |
< viscA4GridMax=1e-2, |
210 |
< viscC4Leith=0.1, |
211 |
fails: |
212 |
< viscA4GridMin=2e-3, |
213 |
< viscA4GridMax=1e-2, |
214 |
< viscC4Leith=0.2, |
215 |
fails: |
216 |
< viscA4GridMin=4e-3, |
217 |
< viscA4GridMax=8e-3, |
218 |
< viscC4Leith=0.2, |
219 |
fails: |
220 |
viscA4GridMin=5e-3, |
221 |
viscA4GridMax=1e-2, |
222 |
viscC4Leith=.2, |
223 |
starting with time step 50112, try: |
224 |
(this mean day 174 becomes day 116) |
225 |
viscA4GridMin=6.e-4, |
226 |
viscA4GridMax=8.e-3, |
227 |
viscC4Leith=0.1, |
228 |
deltat=200 |
229 |
|
230 |
|
231 |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
232 |
|
233 |
1/8 timing tests: |
234 |
|
235 |
columbia20 |
236 |
scsl.1.5.0.0 |
237 |
c18 intel-comp.8.0.070 |
238 |
240-CPU default mpt.1.11-85 |
239 |
------------------------------------------------ |
240 |
MAIN LOOP 7115 5066 |
241 |
EXF_GETFORCING 372 38 |
242 |
DO_OCEANIC_PHYS 1221 1272 |
243 |
STATEKPP 509 532 |
244 |
KPPMIX 370 374 |
245 |
DYNAMICS 1527 1429 |
246 |
SOLVE_FOR_PRESSURE 1804 206 |
247 |
BLOCKING_EXCHANGES 117 115 |
248 |
THERMODYNAMICS 1166 1100 |
249 |
DO_THE_MODEL_IO 173 177 |
250 |
WRITE_CHECKPOINT 341 333 |
251 |
|
252 |
|
253 |
columbia20 |
254 |
scsl.1.5.0.0 |
255 |
c18 intel-comp.8.0.070 |
256 |
480-CPU default mpt.1.11-85 |
257 |
c19 default |
258 |
c17 default, other processes |
259 |
c19 new |
260 |
------------------------------------------------------------------------- |
261 |
MAIN LOOP 4024 2917 5379 8895 3571 |
262 |
EXF_GETFORCING 218 33 880 261 227 |
263 |
DO_OCEANIC_PHYS 622 619 621 913 620 |
264 |
STATEKPP 269 265 268 310 265 |
265 |
KPPMIX 190 188 189 291 189 |
266 |
DYNAMICS 747 718 745 1456 714 |
267 |
SOLVE_FOR_PRESSURE 966 119 1703 2748 900 |
268 |
BLOCKING_EXCHANGES 64 70 392 375 70 |
269 |
THERMODYNAMICS 553 582 566 1510 584 |
270 |
DO_THE_MODEL_IO 178 156 149 184 152 |
271 |
WRITE_CHECKPOINT 364 332 864 420 333 |
272 |
|
273 |
|
274 |
|
275 |
CPU 240 480 960 1440 1920 |
276 |
------------------------------------------------------------- |
277 |
MAIN LOOP 5066 2917 2227 1551 1276 |
278 |
EXF_GETFORCING 38 33 126 60 83 |
279 |
DO_OCEANIC_PHYS 1272 619 320 390 186 |
280 |
STATEKPP 532 265 138 181 77 |
281 |
KPPMIX 374 188 98 112 53 |
282 |
DYNAMICS 1429 718 313 319 173 |
283 |
SOLVE_FOR_PRESSURE 206 119 483 169 310 |
284 |
BLOCKING_EXCHANGES 115 70 53 29 28 |
285 |
THERMODYNAMICS 1100 582 294 279 172 |
286 |
DO_THE_MODEL_IO 177 156 131 202 155 |
287 |
WRITE_CHECKPOINT 333 332 302 524 352 |
288 |
|
289 |
|
290 |
4 clusters timeave |
291 |
Dec 16 08:17 Ttave.0000051840.data |
292 |
Dec 16 08:47 Ttave.0000053136.data |
293 |
Dec 16 09:17 Ttave.0000054432.data |
294 |
Dec 16 09:53 Ttave.0000055728.data |
295 |
Dec 17 14:08 Ttave.0000075168.data |
296 |
Dec 17 14:33 Ttave.0000076464.data |
297 |
Dec 17 14:58 Ttave.0000077760.data |
298 |
Dec 17 15:29 Ttave.0000079056.data |
299 |
|
300 |
1 cluster timeave |
301 |
Dec 16 22:14 Ttave.0000059616.data |
302 |
Dec 16 23:16 Ttave.0000060912.data |
303 |
Dec 17 00:19 Ttave.0000062208.data |
304 |
Dec 17 01:27 Ttave.0000063504.data |
305 |
Dec 17 02:29 Ttave.0000064800.data |
306 |
Dec 17 03:31 Ttave.0000066096.data |
307 |
Dec 17 04:39 Ttave.0000067392.data |
308 |
Dec 17 05:42 Ttave.0000068688.data |
309 |
Dec 17 06:44 Ttave.0000069984.data |
310 |
|
311 |
2 cluster timeave |
312 |
Dec 17 21:19 Ttave.0000086832.data |
313 |
Dec 17 21:55 Ttave.0000088128.data 36 |
314 |
Dec 17 22:30 Ttave.0000089424.data 35 |
315 |
Dec 17 23:12 Ttave.0000090720.data 42 |
316 |
Dec 17 23:47 Ttave.0000092016.data 35 |
317 |
Dec 18 00:23 Ttave.0000093312.data 36 |
318 |
Dec 18 01:04 Ttave.0000094608.data 41 |
319 |
Dec 18 01:39 Ttave.0000095904.data 35 |
320 |
|
321 |
2 cluster diagnostics |
322 |
Dec 18 02:26 THETA.0000086832.data |
323 |
Dec 18 03:05 THETA.0000088128.data 39 |
324 |
Dec 18 03:43 THETA.0000089424.data 38 |
325 |
Dec 18 04:27 THETA.0000090720.data 44 |
326 |
Dec 18 05:06 THETA.0000092016.data 39 |
327 |
Dec 18 05:45 THETA.0000093312.data 39 |
328 |
Dec 18 06:29 THETA.0000094608.data 44 |
329 |
Dec 18 07:08 THETA.0000095904.data 39 |
330 |
|
331 |
|
332 |
%%%%%%%%%%%%%%%%% |
333 |
|
334 |
|
335 |
|
336 |
e1=readbin('test_oldstyle/ETAtave.0000000018.data',[90 40]); |
337 |
e2=readbin('ETAN.0000000018.data',[90 40]); |
338 |
e16=readbin('test_oldstyle/Eta.0000000016.data',[90 40]); |
339 |
e17=readbin('test_oldstyle/Eta.0000000017.data',[90 40]); |
340 |
e18=readbin('test_oldstyle/Eta.0000000018.data',[90 40]); |
341 |
minmax(e1-e2) |
342 |
minmax(e2-(e16+e17+e18)/3) |
343 |
|
344 |
e1=readbin('test_oldstyle/sFluxtave.0000000018.data',[90 40]); |
345 |
e2=readbin('SFLUX.0000000018.data',[90 40]); |
346 |
minmax(e1-e2) |
347 |
|
348 |
e1=readbin('test_oldstyle/tFluxtave.0000000018.data',[90 40]); |
349 |
e2=readbin('TFLUX.0000000018.data',[90 40]); |
350 |
minmax(e1-e2) |
351 |
|
352 |
e1=readbin('TRELAX.0000000018.data',[90 40]); |
353 |
subplot(211),mypcolor(e1');colorbar |
354 |
|
355 |
e1=readbin('SRELAX.0000000018.data',[90 40]); |
356 |
subplot(211),mypcolor(e1');colorbar |
357 |
|
358 |
e1=readbin('TICE.0000000018.data',[90 40]); |
359 |
subplot(211),mypcolor(e1');colorbar |
360 |
|
361 |
|
362 |
e1=readbin('test_oldstyle/wVeltave.0000000018.data',[90 40 15]); |
363 |
e2=readbin('WVEL.0000000018.data',[90 40 15]); |
364 |
e16=readbin('test_oldstyle/W.0000000016.data',[90 40 15]); |
365 |
e17=readbin('test_oldstyle/W.0000000017.data',[90 40 15]); |
366 |
e18=readbin('test_oldstyle/W.0000000018.data',[90 40 15]); |
367 |
minmax(e1-e2) |
368 |
minmax(e2-(e16+e17+e18)/3) |