1 |
|
|
2 |
Instructions for setting up a cube-sphere integration with sea-ice. |
Instructions for setting up a cube-sphere integration with sea-ice. |
3 |
Date last updated: December 13, 2005 |
|
4 |
216-CPU code and input used for integrating cube27. |
Forcing fields and other input files needed for the CS510 |
5 |
|
configuration can be obtained from |
6 |
|
http://ecco2.org/products/model_setups/cs510.html |
7 |
|
|
8 |
|
================= |
9 |
|
UPDATES |
10 |
|
|
11 |
|
February 6, 2006: cube56, same as cube 50 with SEAICEuseFlooding |
12 |
|
|
13 |
|
February 1, 2006: back to cube50 default (bottomDragQuadratic = 0.002) |
14 |
|
removed UVEL and VVEL from diagnostics |
15 |
|
|
16 |
|
January 5, 2006: cube51, like cube50 but with bottomDragQuadratic = 0.0012 |
17 |
|
|
18 |
|
January 5, 2006: cube50, like cube48 but with snow advection |
19 |
|
|
20 |
|
December 30, 2006: updates for cube49, the first |
21 |
|
"blend" integration. |
22 |
|
|
23 |
|
December 11, 2006: preparing for cube47, same as cube43 |
24 |
|
but with Arctic monthly river runoff climatology. |
25 |
|
|
26 |
|
December 5, 2006: 216-CPU code for cube46 |
27 |
|
Changed diagnostics to accomodate JMCs latest and some |
28 |
|
changes to code-mods/CPP_EEOPTIONS.h and W2_OPTIONS.h |
29 |
|
to speed up code. |
30 |
|
|
31 |
|
November 7, 2006: 216-CPU code for cube43 |
32 |
|
Bug fixes in growth.F routine and new set of diagnostics |
33 |
|
|
34 |
|
October 15, 2006: 216-CPU code for cube40 |
35 |
|
|
36 |
|
September 28, 2006: 216-CPU code for cube38 |
37 |
|
|
38 |
|
September 4, 2006 |
39 |
|
Modified 32*32*6*15 experiments to use grid_cs32.* files from |
40 |
|
verification/tutorial_held_suarez_cs/input instead of |
41 |
|
tile* files from verification/global_ocean.cs32x15/input |
42 |
|
|
43 |
|
August 8, 2006: 216-CPU code and input used for cube37 |
44 |
|
using "June 10, 2006" MITgcm code. |
45 |
|
|
46 |
|
for online T/S profile diagnostics all that is needed is |
47 |
|
profiles_init_fixed |
48 |
|
profiles_inloop |
49 |
|
see pkg/ecco/ecco_cost_init_fixed.F |
50 |
|
pkg/ecco/cost_averagesfields.F |
51 |
|
|
52 |
|
########################################################################## |
53 |
|
########################################################################## |
54 |
|
# Get forcing files, etc. |
55 |
|
# Example below is for cube46 on /nobackup2a/menemenl |
56 |
|
|
57 |
|
cd /nobackup2a/menemenl |
58 |
|
mkdir cube47 |
59 |
|
cd cube47 |
60 |
|
ln -sf ../run_template . |
61 |
|
ln -sf ../ncep_rgau . |
62 |
|
|
63 |
|
|
64 |
########################################################################## |
########################################################################## |
65 |
########################################################################## |
########################################################################## |
66 |
# getting the code from anonymous CVS server |
# getting the code from anonymous CVS server |
67 |
|
|
68 |
|
bash or sh shell: |
69 |
|
$ export CVSROOT=':pserver:cvsanon@mitgcm.org:/u/gcmpack' |
70 |
|
$ cvs login |
71 |
|
( enter the CVS password: "cvsanon" ) |
72 |
|
|
73 |
|
tcsh or csh shell: |
74 |
|
$ setenv CVSROOT ':pserver:cvsanon@mitgcm.org:/u/gcmpack' |
75 |
|
$ cvs login |
76 |
|
( enter the CVS password: "cvsanon" ) |
77 |
|
|
78 |
cvs co MITgcm_contrib/high_res_cube/README_ice |
cvs co MITgcm_contrib/high_res_cube/README_ice |
79 |
cvs co MITgcm_contrib/high_res_cube/code-mods |
cvs co MITgcm_contrib/high_res_cube/code-mods |
80 |
cvs co MITgcm_contrib/high_res_cube/input |
cvs co MITgcm_contrib/high_res_cube/input |
81 |
cvs co MITgcm_contrib/high_res_cube/results |
cvs co MITgcm_contrib/high_res_cube/results |
82 |
cvs co MITgcm_code |
cvs co MITgcm_code |
83 |
cvs co MITgcm/verification/global_ocean.cs32x15 |
cvs co MITgcm/verification/global_ocean.cs32x15 |
84 |
|
cvs co MITgcm/verification/tutorial_held_suarez_cs/input |
85 |
|
|
86 |
|
|
87 |
|
########################################################################## |
88 |
|
########################################################################## |
89 |
|
# 510*510*6*50, 216-cpu cube sphere on altix with ifort and pkg/diagnostics |
90 |
|
|
91 |
|
cd MITgcm/verification/global_ocean.cs32x15/code |
92 |
|
\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* . |
93 |
|
\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* . |
94 |
|
\cp SIZE.h_216 SIZE.h |
95 |
|
\cp packages.conf_hr packages.conf |
96 |
|
|
97 |
|
cd ../build |
98 |
|
\rm * |
99 |
|
../../../tools/genmake2 -mods=../code -of ../code/linux_ia64_ifort+mpi_altix_nas |
100 |
|
make depend |
101 |
|
make -j |
102 |
|
|
103 |
|
cd .. |
104 |
|
mkdir run |
105 |
|
cd run |
106 |
|
\rm * |
107 |
|
\cp ../input/* . |
108 |
|
\cp ../../../../MITgcm_contrib/high_res_cube/input/* . |
109 |
|
ln -sf ../../../../run_template/* . |
110 |
|
\rm pickup* eedata.mth data.pkg data.seaice data.exf_clim data.exf data.mnc data pr* |
111 |
|
\rm bathy_Hmin50.bin *192_94* *cs32* *.m lev* tren* job_c* job_altix150 job_altix92_54t |
112 |
|
ln -sf ../../../../run_template/pickup.0000000216.cube49 pickup.0000000216 |
113 |
|
ln -sf ../../../../run_template/pickup_seaice.0000000216.cube48 pickup_seaice.0000000216 |
114 |
|
\cp ../build/mitgcmuv . |
115 |
|
qsub job_altix |
116 |
|
|
117 |
|
|
118 |
########################################################################## |
########################################################################## |
122 |
cd MITgcm/verification/global_ocean.cs32x15/code |
cd MITgcm/verification/global_ocean.cs32x15/code |
123 |
\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* . |
\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* . |
124 |
\cp ../../../utils/exch2/code-mods/s12t_16x32/* . |
\cp ../../../utils/exch2/code-mods/s12t_16x32/* . |
125 |
|
\rm CPP_EEOPTIONS.h |
126 |
cd ../build |
cd ../build |
127 |
\rm * |
\rm * |
128 |
../../../tools/genmake2 -mods=../code |
../../../tools/genmake2 -mods=../code |
132 |
mkdir run |
mkdir run |
133 |
cd run |
cd run |
134 |
\rm * |
\rm * |
135 |
|
\cp ../../tutorial_held_suarez_cs/input/grid* . |
136 |
\cp ../input/* . |
\cp ../input/* . |
137 |
\cp ../../../../MITgcm_contrib/high_res_cube/input/* . |
\cp ../../../../MITgcm_contrib/high_res_cube/input/* . |
138 |
../build/mitgcmuv >& output.txt |
../build/mitgcmuv >& output.txt |
143 |
to use matlab to look at the output |
to use matlab to look at the output |
144 |
matlab |
matlab |
145 |
|
|
146 |
for fld={'VICEtave','UICEtave','HEFFtave','AREAtave','QSWtave', ... |
for fld={'SIuice','SIvice','SIheff','SIarea','oceQsw', ... |
147 |
'QNETtave','FVtave','FUtave','EmPmRtave'} |
'surForcT','oceTAUX','oceTAUY','surForcS'} |
148 |
tmp=permute(readbin([fld{1} '.0000000020.data'],[32,6,32]),[1 3 2]); |
tmp=permute(readbin([fld{1} '.0000000072.data'],[32,6,32]),[1 3 2]); |
149 |
cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause |
cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause |
150 |
end |
end |
151 |
|
|
152 |
heff=permute(readbin(['HEFFtave.0000000020.data'],[32,6,32]),[1 3 2]); |
heff=permute(readbin(['SIheff.0000000072.data'],[32,6,32]),[1 3 2]); |
153 |
mask=0*heff; mask(find(heff>.001))=1; |
mask=0*heff; mask(find(heff>.001))=1; |
154 |
uice=mask.*permute(readbin(['UICEtave.0000000020.data'],[32,6,32]),[1 3 2]); |
uice=mask.*permute(readbin(['SIuice.0000000072.data'],[32,6,32]),[1 3 2]); |
155 |
vice=mask.*permute(readbin(['VICEtave.0000000020.data'],[32,6,32]),[1 3 2]); |
vice=mask.*permute(readbin(['SIvice.0000000072.data'],[32,6,32]),[1 3 2]); |
156 |
clf |
clf |
157 |
subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3') |
subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3') |
158 |
subplot(222), mypcolor(uice(:,:,6)'); colorbar, title('uice, tile 6') |
subplot(222), mypcolor(uice(:,:,6)'); colorbar, title('uice, tile 6') |
176 |
mkdir run |
mkdir run |
177 |
cd run |
cd run |
178 |
\rm * |
\rm * |
179 |
|
\cp ../../tutorial_held_suarez_cs/input/grid* . |
180 |
\cp ../input/* . |
\cp ../input/* . |
181 |
\cp ../../../../MITgcm_contrib/high_res_cube/input/* . |
\cp ../../../../MITgcm_contrib/high_res_cube/input/* . |
182 |
../build/mitgcmuv >& output.txt |
../build/mitgcmuv >& output.txt |
201 |
mkdir run |
mkdir run |
202 |
cd run |
cd run |
203 |
\rm * |
\rm * |
204 |
|
\cp ../../tutorial_held_suarez_cs/input/grid* . |
205 |
\cp ../input/* . |
\cp ../input/* . |
206 |
\cp ../../../../MITgcm_contrib/high_res_cube/input/* . |
\cp ../../../../MITgcm_contrib/high_res_cube/input/* . |
207 |
mpirun -np 2 ../build/mitgcmuv |
mpirun -np 2 ../build/mitgcmuv |
235 |
|
|
236 |
########################################################################## |
########################################################################## |
237 |
########################################################################## |
########################################################################## |
|
# 510*510*6*50, 216-cpu cube sphere on altix with ifort and pkg/diagnostics |
|
|
# requires input files in sudirectories run_template and ncep_rgau |
|
|
# these can be obained from lou:/u/menemenl/cube |
|
|
|
|
|
cd MITgcm/verification/global_ocean.cs32x15/code |
|
|
\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* . |
|
|
\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* . |
|
|
\cp SIZE.h_216 SIZE.h |
|
|
\cp packages.conf_hr packages.conf |
|
|
|
|
|
cd ../build |
|
|
\rm * |
|
|
../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_ifort+mpi_altix_nas |
|
|
make depend |
|
|
make -j |
|
|
|
|
|
cd .. |
|
|
mkdir run |
|
|
cd run |
|
|
\rm * |
|
|
\cp ../input/* . |
|
|
\cp ../../../../MITgcm_contrib/high_res_cube/input/* . |
|
|
ln -sf ../../../../run_template/* . |
|
|
\rm pickup* |
|
|
ln -sf ../../../../run_template/pickup.0000000216.cube27 pickup.0000000216 |
|
|
ln -sf ../../../../run_template/pickup_seaice.0000000216.cube27 pickup_seaice.0000000216 |
|
|
\cp ../build/mitgcmuv . |
|
|
qsub job_altix |
|
|
|
|
|
|
|
|
########################################################################## |
|
|
########################################################################## |
|
238 |
# 510*510*6*50, 1500-tile, 150-cpu, excludes land |
# 510*510*6*50, 1500-tile, 150-cpu, excludes land |
239 |
# uses s1500t_17x51 configuration |
# uses s1500t_17x51 configuration |
240 |
|
|
264 |
qsub job_altix150 |
qsub job_altix150 |
265 |
|
|
266 |
|
|
|
|
|
|
|
|
267 |
########################################################################## |
########################################################################## |
268 |
########################################################################## |
########################################################################## |
269 |
# 510*510*6*50, 1500-tile, 375-cpu, excludes land |
# 510*510*6*50, 1500-tile, 375-cpu, excludes land |
296 |
job_c20_375cpu < /dev/null >&! output.run & |
job_c20_375cpu < /dev/null >&! output.run & |
297 |
|
|
298 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
299 |
########################################################################## |
########################################################################## |
300 |
########################################################################## |
########################################################################## |
301 |
# 510*510*6*50, 1500-tile, 375-cpu, excludes land |
# 510*510*6*50, 1500-tile, 375-cpu, excludes land |
360 |
job_c19_500cpu < /dev/null >&! output.run & |
job_c19_500cpu < /dev/null >&! output.run & |
361 |
|
|
362 |
|
|
|
|
|
|
|
|
|
|
|
363 |
########################################################################## |
########################################################################## |
364 |
########################################################################## |
########################################################################## |
365 |
# 510*510*6*50, 216-cpu cube sphere on cosmos |
# 510*510*6*50, 216-cpu cube sphere on cosmos |
395 |
bsub < job_cosmos |
bsub < job_cosmos |
396 |
|
|
397 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
398 |
########################################################################## |
########################################################################## |
399 |
########################################################################## |
########################################################################## |
400 |
# 510*510*6*50, 216-tile, 54-cpu cube sphere on the altices with efc |
# 510*510*6*50, 216-tile, 54-cpu cube sphere on the altices with efc |