for ff = 605 : 25 : 2805
  fin = ['/asl/s1/sergio/H2008_RUN8_NIRDATABASE/IR_605_2830_H08_CO2/kcomp/'];
  fin = [fin 'cg2v' num2str(ff) '.mat'];
  ee1 = exist(fin);

  fout = ['cg2v' num2str(ff) '.mat'];
  ee2 = exist(fout);
  if ee1 > 0 & ee2 > 0
    a1 = load(fin);
    a2 = load(fout);


    %plot(a1.B ./ a2.B)

    a1_k = squeeze(a1.kcomp(:,:,6));
    a2_k = squeeze(a2.kcomp(:,:,6));

    plot(abs(a1_k'),1:100,'b',abs(a2_k'),1:100,'r')
    whos a1_k a2_k

    title(num2str(ff)); pause(0.1);
  end
end