% Script to edit a file by delete records for which
% a criterion variable falls outside a designated 
% "good" range of values.
%        Fran Hotchkiss 21 August 1998
% Corrected to only update min/max for variables that
% have them.  Fran Hotchkiss 1 Mar 1999

% Use dialog box to choose input file.
disp ('* Select raw data file.')
ding
rawc = netcdf;
% Use dialog box to designate variable that is 
%   basis of truncation, and its range of
%   "good values."
C_range_delete.infile = justname(name(rawc));
C_range_delete.outfile = [justname(name(rawc)) 'e'];
C_range_delete.variable = 'NAV_105';
C_range_delete.good_minimum = 150;
C_range_delete.good_maximum = 250;
uigetparm(C_range_delete);
load C_range_delete
cmin = C_range_delete.good_minimum;
cmax = C_range_delete.good_maximum;

% Find records for which designated 
%   variable is good.
dvar = rawc{C_range_delete.variable}(:);
igood = find(dvar>cmin & dvar<cmax);
ngood = length(igood)%;

% Open output file.
outc = netcdf(C_range_delete.outfile, 'noclobber');

% Global attributes.
%   Copy all global attributes.
rawcatt = att(rawc);
for i = [1:length(rawcatt)]
   copy (rawcatt{i}, outc);
end
% Update global attributes.
history =[' Edited using range_delete.m to delete records unless ' C_range_delete.variable ];  
history =[history ' is in the range ' num2str(cmin) ' to ' num2str(cmax) '.  :' rawc.history(:)];  
outc.history = history;
outc.CREATION_DATE = ncchar(datestr(now,0));

% Dimension definitions and attributes.
rawcdim = dim(rawc);
for i = [1:length(rawcdim)]
   copy (rawcdim{i}, outc);
end
% Variable definitions and attributes.
rawcvar = var(rawc);
for i = [1:length(rawcvar)]
   ivarname = name(rawcvar{i})
   copy (rawcvar{i}, outc, 0, 1);
end
outc{'time'}.type = ncchar('UNEVEN');
outc{'time2'}.type = ncchar('UNEVEN');

% Coordinate variables. 
rawccoord = coord(rawc);
for i = [1:length(rawccoord)]
   ivarname = name(rawccoord{i})
   j = strmatch('time',ivarname,'exact');
   if j
      mstrng = ['Coordinate variable ' ivarname ' found and not filled.']
      disp (mstrng)
   else
      copy (rawccoord{i}, outc, 1, 0);
   end
end

% Record variables.
rawcrecvar = recvar(rawc);
for i = [1:length(rawcrecvar)]
   ivarname = name(rawcrecvar{i})
	data = rawc{ivarname}(:);
  	outc{ivarname}(1:ngood, :) = data(igood);
end

ngood % check agreement in lengths of variables

% Update start time, stop time, variable min/max.
outctime = outc{'time'}(:);
outctime2 = outc{'time2'}(:);
outcstime = ep_datenum([outctime(1) outctime2(1)]);
outcn = length(outctime)%;
outcltime = ep_datenum([outctime(outcn) outctime2(outcn)]);
outc.start_time = datestr(outcstime,0);
outc.stop_time = datestr(outcltime,0);
for i = [1:length(rawcrecvar)]
   ivarname = name(rawcrecvar{i})
	if (~isempty(outc{ivarname}.minimum(:)))
      data = outc{ivarname}(:);
      outc{ivarname}.minimum = ncfloat(min(data));
      outc{ivarname}.maximum = ncfloat(max(data));
   end
end
ncclose