This routine uses the output from combine_waversols, to fit the wavelength residulas from a a simple linaer fit using higher order terms. Additionally, a spectral line regions file is created that can be used to overplot the positions of the final calibrating lines on a ds9 display of the spectral image.

% 
% waversol_refine -h 

usage: waversols_refine 056LL [-v] [-h] [-c] 
arg1 = device designation (056LL)   
Additional options: 
  -h = just show usage message  
  -v = run in verbose/debug mode 
  -c = just run the cleanup function