anchor.diff:  	
   	mag = this%m_mag
   	
  -	start_dist = max(center_point - dist, ib)
  +	start_dist = max(center_point - dist + 1, ib)
   	end_dist = min(center_point + dist, ie)
   		
   	left_start = ib_rel
  
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: 	right_start = max(this%m_middle_point + dist, ib_rel)
  	right_end = ie_rel
  	
  	!print *,mag
  
  	
  	!If we have any humlicek points
  