This configuration is used for the ASOR project.  Please contact Jinlun Zhang
(zhang@apl.washington.edu) before using this configuration.

---------------------------------------------------

Configuring 1-cpu, 40-km regional Arctic model.

cvs co MITgcm_code
cd MITgcm
cvs co -d arctic40km MITgcm_contrib/arctic40km
cd arctic40km/input
wget -nd -r -A"*.bin" http://mit.ecco-group.org/datasets/MITgcm_input_files/MITgcm_contrib/arctic40km/input
wget -nd -r -A"topog.bin_130_102" http://mit.ecco-group.org/datasets/MITgcm_input_files/MITgcm_contrib/arctic40km/input
wget -nd -r -A"*.phc_130_102_25_12" http://mit.ecco-group.org/datasets/MITgcm_input_files/MITgcm_contrib/arctic40km/input
cd ../..
mkdir bin exe
cd bin
\rm *
\cp ../arctic40km/code/* .

>>> nireas
../tools/genmake2

>>> columbia
../tools/genmake2 -of ../tools/build_options/linux_ia64_ifort+mpi_altix_nas

make depend
make -j
cd ../exe
\rm *
\cp ../arctic40km/input/* .
../bin/mitgcmuv >& output.txt &

% looking at output files using matlab
nx=130; ny=102; tm=180;
for fld={'Eta','HEFF','sFluxtave','tFluxtave','uFluxtave','vFluxtave',}
 fn=[fld{1} '.' myint2str(tm,10) '.data'];
 fld=readbin(fn,[nx ny]);
 mypcolor(fld');
 colorbar
 title(fn)
 pause
end

% looking at input/output files using matlab
nx=130; ny=102; tm=36;
fn1='gairx_130_102.1979';
fn2='uFluxtave.0000000036.data';
fn3='gairy_130_102.1979';
fn4='vFluxtave.0000000036.data';
fld1=readbin(fn1,[nx ny]);
fld2=readbin(fn2,[nx ny]);
fld3=readbin(fn3,[nx ny]);
fld4=readbin(fn4,[nx ny]);
clf,cx1=[-1 1]*10;cx2=[-1 1]/4;
subplot(221),mypcolor(fld1');colorbar,title(fn1),caxis(cx1)
subplot(222),mypcolor(fld2');colorbar,title(fn2),caxis(cx2)
subplot(223),mypcolor(fld3');colorbar,title(fn3),caxis(cx1)
subplot(224),mypcolor(fld4');colorbar,title(fn4),caxis(cx2)

nx=130; ny=102; tm=36;
for fn={'topog.bin_130_102','theta.bin_130_102', ...
  'salt.bin_130_102','tair_130_102.1979', ...
  'qa_130_102.1979','prate_130_102', ...
  'fsh_130_102.1979','flo_130_102.1979', ...
  'SSS.bin_130_102','SST.bin_130_102'}
 fld=readbin(fn{1},[nx ny]);
 mypcolor(fld'); colorbar
 title(fn), pause
end