--- MITgcm_contrib/osse/utils/hmaker.m	2004/06/15 13:56:05	1.3
+++ MITgcm_contrib/osse/utils/hmaker.m	2004/06/23 13:28:28	1.4
@@ -3,27 +3,24 @@
 nx=120;
 ny=31;
 nz=29;
-nf=2;
+nf=4; % covers u, v, w, and theta
 h=zeros(nf*nz*ny*nx,1);
-chance=1/3;
 
 vectorsize=nx*ny*nz*nf;
-dim=1;
-obs=0;
+dim=1; % current dimension in loop
+obs=0; % observation counter
 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
-                 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)
+            if ((any(j==[5 10 15 20 25]) & ... % 5 levels...
+                 any(i==[1 2]) & ...  % of u and v...
+                 k>8 & ...  % that clear the center...
+                 mod(l,3)==0 ) | ... % every third "spoke", plus...
+                (i==4 & ... % theta...
+                 any(k==[9 30]))) % in the vert boundary layers
                  h(dim)=dim;
                  fprintf(fid,'%i\n',dim);
                  obs=obs+1;