/[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.52 - (hide annotations) (download)
Sat Feb 19 06:28:28 2005 UTC (19 years, 3 months ago) by dimitri
Branch: MAIN
Changes since 1.51: +5 -36 lines
high_res_cube/*

1 dimitri 1.52 daInstructions for setting up a cube-sphere integration with sea-ice.
2     Date last updated: February 17, 2005
3     Code and input used for integrating cube22.
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 dimitri 1.43 cvs co MITgcm_code
15     cvs co MITgcm/verification/global_ocean.cs32x15
16     cvs co MITgcm/verification/testreport
17 dimitri 1.14
18    
19     ##########################################################################
20     ##########################################################################
21     # 32*32*6*15, 1-cpu cube sphere
22 dimitri 1.3
23 dimitri 1.1 cd MITgcm/verification/global_ocean.cs32x15/code
24 dimitri 1.8 \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
25 dimitri 1.33 \cp ../../../utils/exch2/code-mods/s12t_16x32/* .
26 dimitri 1.37 cd ../build
27     \rm *
28     ../../../tools/genmake2 -mods=../code
29     make depend
30     make
31     cd ..
32     mkdir run
33     cd run
34     \rm *
35     \cp ../input/* .
36     \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
37     ../build/mitgcmuv >& output.txt
38 dimitri 1.3
39 dimitri 1.37 comparison output is in:
40     ../../../../MITgcm_contrib/high_res_cube/results/output.txt
41    
42     to use matlab to look at the output
43     matlab
44    
45     for fld={'Eta','VICE','UICE','HEFF','AREA','Qsw','Qnet','FV','FU','EmPmR'}
46     tmp=permute(readbin([fld{1} '.0000000020.data'],[32,6,32]),[1 3 2]);
47     cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
48     end
49    
50     heff=permute(readbin(['HEFF.0000000020.data'],[32,6,32]),[1 3 2]);
51     mask=0*heff; mask(find(heff>.001))=1;
52     uice=mask.*permute(readbin(['UICE.0000000020.data'],[32,6,32]),[1 3 2]);
53     vice=mask.*permute(readbin(['VICE.0000000020.data'],[32,6,32]),[1 3 2]);
54     clf
55     subplot(221), mypcolor(uice(:,:,3)'), colorbar, title('uice, tile 3')
56     subplot(222), mypcolor(uice(:,:,6)'), colorbar, title('uice, tile 6')
57     subplot(223), mypcolor(vice(:,:,3)'), colorbar, title('vice, tile 3')
58     subplot(224), mypcolor(vice(:,:,6)'), colorbar, title('vice, tile 6')
59    
60    
61     ##########################################################################
62     ##########################################################################
63     # 32*32*6*15, 17-tile, 1-cpu cube sphere experiment
64    
65     cd MITgcm/verification/global_ocean.cs32x15/code
66     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
67     \cp ../code_alt/code.176t_8x4/* .
68 dimitri 1.25 cd ../build
69 dimitri 1.37 \rm *
70 dimitri 1.10 ../../../tools/genmake2 -mods=../code
71 dimitri 1.1 make depend
72     make
73     cd ..
74     mkdir run
75     cd run
76 dimitri 1.37 \rm *
77 dimitri 1.8 \cp ../input/* .
78     \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
79 dimitri 1.3 ../build/mitgcmuv >& output.txt
80    
81     comparison output is in:
82     ../../../../MITgcm_contrib/high_res_cube/results/output.txt
83    
84     to use matlab to look at the output
85     matlab
86 dimitri 1.17
87 dimitri 1.5 for fld={'Eta','VICE','UICE','HEFF','AREA','Qsw','Qnet','FV','FU','EmPmR'}
88 dimitri 1.13 tmp=permute(readbin([fld{1} '.0000000020.data'],[32,6,32]),[1 3 2]);
89 dimitri 1.4 cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
90     end
91 dimitri 1.8
92 dimitri 1.17 heff=permute(readbin(['HEFF.0000000020.data'],[32,6,32]),[1 3 2]);
93     mask=0*heff; mask(find(heff>.001))=1;
94     uice=mask.*permute(readbin(['UICE.0000000020.data'],[32,6,32]),[1 3 2]);
95     vice=mask.*permute(readbin(['VICE.0000000020.data'],[32,6,32]),[1 3 2]);
96     clf
97     subplot(221), mypcolor(uice(:,:,3)'), colorbar, title('uice, tile 3')
98     subplot(222), mypcolor(uice(:,:,6)'), colorbar, title('uice, tile 6')
99     subplot(223), mypcolor(vice(:,:,3)'), colorbar, title('vice, tile 3')
100     subplot(224), mypcolor(vice(:,:,6)'), colorbar, title('vice, tile 6')
101    
102 dimitri 1.8
103     ##########################################################################
104     ##########################################################################
105     # 32*32*6*15, 2-cpu cube sphere on the altices with efc
106    
107     cd MITgcm/verification/global_ocean.cs32x15/code
108     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
109 dimitri 1.33 \cp ../../../utils/exch2/code-mods/s12t_16x32/* .
110 dimitri 1.8 \cp SIZE.h_mpi SIZE.h
111    
112 dimitri 1.25 cd ../build
113 dimitri 1.14 \rm *
114 dimitri 1.41 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi_altix
115 dimitri 1.8 make depend
116 dimitri 1.41 make -j
117 dimitri 1.8
118     cd ..
119     mkdir run
120     cd run
121 dimitri 1.14 \rm *
122 dimitri 1.8 \cp ../input/* .
123     \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
124     mpirun -np 2 ../build/mitgcmuv
125    
126    
127     ##########################################################################
128     ##########################################################################
129 dimitri 1.35 # 510*510*6*50, 54-tile, 54-cpu cube sphere on the altices with efc
130 dimitri 1.36 # starting with pickup file from end of 1992
131 dimitri 1.35 # requires input files in sudirectories run_template and ncep_rgau
132    
133     cd MITgcm/verification/global_ocean.cs32x15/code
134     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
135     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s54t_170x170/* .
136     \cp packages.conf_hr packages.conf
137    
138     cd ../build
139     \rm *
140 dimitri 1.40 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi_altix
141 dimitri 1.35 make depend
142 dimitri 1.41 make -j
143 dimitri 1.35
144     cd ..
145     mkdir run
146     cd run
147     \rm *
148     \cp ../input/* .
149     \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
150 dimitri 1.41 qsub job_altix92_54t
151 dimitri 1.35
152    
153     ##########################################################################
154     ##########################################################################
155 dimitri 1.52 # 510*510*6*50, 216-cpu cube sphere on altix with ifort and pkg/diagnostics
156 dimitri 1.23 # requires input files in sudirectories run_template and ncep_rgau
157 dimitri 1.41 # these can be obained from lou:/u/menemenl/cube
158 dimitri 1.23
159     cd MITgcm/verification/global_ocean.cs32x15/code
160     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
161 dimitri 1.41 \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* .
162     \cp SIZE.h_216 SIZE.h
163 dimitri 1.23 \cp packages.conf_hr packages.conf
164 dimitri 1.8
165 dimitri 1.25 cd ../build
166 dimitri 1.14 \rm *
167 dimitri 1.46 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_ifort+mpi_altix_nas
168 dimitri 1.8 make depend
169 dimitri 1.52 make -j
170 dimitri 1.8
171     cd ..
172     mkdir run
173     cd run
174 dimitri 1.14 \rm *
175 dimitri 1.8 \cp ../input/* .
176     \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
177 dimitri 1.44 ln -sf ../../../../run_template/* .
178 dimitri 1.50 \rm pickup*
179     ln -sf ../../../../run_template/pickup.0000316440.cube20 pickup.0000000216.data
180     ln -sf ../../../../run_template/pickup_seaice.0000316440.cube20 pickup_seaice.0000000216.data
181 dimitri 1.44 \cp ../build/mitgcmuv .
182 dimitri 1.41 qsub job_altix92
183    
184    
185    
186    
187    
188    
189 dimitri 1.51
190    
191    
192    
193    
194 dimitri 1.41
195    
196    
197    
198    
199    
200    
201    
202    
203    
204    
205    
206 dimitri 1.42 ##########################################################################
207     ##########################################################################
208     # 510*510*6*50, 216-tile, 54-cpu cube sphere on the altices with efc
209     # requires input files in sudirectories run_template and ncep_rgau
210     # presently residing under orion:/tmp1/dmenem/cube
211 dimitri 1.41
212 dimitri 1.42 cd MITgcm/verification/global_ocean.cs32x15/code
213     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
214     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* .
215     \cp SIZE.h_54 SIZE.h
216     \cp packages.conf_hr packages.conf
217 dimitri 1.41
218 dimitri 1.42 cd ../build
219     \rm *
220     ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi_altix
221     make depend
222     make -j
223 dimitri 1.41
224 dimitri 1.42 cd ..
225     mkdir run
226     cd run
227     \rm *
228     \cp ../input/* .
229     \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
230     qsub job_orion
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
311 dimitri 1.48
312    
313    
314    
315    
316    
317     #!/bin/csh
318     \cp data.pkg_hr data.pkg
319     \cp data.seaice_hr data.seaice
320     \cp data.exf_clim_hr data.exf_clim
321    
322     \cp data_hr92 data
323     ./modpickup
324     \cp data.exf_hr92 data.exf
325     /opt/sgi/mpt/1.11-85/bin/mpirun.real -v -a a2k c20 216 runon 221-436 mitgcmuv
326     mkdir out.1992
327     \mv STDOUT.0000 *.data pickup*.meta out.1992
328     \cp data* out.1992
329     \rm STD* *.meta

  ViewVC Help
Powered by ViewVC 1.1.22