--- MITgcm_contrib/lab_sea_test/read.me	2004/07/05 20:59:51	1.3
+++ MITgcm_contrib/lab_sea_test/read.me	2004/07/11 06:18:03	1.4
@@ -1,4 +1,3 @@
-
 Test for lab_sea forward model linearity.
 
 cd MITgcm/verification/lab_sea
@@ -10,218 +9,66 @@
 make
 cd ..; mkdir run; cd run; cp ../input/* .
 cp ../../../../MITgcm_contrib/lab_sea_test/* .
-../build/mitgcmuv > output.txt &
+../build/mitgcmuv >! output.txt &
 
 cd ..; mkdir run1; cd run1; cp ../input/* .
 cp ../../../../MITgcm_contrib/lab_sea_test/* .
 mv data.exf tmp
 sed s/tair/tair01/ tmp > data.exf
-../build/mitgcmuv > output.txt &
+../build/mitgcmuv >! output.txt &
 
 cd ..; mkdir run2; cd run2; cp ../input/* .
 cp ../../../../MITgcm_contrib/lab_sea_test/* .
 mv data.exf tmp
 sed s/tair/tair02/ tmp > data.exf
-../build/mitgcmuv > output.txt &
-
-cd ..; mkdir run3; cd run3; cp ../input/* .
-cp ../../../../MITgcm_contrib/lab_sea_test/* .
-mv data.exf tmp
-sed s/tair/tair03/ tmp > data.exf
-../build/mitgcmuv > output.txt &
-
-cd ..; mkdir run4; cd run4; cp ../input/* .
-cp ../../../../MITgcm_contrib/lab_sea_test/* .
-mv data.exf tmp
-sed s/tair/tair04/ tmp > data.exf
-../build/mitgcmuv > output.txt &
-
-cd ..; mkdir run5; cd run5; cp ../input/* .
-cp ../../../../MITgcm_contrib/lab_sea_test/* .
-mv data.exf tmp
-sed s/tair/tair05/ tmp > data.exf
-../build/mitgcmuv > output.txt &
-
-cd ..; mkdir run6; cd run6; cp ../input/* .
-cp ../../../../MITgcm_contrib/lab_sea_test/* .
-mv data.exf tmp
-sed s/tair/tair06/ tmp > data.exf
-../build/mitgcmuv > output.txt &
-
-cd ..; mkdir run7; cd run7; cp ../input/* .
-cp ../../../../MITgcm_contrib/lab_sea_test/* .
-mv data.exf tmp
-sed s/tair/tair07/ tmp > data.exf
-../build/mitgcmuv > output.txt &
-
-cd ..; mkdir run8; cd run8; cp ../input/* .
-cp ../../../../MITgcm_contrib/lab_sea_test/* .
-mv data.exf tmp
-sed s/tair/tair08/ tmp > data.exf
-../build/mitgcmuv > output.txt &
-
-cd ..; mkdir run9; cd run9; cp ../input/* .
-cp ../../../../MITgcm_contrib/lab_sea_test/* .
-mv data.exf tmp
-sed s/tair/tair09/ tmp > data.exf
-../build/mitgcmuv > output.txt &
-
-cd ..; mkdir run10; cd run10; cp ../input/* .
-cp ../../../../MITgcm_contrib/lab_sea_test/* .
-mv data.exf tmp
-sed s/tair/tair10/ tmp > data.exf
-../build/mitgcmuv > output.txt &
-
-cd ..; mkdir run11; cd run11; cp ../input/* .
-cp ../../../../MITgcm_contrib/lab_sea_test/* .
-mv data.exf tmp
-sed s/tair/tair11/ tmp > data.exf
-../build/mitgcmuv > output.txt &
-
-cd ..; mkdir run12; cd run12; cp ../input/* .
-cp ../../../../MITgcm_contrib/lab_sea_test/* .
-mv data.exf tmp
-sed s/tair/tair12/ tmp > data.exf
-../build/mitgcmuv > output.txt &
-
-cd ..; mkdir run13; cd run13; cp ../input/* .
-cp ../../../../MITgcm_contrib/lab_sea_test/* .
-mv data.exf tmp
-sed s/tair/tair13/ tmp > data.exf
-../build/mitgcmuv > output.txt &
-
-cd ..; mkdir run14; cd run14; cp ../input/* .
-cp ../../../../MITgcm_contrib/lab_sea_test/* .
-mv data.exf tmp
-sed s/tair/tair14/ tmp > data.exf
-../build/mitgcmuv > output.txt &
-
-cd ..; mkdir run15; cd run15; cp ../input/* .
-cp ../../../../MITgcm_contrib/lab_sea_test/* .
-mv data.exf tmp
-sed s/tair/tair15/ tmp > data.exf
-../build/mitgcmuv > output.txt &
-
-cd ..; mkdir run16; cd run16; cp ../input/* .
-cp ../../../../MITgcm_contrib/lab_sea_test/* .
-mv data.exf tmp
-sed s/tair/tair16/ tmp > data.exf
-../build/mitgcmuv > output.txt &
-
-cd ..; mkdir run17; cd run17; cp ../input/* .
-cp ../../../../MITgcm_contrib/lab_sea_test/* .
-mv data.exf tmp
-sed s/tair/tair17/ tmp > data.exf
-../build/mitgcmuv > output.txt &
-
-cd ..; mkdir run18; cd run18; cp ../input/* .
-cp ../../../../MITgcm_contrib/lab_sea_test/* .
-mv data.exf tmp
-sed s/tair/tair18/ tmp > data.exf
-../build/mitgcmuv > output.txt &
-
-cd ..; mkdir run19; cd run19; cp ../input/* .
-cp ../../../../MITgcm_contrib/lab_sea_test/* .
-mv data.exf tmp
-sed s/tair/tair19/ tmp > data.exf
-../build/mitgcmuv > output.txt &
+../build/mitgcmuv >! output.txt &
 
-cd ..; mkdir run20; cd run20; cp ../input/* .
-cp ../../../../MITgcm_contrib/lab_sea_test/* .
-mv data.exf tmp
-sed s/tair/tair20/ tmp > data.exf
-../build/mitgcmuv > output.txt &
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
 cd ..
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
 matlab
 
-clear all, clf reset, pn=''; nt=238;
-heff=zeros(20,16,nt,21);
-area=zeros(20,16,nt,21);
-for e=0:20
- for t=0:nt, disp([e t])
-  if e==0, fn=[pn 'run/HEFF.'];
-  else fn=[pn 'run' int2str(e) '/HEFF.']; end
+clear all, clf reset, pn=''; ne=5; nt=120;
+heff=zeros(20,16,nt,ne+1);
+area=zeros(20,16,nt,ne+1);
+for e=0:ne
+ for t=1:nt, disp([e t])
+  if e==0, fn=[pn 'run/HEFFtave.'];
+  else fn=[pn 'run' int2str(e) '/HEFFtave.']; end
   fn=[fn myint2str(t*240,10) '.data'];
   heff(:,:,t+1,e+1)=readbin(fn,[20 16],1);
-  if e==0, fn=[pn 'run/AREA.'];
-  else fn=[pn 'run' int2str(e) '/AREA.']; end
+  if e==0, fn=[pn 'run/AREAtave.'];
+  else fn=[pn 'run' int2str(e) '/AREAtave.']; end
   fn=[fn myint2str(t*240,10) '.data'];
   area(:,:,t+1,e+1)=readbin(fn,[20 16],1);
  end
 end
-for t=1:nt
-tmp=squeeze(sum(sum(heff(:,:,t,:))));
+
+tmp=squeeze(sum(sum(sum(heff))));
 clf,subplot(211),plot(tmp-mean(tmp))
 grid,title([int2str((t-1)*10) ' days'])
-tmp=squeeze(sum(sum(area(:,:,t,:))));
+tmp=squeeze(sum(sum(sum(area))));
+subplot(212),plot(tmp-mean(tmp))
+grid,title([int2str((t-1)*10) ' days'])
+tmp=squeeze(sum(sum(heff(:,:,:,1))));
+clf,subplot(211),plot(tmp-mean(tmp))
+grid,title([int2str((t-1)*10) ' days'])
+tmp=squeeze(sum(sum(area(:,:,:,1))));
 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 &
+tmp1=squeeze(sum(sum(heff(:,:,:,1))));
+tmp2=squeeze(sum(sum(heff(:,:,:,2))));
+tmp3=squeeze(sum(sum(heff(:,:,:,3))));
+t=1:length(tmp1); t=(t-1)*10
+clf,subplot(211),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2)
 
-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 &
 
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
-matlab
 
-heff=zeros(20,16,121,6);
-area=zeros(20,16,121,6);
-for e=10:15, disp(e)
- for t=144:264
-  fn=['run' int2str(e) '/HEFF.'];
-  fn=[fn myint2str(t,10) '.data'];
-  heff(:,:,t-143,e-9)=readbin(fn,[20 16],1);
-  fn=['run' int2str(e) '/AREA.'];
-  fn=[fn myint2str(t,10) '.data'];
-  area(:,:,t-143,e-9)=readbin(fn,[20 16],1);
- end
-end
 
-for t=144:264
-tmp=squeeze(sum(sum(heff(:,:,t-143,:))));
-clf,subplot(211), plot(tmp-mean(tmp)), grid, title(t/24)
-tmp=squeeze(sum(sum(area(:,:,t-143,:))));
-subplot(212), plot(tmp-mean(tmp)), grid, title(t/24)
-pause, end
+grid,title([int2str((t-1)*10) ' days'])
+tmp=squeeze(sum(sum(area(:,:,:,1))));
+subplot(212),plot(tmp-mean(tmp))
+grid,title([int2str((t-1)*10) ' days'])