%======================================================================
%                    M M P 2 A N T S . M 
%                    doc: Tue Apr 21 17:31:52 2009
%                    dlm: Tue Oct 13 10:04:17 2009
%                    (c) 2009 A.M. Thurnherr
%                    uE-Info: 16 54 NIL 0 0 72 2 2 4 NIL ofnI
%======================================================================

% HISTORY:
%	Apr 21, 2009: - created
%	May 15, 2009: - added meta data (requires modified [mmp_pgrid_main.m])
%				  - added ancillary data (hdg, etc)
%				  - added depth calc
%	May 18, 2009: - added export of raw ACM file
%				  - BUG: hdg actually was vel_dir
%	Oct 13, 2009: - adapted to new struct2ANTS (V4.0: dependencies)

function mmp2ANTS(plist,epoch)

if nargin ~= 2											% default epoch could be dangerous
	error('Usage: mmp2ANTS(plist,epoch)');				% for year-end-spanning data sets
end

for p = plist

	%--------------------------------
	% EXPORT PRESSURE-GRIDDED PROFILE
	%--------------------------------

	okay = 1;												% attempt to load data
	eval(sprintf('load grd%04d',p),'okay = 0;');

	if okay

		disp(sprintf('exporting gridded profile %04d...',p));		% process
		
		%----------
		% META DATA
		%----------

		prof.MMP_id = MMP_id;
		prof.ACM_id = ACM_id;
		prof.experiment = experiment_name;

		prof.min_press = gmin;
		prof.max_press = gmax;
		prof.press_step = gstep;

		prof.ACM_calfile = acm_cal_filename;

		prof.vel_l_bound = vel_bounds(1) / 100;
		prof.vel_u_bound = vel_bounds(2) / 100;

		pos = sscanf(mooring_position,'[%d %f,%d %f]');
		if pos(1) >= 0, prof.lat = pos(1) + pos(2)/60;
		else, 			prof.lat = pos(1) - pos(2)/60;
		end
		if pos(3) >= 0, prof.lon = pos(3) + pos(4)/60;
		else, 			prof.lon = pos(3) - pos(4)/60;
		end

		prof.magdecl = 180 * mag_dev/pi;

		eval(sprintf('prof.start_dn%02d = %f;',epoch-2000,...	% start/stop times in ANTS dn format
			startdaytime - datenum(sprintf('1-1-%d 00:00:00',epoch)) + 1));
		eval(sprintf('prof.stop_dn%02d = %f;',epoch-2000,...
			stopdaytime - datenum(sprintf('1-1-%d 00:00:00',epoch)) + 1));

		%-----
		% DATA	
		%-----

		stimes = ctimave - datenum(sprintf('1-1-%d 00:00:00',epoch)) + 1;	% bin-averaged time in dn format
		igood = isfinite(stimes);
		eval(sprintf('prof.dn%02d = stimes(igood);',epoch-2000));
	
		prof.press = pgrid(igood);									% center of bins used for gridding
		if exist('sw_dpth','file')
			prof.depth = sw_dpth(prof.press',prof.lat);
		end
	    
		prof.binned_press = pave(igood);							% bin-averaged CTD data
		prof.temp	= tave(igood);
		prof.salin	= s_ave(igood);
		prof.theta0 = thetave(igood);
		prof.sigma0 = sigthave(igood);
		prof.cond	= cave(igood);
		prof.CTD_Nsamp = cscan2(igood) - cscan1(igood);
	    
		prof.u		= uave(igood) / 100;							% bin-averaged ACM data in [m/s]
		prof.v		= vave(igood) / 100;
		prof.w		= wave(igood) / 100;

		tmp = 180 * mvpdirave(igood)/pi;
		prof.hdg	= (270-tmp)-(fix((270-tmp)/360)*360);			% standard heading btw 0 and 360
		prof.pitchx	= aTxave(igood);
		prof.pitchy	= aTyave(igood);

		prof.dpdt	= dpdtave(igood);
		prof.ACM_Nsamp = ascan2(igood) - ascan1(igood);
	    
	    struct2ANTS(prof,sprintf('grd%04d.mat',p),sprintf('%04d.prof',p));
	    
	end % if grd file loaded

	%--------------------
	% EXPORT RAW PROFILES
	%--------------------

	okay = 1;														% attempt to load data
	eval(sprintf('load raw%04d',p),'okay = 0;');

	if okay

		disp(sprintf('exporting raw data %04d...',p));				% process
		
		ACM.Vab = Vab;
		ACM.Vcd = Vcd;
		ACM.Vef = Vef;
		ACM.Vgh = Vgh;

		ACM.aHx = aHx;
		ACM.aHy = aHy;
		ACM.aHz = aHz;
	    
		ACM.aTx = aTx;
		ACM.aTy = aTy;

	    struct2ANTS(ACM,sprintf('raw%04d.mat',p),sprintf('%04d.ACM',p));
	    
	end % if grd file loaded

end % for p