Instructions for setting up a cube-sphere integration
with sea-ice on the cs32 grid.  Is meant to be used as
a preliminary sanity check for the cs510 configuration.

##########################################################################
# getting the code from anonymous CVS server

cvs co -P MITgcm_contrib/high_res_cube/README.cs32
cvs co -P MITgcm_contrib/high_res_cube/code-mods
cvs co -P MITgcm_contrib/high_res_cube/input
cvs co -P MITgcm_contrib/high_res_cube/results
cvs co -P MITgcm_code
cvs co -P MITgcm/verification/global_ocean.cs32x15/input
cvs co -P MITgcm/verification/tutorial_held_suarez_cs/input

##########################################################################
# 32*32*6*15, 1-cpu cube sphere

cd MITgcm
mkdir build run
cd build
\rm *
\cp ../../MITgcm_contrib/high_res_cube/code-mods/* .
\rm CPP_EEOPTIONS.h
../tools/genmake2
make depend
make -j
cd ../run
\rm *
\cp ../verification/tutorial_held_suarez_cs/input/grid* .
\cp ../verification/global_ocean.cs32x15/input/lev_*_cs_15k.bin .
\cp ../verification/global_ocean.cs32x15/input/bathy_Hmin50.bin .
\cp ../../MITgcm_contrib/high_res_cube/input/* .
../build/mitgcmuv >& output.txt

emacs output.txt ../../MITgcm_contrib/high_res_cube/results/output.txt

to use matlab to look at the output:
matlab
for fld={'SIuice','SIvice','SIheff','SIarea','oceQsw', ...
         'surForcT','oceTAUX','oceTAUY','surForcS'}
    tmp=permute(readbin([fld{1} '.0000000072.data'],[32,6,32]),[1 3 2]);
    cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
end

heff=permute(readbin(['SIheff.0000000072.data'],[32,6,32]),[1 3 2]);
mask=0*heff; mask(find(heff>.001))=1;
uice=mask.*permute(readbin(['SIuice.0000000072.data'],[32,6,32]),[1 3 2]);
vice=mask.*permute(readbin(['SIvice.0000000072.data'],[32,6,32]),[1 3 2]);
clf
subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3')
subplot(222), mypcolor(uice(:,:,6)'); colorbar, title('uice, tile 6')
subplot(223), mypcolor(vice(:,:,3)'); colorbar, title('vice, tile 3')
subplot(224), mypcolor(vice(:,:,6)'); colorbar, title('vice, tile 6')


##########################################################################
# 32*32*6*15, 176-tile, 1-cpu cube sphere experiment
# example that has 16 blank tiles

cd ../build
\rm *
\cp ../../MITgcm_contrib/high_res_cube/code-mods/* .
\cp ../utils/exch2/code-mods/s176t_8x4/SIZE.h .
ex - SIZE.h >> /dev/null <<EOF
/OLx =
c
     &           OLx =   5,
.
/OLy =
c
     &           OLy =   5,
.
w
q
EOF
\rm CPP_EEOPTIONS.h
../tools/genmake2
make depend
make -j
cd ..
mv run run_orig
mkdir run
cd run
\cp ../verification/tutorial_held_suarez_cs/input/grid* .
\cp ../verification/global_ocean.cs32x15/input/lev_*_cs_15k.bin .
\cp ../verification/global_ocean.cs32x15/input/bathy_Hmin50.bin .
\cp ../utils/exch2/code-mods/data.exch2.16_blk data.exch2
\cp ../../MITgcm_contrib/high_res_cube/input/* .
../build/mitgcmuv >& output.txt

to use matlab to compare to previous solution:
matlab
for f={'ETAN','oceTAUY','SIhsnow','surForcS','KPPhbl','PHIBOT', ...
       'SIuice','surForcT','MXLDEPTH','SIarea','SIvice','UVEL_k1', ...
       'oceQsw','SIheff','SSS','VVEL_k1','oceTAUX','SIhsalt','SST'}
 t1=readbin(['../run_orig/' f{1} '.0000000072.data'],[32*6 32]);
 t2=readbin([f{1} '.0000000072.data'],[32*6 32]);
 clf, subplot(311), mypcolor(t1'); thincolorbar, title(f{1})
 subplot(312), mypcolor(t2'); thincolorbar
 subplot(313), mypcolor(t2'-t1'); thincolorbar, pause
end
for f={'YG','YC','XG','XC','RAZ','RAW','RAS','RAC','hFacW','hFacS', ...
     'hFacC','Depth','DYG','DYC','DXG','DXC','AngleSN','AngleCS'}
 t1=readbin(['../run_orig/' f{1} '.data'],[32*6 32]);
 t2=readbin([f{1} '.data'],[32*6 32]);
 clf, subplot(311), mypcolor(t1'); thincolorbar, title(f{1})
 subplot(312), mypcolor(t2'); thincolorbar
 subplot(313), mypcolor(t2'-t1'); thincolorbar, pause
end


##########################################################################
# 32*32*6*15, 2-cpu cube sphere on the NAS altices with ifort

cd MITgcm
mkdir build run
cd build
\rm *
\cp ../../MITgcm_contrib/high_res_cube/code-mods/* .
\mv SIZE.h_mpi SIZE.h
../tools/genmake2 -of ../../MITgcm_contrib/high_res_cube/code-mods/linux_ia64_ifort+mpi_altix_nas
make depend
make -j 8
cd ../run
\rm *
\cp ../verification/tutorial_held_suarez_cs/input/grid* .
\cp ../verification/global_ocean.cs32x15/input/lev_*_cs_15k.bin .
\cp ../verification/global_ocean.cs32x15/input/bathy_Hmin50.bin .
\cp ../../MITgcm_contrib/high_res_cube/input/* .
mpirun -np 2 ../build/mitgcmuv


##########################################################################
# 32*32*6*15, 2-cpu cube sphere on the JPL altices with ifort

cd MITgcm
mkdir build run
cd build
\rm *
\cp ../../MITgcm_contrib/high_res_cube/code-mods/* .
\mv SIZE.h_mpi SIZE.h
../tools/genmake2 -of ../tools/build_options/linux_ia64_ifort+mpi_altix_jpl
make depend
make -j
cd ../run
\rm *
\cp ../verification/tutorial_held_suarez_cs/input/grid* .
\cp ../verification/global_ocean.cs32x15/input/lev_*_cs_15k.bin .
\cp ../verification/global_ocean.cs32x15/input/bathy_Hmin50.bin .
\cp ../../MITgcm_contrib/high_res_cube/input/* .
mpirun -np 2 ../build/mitgcmuv