/[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.2 by dimitri, Mon Jul 5 14:15:44 2004 UTC revision 1.3 by dimitri, Mon Jul 5 20:59:51 2004 UTC
# Line 4  Test for lab_sea forward model linearity Line 4  Test for lab_sea forward model linearity
4  cd MITgcm/verification/lab_sea  cd MITgcm/verification/lab_sea
5  mkdir build  mkdir build
6  cd build  cd build
7  cp ../code/*.h ../code/packages.conf .  cp ../code/*.h ../code/*.F ../code/packages.conf .
8  ../../../tools/genmake2  ../../../tools/genmake2
9  make depend  make depend
10  make  make
# Line 138  cd .. Line 138  cd ..
138    
139  matlab  matlab
140    
141  pn='nokpp_nogm_withseaicedyn/';  clear all, clf reset, pn=''; nt=238;
142    heff=zeros(20,16,nt,21);
143  heff=zeros(20,16,21,21);  area=zeros(20,16,nt,21);
 area=zeros(20,16,21,21);  
144  for e=0:20  for e=0:20
145   for t=0:4, disp([e t])   for t=0:nt, disp([e t])
146    if e==0, fn=[pn 'run/HEFF.'];    if e==0, fn=[pn 'run/HEFF.'];
147    else fn=[pn 'run' int2str(e) '/HEFF.']; end    else fn=[pn 'run' int2str(e) '/HEFF.']; end
148    fn=[fn myint2str(t*240,10) '.data'];    fn=[fn myint2str(t*240,10) '.data'];
# Line 154  for e=0:20 Line 153  for e=0:20
153    area(:,:,t+1,e+1)=readbin(fn,[20 16],1);    area(:,:,t+1,e+1)=readbin(fn,[20 16],1);
154   end   end
155  end  end
156    for t=1:nt
 for t=1:21  
157  tmp=squeeze(sum(sum(heff(:,:,t,:))));  tmp=squeeze(sum(sum(heff(:,:,t,:))));
158  clf,subplot(211),plot(tmp-mean(tmp))  clf,subplot(211),plot(tmp-mean(tmp))
159  grid,title([int2str((t-1)*10) ' days'])  grid,title([int2str((t-1)*10) ' days'])
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22