/[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.3 by dimitri, Mon Jul 5 20:59:51 2004 UTC revision 1.5 by dimitri, Mon Jul 12 01:00:20 2004 UTC
# Line 1  Line 1 
1    Tests for lab_sea forward model linearity.
2    
3  Test for lab_sea forward model linearity.  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
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 ..; mkdir run3; cd run3; cp ../input/* .  cd ..
31  cp ../../../../MITgcm_contrib/lab_sea_test/* .  matlab
 mv data.exf tmp  
 sed s/tair/tair03/ tmp > data.exf  
 ../build/mitgcmuv > output.txt &  
32    
33  cd ..; mkdir run4; cd run4; cp ../input/* .  clear all, clf reset, pn=''; ne=2; nt=53;
34  cp ../../../../MITgcm_contrib/lab_sea_test/* .  heff=zeros(20,16,nt,ne+1);
35  mv data.exf tmp  area=zeros(20,16,nt,ne+1);
36  sed s/tair/tair04/ tmp > data.exf  for e=0:ne
37  ../build/mitgcmuv > output.txt &   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  cd ..; mkdir run5; cd run5; cp ../input/* .  clf reset, orient tall, wysiwyg
50  cp ../../../../MITgcm_contrib/lab_sea_test/* .  tmp1=squeeze(sum(sum(heff(:,:,:,1))));
51  mv data.exf tmp  tmp2=squeeze(sum(sum(heff(:,:,:,2))));
52  sed s/tair/tair05/ tmp > data.exf  tmp3=squeeze(sum(sum(heff(:,:,:,3))));
53  ../build/mitgcmuv > output.txt &  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  cd ..; mkdir run6; cd run6; cp ../input/* .  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
78  cp ../../../../MITgcm_contrib/lab_sea_test/* .  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 mv data.exf tmp  
 sed s/tair/tair06/ tmp > data.exf  
 ../build/mitgcmuv > output.txt &  
79    
80  cd ..; mkdir run7; cd run7; cp ../input/* .  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/* .  cp ../../../../MITgcm_contrib/lab_sea_test/* .
90  mv data.exf tmp  mv data.pkg tmp
91  sed s/tair/tair07/ tmp > data.exf  sed s/"GMRedi = .TRUE"/"GMRedi = .FALSE"/ tmp > tmp1
92  ../build/mitgcmuv > output.txt &  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 run8; cd run8; cp ../input/* .  cd ..; mkdir run1; cd run1; cp ../run/* .
 cp ../../../../MITgcm_contrib/lab_sea_test/* .  
98  mv data.exf tmp  mv data.exf tmp
99  sed s/tair/tair08/ tmp > data.exf  sed s/tair/tair01/ tmp > data.exf
100  ../build/mitgcmuv > output.txt &  ../build/mitgcmuv >& output.txt &
101    
102  cd ..; mkdir run9; cd run9; cp ../input/* .  cd ..; mkdir run2; cd run2; cp ../run/* .
 cp ../../../../MITgcm_contrib/lab_sea_test/* .  
103  mv data.exf tmp  mv data.exf tmp
104  sed s/tair/tair09/ tmp > data.exf  sed s/tair/tair02/ tmp > data.exf
105  ../build/mitgcmuv > output.txt &  ../build/mitgcmuv >& output.txt &
106    
107  cd ..; mkdir run10; cd run10; cp ../input/* .  cd ..
108  cp ../../../../MITgcm_contrib/lab_sea_test/* .  matlab
 mv data.exf tmp  
 sed s/tair/tair10/ tmp > data.exf  
 ../build/mitgcmuv > output.txt &  
109    
110  cd ..; mkdir run11; cd run11; cp ../input/* .  clear all, pn=''; ne=2; nt=53;
111  cp ../../../../MITgcm_contrib/lab_sea_test/* .  heff=zeros(20,16,nt,ne+1);
112  mv data.exf tmp  area=zeros(20,16,nt,ne+1);
113  sed s/tair/tair11/ tmp > data.exf  for e=0:ne
114  ../build/mitgcmuv > output.txt &   for t=1:nt, disp([e t])
115      if e==0, fn=[pn 'run/HEFFtave.'];
116      else fn=[pn 'run' int2str(e) '/HEFFtave.']; end
117      fn=[fn myint2str(t*240,10) '.data'];
118      heff(:,:,t+1,e+1)=readbin(fn,[20 16],1);
119      if e==0, fn=[pn 'run/AREAtave.'];
120      else fn=[pn 'run' int2str(e) '/AREAtave.']; end
121      fn=[fn myint2str(t*240,10) '.data'];
122      area(:,:,t+1,e+1)=readbin(fn,[20 16],1);
123     end
124    end
125    
126  cd ..; mkdir run12; cd run12; cp ../input/* .  clf reset, orient tall, wysiwyg
127  cp ../../../../MITgcm_contrib/lab_sea_test/* .  tmp1=squeeze(sum(sum(heff(:,:,:,1))));
128  mv data.exf tmp  tmp2=squeeze(sum(sum(heff(:,:,:,2))));
129  sed s/tair/tair12/ tmp > data.exf  tmp3=squeeze(sum(sum(heff(:,:,:,3))));
130  ../build/mitgcmuv > output.txt &  t=1:length(tmp1); t=(t-1)*10
131    subplot(411),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
132    axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ...
133            max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid
134    title('sum(heff), original config')
135    legend('.01^oC','.02^oC',0)
136    subplot(412),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
137    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  cd ..; mkdir run13; cd run13; cp ../input/* .  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
155  cp ../../../../MITgcm_contrib/lab_sea_test/* .  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 mv data.exf tmp  
 sed s/tair/tair13/ tmp > data.exf  
 ../build/mitgcmuv > output.txt &  
156    
157  cd ..; mkdir run14; cd run14; cp ../input/* .  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/* .  cp ../../../../MITgcm_contrib/lab_sea_test/* .
167  mv data.exf tmp  ../build/mitgcmuv >& output.txt &
 sed s/tair/tair14/ tmp > data.exf  
 ../build/mitgcmuv > output.txt &  
168    
169  cd ..; mkdir run15; cd run15; cp ../input/* .  cd ..; mkdir run1; cd run1; cp ../input/* .
170  cp ../../../../MITgcm_contrib/lab_sea_test/* .  cp ../../../../MITgcm_contrib/lab_sea_test/* .
171  mv data.exf tmp  mv data.exf tmp
172  sed s/tair/tair15/ tmp > data.exf  sed s/tair/tair01/ tmp > data.exf
173  ../build/mitgcmuv > output.txt &  ../build/mitgcmuv >& output.txt &
174    
175  cd ..; mkdir run16; cd run16; cp ../input/* .  cd ..; mkdir run2; cd run2; cp ../input/* .
176  cp ../../../../MITgcm_contrib/lab_sea_test/* .  cp ../../../../MITgcm_contrib/lab_sea_test/* .
177  mv data.exf tmp  mv data.exf tmp
178  sed s/tair/tair16/ tmp > data.exf  sed s/tair/tair02/ tmp > data.exf
179  ../build/mitgcmuv > output.txt &  ../build/mitgcmuv >& output.txt &
180    
181  cd ..; mkdir run17; cd run17; cp ../input/* .  cd ..
182  cp ../../../../MITgcm_contrib/lab_sea_test/* .  matlab
 mv data.exf tmp  
 sed s/tair/tair17/ tmp > data.exf  
 ../build/mitgcmuv > output.txt &  
183    
184  cd ..; mkdir run18; cd run18; cp ../input/* .  clear all, clf reset, pn=''; ne=2; nt=53;
185  cp ../../../../MITgcm_contrib/lab_sea_test/* .  heff=zeros(20,16,nt,ne+1);
186  mv data.exf tmp  area=zeros(20,16,nt,ne+1);
187  sed s/tair/tair18/ tmp > data.exf  for e=0:ne
188  ../build/mitgcmuv > output.txt &   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))));
202    tmp2=squeeze(sum(sum(heff(:,:,:,2))));
203    tmp3=squeeze(sum(sum(heff(:,:,:,3))));
204    t=1:length(tmp1); t=(t-1)*10
205    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  cd ..; mkdir run19; cd run19; cp ../input/* .  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
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/* .  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  mv data.exf tmp
249  sed s/tair/tair19/ tmp > data.exf  sed s/tair/tair01/ tmp > data.exf
250  ../build/mitgcmuv > output.txt &  ../build/mitgcmuv >& output.txt &
251    
252  cd ..; mkdir run20; cd run20; cp ../input/* .  cd ..; mkdir run2; cd run2; cp ../run/* .
 cp ../../../../MITgcm_contrib/lab_sea_test/* .  
253  mv data.exf tmp  mv data.exf tmp
254  sed s/tair/tair20/ tmp > data.exf  sed s/tair/tair02/ tmp > data.exf
255  ../build/mitgcmuv > output.txt &  ../build/mitgcmuv >& output.txt &
256    
257  cd ..  cd ..
   
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  
   
258  matlab  matlab
259    
260  clear all, clf reset, pn=''; nt=238;  clear all, clf reset, pn=''; ne=2; nt=53;
261  heff=zeros(20,16,nt,21);  heff=zeros(20,16,nt,ne+1);
262  area=zeros(20,16,nt,21);  area=zeros(20,16,nt,ne+1);
263  for e=0:20  for e=0:ne
264   for t=0:nt, disp([e t])   for t=1:nt, disp([e t])
265    if e==0, fn=[pn 'run/HEFF.'];    if e==0, fn=[pn 'run/HEFFtave.'];
266    else fn=[pn 'run' int2str(e) '/HEFF.']; end    else fn=[pn 'run' int2str(e) '/HEFFtave.']; end
267    fn=[fn myint2str(t*240,10) '.data'];    fn=[fn myint2str(t*240,10) '.data'];
268    heff(:,:,t+1,e+1)=readbin(fn,[20 16],1);    heff(:,:,t+1,e+1)=readbin(fn,[20 16],1);
269    if e==0, fn=[pn 'run/AREA.'];    if e==0, fn=[pn 'run/AREAtave.'];
270    else fn=[pn 'run' int2str(e) '/AREA.']; end    else fn=[pn 'run' int2str(e) '/AREAtave.']; end
271    fn=[fn myint2str(t*240,10) '.data'];    fn=[fn myint2str(t*240,10) '.data'];
272    area(:,:,t+1,e+1)=readbin(fn,[20 16],1);    area(:,:,t+1,e+1)=readbin(fn,[20 16],1);
273   end   end
274  end  end
 for t=1:nt  
 tmp=squeeze(sum(sum(heff(:,:,t,:))));  
 clf,subplot(211),plot(tmp-mean(tmp))  
 grid,title([int2str((t-1)*10) ' days'])  
 tmp=squeeze(sum(sum(area(:,:,t,:))));  
 subplot(212),plot(tmp-mean(tmp))  
 grid,title([int2str((t-1)*10) ' days'])  
 pause, end  
   
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  
   
 repeat with daily dumps for experiments 10-15, days 6-10  
   
 cd run10; mv data tmp  
 sed s/startTime=0.0,/startTime=518400.0,/ tmp > tmp2  
 sed s/endTime=1728000.,/endTime=950400.0,/ tmp2 > tmp3  
 sed s/"dumpFreq = 86400.,"/dumpFreq=3600.,/ tmp3 > data  
 ../build/mitgcmuv >! output.txt &  
   
 cd ../run11; mv data tmp  
 sed s/startTime=0.0,/startTime=518400.0,/ tmp > tmp2  
 sed s/endTime=1728000.,/endTime=950400.0,/ tmp2 > tmp3  
 sed s/"dumpFreq = 86400.,"/dumpFreq=3600.,/ tmp3 > data  
 ../build/mitgcmuv >! output.txt &  
   
 cd ../run12; mv data tmp  
 sed s/startTime=0.0,/startTime=518400.0,/ tmp > tmp2  
 sed s/endTime=1728000.,/endTime=950400.0,/ tmp2 > tmp3  
 sed s/"dumpFreq = 86400.,"/dumpFreq=3600.,/ tmp3 > data  
 ../build/mitgcmuv >! output.txt &  
   
 cd ../run13; mv data tmp  
 sed s/startTime=0.0,/startTime=518400.0,/ tmp > tmp2  
 sed s/endTime=1728000.,/endTime=950400.0,/ tmp2 > tmp3  
 sed s/"dumpFreq = 86400.,"/dumpFreq=3600.,/ tmp3 > data  
 ../build/mitgcmuv >! output.txt &  
   
 cd ../run14; mv data tmp  
 sed s/startTime=0.0,/startTime=518400.0,/ tmp > tmp2  
 sed s/endTime=1728000.,/endTime=950400.0,/ tmp2 > tmp3  
 sed s/"dumpFreq = 86400.,"/dumpFreq=3600.,/ tmp3 > data  
 ../build/mitgcmuv >! output.txt &  
   
 cd ../run15; mv data tmp  
 sed s/startTime=0.0,/startTime=518400.0,/ tmp > tmp2  
 sed s/endTime=1728000.,/endTime=950400.0,/ tmp2 > tmp3  
 sed s/"dumpFreq = 86400.,"/dumpFreq=3600.,/ tmp3 > data  
 ../build/mitgcmuv >! output.txt &  
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  matlab
332    
333  heff=zeros(20,16,121,6);  clear all, clf reset, pn=''; ne=2; nt=53;
334  area=zeros(20,16,121,6);  heff=zeros(20,16,nt,ne+1);
335  for e=10:15, disp(e)  area=zeros(20,16,nt,ne+1);
336   for t=144:264  for e=0:ne
337    fn=['run' int2str(e) '/HEFF.'];   for t=1:nt, disp([e t])
338    fn=[fn myint2str(t,10) '.data'];    if e==0, fn=[pn 'run/HEFFtave.'];
339    heff(:,:,t-143,e-9)=readbin(fn,[20 16],1);    else fn=[pn 'run' int2str(e) '/HEFFtave.']; end
340    fn=['run' int2str(e) '/AREA.'];    fn=[fn myint2str(t*240,10) '.data'];
341    fn=[fn myint2str(t,10) '.data'];    heff(:,:,t+1,e+1)=readbin(fn,[20 16],1);
342    area(:,:,t-143,e-9)=readbin(fn,[20 16],1);    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   end
347  end  end
348    
349  for t=144:264  clf reset, orient tall, wysiwyg
350  tmp=squeeze(sum(sum(heff(:,:,t-143,:))));  tmp1=squeeze(sum(sum(heff(:,:,:,1))));
351  clf,subplot(211), plot(tmp-mean(tmp)), grid, title(t/24)  tmp2=squeeze(sum(sum(heff(:,:,:,2))));
352  tmp=squeeze(sum(sum(area(:,:,t-143,:))));  tmp3=squeeze(sum(sum(heff(:,:,:,3))));
353  subplot(212), plot(tmp-mean(tmp)), grid, title(t/24)  t=1:length(tmp1); t=(t-1)*10
354  pause, end  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    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
378    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22