1 |
Instructions for setting up a cube-sphere integration |
2 |
with sea-ice on the cs32 grid. Is meant to be used as |
3 |
a preliminary sanity check for the cs510 configuration. |
4 |
|
5 |
########################################################################## |
6 |
# getting the code from anonymous CVS server |
7 |
|
8 |
bash or sh shell: |
9 |
$ export CVSROOT=':pserver:cvsanon@mitgcm.org:/u/gcmpack' |
10 |
$ cvs login |
11 |
( enter the CVS password: "cvsanon" ) |
12 |
|
13 |
tcsh or csh shell: |
14 |
$ setenv CVSROOT ':pserver:cvsanon@mitgcm.org:/u/gcmpack' |
15 |
$ cvs login |
16 |
( enter the CVS password: "cvsanon" ) |
17 |
|
18 |
cvs co MITgcm_contrib/high_res_cube/README.cs32 |
19 |
cvs co MITgcm_contrib/high_res_cube/code-mods |
20 |
cvs co MITgcm_contrib/high_res_cube/input |
21 |
cvs co MITgcm_contrib/high_res_cube/results |
22 |
cvs co MITgcm_code |
23 |
cvs co MITgcm/verification/global_ocean.cs32x15 |
24 |
cvs co MITgcm/verification/tutorial_held_suarez_cs/input |
25 |
|
26 |
########################################################################## |
27 |
# 32*32*6*15, 1-cpu cube sphere |
28 |
|
29 |
cd MITgcm |
30 |
mkdir build run |
31 |
cd build |
32 |
\rm * |
33 |
\cp ../../MITgcm_contrib/high_res_cube/code-mods/* . |
34 |
\cp ../utils/exch2/code-mods/s12t_16x32/* . |
35 |
ex - SIZE.h >> /dev/null <<EOF |
36 |
/OLx = |
37 |
c |
38 |
& OLx = 5, |
39 |
. |
40 |
/OLy = |
41 |
c |
42 |
& OLy = 5, |
43 |
. |
44 |
w |
45 |
q |
46 |
EOF |
47 |
\rm CPP_EEOPTIONS.h |
48 |
../tools/genmake2 |
49 |
make depend |
50 |
make -j |
51 |
cd ../run |
52 |
\rm * |
53 |
\cp ../verification/tutorial_held_suarez_cs/input/grid* . |
54 |
\cp ../verification/global_ocean.cs32x15/input/lev_*_cs_15k.bin . |
55 |
\cp ../verification/global_ocean.cs32x15/input/bathy_Hmin50.bin . |
56 |
\cp ../../MITgcm_contrib/high_res_cube/input/* . |
57 |
../build/mitgcmuv >& output.txt |
58 |
|
59 |
emacs output.txt ../../MITgcm_contrib/high_res_cube/results/output.txt |
60 |
|
61 |
to use matlab to look at the output |
62 |
matlab |
63 |
|
64 |
for fld={'SIuice','SIvice','SIheff','SIarea','oceQsw', ... |
65 |
'surForcT','oceTAUX','oceTAUY','surForcS'} |
66 |
tmp=permute(readbin([fld{1} '.0000000072.data'],[32,6,32]),[1 3 2]); |
67 |
cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause |
68 |
end |
69 |
|
70 |
heff=permute(readbin(['SIheff.0000000072.data'],[32,6,32]),[1 3 2]); |
71 |
mask=0*heff; mask(find(heff>.001))=1; |
72 |
uice=mask.*permute(readbin(['SIuice.0000000072.data'],[32,6,32]),[1 3 2]); |
73 |
vice=mask.*permute(readbin(['SIvice.0000000072.data'],[32,6,32]),[1 3 2]); |
74 |
clf |
75 |
subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3') |
76 |
subplot(222), mypcolor(uice(:,:,6)'); colorbar, title('uice, tile 6') |
77 |
subplot(223), mypcolor(vice(:,:,3)'); colorbar, title('vice, tile 3') |
78 |
subplot(224), mypcolor(vice(:,:,6)'); colorbar, title('vice, tile 6') |
79 |
|
80 |
|
81 |
########################################################################## |
82 |
# 32*32*6*15, 17-tile, 1-cpu cube sphere experiment |
83 |
# needs fixing; does not work as of checkpoint58x_post |
84 |
# see discussion on MITgcm-devel |
85 |
|
86 |
cd ../build |
87 |
\rm * |
88 |
\cp ../../MITgcm_contrib/high_res_cube/code-mods/* . |
89 |
\cp ../verification/global_ocean.cs32x15/code_alt/code.176t_8x4/* . |
90 |
\rm CPP_EEOPTIONS.h |
91 |
../tools/genmake2 |
92 |
make depend |
93 |
make -j |
94 |
cd ../run |
95 |
\rm * |
96 |
\cp ../verification/tutorial_held_suarez_cs/input/grid* . |
97 |
\cp ../verification/global_ocean.cs32x15/input/lev_*_cs_15k.bin . |
98 |
\cp ../verification/global_ocean.cs32x15/input/bathy_Hmin50.bin . |
99 |
\cp ../../MITgcm_contrib/high_res_cube/input/* . |
100 |
../build/mitgcmuv >& output.txt |
101 |
|
102 |
|
103 |
########################################################################## |
104 |
# 32*32*6*15, 2-cpu cube sphere on the NAS altices with ifort |
105 |
|
106 |
cd MITgcm |
107 |
mkdir build run |
108 |
cd build |
109 |
\rm * |
110 |
\cp ../../MITgcm_contrib/high_res_cube/code-mods/* . |
111 |
\cp ../utils/exch2/code-mods/s12t_16x32/* . |
112 |
\mv SIZE.h_mpi SIZE.h |
113 |
../tools/genmake2 -of ../../MITgcm_contrib/high_res_cube/code-mods/linux_ia64_ifort+mpi_altix_nas |
114 |
make depend |
115 |
make -j |
116 |
cd ../run |
117 |
\rm * |
118 |
\cp ../verification/tutorial_held_suarez_cs/input/grid* . |
119 |
\cp ../verification/global_ocean.cs32x15/input/lev_*_cs_15k.bin . |
120 |
\cp ../verification/global_ocean.cs32x15/input/bathy_Hmin50.bin . |
121 |
\cp ../../MITgcm_contrib/high_res_cube/input/* . |
122 |
mpirun -np 2 ../build/mitgcmuv |
123 |
|
124 |
|
125 |
########################################################################## |
126 |
# 32*32*6*15, 2-cpu cube sphere on the JPL altices with ifort |
127 |
|
128 |
cd MITgcm |
129 |
mkdir build run |
130 |
cd build |
131 |
\rm * |
132 |
\cp ../../MITgcm_contrib/high_res_cube/code-mods/* . |
133 |
\cp ../utils/exch2/code-mods/s12t_16x32/* . |
134 |
\mv SIZE.h_mpi SIZE.h |
135 |
../tools/genmake2 -of ../tools/build_options/linux_ia64_ifort+mpi_altix_jpl |
136 |
make depend |
137 |
make -j |
138 |
cd ../run |
139 |
\rm * |
140 |
\cp ../verification/tutorial_held_suarez_cs/input/grid* . |
141 |
\cp ../verification/global_ocean.cs32x15/input/lev_*_cs_15k.bin . |
142 |
\cp ../verification/global_ocean.cs32x15/input/bathy_Hmin50.bin . |
143 |
\cp ../../MITgcm_contrib/high_res_cube/input/* . |
144 |
mpirun -np 2 ../build/mitgcmuv |