--- MITgcm_contrib/lab_sea_test/read.me 2004/07/05 20:59:51 1.3 +++ MITgcm_contrib/lab_sea_test/read.me 2004/07/12 01:00:20 1.5 @@ -1,227 +1,378 @@ +Tests for lab_sea forward model linearity. -Test for lab_sea forward model linearity. +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% cd MITgcm/verification/lab_sea mkdir build cd build -cp ../code/*.h ../code/*.F ../code/packages.conf . +cp ../code/*.h ../code/packages.conf . +cp ../../../../MITgcm_contrib/lab_sea_test/*.F . ../../../tools/genmake2 make depend 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 & +../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 .. +matlab -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 & +clear all, clf reset, pn=''; ne=2; nt=53; +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/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 -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 & +clf reset, orient tall, wysiwyg +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 +subplot(411),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ... + max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid +title('sum(heff), original config') +legend('.01^oC','.02^oC',0) +subplot(412),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) -2e-5 1e-5]), grid +title('sum(heff), original config') +legend('.01^oC','.02^oC',0) +tmp1=squeeze(sum(sum(area(:,:,:,1)))); +tmp2=squeeze(sum(sum(area(:,:,:,2)))); +tmp3=squeeze(sum(sum(area(:,:,:,3)))); +subplot(413),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ... + max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid +title('sum(area), original config') +legend('.01^oC','.02^oC',0) +subplot(414),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) -2e-5 1e-5]), grid, xlabel('days') +title('sum(area), original config') +legend('.01^oC','.02^oC',0) +print -djpeg FIG_ynegminmaxrelax -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/* . +cd MITgcm/verification/lab_sea +mkdir build +cd build +cp ../code/*.h ../code/packages.conf . +cp ../../../../MITgcm_contrib/lab_sea_test/growth.F . +../../../tools/genmake2 +make depend +make +cd ..; mkdir run; cd run; cp ../input/* . cp ../../../../MITgcm_contrib/lab_sea_test/* . -mv data.exf tmp -sed s/tair/tair07/ tmp > data.exf -../build/mitgcmuv > output.txt & +mv data.pkg tmp +sed s/"GMRedi = .TRUE"/"GMRedi = .FALSE"/ tmp > tmp1 +sed s/"KPP = .TRUE"/"KPP = .FALSE"/ tmp1 > data.pkg +mv data.seaice tmp +sed s/"DYNAMICS = .TRUE"/"DYNAMICS = .FALSE"/ tmp > data.seaice +../build/mitgcmuv >& output.txt & -cd ..; mkdir run8; cd run8; cp ../input/* . -cp ../../../../MITgcm_contrib/lab_sea_test/* . +cd ..; mkdir run1; cd run1; cp ../run/* . mv data.exf tmp -sed s/tair/tair08/ tmp > data.exf -../build/mitgcmuv > output.txt & +sed s/tair/tair01/ tmp > data.exf +../build/mitgcmuv >& output.txt & -cd ..; mkdir run9; cd run9; cp ../input/* . -cp ../../../../MITgcm_contrib/lab_sea_test/* . +cd ..; mkdir run2; cd run2; cp ../run/* . mv data.exf tmp -sed s/tair/tair09/ tmp > data.exf -../build/mitgcmuv > output.txt & +sed s/tair/tair02/ 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 .. +matlab -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 & +clear all, pn=''; ne=2; nt=53; +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/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 -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 & +clf reset, orient tall, wysiwyg +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 +subplot(411),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ... + max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid +title('sum(heff), original config') +legend('.01^oC','.02^oC',0) +subplot(412),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) -2e-5 1e-5]), grid +title('sum(heff), original config') +legend('.01^oC','.02^oC',0) +tmp1=squeeze(sum(sum(area(:,:,:,1)))); +tmp2=squeeze(sum(sum(area(:,:,:,2)))); +tmp3=squeeze(sum(sum(area(:,:,:,3)))); +subplot(413),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ... + max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid +title('sum(area), original config') +legend('.01^oC','.02^oC',0) +subplot(414),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) -2e-5 1e-5]), grid, xlabel('days') +title('sum(area), original config') +legend('.01^oC','.02^oC',0) +print -djpeg FIG_ynegrelax_no_kpp_gm_dyn -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/* . +cd MITgcm/verification/lab_sea +mkdir build +cd build +cp ../code/*.h ../code/packages.conf . +cp ../../../../MITgcm_contrib/lab_sea_test/growth.F . +../../../tools/genmake2 +make depend +make +cd ..; mkdir run; cd run; cp ../input/* . cp ../../../../MITgcm_contrib/lab_sea_test/* . -mv data.exf tmp -sed s/tair/tair14/ tmp > data.exf -../build/mitgcmuv > output.txt & +../build/mitgcmuv >& output.txt & -cd ..; mkdir run15; cd run15; cp ../input/* . +cd ..; mkdir run1; cd run1; cp ../input/* . cp ../../../../MITgcm_contrib/lab_sea_test/* . mv data.exf tmp -sed s/tair/tair15/ tmp > data.exf -../build/mitgcmuv > output.txt & +sed s/tair/tair01/ tmp > data.exf +../build/mitgcmuv >& output.txt & -cd ..; mkdir run16; cd run16; cp ../input/* . +cd ..; mkdir run2; cd run2; cp ../input/* . cp ../../../../MITgcm_contrib/lab_sea_test/* . mv data.exf tmp -sed s/tair/tair16/ tmp > data.exf -../build/mitgcmuv > output.txt & +sed s/tair/tair02/ 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 .. +matlab -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 & +clear all, clf reset, pn=''; ne=2; nt=53; +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/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 + +clf reset, orient tall, wysiwyg +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 +subplot(411),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ... + max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid +title('sum(heff), original config') +legend('.01^oC','.02^oC',0) +subplot(412),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) -2e-5 1e-5]), grid +title('sum(heff), original config') +legend('.01^oC','.02^oC',0) +tmp1=squeeze(sum(sum(area(:,:,:,1)))); +tmp2=squeeze(sum(sum(area(:,:,:,2)))); +tmp3=squeeze(sum(sum(area(:,:,:,3)))); +subplot(413),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ... + max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid +title('sum(area), original config') +legend('.01^oC','.02^oC',0) +subplot(414),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) -2e-5 1e-5]), grid, xlabel('days') +title('sum(area), original config') +legend('.01^oC','.02^oC',0) +print -djpeg FIG_ynegrelax -cd ..; mkdir run19; cd run19; cp ../input/* . +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + +cd MITgcm/verification/lab_sea +mkdir build +cd build +cp ../code/*.h ../code/packages.conf . +../../../tools/genmake2 +make depend +make +cd ..; mkdir run; cd run; cp ../input/* . cp ../../../../MITgcm_contrib/lab_sea_test/* . +mv data.pkg tmp +sed s/"GMRedi = .TRUE"/"GMRedi = .FALSE"/ tmp > tmp1 +sed s/"KPP = .TRUE"/"KPP = .FALSE"/ tmp1 > data.pkg +mv data.seaice tmp +sed s/"DYNAMICS = .TRUE"/"DYNAMICS = .FALSE"/ tmp > data.seaice +../build/mitgcmuv >& output.txt & + +cd ..; mkdir run1; cd run1; cp ../run/* . mv data.exf tmp -sed s/tair/tair19/ tmp > data.exf -../build/mitgcmuv > output.txt & +sed s/tair/tair01/ tmp > data.exf +../build/mitgcmuv >& output.txt & -cd ..; mkdir run20; cd run20; cp ../input/* . -cp ../../../../MITgcm_contrib/lab_sea_test/* . +cd ..; mkdir run2; cd run2; cp ../run/* . mv data.exf tmp -sed s/tair/tair20/ tmp > data.exf -../build/mitgcmuv > output.txt & +sed s/tair/tair02/ 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=2; nt=53; +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,:)))); -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 & -%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +clf reset, orient tall, wysiwyg +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 +subplot(411),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ... + max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid +title('sum(heff), original config') +legend('.01^oC','.02^oC',0) +subplot(412),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) -2e-5 1e-5]), grid +title('sum(heff), original config') +legend('.01^oC','.02^oC',0) +tmp1=squeeze(sum(sum(area(:,:,:,1)))); +tmp2=squeeze(sum(sum(area(:,:,:,2)))); +tmp3=squeeze(sum(sum(area(:,:,:,3)))); +subplot(413),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ... + max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid +title('sum(area), original config') +legend('.01^oC','.02^oC',0) +subplot(414),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) -2e-5 1e-5]), grid, xlabel('days') +title('sum(area), original config') +legend('.01^oC','.02^oC',0) +print -djpeg FIG_no_kpp_gm_dyn + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + +cd MITgcm/verification/lab_sea +mkdir build +cd build +cp ../code/*.h ../code/packages.conf . +../../../tools/genmake2 +make depend +make +cd ..; mkdir run; cd run; cp ../input/* . +cp ../../../../MITgcm_contrib/lab_sea_test/* . +../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 & + +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 .. 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); +clear all, clf reset, pn=''; ne=2; nt=53; +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/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=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 +clf reset, orient tall, wysiwyg +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 +subplot(411),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ... + max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid +title('sum(heff), original config') +legend('.01^oC','.02^oC',0) +subplot(412),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) -2e-5 1e-5]), grid +title('sum(heff), original config') +legend('.01^oC','.02^oC',0) +tmp1=squeeze(sum(sum(area(:,:,:,1)))); +tmp2=squeeze(sum(sum(area(:,:,:,2)))); +tmp3=squeeze(sum(sum(area(:,:,:,3)))); +subplot(413),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) min(min(tmp3-tmp1),min(tmp2-tmp1)) ... + max(max(tmp3-tmp1),max(tmp2-tmp1))]), grid +title('sum(area), original config') +legend('.01^oC','.02^oC',0) +subplot(414),plot(t,tmp2-tmp1,t,tmp3-tmp1,'linewidth',2) +axis([0 max(t) -2e-5 1e-5]), grid, xlabel('days') +title('sum(area), original config') +legend('.01^oC','.02^oC',0) +print -djpeg FIG_orig + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%