/[MITgcm]/MITgcm_contrib/afe/osse_MkII/utils/ensavg.m
ViewVC logotype

Diff of /MITgcm_contrib/afe/osse_MkII/utils/ensavg.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by afe, Thu Feb 17 22:42:13 2005 UTC revision 1.2 by afe, Fri Mar 11 19:29:00 2005 UTC
# Line 1  Line 1 
1  nx=120;ny=31;nz=29;n=nx*ny*nz;nens=30;  nx=120;ny=31;nz=29;n=nx*ny*nz;nens=30;
2    %xdata=nx;ydata=ny;zdata=nz;ndata=xdata*ydata*zdata
3    xdata=90;ydata=21;zdata=29;ndata=xdata*ydata*zdata
4  %iter=68;  %iter=68;
5  z=15;  z=15;
6  disp=1;  disp=1;
# Line 12  P=find(vars=='P'); Line 14  P=find(vars=='P');
14  T=find(vars=='T');  T=find(vars=='T');
15  varnum= size(vars,2);  varnum= size(vars,2);
16    
17    if(1)
18  %filename=sprintf('%s','../da/00/assimilate/pickup.in')  %filename=sprintf('%s','../da/00/assimilate/pickup.in')
19  filename=sprintf('%s%03i%s','../da/00/assimilate/pickup.',iter,'.out');  %filename=sprintf('%s%03i%s','../da/00/assimilate/pickup.',iter,'.out');
20  %filename=sprintf('%s','../da/i90/00/assimilate/pickup.0000016820.001.001.data')  filename=sprintf('%s','../da/00/assimilate/pickup.0000016820.001.001.data')
21   fid=fopen(filename,'r','ieee-be');   fid=fopen(filename,'r','ieee-be');
22  foo=fread(fid,13*n+ny*nx,'float64');  foo=fread(fid,13*ndata+ydata*xdata,'float64');
23  i=0;tru(:,:,:,U)=reshape(foo(i*n+1:(i+1)*n),[nx ny nz]);  i=0;tru(:,:,:,U)=reshape(foo(i*ndata+1:(i+1)*ndata),[xdata ydata zdata]);
24  i=3;tru(:,:,:,V)=reshape(foo(i*n+1:(i+1)*n),[nx ny nz]);  i=3;tru(:,:,:,V)=reshape(foo(i*ndata+1:(i+1)*ndata),[xdata ydata zdata]);
25  i=6;tru(:,:,:,W)=reshape(foo(i*n+1:(i+1)*n),[nx ny nz]);  i=6;tru(:,:,:,W)=reshape(foo(i*ndata+1:(i+1)*ndata),[xdata ydata zdata]);
26  i=7;tru(:,:,:,T)=reshape(foo(i*n+1:(i+1)*n),[nx ny nz]);  i=7;tru(:,:,:,T)=reshape(foo(i*ndata+1:(i+1)*ndata),[xdata ydata zdata]);
27  fclose(fid);  fclose(fid);
28    end
29    
30  for k=1:nens  for k=1:nens
31  %k=0;  %k=0;
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22