anchor.diff:  	right_start = max(this%m_middle_point + dist, ib_rel)
   	right_end = ie_rel
   	
  +	!print *,mag
   
   	
   	!If we have any humlicek points
  
positive.code:   	this%m_lorentz_cutoff = pLorentzCutoff
   	this%normalize = pNormalize
    	this%m_Npoints = int(2.0*this%m_lorentz_cutoff/this%m_res,ik);
    	this%m_middle_point = (this%m_Npoints/2)
    	!Construct the voigt grid
    	
    	this%m_mag = 0.0
  
negative.code:   	real(rk),allocatable		::	temp_humlicek(:)
    	
    	
    	center_point = (nu0-start_nu)/this%m_res + 1
    	middle_shift = (center_point - this%m_middle_point) + 1
    	
    	
    	
    	dist = max(DISTANCE_MAGIC_NUMBER/this%m_res,3.0)
    	
  