| 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 |
|
January 5, 2006: cube51, like cube50 but with bottomDragQuadratic = 0.0012 |
| 12 |
|
|
| 13 |
|
January 5, 2006: cube50, like cube48 but with snow advection |
| 14 |
|
|
| 15 |
|
December 30, 2006: updates for cube49, the first |
| 16 |
|
"blend" integration. |
| 17 |
|
|
| 18 |
|
December 11, 2006: preparing for cube47, same as cube43 |
| 19 |
|
but with Arctic monthly river runoff climatology. |
| 20 |
|
|
| 21 |
|
December 5, 2006: 216-CPU code for cube46 |
| 22 |
|
Changed diagnostics to accomodate JMCs latest and some |
| 23 |
|
changes to code-mods/CPP_EEOPTIONS.h and W2_OPTIONS.h |
| 24 |
|
to speed up code. |
| 25 |
|
|
| 26 |
|
November 7, 2006: 216-CPU code for cube43 |
| 27 |
|
Bug fixes in growth.F routine and new set of diagnostics |
| 28 |
|
|
| 29 |
|
October 15, 2006: 216-CPU code for cube40 |
| 30 |
|
|
| 31 |
September 28, 2006: 216-CPU code for cube38 |
September 28, 2006: 216-CPU code for cube38 |
| 32 |
|
|
| 33 |
September 4, 2006 |
September 4, 2006 |
| 35 |
verification/tutorial_held_suarez_cs/input instead of |
verification/tutorial_held_suarez_cs/input instead of |
| 36 |
tile* files from verification/global_ocean.cs32x15/input |
tile* files from verification/global_ocean.cs32x15/input |
| 37 |
|
|
| 38 |
August 8, 2006 |
August 8, 2006: 216-CPU code and input used for cube37 |
|
216-CPU code and input used for integrating cube37 |
|
| 39 |
using "June 10, 2006" MITgcm code. |
using "June 10, 2006" MITgcm code. |
| 40 |
|
|
| 41 |
|
for online T/S profile diagnostics all that is needed is |
| 42 |
|
profiles_init_fixed |
| 43 |
|
profiles_inloop |
| 44 |
|
see pkg/ecco/ecco_cost_init_fixed.F |
| 45 |
|
pkg/ecco/cost_averagesfields.F |
| 46 |
|
|
| 47 |
|
########################################################################## |
| 48 |
|
########################################################################## |
| 49 |
|
# Get forcing files, etc. |
| 50 |
|
# Example below is for cube46 on /nobackup2a/menemenl |
| 51 |
|
|
| 52 |
|
cd /nobackup2a/menemenl |
| 53 |
|
mkdir cube47 |
| 54 |
|
cd cube47 |
| 55 |
|
ln -sf ../run_template . |
| 56 |
|
ln -sf ../ncep_rgau . |
| 57 |
|
|
| 58 |
|
|
| 59 |
########################################################################## |
########################################################################## |
| 60 |
########################################################################## |
########################################################################## |
| 61 |
# getting the code from anonymous CVS server |
# getting the code from anonymous CVS server |
| 62 |
|
|
| 63 |
|
bash or sh shell: |
| 64 |
|
$ export CVSROOT=':pserver:cvsanon@mitgcm.org:/u/gcmpack' |
| 65 |
|
$ cvs login |
| 66 |
|
( enter the CVS password: "cvsanon" ) |
| 67 |
|
|
| 68 |
|
tcsh or csh shell: |
| 69 |
|
$ setenv CVSROOT ':pserver:cvsanon@mitgcm.org:/u/gcmpack' |
| 70 |
|
$ cvs login |
| 71 |
|
( enter the CVS password: "cvsanon" ) |
| 72 |
|
|
| 73 |
cvs co MITgcm_contrib/high_res_cube/README_ice |
cvs co MITgcm_contrib/high_res_cube/README_ice |
| 74 |
cvs co MITgcm_contrib/high_res_cube/code-mods |
cvs co MITgcm_contrib/high_res_cube/code-mods |
| 75 |
cvs co MITgcm_contrib/high_res_cube/input |
cvs co MITgcm_contrib/high_res_cube/input |
| 81 |
|
|
| 82 |
########################################################################## |
########################################################################## |
| 83 |
########################################################################## |
########################################################################## |
| 84 |
|
# 510*510*6*50, 216-cpu cube sphere on altix with ifort and pkg/diagnostics |
| 85 |
|
|
| 86 |
|
cd MITgcm/verification/global_ocean.cs32x15/code |
| 87 |
|
\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* . |
| 88 |
|
\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s216t_85x85/* . |
| 89 |
|
\cp SIZE.h_216 SIZE.h |
| 90 |
|
\cp packages.conf_hr packages.conf |
| 91 |
|
|
| 92 |
|
cd ../build |
| 93 |
|
\rm * |
| 94 |
|
../../../tools/genmake2 -mods=../code -of ../code/linux_ia64_ifort+mpi_altix_nas |
| 95 |
|
make depend |
| 96 |
|
make -j |
| 97 |
|
|
| 98 |
|
cd .. |
| 99 |
|
mkdir run |
| 100 |
|
cd run |
| 101 |
|
\rm * |
| 102 |
|
\cp ../input/* . |
| 103 |
|
\cp ../../../../MITgcm_contrib/high_res_cube/input/* . |
| 104 |
|
ln -sf ../../../../run_template/* . |
| 105 |
|
\rm pickup* eedata.mth data.pkg data.seaice data.exf_clim data.exf data.mnc data pr* |
| 106 |
|
\rm bathy_Hmin50.bin *192_94* *cs32* *.m lev* tren* job_c* job_altix150 job_altix92_54t |
| 107 |
|
ln -sf ../../../../run_template/pickup.0000000216.cube49 pickup.0000000216 |
| 108 |
|
ln -sf ../../../../run_template/pickup_seaice.0000000216.cube48 pickup_seaice.0000000216 |
| 109 |
|
\cp ../build/mitgcmuv . |
| 110 |
|
qsub job_altix |
| 111 |
|
|
| 112 |
|
|
| 113 |
|
########################################################################## |
| 114 |
|
########################################################################## |
| 115 |
# 32*32*6*15, 1-cpu cube sphere |
# 32*32*6*15, 1-cpu cube sphere |
| 116 |
|
|
| 117 |
cd MITgcm/verification/global_ocean.cs32x15/code |
cd MITgcm/verification/global_ocean.cs32x15/code |
| 118 |
\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* . |
\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* . |
| 119 |
\cp ../../../utils/exch2/code-mods/s12t_16x32/* . |
\cp ../../../utils/exch2/code-mods/s12t_16x32/* . |
| 120 |
|
\rm CPP_EEOPTIONS.h |
| 121 |
cd ../build |
cd ../build |
| 122 |
\rm * |
\rm * |
| 123 |
../../../tools/genmake2 -mods=../code |
../../../tools/genmake2 -mods=../code |
| 138 |
to use matlab to look at the output |
to use matlab to look at the output |
| 139 |
matlab |
matlab |
| 140 |
|
|
| 141 |
for fld={'VICEtave','UICEtave','HEFFtave','AREAtave','QSWtave', ... |
for fld={'SIuice','SIvice','SIheff','SIarea','oceQsw', ... |
| 142 |
'QNETtave','FVtave','FUtave','EmPmRtave'} |
'surForcT','oceTAUX','oceTAUY','surForcS'} |
| 143 |
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]); |
| 144 |
cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause |
cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause |
| 145 |
end |
end |
| 146 |
|
|
| 147 |
heff=permute(readbin(['HEFFtave.0000000020.data'],[32,6,32]),[1 3 2]); |
heff=permute(readbin(['SIheff.0000000072.data'],[32,6,32]),[1 3 2]); |
| 148 |
mask=0*heff; mask(find(heff>.001))=1; |
mask=0*heff; mask(find(heff>.001))=1; |
| 149 |
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]); |
| 150 |
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]); |
| 151 |
clf |
clf |
| 152 |
subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3') |
subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3') |
| 153 |
subplot(222), mypcolor(uice(:,:,6)'); colorbar, title('uice, tile 6') |
subplot(222), mypcolor(uice(:,:,6)'); colorbar, title('uice, tile 6') |
| 230 |
|
|
| 231 |
########################################################################## |
########################################################################## |
| 232 |
########################################################################## |
########################################################################## |
|
# 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* 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 |
|
|
|
|
|
|
|
|
########################################################################## |
|
|
########################################################################## |
|
| 233 |
# 510*510*6*50, 1500-tile, 150-cpu, excludes land |
# 510*510*6*50, 1500-tile, 150-cpu, excludes land |
| 234 |
# uses s1500t_17x51 configuration |
# uses s1500t_17x51 configuration |
| 235 |
|
|
| 259 |
qsub job_altix150 |
qsub job_altix150 |
| 260 |
|
|
| 261 |
|
|
|
|
|
|
|
|
| 262 |
########################################################################## |
########################################################################## |
| 263 |
########################################################################## |
########################################################################## |
| 264 |
# 510*510*6*50, 1500-tile, 375-cpu, excludes land |
# 510*510*6*50, 1500-tile, 375-cpu, excludes land |
| 291 |
job_c20_375cpu < /dev/null >&! output.run & |
job_c20_375cpu < /dev/null >&! output.run & |
| 292 |
|
|
| 293 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| 294 |
########################################################################## |
########################################################################## |
| 295 |
########################################################################## |
########################################################################## |
| 296 |
# 510*510*6*50, 1500-tile, 375-cpu, excludes land |
# 510*510*6*50, 1500-tile, 375-cpu, excludes land |
| 355 |
job_c19_500cpu < /dev/null >&! output.run & |
job_c19_500cpu < /dev/null >&! output.run & |
| 356 |
|
|
| 357 |
|
|
|
|
|
|
|
|
|
|
|
| 358 |
########################################################################## |
########################################################################## |
| 359 |
########################################################################## |
########################################################################## |
| 360 |
# 510*510*6*50, 216-cpu cube sphere on cosmos |
# 510*510*6*50, 216-cpu cube sphere on cosmos |
| 390 |
bsub < job_cosmos |
bsub < job_cosmos |
| 391 |
|
|
| 392 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| 393 |
########################################################################## |
########################################################################## |
| 394 |
########################################################################## |
########################################################################## |
| 395 |
# 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 |