/[MITgcm]/MITgcm_contrib/lab_sea_test/read.me
ViewVC logotype

Diff of /MITgcm_contrib/lab_sea_test/read.me

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

revision 1.1 by dimitri, Tue Jun 29 20:51:16 2004 UTC revision 1.2 by dimitri, Mon Jul 5 14:15:44 2004 UTC
# Line 1  Line 1 
1  test for lab_sea forward model linearity  
2    Test for lab_sea forward model linearity.
3    
4    cd MITgcm/verification/lab_sea
5    mkdir build
6    cd build
7    cp ../code/*.h ../code/packages.conf .
8    ../../../tools/genmake2
9    make depend
10    make
11    cd ..; mkdir run; cd run; cp ../input/* .
12    cp ../../../../MITgcm_contrib/lab_sea_test/* .
13    ../build/mitgcmuv > output.txt &
14    
15    cd ..; mkdir run1; cd run1; cp ../input/* .
16    cp ../../../../MITgcm_contrib/lab_sea_test/* .
17    mv data.exf tmp
18    sed s/tair/tair01/ tmp > data.exf
19    ../build/mitgcmuv > output.txt &
20    
21    cd ..; mkdir run2; cd run2; cp ../input/* .
22    cp ../../../../MITgcm_contrib/lab_sea_test/* .
23    mv data.exf tmp
24    sed s/tair/tair02/ tmp > data.exf
25    ../build/mitgcmuv > output.txt &
26    
27    cd ..; mkdir run3; cd run3; cp ../input/* .
28    cp ../../../../MITgcm_contrib/lab_sea_test/* .
29    mv data.exf tmp
30    sed s/tair/tair03/ tmp > data.exf
31    ../build/mitgcmuv > output.txt &
32    
33    cd ..; mkdir run4; cd run4; cp ../input/* .
34    cp ../../../../MITgcm_contrib/lab_sea_test/* .
35    mv data.exf tmp
36    sed s/tair/tair04/ tmp > data.exf
37    ../build/mitgcmuv > output.txt &
38    
39    cd ..; mkdir run5; cd run5; cp ../input/* .
40    cp ../../../../MITgcm_contrib/lab_sea_test/* .
41    mv data.exf tmp
42    sed s/tair/tair05/ tmp > data.exf
43    ../build/mitgcmuv > output.txt &
44    
45    cd ..; mkdir run6; cd run6; cp ../input/* .
46    cp ../../../../MITgcm_contrib/lab_sea_test/* .
47    mv data.exf tmp
48    sed s/tair/tair06/ tmp > data.exf
49    ../build/mitgcmuv > output.txt &
50    
51    cd ..; mkdir run7; cd run7; cp ../input/* .
52    cp ../../../../MITgcm_contrib/lab_sea_test/* .
53    mv data.exf tmp
54    sed s/tair/tair07/ tmp > data.exf
55    ../build/mitgcmuv > output.txt &
56    
57    cd ..; mkdir run8; cd run8; cp ../input/* .
58    cp ../../../../MITgcm_contrib/lab_sea_test/* .
59    mv data.exf tmp
60    sed s/tair/tair08/ tmp > data.exf
61    ../build/mitgcmuv > output.txt &
62    
63    cd ..; mkdir run9; cd run9; cp ../input/* .
64    cp ../../../../MITgcm_contrib/lab_sea_test/* .
65    mv data.exf tmp
66    sed s/tair/tair09/ tmp > data.exf
67    ../build/mitgcmuv > output.txt &
68    
69    cd ..; mkdir run10; cd run10; cp ../input/* .
70    cp ../../../../MITgcm_contrib/lab_sea_test/* .
71    mv data.exf tmp
72    sed s/tair/tair10/ tmp > data.exf
73    ../build/mitgcmuv > output.txt &
74    
75    cd ..; mkdir run11; cd run11; cp ../input/* .
76    cp ../../../../MITgcm_contrib/lab_sea_test/* .
77    mv data.exf tmp
78    sed s/tair/tair11/ tmp > data.exf
79    ../build/mitgcmuv > output.txt &
80    
81    cd ..; mkdir run12; cd run12; cp ../input/* .
82    cp ../../../../MITgcm_contrib/lab_sea_test/* .
83    mv data.exf tmp
84    sed s/tair/tair12/ tmp > data.exf
85    ../build/mitgcmuv > output.txt &
86    
87    cd ..; mkdir run13; cd run13; cp ../input/* .
88    cp ../../../../MITgcm_contrib/lab_sea_test/* .
89    mv data.exf tmp
90    sed s/tair/tair13/ tmp > data.exf
91    ../build/mitgcmuv > output.txt &
92    
93    cd ..; mkdir run14; cd run14; cp ../input/* .
94    cp ../../../../MITgcm_contrib/lab_sea_test/* .
95    mv data.exf tmp
96    sed s/tair/tair14/ tmp > data.exf
97    ../build/mitgcmuv > output.txt &
98    
99    cd ..; mkdir run15; cd run15; cp ../input/* .
100    cp ../../../../MITgcm_contrib/lab_sea_test/* .
101    mv data.exf tmp
102    sed s/tair/tair15/ tmp > data.exf
103    ../build/mitgcmuv > output.txt &
104    
105    cd ..; mkdir run16; cd run16; cp ../input/* .
106    cp ../../../../MITgcm_contrib/lab_sea_test/* .
107    mv data.exf tmp
108    sed s/tair/tair16/ tmp > data.exf
109    ../build/mitgcmuv > output.txt &
110    
111    cd ..; mkdir run17; cd run17; cp ../input/* .
112    cp ../../../../MITgcm_contrib/lab_sea_test/* .
113    mv data.exf tmp
114    sed s/tair/tair17/ tmp > data.exf
115    ../build/mitgcmuv > output.txt &
116    
117    cd ..; mkdir run18; cd run18; cp ../input/* .
118    cp ../../../../MITgcm_contrib/lab_sea_test/* .
119    mv data.exf tmp
120    sed s/tair/tair18/ tmp > data.exf
121    ../build/mitgcmuv > output.txt &
122    
123    cd ..; mkdir run19; cd run19; cp ../input/* .
124    cp ../../../../MITgcm_contrib/lab_sea_test/* .
125    mv data.exf tmp
126    sed s/tair/tair19/ tmp > data.exf
127    ../build/mitgcmuv > output.txt &
128    
129    cd ..; mkdir run20; cd run20; cp ../input/* .
130    cp ../../../../MITgcm_contrib/lab_sea_test/* .
131    mv data.exf tmp
132    sed s/tair/tair20/ tmp > data.exf
133    ../build/mitgcmuv > output.txt &
134    
135    cd ..
136    
137    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
138    
139    matlab
140    
141    pn='nokpp_nogm_withseaicedyn/';
142    
143    heff=zeros(20,16,21,21);
144    area=zeros(20,16,21,21);
145    for e=0:20
146     for t=0:4, disp([e t])
147      if e==0, fn=[pn 'run/HEFF.'];
148      else fn=[pn 'run' int2str(e) '/HEFF.']; end
149      fn=[fn myint2str(t*240,10) '.data'];
150      heff(:,:,t+1,e+1)=readbin(fn,[20 16],1);
151      if e==0, fn=[pn 'run/AREA.'];
152      else fn=[pn 'run' int2str(e) '/AREA.']; end
153      fn=[fn myint2str(t*240,10) '.data'];
154      area(:,:,t+1,e+1)=readbin(fn,[20 16],1);
155     end
156    end
157    
158    for t=1:21
159    tmp=squeeze(sum(sum(heff(:,:,t,:))));
160    clf,subplot(211),plot(tmp-mean(tmp))
161    grid,title([int2str((t-1)*10) ' days'])
162    tmp=squeeze(sum(sum(area(:,:,t,:))));
163    subplot(212),plot(tmp-mean(tmp))
164    grid,title([int2str((t-1)*10) ' days'])
165    pause, end
166    
167    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
168    
169    repeat with daily dumps for experiments 10-15, days 6-10
170    
171    cd run10; mv data tmp
172    sed s/startTime=0.0,/startTime=518400.0,/ tmp > tmp2
173    sed s/endTime=1728000.,/endTime=950400.0,/ tmp2 > tmp3
174    sed s/"dumpFreq = 86400.,"/dumpFreq=3600.,/ tmp3 > data
175    ../build/mitgcmuv >! output.txt &
176    
177    cd ../run11; mv data tmp
178    sed s/startTime=0.0,/startTime=518400.0,/ tmp > tmp2
179    sed s/endTime=1728000.,/endTime=950400.0,/ tmp2 > tmp3
180    sed s/"dumpFreq = 86400.,"/dumpFreq=3600.,/ tmp3 > data
181    ../build/mitgcmuv >! output.txt &
182    
183    cd ../run12; mv data tmp
184    sed s/startTime=0.0,/startTime=518400.0,/ tmp > tmp2
185    sed s/endTime=1728000.,/endTime=950400.0,/ tmp2 > tmp3
186    sed s/"dumpFreq = 86400.,"/dumpFreq=3600.,/ tmp3 > data
187    ../build/mitgcmuv >! output.txt &
188    
189    cd ../run13; mv data tmp
190    sed s/startTime=0.0,/startTime=518400.0,/ tmp > tmp2
191    sed s/endTime=1728000.,/endTime=950400.0,/ tmp2 > tmp3
192    sed s/"dumpFreq = 86400.,"/dumpFreq=3600.,/ tmp3 > data
193    ../build/mitgcmuv >! output.txt &
194    
195    cd ../run14; mv data tmp
196    sed s/startTime=0.0,/startTime=518400.0,/ tmp > tmp2
197    sed s/endTime=1728000.,/endTime=950400.0,/ tmp2 > tmp3
198    sed s/"dumpFreq = 86400.,"/dumpFreq=3600.,/ tmp3 > data
199    ../build/mitgcmuv >! output.txt &
200    
201    cd ../run15; mv data tmp
202    sed s/startTime=0.0,/startTime=518400.0,/ tmp > tmp2
203    sed s/endTime=1728000.,/endTime=950400.0,/ tmp2 > tmp3
204    sed s/"dumpFreq = 86400.,"/dumpFreq=3600.,/ tmp3 > data
205    ../build/mitgcmuv >! output.txt &
206    
207    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
208    
209    matlab
210    
211    heff=zeros(20,16,121,6);
212    area=zeros(20,16,121,6);
213    for e=10:15, disp(e)
214     for t=144:264
215      fn=['run' int2str(e) '/HEFF.'];
216      fn=[fn myint2str(t,10) '.data'];
217      heff(:,:,t-143,e-9)=readbin(fn,[20 16],1);
218      fn=['run' int2str(e) '/AREA.'];
219      fn=[fn myint2str(t,10) '.data'];
220      area(:,:,t-143,e-9)=readbin(fn,[20 16],1);
221     end
222    end
223    
224    for t=144:264
225    tmp=squeeze(sum(sum(heff(:,:,t-143,:))));
226    clf,subplot(211), plot(tmp-mean(tmp)), grid, title(t/24)
227    tmp=squeeze(sum(sum(area(:,:,t-143,:))));
228    subplot(212), plot(tmp-mean(tmp)), grid, title(t/24)
229    pause, end
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22