
toptsRMV.iKCBand = [2310 2311 2320 2321 2350 2351 2352 2353 2354 2355];
toptsRMV.iISO =    [   1    2    1    2    1    2    3    1    1    4];  
toptsRMV.iL =      [   4    4    2    2    1    1    1    3    5    1];
toptsRMV.iU =      [  24   24   16   16    9    9    9   23   25    9];
topts.HITRAN = '/asl/data/hitran/h2k.oldiso/';

idgas = 2;
f1 = 2000;
f2 = 2600;
fnameout1 = '/carrot/s1/sergio/AIRSCO2/CO2_BANDS_PARAM/co2_allpqr_';
fnameout2 = '/carrot/s1/sergio/AIRSCO2/CO2_BANDS_PARAM/co2_allpqrH96_';
for ii = 1 : length(toptsRMV.iU)
  fprintf(1,'band = %6i \n\n',toptsRMV.iKCBand(ii));

  file1 = ['load ' fnameout1 num2str(toptsRMV.iKCBand(ii)) '.mat'];
  eval(file1);
  freq2000 = freq; stren2000 = stren; w2000 = w; 
  ps2000 = p_shift;                   w2000_s = w_s;
                                      w2000_temp = w_temp;
 

  file2 = ['load ' fnameout2 num2str(toptsRMV.iKCBand(ii)) '.mat'];
  eval(file2);
  freq1996 = freq; stren1996 = stren; w1996 = w; 
  ps1996 = p_shift;                   w1996_s = w_s;
                                      w1996_temp = w_temp;
 
  plot(freq2000./freq1996,'.-'); title('freq'); pause(1);
  plot(stren2000./stren1996,'.-'); title('stren'); pause(1);
  plot(w2000./w1996,'.-'); title('w'); pause(1);
  plot(w2000_s./w1996_s,'.-'); title('w_s'); pause(1);
  plot(w2000_temp./w1996_temp,'.-'); title('w_temp'); pause(1);
  semilogy(freq2000,stren1996,'.-')
  pause;
  end
  