| 1 | 
Configuring 1-cpu, 40-km regional Arctic model. | 
| 2 | 
 | 
| 3 | 
cvs co MITgcm_code | 
| 4 | 
cd MITgcm | 
| 5 | 
cvs co -d arctic40km MITgcm_contrib/arctic40km | 
| 6 | 
cd arctic40km/input | 
| 7 | 
wget -nd -r -A"*.bin" http://channel.mit.edu/datasets/MITgcm_input_files/MITgcm_contrib/arctic40km/input | 
| 8 | 
wget -nd -r -A"*130_102*" http://channel.mit.edu/datasets/MITgcm_input_files/MITgcm_contrib/arctic40km/input | 
| 9 | 
cd ../.. | 
| 10 | 
mkdir bin exe | 
| 11 | 
cd bin | 
| 12 | 
\rm * | 
| 13 | 
\cp ../arctic40km/code/* . | 
| 14 | 
../tools/genmake2 | 
| 15 | 
make depend | 
| 16 | 
make | 
| 17 | 
cd ../exe | 
| 18 | 
\rm * | 
| 19 | 
\cp ../arctic40km/input/* . | 
| 20 | 
../bin/mitgcmuv >& output.txt & | 
| 21 | 
 | 
| 22 | 
% looking at output files using matlab | 
| 23 | 
nx=130; ny=102; tm=180; | 
| 24 | 
for fld={'Eta','HEFF','sFluxtave','tFluxtave','uFluxtave','vFluxtave',} | 
| 25 | 
 fn=[fld{1} '.' myint2str(tm,10) '.data']; | 
| 26 | 
 fld=readbin(fn,[nx ny]); | 
| 27 | 
 mypcolor(fld'); | 
| 28 | 
 colorbar | 
| 29 | 
 title(fn) | 
| 30 | 
 pause | 
| 31 | 
end | 
| 32 | 
 | 
| 33 | 
% looking at input/output files using matlab | 
| 34 | 
nx=130; ny=102; tm=36; | 
| 35 | 
fn1='gairx_130_102.1979'; | 
| 36 | 
fn2='uFluxtave.0000000036.data'; | 
| 37 | 
fn3='gairy_130_102.1979'; | 
| 38 | 
fn4='vFluxtave.0000000036.data'; | 
| 39 | 
fld1=readbin(fn1,[nx ny]); | 
| 40 | 
fld2=readbin(fn2,[nx ny]); | 
| 41 | 
fld3=readbin(fn3,[nx ny]); | 
| 42 | 
fld4=readbin(fn4,[nx ny]); | 
| 43 | 
clf,cx1=[-1 1]*10;cx2=[-1 1]/4; | 
| 44 | 
subplot(221),mypcolor(fld1');colorbar,title(fn1),caxis(cx1) | 
| 45 | 
subplot(222),mypcolor(fld2');colorbar,title(fn2),caxis(cx2) | 
| 46 | 
subplot(223),mypcolor(fld3');colorbar,title(fn3),caxis(cx1) | 
| 47 | 
subplot(224),mypcolor(fld4');colorbar,title(fn4),caxis(cx2) | 
| 48 | 
 | 
| 49 | 
nx=130; ny=102; tm=36; | 
| 50 | 
for fn={'topog.bin_130_102','theta.bin_130_102', ... | 
| 51 | 
  'salt.bin_130_102','tair_130_102.1979', ... | 
| 52 | 
  'qa_130_102.1979','prate_130_102', ... | 
| 53 | 
  'fsh_130_102.1979','flo_130_102.1979', ... | 
| 54 | 
  'SSS.bin_130_102','SST.bin_130_102'} | 
| 55 | 
 fld=readbin(fn{1},[nx ny]); | 
| 56 | 
 mypcolor(fld'); colorbar | 
| 57 | 
 title(fn), pause | 
| 58 | 
end |