CPSbt E

% AC9BT
% FEB 98,FEB 99
% Stephen Wolfe
% Input:	none
% Look for:	list_prr
% Output: 	z*.lcd.?, z*.lcd.strip
% Assumptions:	fixed for either single casts or yo-yo's. Have to
%               pick something for at least 1 profile of each yo-yo
%               for the profile numbers to work correctly. Can't tell
%               if you skip a down-up pair.
%		
% Purpose:	pick tops and bottoms 
%		
% Calls:	bt_ac9
% Maintain: 	{filename, columnarraylist, disp_text}
% Uses:     	perfect for prr processing, ref: prr.ins
% Example:  	prrbt
% Notes:        This will run on NT if the call to bbopstrip is
%               commented out and the files are already stripped
%               Should label direction in castid table as well as
%               profile number. Work from left to right. If you
%               choose not to save points for a profile, you can
%               repick them and choose save and then go on.
% TO DO:        The archive scheme does not work due to changes in
%               naming convention. Doesn't create the log directory
%               if it's not there either.
%Erik Stassinos Apr 2013
clf;
list_prr=[];
prr_files=[];
not_prr_files=[];
tmpfile=fullfile(pwd,'goprrbt_tmp');
                                                % check for list_prr in pwd
[s,m]=unix('ls list_prr');
if isempty(m)
 msg=sprintf('Check your list_prr.');
 disp(msg)
 return
end
clear s m
						% load list_prr
fid=fopen('list_prr.txt','r');
while fid,
  check=fgetl(fid);
  if check == -1, break, end
  list_prr=[list_prr; check];
end
list_prr;
clear fid check
						% check presence of prr files
badfilectr=0;
for i=1:size(list_prr,1)
 [s,m]=unix(['ls ',list_prr(i,:)]);
 if s
  badfilectr=badfilectr+1;
 end
end
clear i
if badfilectr
 if badfilectr==1
  msg=sprintf('%i missing prr file.',badfilectr);
  disp(msg)
  error('Your prrlist includes a file not in pwd.')
 else
  msg=sprintf('%i missing prr files.',badfilectr);
  disp(msg)
  error('Your prrlist includes files not in pwd.')
 end
end
clear s m badfilectr msg
						% loop through 
for i=1:size(list_prr,1)
  prr_lcd=list_prr(i,:);
						% strip em
  prr_strip=[prr_lcd '.strip'];
  %eval(['!sbtprr_mob ',prr_lcd])
  %eval(['!bbopstrip -sz  1mer_time 1depth ' prr_lcd ' 'prr_strip]);
  %prrstrip(1,prr_lcd,'1mer_time','depth_ed') %changed 9-18-2012
  %prrstrip(1,prr_lcd,'1mer_time','1Depth') %need to change to accomodate screwed up LCD conversion 12-3-12
 % prrstrip(1,prr_lcd,'1mer_time','1Depth')
   prrstrip(1,prr_lcd,'FrameCount','LuZDepth')
  [fileA,headA]=read_strip(prr_strip);
  %eval(['!rm ''prr_strip']);
  
  [bt_pts]=bt_ac9(fileA);
  
  if ~isempty(bt_pts),
    %create header, castid matrix
    %castid  index   1mer_time       1depth
    %zub11008adt1     7.5000000e+01   2.1710000e+04   2.0443990e+02
    %zub11008adb1     1.2810000e+03   3.7655000e+05   1.3007170e+00 
    %castdir={ 'dt','db','ub','ut' };
    %castdir=[castdir castdir castdir castdir]; %  good for 4 yoyo's
    %if size(bt_pts,1)>size(castdir,2),
    %  error('prrbt: increase size of castdir.');
    %end
    dd=[];
    idir=[];
    dd=diff(bt_pts(:,3));
    dd=[dd; 0];
    id=1;
    fid=fopen( 'bt.pts', 'w' );
    fprintf(fid,'castid	index	1mer_time	1depth\n');
    for i=1:2:size(bt_pts,1),
      if dd(i)>0,
	if idir==1,
	  id=id+1;
	end
	idir=1;
	castdir1(i:i+1,:)=['dt'; 'db'];
      else
	idir=2;
	castdir1(i:i+1,:)=['ub'; 'ut'];
      end
      bt_castname=[ prr_lcd(1:9),castdir1(i,:),num2str(id) ];
      fprintf(fid,'%s %1.7e %1.7e %1.7e \n',bt_castname,bt_pts(i,:));
      bt_castname=[ prr_lcd(1:9),castdir1(i+1,:),num2str(id) ];
      fprintf(fid,'%s %1.7e %1.7e %1.7e \n',bt_castname,bt_pts(i+1,:));
      if idir==2,
	id=id+1;
      end
    end
    fclose(fid);
    
    %do this by calling a script that takes care of the awk and archive
    eval([ '!awk < ' prr_lcd ' > ' tmpfile...
	  ' ''$1=="<sampled_parameters>"{system("cat bt.pts")}{print $0}'''])
    eval([ '!mv ' tmpfile ' ' prr_lcd ]);

    %archive
    %cruisename=eval(['echo ' prr_lcd ' | cut -b 3-6'])
    %cruisename=prr_lcd(3:6);
    %%cruisename=prr_lcd(11:17); commented 9-18-2012
    %archive=[ pwd '/ARCHIVE/' ];
    %=>archive=[ '/home/data92/bbop/ARCHIVE/' ];
    %cp $lcd_file /home/data92/bbop/ARCHIVE/$cruisename/prr
   % eval(['!cp ' prr_lcd ' ' archive cruisename '/prr']);
    %echo $lcd_file has been archived to /home/data92/bbop/ARCHIVE/$cruisename/prr
    %=>eval(['!echo ' prr_lcd ' has been archived to 'archive cruisename '/prr']);

    prr_files=[prr_files; prr_lcd];
  else
    not_prr_files=[not_prr_files; prr_lcd];
  end
      
  %cleanup
  %clear prr_lcd bt_pts fileA headA prr_strip;

  %eval('!rm bt.pts');
end
prr_files
not_prr_files

back to C-OPS/CERBERUS Mfiles