--- MITgcm_contrib/osse/utils/hmaker.m	2004/06/10 16:21:40	1.2
+++ MITgcm_contrib/osse/utils/hmaker.m	2004/06/15 13:56:05	1.3
@@ -1,26 +1,39 @@
 
+
 nx=120;
 ny=31;
 nz=29;
 nf=2;
+h=zeros(nf*nz*ny*nx,1);
+chance=1/3;
+
 vectorsize=nx*ny*nz*nf;
 dim=1;
-fid = fopen('hfile.txt','w')
+obs=0;
+fid = fopen('hfile.txt','w');
 
 for i=1:nf
    for j=1:nz
       for k=1:ny 
         for l=1:nx
-            if (any(j==[5 10 15 20 25]) && ... % 5 levels
-                any(i==[1 2]) && ... % u, v
-                mod(k,2)==0)
+            if (any(j==[5 10 15 20 25]) & ... % 5 levels
+                 k>8 )
+%                 mod(l,3)==0 )
+%                any(i==[1 2]) & ... % u, v
+%                mod(k+1,2)==0)
+%                mod(l,3)==0)
 %	   if mod(dim,359)==0
+%            if (mod(dim,3)==0)
+                 h(dim)=dim;
                  fprintf(fid,'%i\n',dim);
+                 obs=obs+1;
             end
             dim=dim+1;
          end
       end
    end
 end
-
+fclose(fid);
+obs
+h3=reshape(h,[nx,ny,nz,nf]);