anchor.diff:    	real(rk),allocatable		::	temp_humlicek(:)
     	
     	
  -  	center_point = (nu0-start_nu)/this%m_res
  -  	middle_shift = (center_point - this%m_middle_point)
  +  	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)
     	
  
positive.code: 	right_start = max(this%m_middle_point + dist, ib_rel)
  	right_end = ie_rel
  	
  
  	
  	!If we have any humlicek points
  
negative.code:   		this%m_mag = this%m_mag + this%m_voigt_grid(ido)
    		nu = nu+ this%m_res
    	enddo
    	
    	this%m_mag = this%m_mag*this%m_res
    	
    	this%constructed = .true.
    	
  