r1 = rfits('r312_bir.fit');
i1 = rfits('r312_bii.fit');
x1 = r1+j*i1;
figure(1);
clf reset;
hs_bs(x1.', 'sqrt', 1E8, 'Baseline (wavelengths)', 'Baseline (wavelengths)', '', [0:112]*59027, [0:56]*59027);
print -dpsc r312.ps
r2 = rfits('r317_bir.fit');
i2 = rfits('r317_bii.fit');
x2 = r2+j*i2;
figure(2);
clf reset;
hs_bs(x2.', 'sqrt', 5E8, 'Baseline (wavelengths)', 'Baseline (wavelengths)', '', [0:112]*59027, [0:56]*59027);
print -dpsc r317.ps