%======================================================================
%                    I F M _ L A D C P 2 A N T S . M 
%                    doc: Sun Jan 22 15:19:00 2006
%                    dlm: Mon Oct 12 22:58:20 2009
%                    (c) 2006 A.M. Thurnherr
%                    uE-Info: 68 0 NIL 0 0 72 2 2 4 NIL ofnI
%======================================================================
%
% export IfM-GEOMAR LADCP output to ANTS file
%
% USAGE: IfM_LADCP2ANTS(dr,p,outBaseName)
%

% HISTORY:
%	Jul 21, 2008: - created from [LDEO_LADCP2ANTS.m]
%	Jul 22, 2008: - cosmetics
%	Aug 26, 2008: - cosmetics
%	Oct 12, 2009: - adapted to new struct2ANTS

function [] = IfM_LADCP2ANTS(dr,p,obn)

	%----------------------------------------------------------------------
	% INVERSE SOLUTION
	%----------------------------------------------------------------------

	prof.name = dr.name;
	prof.cruise = p.cruise_id;
	prof.software = p.software;
	prof.magdecl = p.magdev;	% always zero in V10.6
	prof.procdir = pwd;
	
	prof.start_date  = sprintf('%d/%02d/%02d',p.time_start(1),p.time_start(2),p.time_start(3));
	prof.start_time	 = sprintf('%02d:%02d:%02d',p.time_start(4),p.time_start(5),p.time_start(6));

	prof.end_date    = sprintf('%d/%02d/%02d',p.time_end(1),p.time_end(2),p.time_end(3));
	prof.end_time	 = sprintf('%02d:%02d:%02d',p.time_end(4),p.time_end(5),p.time_end(6));

	prof.median_time = sprintf('%02d:%02d:%02d',dr.date(4),dr.date(5),dr.date(6));

	prof.lat 		= dr.lat; 			  prof.lon  	  = dr.lon;				% (start+end)/2
	prof.mean_lat 	= mean(dr.shiplat);   prof.mean_lon   = mean(dr.shiplon);
	prof.median_lat = median(dr.shiplat); prof.median_lon = median(dr.shiplon);
	prof.start_lat	= dr.shiplat(1);	  prof.start_lon  = dr.shiplon(1);
	prof.end_lat	= dr.shiplat(end);	  prof.end_lon	  = dr.shiplon(end);
	i_bot = find(dr.zctd==min(dr.zctd));
	prof.bot_lat	= dr.shiplat(i_bot);  prof.bot_lon    = dr.shiplon(i_bot);

	prof.depth = dr.z;
	prof.max_depth = max(prof.depth);
	
	prof.u	   = dr.u;
	prof.dn_u  = dr.u_do;
	prof.up_u  = dr.u_up;
	prof.v     = dr.v;
	prof.dn_v  = dr.v_do;
	prof.up_v  = dr.v_up;
	
	prof.u_fromshear = dr.u_shear_method;
	prof.v_fromshear = dr.v_shear_method;

	prof.samp  = dr.nvel;
	prof.err   = dr.uerr;
	prof.range = dr.range;
	if existf(dr,'ensemble_vel_err')
		prof.ensemble_vel_err = dr.ensemble_vel_err;
	end

	struct2ANTS(prof,sprintf('%s.mat',f.res),sprintf('%s.prof',obn));

	%----------------------------------------------------------------------
	% SADCP
	%----------------------------------------------------------------------

	if existf(dr,'u_sadcp')

		SADCP.name = prof.name;
		SADCP.cruise = prof.cruise;
		SADCP.software = prof.software;
	    SADCP.magdecl = prof.magdecl;
		SADCP.procdir = prof.procdir;
		SADCP.date = sprintf('%d/%02d/%02d',dr.date(1),dr.date(2),dr.date(3)); % median
		SADCP.time = prof.median_time;
		SADCP.lat  = prof.lat; SADCP.lon  = prof.lon;
		SADCP.start_lat  = prof.start_lat; SADCP.start_lon  = prof.start_lon;
		SADCP.end_lat  = prof.end_lat; SADCP.end_lon  = prof.end_lon;
		SADCP.mean_lat  = prof.mean_lat; SADCP.mean_lon  = prof.mean_lon;
		SADCP.median_lat  = prof.median_lat; SADCP.median_lon  = prof.median_lon;

		SADCP.depth = dr.z_sadcp;
		SADCP.max_depth = max(dr.z_sadcp);

		SADCP.u		= dr.u_sadcp;
		SADCP.v		= dr.v_sadcp;
		if existf(dr,'uerr_sadcp')			% V7 does not have this
			SADCP.err	= dr.uerr_sadcp;
		end
		
		struct2ANTS(SADCP,sprintf('%s.mat',f.res),sprintf('%s.SADCP',obn));
		
	end

	%----------------------------------------------------------------------
	% BT
	%----------------------------------------------------------------------

	if existf(dr,'ubot')
		BT.name = prof.name;
		BT.cruise = prof.cruise;
		BT.software = prof.software;
		BT.magdecl = prof.magdecl;
		BT.procdir = prof.procdir;
		BT.date = sprintf('%d/%02d/%02d',dr.date(1),dr.date(2),dr.date(3)); % median
		BT.time = prof.median_time;
		BT.lat  = prof.bot_lat;
		BT.lon  = prof.bot_lon;

		BT.depth = dr.zbot;
		BT.max_depth = max(dr.zbot);

		BT.u		= dr.ubot;
		BT.v		= dr.vbot;
		BT.err		= dr.uerrbot;
		
		struct2ANTS(BT,sprintf('%s.mat',f.res),sprintf('%s.BT',obn));
	end