/[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.4 by dimitri, Sun Jul 11 06:18:03 2004 UTC revision 1.5 by dimitri, Mon Jul 12 01:00:20 2004 UTC
# Line 1  Line 1 
1  Test for lab_sea forward model linearity.  Tests for lab_sea forward model linearity.
2    
3    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5    
6  cd MITgcm/verification/lab_sea  cd MITgcm/verification/lab_sea
7  mkdir build  mkdir build
8  cd build  cd build
9  cp ../code/*.h ../code/*.F ../code/packages.conf .  cp ../code/*.h ../code/packages.conf .
10    cp ../../../../MITgcm_contrib/lab_sea_test/*.F .
11  ../../../tools/genmake2  ../../../tools/genmake2
12  make depend  make depend
13  make  make
14  cd ..; mkdir run; cd run; cp ../input/* .  cd ..; mkdir run; cd run; cp ../input/* .
15  cp ../../../../MITgcm_contrib/lab_sea_test/* .  cp ../../../../MITgcm_contrib/lab_sea_test/* .
16  ../build/mitgcmuv >! output.txt &  ../build/mitgcmuv >& output.txt &
17    
18  cd ..; mkdir run1; cd run1; cp ../input/* .  cd ..; mkdir run1; cd run1; cp ../input/* .
19  cp ../../../../MITgcm_contrib/lab_sea_test/* .  cp ../../../../MITgcm_contrib/lab_sea_test/* .
20  mv data.exf tmp  mv data.exf tmp
21  sed s/tair/tair01/ tmp > data.exf  sed s/tair/tair01/ tmp > data.exf
22  ../build/mitgcmuv >! output.txt &  ../build/mitgcmuv >& output.txt &
23    
24  cd ..; mkdir run2; cd run2; cp ../input/* .  cd ..; mkdir run2; cd run2; cp ../input/* .
25  cp ../../../../MITgcm_contrib/lab_sea_test/* .  cp ../../../../MITgcm_contrib/lab_sea_test/* .
26  mv data.exf tmp  mv data.exf tmp
27  sed s/tair/tair02/ tmp > data.exf  sed s/tair/tair02/ tmp > data.exf
28  ../build/mitgcmuv >! output.txt &  ../build/mitgcmuv >& output.txt &
29    
30    cd ..
31    matlab
32    
33    clear all, clf reset, pn=''; ne=2; nt=53;
34    heff=zeros(20,16,nt,ne+1);
35    area=zeros(20,16,nt,ne+1);
36    for e=0:ne
37     for t=1:nt, disp([e t])
38      if e==0, fn=[pn 'run/HEFFtave.'];
39      else fn=[pn 'run' int2str(e) '/HEFFtave.']; end
40      fn=[fn myint2str(t*240,10) '.data'];
41      heff(:,:,t+1,e+1)=readbin(fn,[20 16],1);
42      if e==0, fn=[pn 'run/AREAtave.'];
43      else fn=[pn 'run' int2str(e) '/AREAtave.']; end
44      fn=[fn myint2str(t*240,10) '.data'];
45      area(:,:,t+1,e+1)=readbin(fn,[20 16],1);
46     end
47    end
48    
49    clf reset, orient tall, wysiwyg
50    tmp1=squeeze(sum(sum(heff(:,:,:,1))));
51    tmp2=squeeze(sum(sum(heff(:,:,:,2))));
52    tmp3=squeeze(sum(sum(heff(:,:,:,3))));
53    t=1:length(tmp1); t=(t-1)*10
54    subplot(411),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
55    axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ...
56            max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid
57    title('sum(heff), original config')
58    legend('.01^oC','.02^oC',0)
59    subplot(412),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
60    axis([0 max(t) -2e-5 1e-5]), grid
61    title('sum(heff), original config')
62    legend('.01^oC','.02^oC',0)
63    tmp1=squeeze(sum(sum(area(:,:,:,1))));
64    tmp2=squeeze(sum(sum(area(:,:,:,2))));
65    tmp3=squeeze(sum(sum(area(:,:,:,3))));
66    subplot(413),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
67    axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ...
68            max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid
69    title('sum(area), original config')
70    legend('.01^oC','.02^oC',0)
71    subplot(414),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
72    axis([0 max(t) -2e-5 1e-5]), grid, xlabel('days')
73    title('sum(area), original config')
74    legend('.01^oC','.02^oC',0)
75    print -djpeg FIG_ynegminmaxrelax
76    
77  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
78  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
79    
80    cd MITgcm/verification/lab_sea
81    mkdir build
82    cd build
83    cp ../code/*.h ../code/packages.conf .
84    cp ../../../../MITgcm_contrib/lab_sea_test/growth.F .
85    ../../../tools/genmake2
86    make depend
87    make
88    cd ..; mkdir run; cd run; cp ../input/* .
89    cp ../../../../MITgcm_contrib/lab_sea_test/* .
90    mv data.pkg tmp
91    sed s/"GMRedi = .TRUE"/"GMRedi = .FALSE"/ tmp > tmp1
92    sed s/"KPP    = .TRUE"/"KPP    = .FALSE"/ tmp1 > data.pkg
93    mv data.seaice tmp
94    sed s/"DYNAMICS  = .TRUE"/"DYNAMICS  = .FALSE"/ tmp > data.seaice
95    ../build/mitgcmuv >& output.txt &
96    
97    cd ..; mkdir run1; cd run1; cp ../run/* .
98    mv data.exf tmp
99    sed s/tair/tair01/ tmp > data.exf
100    ../build/mitgcmuv >& output.txt &
101    
102    cd ..; mkdir run2; cd run2; cp ../run/* .
103    mv data.exf tmp
104    sed s/tair/tair02/ tmp > data.exf
105    ../build/mitgcmuv >& output.txt &
106    
107  cd ..  cd ..
108  matlab  matlab
109    
110  clear all, clf reset, pn=''; ne=5; nt=120;  clear all, pn=''; ne=2; nt=53;
111  heff=zeros(20,16,nt,ne+1);  heff=zeros(20,16,nt,ne+1);
112  area=zeros(20,16,nt,ne+1);  area=zeros(20,16,nt,ne+1);
113  for e=0:ne  for e=0:ne
# Line 45  for e=0:ne Line 123  for e=0:ne
123   end   end
124  end  end
125    
126  tmp=squeeze(sum(sum(sum(heff))));  clf reset, orient tall, wysiwyg
127  clf,subplot(211),plot(tmp-mean(tmp))  tmp1=squeeze(sum(sum(heff(:,:,:,1))));
128  grid,title([int2str((t-1)*10) ' days'])  tmp2=squeeze(sum(sum(heff(:,:,:,2))));
129  tmp=squeeze(sum(sum(sum(area))));  tmp3=squeeze(sum(sum(heff(:,:,:,3))));
130  subplot(212),plot(tmp-mean(tmp))  t=1:length(tmp1); t=(t-1)*10
131  grid,title([int2str((t-1)*10) ' days'])  subplot(411),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
132  tmp=squeeze(sum(sum(heff(:,:,:,1))));  axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ...
133  clf,subplot(211),plot(tmp-mean(tmp))          max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid
134  grid,title([int2str((t-1)*10) ' days'])  title('sum(heff), original config')
135  tmp=squeeze(sum(sum(area(:,:,:,1))));  legend('.01^oC','.02^oC',0)
136  subplot(212),plot(tmp-mean(tmp))  subplot(412),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
137  grid,title([int2str((t-1)*10) ' days'])  axis([0 max(t) -2e-5 1e-5]), grid
138    title('sum(heff), original config')
139    legend('.01^oC','.02^oC',0)
140    tmp1=squeeze(sum(sum(area(:,:,:,1))));
141    tmp2=squeeze(sum(sum(area(:,:,:,2))));
142    tmp3=squeeze(sum(sum(area(:,:,:,3))));
143    subplot(413),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
144    axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ...
145            max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid
146    title('sum(area), original config')
147    legend('.01^oC','.02^oC',0)
148    subplot(414),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
149    axis([0 max(t) -2e-5 1e-5]), grid, xlabel('days')
150    title('sum(area), original config')
151    legend('.01^oC','.02^oC',0)
152    print -djpeg FIG_ynegrelax_no_kpp_gm_dyn
153    
154    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
155    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
156    
157    cd MITgcm/verification/lab_sea
158    mkdir build
159    cd build
160    cp ../code/*.h ../code/packages.conf .
161    cp ../../../../MITgcm_contrib/lab_sea_test/growth.F .
162    ../../../tools/genmake2
163    make depend
164    make
165    cd ..; mkdir run; cd run; cp ../input/* .
166    cp ../../../../MITgcm_contrib/lab_sea_test/* .
167    ../build/mitgcmuv >& output.txt &
168    
169    cd ..; mkdir run1; cd run1; cp ../input/* .
170    cp ../../../../MITgcm_contrib/lab_sea_test/* .
171    mv data.exf tmp
172    sed s/tair/tair01/ tmp > data.exf
173    ../build/mitgcmuv >& output.txt &
174    
175    cd ..; mkdir run2; cd run2; cp ../input/* .
176    cp ../../../../MITgcm_contrib/lab_sea_test/* .
177    mv data.exf tmp
178    sed s/tair/tair02/ tmp > data.exf
179    ../build/mitgcmuv >& output.txt &
180    
181    cd ..
182    matlab
183    
184    clear all, clf reset, pn=''; ne=2; nt=53;
185    heff=zeros(20,16,nt,ne+1);
186    area=zeros(20,16,nt,ne+1);
187    for e=0:ne
188     for t=1:nt, disp([e t])
189      if e==0, fn=[pn 'run/HEFFtave.'];
190      else fn=[pn 'run' int2str(e) '/HEFFtave.']; end
191      fn=[fn myint2str(t*240,10) '.data'];
192      heff(:,:,t+1,e+1)=readbin(fn,[20 16],1);
193      if e==0, fn=[pn 'run/AREAtave.'];
194      else fn=[pn 'run' int2str(e) '/AREAtave.']; end
195      fn=[fn myint2str(t*240,10) '.data'];
196      area(:,:,t+1,e+1)=readbin(fn,[20 16],1);
197     end
198    end
199    
200    clf reset, orient tall, wysiwyg
201  tmp1=squeeze(sum(sum(heff(:,:,:,1))));  tmp1=squeeze(sum(sum(heff(:,:,:,1))));
202  tmp2=squeeze(sum(sum(heff(:,:,:,2))));  tmp2=squeeze(sum(sum(heff(:,:,:,2))));
203  tmp3=squeeze(sum(sum(heff(:,:,:,3))));  tmp3=squeeze(sum(sum(heff(:,:,:,3))));
204  t=1:length(tmp1); t=(t-1)*10  t=1:length(tmp1); t=(t-1)*10
205  clf,subplot(211),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)  subplot(411),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
206    axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ...
207            max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid
208    title('sum(heff), original config')
209    legend('.01^oC','.02^oC',0)
210    subplot(412),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
211    axis([0 max(t) -2e-5 1e-5]), grid
212    title('sum(heff), original config')
213    legend('.01^oC','.02^oC',0)
214    tmp1=squeeze(sum(sum(area(:,:,:,1))));
215    tmp2=squeeze(sum(sum(area(:,:,:,2))));
216    tmp3=squeeze(sum(sum(area(:,:,:,3))));
217    subplot(413),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
218    axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ...
219            max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid
220    title('sum(area), original config')
221    legend('.01^oC','.02^oC',0)
222    subplot(414),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
223    axis([0 max(t) -2e-5 1e-5]), grid, xlabel('days')
224    title('sum(area), original config')
225    legend('.01^oC','.02^oC',0)
226    print -djpeg FIG_ynegrelax
227    
228    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
229    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
230    
231    cd MITgcm/verification/lab_sea
232    mkdir build
233    cd build
234    cp ../code/*.h ../code/packages.conf .
235    ../../../tools/genmake2
236    make depend
237    make
238    cd ..; mkdir run; cd run; cp ../input/* .
239    cp ../../../../MITgcm_contrib/lab_sea_test/* .
240    mv data.pkg tmp
241    sed s/"GMRedi = .TRUE"/"GMRedi = .FALSE"/ tmp > tmp1
242    sed s/"KPP    = .TRUE"/"KPP    = .FALSE"/ tmp1 > data.pkg
243    mv data.seaice tmp
244    sed s/"DYNAMICS  = .TRUE"/"DYNAMICS  = .FALSE"/ tmp > data.seaice
245    ../build/mitgcmuv >& output.txt &
246    
247    cd ..; mkdir run1; cd run1; cp ../run/* .
248    mv data.exf tmp
249    sed s/tair/tair01/ tmp > data.exf
250    ../build/mitgcmuv >& output.txt &
251    
252    cd ..; mkdir run2; cd run2; cp ../run/* .
253    mv data.exf tmp
254    sed s/tair/tair02/ tmp > data.exf
255    ../build/mitgcmuv >& output.txt &
256    
257    cd ..
258    matlab
259    
260    clear all, clf reset, pn=''; ne=2; nt=53;
261    heff=zeros(20,16,nt,ne+1);
262    area=zeros(20,16,nt,ne+1);
263    for e=0:ne
264     for t=1:nt, disp([e t])
265      if e==0, fn=[pn 'run/HEFFtave.'];
266      else fn=[pn 'run' int2str(e) '/HEFFtave.']; end
267      fn=[fn myint2str(t*240,10) '.data'];
268      heff(:,:,t+1,e+1)=readbin(fn,[20 16],1);
269      if e==0, fn=[pn 'run/AREAtave.'];
270      else fn=[pn 'run' int2str(e) '/AREAtave.']; end
271      fn=[fn myint2str(t*240,10) '.data'];
272      area(:,:,t+1,e+1)=readbin(fn,[20 16],1);
273     end
274    end
275    
276    clf reset, orient tall, wysiwyg
277    tmp1=squeeze(sum(sum(heff(:,:,:,1))));
278    tmp2=squeeze(sum(sum(heff(:,:,:,2))));
279    tmp3=squeeze(sum(sum(heff(:,:,:,3))));
280    t=1:length(tmp1); t=(t-1)*10
281    subplot(411),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
282    axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ...
283            max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid
284    title('sum(heff), original config')
285    legend('.01^oC','.02^oC',0)
286    subplot(412),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
287    axis([0 max(t) -2e-5 1e-5]), grid
288    title('sum(heff), original config')
289    legend('.01^oC','.02^oC',0)
290    tmp1=squeeze(sum(sum(area(:,:,:,1))));
291    tmp2=squeeze(sum(sum(area(:,:,:,2))));
292    tmp3=squeeze(sum(sum(area(:,:,:,3))));
293    subplot(413),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
294    axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ...
295            max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid
296    title('sum(area), original config')
297    legend('.01^oC','.02^oC',0)
298    subplot(414),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
299    axis([0 max(t) -2e-5 1e-5]), grid, xlabel('days')
300    title('sum(area), original config')
301    legend('.01^oC','.02^oC',0)
302    print -djpeg FIG_no_kpp_gm_dyn
303    
304    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
305    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
306    
307    cd MITgcm/verification/lab_sea
308    mkdir build
309    cd build
310    cp ../code/*.h ../code/packages.conf .
311    ../../../tools/genmake2
312    make depend
313    make
314    cd ..; mkdir run; cd run; cp ../input/* .
315    cp ../../../../MITgcm_contrib/lab_sea_test/* .
316    ../build/mitgcmuv >& output.txt &
317    
318    cd ..; mkdir run1; cd run1; cp ../input/* .
319    cp ../../../../MITgcm_contrib/lab_sea_test/* .
320    mv data.exf tmp
321    sed s/tair/tair01/ tmp > data.exf
322    ../build/mitgcmuv >& output.txt &
323    
324    cd ..; mkdir run2; cd run2; cp ../input/* .
325    cp ../../../../MITgcm_contrib/lab_sea_test/* .
326    mv data.exf tmp
327    sed s/tair/tair02/ tmp > data.exf
328    ../build/mitgcmuv >& output.txt &
329    
330    cd ..
331    matlab
332    
333    clear all, clf reset, pn=''; ne=2; nt=53;
334    heff=zeros(20,16,nt,ne+1);
335    area=zeros(20,16,nt,ne+1);
336    for e=0:ne
337     for t=1:nt, disp([e t])
338      if e==0, fn=[pn 'run/HEFFtave.'];
339      else fn=[pn 'run' int2str(e) '/HEFFtave.']; end
340      fn=[fn myint2str(t*240,10) '.data'];
341      heff(:,:,t+1,e+1)=readbin(fn,[20 16],1);
342      if e==0, fn=[pn 'run/AREAtave.'];
343      else fn=[pn 'run' int2str(e) '/AREAtave.']; end
344      fn=[fn myint2str(t*240,10) '.data'];
345      area(:,:,t+1,e+1)=readbin(fn,[20 16],1);
346     end
347    end
348    
349    clf reset, orient tall, wysiwyg
350    tmp1=squeeze(sum(sum(heff(:,:,:,1))));
351    tmp2=squeeze(sum(sum(heff(:,:,:,2))));
352    tmp3=squeeze(sum(sum(heff(:,:,:,3))));
353    t=1:length(tmp1); t=(t-1)*10
354    subplot(411),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
355    axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ...
356            max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid
357    title('sum(heff), original config')
358    legend('.01^oC','.02^oC',0)
359    subplot(412),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
360    axis([0 max(t) -2e-5 1e-5]), grid
361    title('sum(heff), original config')
362    legend('.01^oC','.02^oC',0)
363    tmp1=squeeze(sum(sum(area(:,:,:,1))));
364    tmp2=squeeze(sum(sum(area(:,:,:,2))));
365    tmp3=squeeze(sum(sum(area(:,:,:,3))));
366    subplot(413),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
367    axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ...
368            max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid
369    title('sum(area), original config')
370    legend('.01^oC','.02^oC',0)
371    subplot(414),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
372    axis([0 max(t) -2e-5 1e-5]), grid, xlabel('days')
373    title('sum(area), original config')
374    legend('.01^oC','.02^oC',0)
375    print -djpeg FIG_orig
376    
377  grid,title([int2str((t-1)*10) ' days'])  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
378  tmp=squeeze(sum(sum(area(:,:,:,1))));  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 subplot(212),plot(tmp-mean(tmp))  
 grid,title([int2str((t-1)*10) ' days'])  
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22