function [] = mkSADCP(CODAS_CONTOUR_DIR,SADCP_OUTPUT_FILE)
% mkSADCP(CODAS_CONTOUR_DIR,SADCP_OUTPUT_FILE)
%
% CODAS_CONTOUR_DIR path with CODAS output files
% [contour_xy.mat] and [contour_uv.mat]
% SADCP_OUTPUT_FILE output file (including path)
%======================================================================
% M K S A D C P . M
% doc: Mon Mar 15 06:00:43 2004
% dlm: Fri Jan 6 12:28:45 2017
% (c) 2004 A.M. Thurnherr
% uE-Info: 31 61 NIL 0 0 72 0 2 8 NIL ofnI
%======================================================================
%
% pre-processing script to get UH/CODAS-processed SADCP data (usually
% available in near-real time) into a format that can be loaded with
% [loadsadcp.m].
%
% NOTES:
% - this script needs only be called after re-processing the SADCP
% data with the UH/CODAS software
% - during ship-board processing, this script will be called near the
% beginning of [set_cast_params.m] to update the SADCP matrix
% - during post-cruise processing this script will only have to be
% called once, when the finalized SADCP data have become available
%
% HISTORY:
% Mar 15, 2004: - created
% Jan 5, 2007: - documented and added to IX_4 distribution
% Jan 6, 2016: - added time of first SADCP record to output
%======================================================================
eval(['load ' CODAS_CONTOUR_DIR '/contour_xy']);
eval(['load ' CODAS_CONTOUR_DIR '/contour_uv']);
u_sadcp = uv(:,1:2:end);
v_sadcp = uv(:,2:2:end);
lat_sadcp = xyt(2,:);
lon_sadcp = xyt(1,:)-360;
dday = xyt(3,:);
z_sadcp = zc;
badi = find(isnan(dday) | isnan(lat_sadcp) | isnan(lon_sadcp));
dday(badi) = [];
lat_sadcp(badi) = [];
lon_sadcp(badi) = [];
u_sadcp(:, badi) = [];
v_sadcp(:, badi) = [];
tim_sadcp = julian([year_base 1 1 0 0 0])+dday;
lat_sadcp = lat_sadcp(:);
lon_sadcp = lon_sadcp(:);
ii=find(lon_sadcp<-180);
lon_sadcp(ii)=lon_sadcp(ii)+360;
z_sadcp = z_sadcp(:);
tim_sadcp = tim_sadcp(:);
disp(sprintf('first SADCP processing date %d-%d-%d %d:%d:%d',...
gregoria(tim_sadcp(1))));
disp(sprintf(' last SADCP processing date %d-%d-%d %d:%d:%d',...
gregoria(tim_sadcp(end))));
disp([' Time is now: ',datestr(now)])
eval(['save ' SADCP_OUTPUT_FILE ' tim_sadcp lat_sadcp lon_sadcp u_sadcp v_sadcp z_sadcp']);