/[MITgcm]/MITgcm_contrib/osse/utils/hmaker.m
ViewVC logotype

Diff of /MITgcm_contrib/osse/utils/hmaker.m

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

revision 1.1 by afe, Tue Jun 8 19:38:54 2004 UTC revision 1.3 by afe, Tue Jun 15 13:56:05 2004 UTC
# Line 1  Line 1 
1    
2    
3  nx=120;  nx=120;
4  ny=31;  ny=31;
5  nz=29;  nz=29;
6  nf=2;  nf=2;
7    h=zeros(nf*nz*ny*nx,1);
8    chance=1/3;
9    
10  vectorsize=nx*ny*nz*nf;  vectorsize=nx*ny*nz*nf;
11  dim=1;  dim=1;
12  fid = fopen('hfile.txt','w')  obs=0;
13    fid = fopen('hfile.txt','w');
14    
15  for i=1:nf  for i=1:nf
16     for j=1:nz     for j=1:nz
17        for k=1:ny % get every other ring        for k=1:ny
18          for l=1:nx          for l=1:nx
19              if (any(j==[5 10 15 20 25]) && ... % 5 levels              if (any(j==[5 10 15 20 25]) & ... % 5 levels
20                  any(i==[1 2]) && ... % u, v                   k>8 )
21                  mod(k,2)==0)  %                 mod(l,3)==0 )
22    %                any(i==[1 2]) & ... % u, v
23    %                mod(k+1,2)==0)
24    %                mod(l,3)==0)
25  %          if mod(dim,359)==0  %          if mod(dim,359)==0
26    %            if (mod(dim,3)==0)
27                     h(dim)=dim;
28                   fprintf(fid,'%i\n',dim);                   fprintf(fid,'%i\n',dim);
29                     obs=obs+1;
30              end              end
31              dim=dim+1;              dim=dim+1;
32           end           end
33        end        end
34     end     end
35  end  end
36    fclose(fid);
37    obs
38    h3=reshape(h,[nx,ny,nz,nf]);
39    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22