Instructions for setting up a cube-sphere integration with sea-ice.
Date last tested: November 21, 2003.

setenv CVSROOT :pserver:cvsanon@mitgcm.org:/u/u0/gcmpack
cvs login ( CVS password: cvsanon )

cvs co MITgcm_contrib/high_res_cube/README_ice
cvs co MITgcm_contrib/high_res_cube/code-mods
cvs co MITgcm_contrib/high_res_cube/matlab-grid-converter
cvs co MITgcm_contrib/high_res_cube/input
cvs co MITgcm_contrib/high_res_cube/results

cvs co -r checkpoint52b_post MITgcm_code
cvs co -r checkpoint52b_post MITgcm/verification/global_ocean.cs32x15

cd MITgcm/verification/global_ocean.cs32x15/code
cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s12t_16x32/* .
\rm mdsio_readfield.F

cd ..
mkdir build
cd build
../code/genmake2 -mods=../code
make depend
make

cd ..
mkdir run
cd run
cp ../inp_thsice/* .
cp ../input/* .
cp ../../../../MITgcm_contrib/high_res_cube/matlab-grid-converter/tile00* .
cp ../../../../MITgcm_contrib/high_res_cube/input/* .
../build/mitgcmuv >& output.txt

comparison output is in:
../../../../MITgcm_contrib/high_res_cube/results/output.txt


to use matlab to look at the output
matlab
for fld={'Eta','VICE','UICE','HEFF','AREA','Qsw','Qnet','FV','FU','EmPmR'}
    tmp=permute(readbin([fld{1} '.0000000020.data'],[32,6,32]),[1 3 2]);
    cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
end