% make a complete spectroscopic error calculation;
% call structure 
Q=qpack('spectro_sample');
% calculate H, Se, and Se
qp_H(Q); 
qp_Se(Q);
qp_Sx(Q);

% calculate Kb, spectroscopic error and plots the results.
qpcls_invchar(Q);