1 |
dimitri |
1.2 |
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 |
dimitri |
1.1 |
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 |
dimitri |
1.4 |
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 |
dimitri |
1.1 |
cd ../.. |
16 |
|
|
mkdir bin exe |
17 |
|
|
cd bin |
18 |
|
|
\rm * |
19 |
|
|
\cp ../arctic40km/code/* . |
20 |
dimitri |
1.3 |
|
21 |
|
|
>>> nireas |
22 |
dimitri |
1.1 |
../tools/genmake2 |
23 |
dimitri |
1.3 |
|
24 |
|
|
>>> columbia |
25 |
|
|
../tools/genmake2 -of ../tools/build_options/linux_ia64_ifort+mpi_altix_nas |
26 |
|
|
|
27 |
dimitri |
1.1 |
make depend |
28 |
dimitri |
1.3 |
make -j |
29 |
dimitri |
1.1 |
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 |