/[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.31 - (hide annotations) (download)
Wed Mar 3 06:56:37 2004 UTC (21 years, 4 months ago) by dimitri
Branch: MAIN
Changes since 1.30: +5 -0 lines 
modified high_res_cube/*
1 dimitri 1.1 Instructions for setting up a cube-sphere integration with sea-ice.
2 dimitri 1.30 Date last tested: March 2, 2004.
3 dimitri 1.28 This is the code used for cube4, the third 1992-2002 510x510x6
4 dimitri 1.14 cube-sphere integration.
5 dimitri 1.8
6    
7     ##########################################################################
8     ##########################################################################
9 dimitri 1.14 # getting the code from anonymous CVS server
10 dimitri 1.1
11     setenv CVSROOT :pserver:cvsanon@mitgcm.org:/u/u0/gcmpack
12     cvs login ( CVS password: cvsanon )
13 dimitri 1.30 cvs co -r hrcube4a MITgcm_contrib/high_res_cube/README_ice
14     cvs co -r hrcube4a MITgcm_contrib/high_res_cube/code-mods
15     cvs co -r hrcube4a MITgcm_contrib/high_res_cube/matlab-grid-converter
16     cvs co -r hrcube4a MITgcm_contrib/high_res_cube/input
17     cvs co -r hrcube4a MITgcm_contrib/high_res_cube/results
18     cvs co -r hrcube4a MITgcm_code
19     cvs co -r hrcube4a MITgcm/verification/global_ocean.cs32x15
20 dimitri 1.14
21    
22     ##########################################################################
23     ##########################################################################
24     # 32*32*6*15, 1-cpu cube sphere
25 dimitri 1.3
26 dimitri 1.1 cd MITgcm/verification/global_ocean.cs32x15/code
27 dimitri 1.8 \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
28     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s12t_16x32/* .
29 dimitri 1.3
30 dimitri 1.25 cd ../build
31 dimitri 1.10 ../../../tools/genmake2 -mods=../code
32 dimitri 1.1 make depend
33 dimitri 1.31 \rm getcon.F
34 dimitri 1.1 make
35 dimitri 1.3
36 dimitri 1.1 cd ..
37     mkdir run
38     cd run
39 dimitri 1.8 \cp ../input/* .
40     \cp ../../../../MITgcm_contrib/high_res_cube/matlab-grid-converter/tile00* .
41     \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
42 dimitri 1.3 ../build/mitgcmuv >& output.txt
43    
44     comparison output is in:
45     ../../../../MITgcm_contrib/high_res_cube/results/output.txt
46    
47     to use matlab to look at the output
48     matlab
49 dimitri 1.17
50 dimitri 1.5 for fld={'Eta','VICE','UICE','HEFF','AREA','Qsw','Qnet','FV','FU','EmPmR'}
51 dimitri 1.13 tmp=permute(readbin([fld{1} '.0000000020.data'],[32,6,32]),[1 3 2]);
52 dimitri 1.4 cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
53     end
54 dimitri 1.8
55 dimitri 1.17 heff=permute(readbin(['HEFF.0000000020.data'],[32,6,32]),[1 3 2]);
56     mask=0*heff; mask(find(heff>.001))=1;
57     uice=mask.*permute(readbin(['UICE.0000000020.data'],[32,6,32]),[1 3 2]);
58     vice=mask.*permute(readbin(['VICE.0000000020.data'],[32,6,32]),[1 3 2]);
59     clf
60     subplot(221), mypcolor(uice(:,:,3)'), colorbar, title('uice, tile 3')
61     subplot(222), mypcolor(uice(:,:,6)'), colorbar, title('uice, tile 6')
62     subplot(223), mypcolor(vice(:,:,3)'), colorbar, title('vice, tile 3')
63     subplot(224), mypcolor(vice(:,:,6)'), colorbar, title('vice, tile 6')
64    
65 dimitri 1.8
66     ##########################################################################
67     ##########################################################################
68     # 32*32*6*15, 2-cpu cube sphere on the altices with efc
69    
70     cd MITgcm/verification/global_ocean.cs32x15/code
71     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
72     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s12t_16x32/* .
73     \cp SIZE.h_mpi SIZE.h
74    
75 dimitri 1.25 cd ../build
76 dimitri 1.14 \rm *
77 dimitri 1.10 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi
78 dimitri 1.8 make depend
79 dimitri 1.31 \rm getcon.F
80 dimitri 1.8 make
81    
82     cd ..
83     mkdir run
84     cd run
85 dimitri 1.14 \rm *
86 dimitri 1.8 \cp ../input/* .
87     \cp ../../../../MITgcm_contrib/high_res_cube/matlab-grid-converter/tile00* .
88     \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
89     mpirun -np 2 ../build/mitgcmuv
90    
91    
92     ##########################################################################
93     ##########################################################################
94 dimitri 1.9 # 510*510*6*50, 54-cpu cube sphere on the altices with efc
95 dimitri 1.14 # requires input files in sudirectories run_template and ncep_rgau
96     # presently residing under orion:/tmp1/dmenem/cube
97 dimitri 1.8
98     cd MITgcm/verification/global_ocean.cs32x15/code
99     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
100     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* .
101     \cp SIZE.h_54 SIZE.h
102 dimitri 1.11 \cp packages.conf_hr packages.conf
103 dimitri 1.23
104 dimitri 1.25 cd ../build
105 dimitri 1.23 \rm *
106     ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi
107     make depend
108 dimitri 1.31 \rm getcon.F
109 dimitri 1.23 make
110    
111     cd ..
112     mkdir run
113     cd run
114     \rm *
115     \cp ../input/* .
116     \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
117     qsub job_orion
118    
119    
120     ##########################################################################
121     ##########################################################################
122     # 510*510*6*50, 205-tile, 41-cpu cube sphere on the altices with efc
123     # requires input files in sudirectories run_template and ncep_rgau
124     # presently residing under orion:/tmp1/dmenem/cube
125    
126     cd MITgcm/verification/global_ocean.cs32x15/code
127     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
128     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s205t_85x85/* .
129     \cp SIZE.h_41 SIZE.h
130     \cp packages.conf_hr packages.conf
131 dimitri 1.8
132 dimitri 1.25 cd ../build
133 dimitri 1.14 \rm *
134 dimitri 1.10 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi
135 dimitri 1.8 make depend
136 dimitri 1.31 \rm getcon.F
137 dimitri 1.8 make
138    
139     cd ..
140     mkdir run
141     cd run
142 dimitri 1.14 \rm *
143 dimitri 1.8 \cp ../input/* .
144     \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
145 dimitri 1.24 qsub job_orion41
146 dimitri 1.9
147    
148     ##########################################################################
149     ##########################################################################
150 dimitri 1.13 # 510*510*6*50, 216-cpu cube sphere on altix with efc
151 dimitri 1.17 # requires input files in sudirectories run_template and ncep_rgau
152     # these can be obained from lou:/u/menemenl/cube
153 dimitri 1.9
154     cd MITgcm/verification/global_ocean.cs32x15/code
155     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
156     \cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* .
157     \cp SIZE.h_216 SIZE.h
158 dimitri 1.12 \cp packages.conf_hr packages.conf
159 dimitri 1.9
160 dimitri 1.25 cd ../build
161 dimitri 1.15 \rm *
162 dimitri 1.13 ../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_efc+mpi_altix
163 dimitri 1.9 make depend
164 dimitri 1.31 \rm getcon.F
165 dimitri 1.9 make
166    
167     cd ..
168     mkdir run
169     cd run
170 dimitri 1.15 \rm *
171     \cp ../input/* .
172     \cp ../../../../MITgcm_contrib/high_res_cube/input/* .
173 dimitri 1.29 qsub job_altix
174 dimitri 1.17
175    
176     ##########################################################################
177     ##########################################################################
178     # matlab script for generating pChkptFreq
179    
180     nTimeSteps=26352;
181     for niter0=216:26352:(26352*11)
182     pChkptFreq=(niter0+nTimeSteps)*1200;
183     day=pChkptFreq/24/60/60;
184     disp(['niter0=' int2str(niter0) ',day=' int2str(day) ', ' ...
185     datestr(datenum(1992,1,1)+day) ', pChkptFreq=' int2str(pChkptFreq)])
186     end
187    
188     niter0 = 216 , day = 369 , 04-Jan-1993, pChkptFreq = 31881600
189     niter0 = 26568 , day = 735 , 05-Jan-1994, pChkptFreq = 63504000
190     niter0 = 52920 , day = 1101, 06-Jan-1995, pChkptFreq = 95126400
191     niter0 = 79272 , day = 1467, 07-Jan-1996, pChkptFreq = 126748800
192     niter0 = 105624, day = 1833, 07-Jan-1997, pChkptFreq = 158371200
193     niter0 = 131976, day = 2199, 08-Jan-1998, pChkptFreq = 189993600
194     niter0 = 158328, day = 2565, 09-Jan-1999, pChkptFreq = 221616000
195     niter0 = 184680, day = 2931, 10-Jan-2000, pChkptFreq = 253238400
196     niter0 = 211032, day = 3297, 10-Jan-2001, pChkptFreq = 284860800
197     niter0 = 237384, day = 3663, 11-Jan-2002, pChkptFreq = 316483200
198     niter0 = 263736, day = 4029, 12-Jan-2003, pChkptFreq = 348105600
  ViewVC Help
Powered by ViewVC 1.1.22