--- MITgcm_contrib/lab_sea_test/read.me 2004/07/05 14:15:44 1.2 +++ MITgcm_contrib/lab_sea_test/read.me 2004/07/05 20:59:51 1.3 @@ -4,7 +4,7 @@ cd MITgcm/verification/lab_sea mkdir build cd build -cp ../code/*.h ../code/packages.conf . +cp ../code/*.h ../code/*.F ../code/packages.conf . ../../../tools/genmake2 make depend make @@ -138,12 +138,11 @@ matlab -pn='nokpp_nogm_withseaicedyn/'; - -heff=zeros(20,16,21,21); -area=zeros(20,16,21,21); +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:4, disp([e t]) + for t=0:nt, disp([e t]) if e==0, fn=[pn 'run/HEFF.']; else fn=[pn 'run' int2str(e) '/HEFF.']; end fn=[fn myint2str(t*240,10) '.data']; @@ -154,8 +153,7 @@ area(:,:,t+1,e+1)=readbin(fn,[20 16],1); end end - -for t=1:21 +for t=1:nt tmp=squeeze(sum(sum(heff(:,:,t,:)))); clf,subplot(211),plot(tmp-mean(tmp)) grid,title([int2str((t-1)*10) ' days'])