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 |
|
|
4 |
|
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 1, 2006: back to cube50 default (bottomDragQuadratic = 0.0012) |
12 |
|
removed UVEL and VVEL from diagnostics |
13 |
|
|
14 |
|
January 5, 2006: cube51, like cube50 but with bottomDragQuadratic = 0.0012 |
15 |
|
|
16 |
|
January 5, 2006: cube50, like cube48 but with snow advection |
17 |
|
|
18 |
|
December 30, 2006: updates for cube49, the first |
19 |
|
"blend" integration. |
20 |
|
|
21 |
|
December 11, 2006: preparing for cube47, same as cube43 |
22 |
|
but with Arctic monthly river runoff climatology. |
23 |
|
|
24 |
|
December 5, 2006: 216-CPU code for cube46 |
25 |
|
Changed diagnostics to accomodate JMCs latest and some |
26 |
|
changes to code-mods/CPP_EEOPTIONS.h and W2_OPTIONS.h |
27 |
|
to speed up code. |
28 |
|
|
29 |
|
November 7, 2006: 216-CPU code for cube43 |
30 |
|
Bug fixes in growth.F routine and new set of diagnostics |
31 |
|
|
32 |
|
October 15, 2006: 216-CPU code for cube40 |
33 |
|
|
34 |
September 28, 2006: 216-CPU code for cube38 |
September 28, 2006: 216-CPU code for cube38 |
35 |
|
|
36 |
September 4, 2006 |
September 4, 2006 |
38 |
verification/tutorial_held_suarez_cs/input instead of |
verification/tutorial_held_suarez_cs/input instead of |
39 |
tile* files from verification/global_ocean.cs32x15/input |
tile* files from verification/global_ocean.cs32x15/input |
40 |
|
|
41 |
August 8, 2006 |
August 8, 2006: 216-CPU code and input used for cube37 |
|
216-CPU code and input used for integrating cube37 |
|
42 |
using "June 10, 2006" MITgcm code. |
using "June 10, 2006" MITgcm code. |
43 |
|
|
44 |
for online T/S profile diagnostics all that is needed is |
for online T/S profile diagnostics all that is needed is |
47 |
see pkg/ecco/ecco_cost_init_fixed.F |
see pkg/ecco/ecco_cost_init_fixed.F |
48 |
pkg/ecco/cost_averagesfields.F |
pkg/ecco/cost_averagesfields.F |
49 |
|
|
50 |
|
########################################################################## |
51 |
|
########################################################################## |
52 |
|
# Get forcing files, etc. |
53 |
|
# Example below is for cube46 on /nobackup2a/menemenl |
54 |
|
|
55 |
|
cd /nobackup2a/menemenl |
56 |
|
mkdir cube47 |
57 |
|
cd cube47 |
58 |
|
ln -sf ../run_template . |
59 |
|
ln -sf ../ncep_rgau . |
60 |
|
|
61 |
|
|
62 |
########################################################################## |
########################################################################## |
63 |
########################################################################## |
########################################################################## |
64 |
# getting the code from anonymous CVS server |
# getting the code from anonymous CVS server |
|
# |
|
65 |
|
|
66 |
bash or sh shell: |
bash or sh shell: |
67 |
$ export CVSROOT=':pserver:cvsanon@mitgcm.org:/u/gcmpack' |
$ export CVSROOT=':pserver:cvsanon@mitgcm.org:/u/gcmpack' |
84 |
|
|
85 |
########################################################################## |
########################################################################## |
86 |
########################################################################## |
########################################################################## |
87 |
|
# 510*510*6*50, 216-cpu cube sphere on altix with ifort and pkg/diagnostics |
88 |
|
|
89 |
|
cd MITgcm/verification/global_ocean.cs32x15/code |
90 |
|
\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* . |
91 |
|
\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* . |
92 |
|
\cp SIZE.h_216 SIZE.h |
93 |
|
\cp packages.conf_hr packages.conf |
94 |
|
|
95 |
|
cd ../build |
96 |
|
\rm * |
97 |
|
../../../tools/genmake2 -mods=../code -of ../code/linux_ia64_ifort+mpi_altix_nas |
98 |
|
make depend |
99 |
|
make -j |
100 |
|
|
101 |
|
cd .. |
102 |
|
mkdir run |
103 |
|
cd run |
104 |
|
\rm * |
105 |
|
\cp ../input/* . |
106 |
|
\cp ../../../../MITgcm_contrib/high_res_cube/input/* . |
107 |
|
ln -sf ../../../../run_template/* . |
108 |
|
\rm pickup* eedata.mth data.pkg data.seaice data.exf_clim data.exf data.mnc data pr* |
109 |
|
\rm bathy_Hmin50.bin *192_94* *cs32* *.m lev* tren* job_c* job_altix150 job_altix92_54t |
110 |
|
ln -sf ../../../../run_template/pickup.0000000216.cube49 pickup.0000000216 |
111 |
|
ln -sf ../../../../run_template/pickup_seaice.0000000216.cube48 pickup_seaice.0000000216 |
112 |
|
\cp ../build/mitgcmuv . |
113 |
|
qsub job_altix |
114 |
|
|
115 |
|
|
116 |
|
########################################################################## |
117 |
|
########################################################################## |
118 |
# 32*32*6*15, 1-cpu cube sphere |
# 32*32*6*15, 1-cpu cube sphere |
119 |
|
|
120 |
cd MITgcm/verification/global_ocean.cs32x15/code |
cd MITgcm/verification/global_ocean.cs32x15/code |
121 |
\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* . |
\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* . |
122 |
\cp ../../../utils/exch2/code-mods/s12t_16x32/* . |
\cp ../../../utils/exch2/code-mods/s12t_16x32/* . |
123 |
|
\rm CPP_EEOPTIONS.h |
124 |
cd ../build |
cd ../build |
125 |
\rm * |
\rm * |
126 |
../../../tools/genmake2 -mods=../code |
../../../tools/genmake2 -mods=../code |
141 |
to use matlab to look at the output |
to use matlab to look at the output |
142 |
matlab |
matlab |
143 |
|
|
144 |
for fld={'VICEtave','UICEtave','HEFFtave','AREAtave','QSWtave', ... |
for fld={'SIuice','SIvice','SIheff','SIarea','oceQsw', ... |
145 |
'QNETtave','FVtave','FUtave','EmPmRtave'} |
'surForcT','oceTAUX','oceTAUY','surForcS'} |
146 |
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]); |
147 |
cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause |
cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause |
148 |
end |
end |
149 |
|
|
150 |
heff=permute(readbin(['HEFFtave.0000000020.data'],[32,6,32]),[1 3 2]); |
heff=permute(readbin(['SIheff.0000000072.data'],[32,6,32]),[1 3 2]); |
151 |
mask=0*heff; mask(find(heff>.001))=1; |
mask=0*heff; mask(find(heff>.001))=1; |
152 |
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]); |
153 |
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]); |
154 |
clf |
clf |
155 |
subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3') |
subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3') |
156 |
subplot(222), mypcolor(uice(:,:,6)'); colorbar, title('uice, tile 6') |
subplot(222), mypcolor(uice(:,:,6)'); colorbar, title('uice, tile 6') |
159 |
|
|
160 |
|
|
161 |
########################################################################## |
########################################################################## |
|
########################################################################## |
|
|
# 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_ia\ |
|
|
64_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* data.pkg data.seaice data.exf_clim data.exf data bathy_Hmin50.bin |
|
|
\rm *192_94* *cs32* *.m lev* tren* job_c* job_altix150 job_altix92_54t |
|
|
ln -sf ../../../../run_template/pickup.0000000216.cube38 pickup.0000000216 |
|
|
ln -sf ../../../../run_template/pickup_seaice.0000000216.cube38 pickup_seaice.0000000216 |
|
|
\cp ../build/mitgcmuv . |
|
|
qsub job_altix |
|
|
|
|
|
|
|
|
########################################################################## |
|
162 |
########################################################################## |
########################################################################## |
163 |
# 32*32*6*15, 17-tile, 1-cpu cube sphere experiment |
# 32*32*6*15, 17-tile, 1-cpu cube sphere experiment |
164 |
|
|