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

Annotation of /MITgcm_contrib/arctic40km/README

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


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

1 dimitri 1.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