Createdepth ac9 pnb.m

function out=createdepth_ac9_pnb(ctd_strip,ac9_strip,pb)
  %have olga make unbinned file in stripped form
  %createdepth_ac9_pnb('X990614D.strip','a990614d.lcd')
                         % strip the ac9 lcd
  %cmd1=sprintf('!bboprecal -r 1depth 1 0 %s ac9.raw',ac9_lcd);
  %eval([cmd1]);
                         % strip the ctd lcd, strip by hand for now
  %eval(['!bbopstrip -sz 1mer_time 1depth ac9.raw ac9.strip']);
                         % load stripped ac9 pb = 'pb218'
pb;
  ac9=read_strip(ac9_strip);
                         % load stripped ctd  
  subplot(211)
  plot(ac9(:,1),ac9(:,2))
  title('AC9 Depth v. Time');
  %cd ..%/..
  sprintf('second time to read_strip')
  
  ctdval = read_strip(ctd_strip);
  %cd ..%/..
  pwd
  ctd=excise(ctdval);
  subplot(212)
  plot(ctd(:,1),ctd(:,2))
  title('Pressure v. CTD scan line');
  disp('calling fixdepth_ac9')
  %cd ..
  [out.scale,out.offset,out.matrix]=fixdepth_ac9_pnb(ctd,ac9);
  disp('createdepth_ac9 finished')
  return

createdepth_ac9_pnb