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