3 |
================= |
================= |
4 |
UPDATES |
UPDATES |
5 |
|
|
6 |
|
December 11, 2006: preparing for cube47, same as cube43 |
7 |
|
but with Arctic monthly river runoff climatology |
8 |
|
|
9 |
|
December 5, 2006: 216-CPU code for cube46 |
10 |
|
Changed diagnostics to accomodate JMCs latest and some |
11 |
|
changes to code-mods/CPP_EEOPTIONS.h and W2_OPTIONS.h |
12 |
|
to speed up code. |
13 |
|
|
14 |
November 7, 2006: 216-CPU code for cube43 |
November 7, 2006: 216-CPU code for cube43 |
15 |
Bug fixes in growth.F routine and new set of diagnostics |
Bug fixes in growth.F routine and new set of diagnostics |
16 |
|
|
32 |
see pkg/ecco/ecco_cost_init_fixed.F |
see pkg/ecco/ecco_cost_init_fixed.F |
33 |
pkg/ecco/cost_averagesfields.F |
pkg/ecco/cost_averagesfields.F |
34 |
|
|
35 |
|
########################################################################## |
36 |
|
########################################################################## |
37 |
|
# Get forcing files, etc. |
38 |
|
# Example below is for cube46 on /nobackup2a/menemenl |
39 |
|
|
40 |
|
cd /nobackup2a/menemenl |
41 |
|
mkdir cube46 |
42 |
|
cd cube46 |
43 |
|
ln -sf ../run_template . |
44 |
|
ln -sf ../ncep_rgau . |
45 |
|
|
46 |
|
|
47 |
########################################################################## |
########################################################################## |
48 |
########################################################################## |
########################################################################## |
49 |
# getting the code from anonymous CVS server |
# getting the code from anonymous CVS server |
|
# |
|
50 |
|
|
51 |
bash or sh shell: |
bash or sh shell: |
52 |
$ export CVSROOT=':pserver:cvsanon@mitgcm.org:/u/gcmpack' |
$ export CVSROOT=':pserver:cvsanon@mitgcm.org:/u/gcmpack' |
70 |
########################################################################## |
########################################################################## |
71 |
########################################################################## |
########################################################################## |
72 |
# 510*510*6*50, 216-cpu cube sphere on altix with ifort and pkg/diagnostics |
# 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 |
|
73 |
|
|
74 |
cd MITgcm/verification/global_ocean.cs32x15/code |
cd MITgcm/verification/global_ocean.cs32x15/code |
75 |
\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* . |
\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* . |
79 |
|
|
80 |
cd ../build |
cd ../build |
81 |
\rm * |
\rm * |
82 |
../../../tools/genmake2 -mods=../code -of ../../../tools/build_options/linux_ia64_ifort+mpi_altix_nas |
../../../tools/genmake2 -mods=../code -of ../code/linux_ia64_ifort+mpi_altix_nas |
83 |
make depend |
make depend |
84 |
make -j |
make -j |
85 |
|
|
125 |
to use matlab to look at the output |
to use matlab to look at the output |
126 |
matlab |
matlab |
127 |
|
|
128 |
for fld={'VICEtave','UICEtave','HEFFtave','AREAtave','QSWtave', ... |
for fld={'SIuice','SIvice','SIheff','SIarea','oceQsw', ... |
129 |
'QNETtave','FVtave','FUtave','EmPmRtave'} |
'surForcT','oceTAUX','oceTAUY','surForcS'} |
130 |
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]); |
131 |
cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause |
cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause |
132 |
end |
end |
133 |
|
|
134 |
heff=permute(readbin(['HEFFtave.0000000020.data'],[32,6,32]),[1 3 2]); |
heff=permute(readbin(['SIheff.0000000072.data'],[32,6,32]),[1 3 2]); |
135 |
mask=0*heff; mask(find(heff>.001))=1; |
mask=0*heff; mask(find(heff>.001))=1; |
136 |
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]); |
137 |
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]); |
138 |
clf |
clf |
139 |
subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3') |
subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3') |
140 |
subplot(222), mypcolor(uice(:,:,6)'); colorbar, title('uice, tile 6') |
subplot(222), mypcolor(uice(:,:,6)'); colorbar, title('uice, tile 6') |