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