/[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

--- 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'])

 

  ViewVC Help
Powered by ViewVC 1.1.22