/[MITgcm]/MITgcm_contrib/arctic40km/README
ViewVC logotype

Contents of /MITgcm_contrib/arctic40km/README

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.1 - (show annotations) (download)
Thu Sep 1 14:16:40 2005 UTC (19 years, 10 months ago) by dimitri
Branch: MAIN
adding MITgcm_contrib/arctic40km

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

  ViewVC Help
Powered by ViewVC 1.1.22