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  5 15:06:38 2007

%                    (c) 2004 A.M. Thurnherr

%                    uE-Info: 63 33 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

%======================================================================

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(' 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']);