depth_control_no_nsa/input/topbot.m

function [tfld,bfld]=topbot(g,fld);
%function [tfld,bfld] = mit_bottomfield(g,fld); or
%function [tfld,bfld] = mit_bottomfield(cmask,fld);

  if isstruct(g)
    cm = change(g.cmask,'==',NaN,0);
  else
    cm = change(g,'==',NaN,0);
  end
  tfld = repmat(NaN,size(fld(:,:,1,:)));
  bfld = repmat(NaN,size(fld(:,:,1,:)));
  for kx=1:size(fld,1);
    for ky=1:size(fld,2);
      kt=min(find(cm(kx,ky,:)>0));
      kb=max(find(cm(kx,ky,:)>0));
      if kt > 0;
        tfld(kx,ky,1,:) = fld(kx,ky,kt,:);
        bfld(kx,ky,1,:) = fld(kx,ky,kb,:);
      end
    end
  end
  tfld = squeeze(tfld);
  bfld = squeeze(bfld);
  
  return
  
1	dgoldberg	1.1	function [tfld,bfld]=topbot(g,fld);
2			%function [tfld,bfld] = mit_bottomfield(g,fld); or
3			%function [tfld,bfld] = mit_bottomfield(cmask,fld);
4
5			if isstruct(g)
6			cm = change(g.cmask,'==',NaN,0);
7			else
8			cm = change(g,'==',NaN,0);
9			end
10			tfld = repmat(NaN,size(fld(:,:,1,:)));
11			bfld = repmat(NaN,size(fld(:,:,1,:)));
12			for kx=1:size(fld,1);
13			for ky=1:size(fld,2);
14			kt=min(find(cm(kx,ky,:)>0));
15			kb=max(find(cm(kx,ky,:)>0));
16			if kt > 0;
17			tfld(kx,ky,1,:) = fld(kx,ky,kt,:);
18			bfld(kx,ky,1,:) = fld(kx,ky,kb,:);
19			end
20			end
21			end
22			tfld = squeeze(tfld);
23			bfld = squeeze(bfld);
24
25			return
26