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

%                    S A V E P R O T . M 

%                    doc: Wed Jan  7 16:52:50 2009

%                    dlm: Wed Jan  7 16:53:04 2009

%                    (c) 2009 A.M. Thurnherr

%                    uE-Info: 11 0 NIL 0 0 72 0 2 8 NIL ofnI

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

% CHANGES BY ANT:

%   Jan  7, 2009: - tightened use of exist()

% save a processing protocol sheet for LADCP profiles

disp(' ')

disp(['Writing processing protocol ',f.res,'.txt'])

% open file

fid = fopen([f.res,'.txt'],'wt');

fprintf(fid,'    LADCP processing protocol \n',[]);

fprintf(fid,' \n',[]);

fprintf(fid,'Station name         : %s \n',p.name);

fprintf(fid,' \n',[]);

fprintf(fid,'LADCP input files    : \n',[]);

if isfield(f,'ladcpdo')

  fprintf(fid,'down                 : %s \n',f.ladcpdo);

else

  fprintf(fid,'down                 : %s \n','none');

end

if isfield(f,'ladcpup')

  fprintf(fid,'up                   : %s \n',f.ladcpup);

else

  fprintf(fid,'up                   : %s \n','none');

end

fprintf(fid,' \n',[]);

fprintf(fid,'Command file info    : %s \n',[]);

fprintf(fid,' \n',[]);

if exist('pinfo','var')

  fprintf(fid,'Start time (CTD)     : %4d %2d %2d %2d %2d %2d \n',...

	pinfo.ladcpprofstart);

  fprintf(fid,'End time (CTD)       : %4d %2d %2d %2d %2d %2d \n',...

	pinfo.ladcpprofend);

  fprintf(fid,'Rec start (LADCP)    : %4d %2d %2d %2d %2d %2d \n',...

	pinfo.ladcpstart);

  fprintf(fid,'Rec end (LADCP)      : %4d %2d %2d %2d %2d %2d \n',...

	pinfo.ladcpend);

end

fprintf(fid,'Prof start (LADCP)   : %4d %2d %2d %2d %2d %2d \n',...

	p.time_start);

fprintf(fid,'Prof end (LADCP)     : %4d %2d %2d %2d %2d %2d \n',...

	p.time_end);

[strlat,strlon] = pos2str([p.poss(1)+p.poss(2)/60,p.poss(3)+p.poss(4)/60]);

fprintf(fid,'Start position       : %s   ',strlat);

fprintf(fid,'%s \n',strlon);

[strlat,strlon] = pos2str([p.pose(1)+p.pose(2)/60,p.pose(3)+p.pose(4)/60]);

fprintf(fid,'End position         : %s   ',strlat);

fprintf(fid,'%s \n',strlon);

fprintf(fid,'Magnetic deviation   : %f \n',p.drot);

if isfield(p,'timoff')

  fprintf(fid,'Time offset          : %f \n',p.timoff);

else

  fprintf(fid,'Time offset          : 0.0 \n',[]);

end

fprintf(fid,'Depths               : %d %d %d \n',p.zpar);

fprintf(fid,' \n',[]);

fprintf(fid,'Processing thresholds  \n',[]);

fprintf(fid,' \n',[]);

fprintf(fid,'min percent good     : %f \n',p.pglim);

fprintf(fid,'max error velocity   : %f \n',p.elim);

fprintf(fid,'max w variability    : %f \n',p.wlim);

fprintf(fid,' \n',[]);

fprintf(fid,'Processing results  \n',[]);

fprintf(fid,' \n',[]);

fprintf(fid,'u ship speed   (m/s) : %f \n',p.uship);

fprintf(fid,'v ship speed   (m/s) : %f \n',p.vship);

fprintf(fid,'x displacement   (m) : %f \n',p.xdisp);

fprintf(fid,'y displacement   (m) : %f \n',p.ydisp);

fprintf(fid,'bottom found at  (m) : %f \n',p.zbottom);

fprintf(fid,'maximum instr dep(m) : %f \n',p.maxdepth);

fprintf(fid,'maximum prof dep (m) : %f \n',max(dr.z));

fprintf(fid,'bottom track mode    : %f \n',p.btrk_mode);

fprintf(fid,'bottom track thrsd   : %f \n',p.btrk_ts);

fprintf(fid,'bottom track range   : %f \n',p.btrk_range);

if p.btrk_used == 1

  fprintf(fid,'bottom track source  : RDI \n',[]);

elseif p.btrk_used == 2

  fprintf(fid,'bottom track source  : normal data \n',[]);

elseif p.btrk_used == -1

  fprintf(fid,'bottom track source  : not enough data \n',[]);

elseif p.btrk_used == 0

  fprintf(fid,'bottom track source  : none \n',[]);

end

fprintf(fid,'profile ensembles    : %d to %d \n',d.firstlastindx);

fprintf(fid,' \n',[]);

fprintf(fid,'Computer : %s \n',computer);

fprintf(fid,'Total processing time: %4.0f seconds \n',toc);

fclose(fid);