anchor.diff:  			nu=freq(ido)
   			hum_res = voigt_humlicek(nu,nu0,gammaD,this%m_gammaL)
   			temp_humlicek(ido-start_dist+1) = hum_res
  -			mag = mag - this%m_voigt_grid(ido -middle_shift)
  -			mag = mag + hum_res
  +			mag = mag - this%m_voigt_grid(ido -middle_shift+1)*this%m_res
  +			mag = mag + hum_res*this%m_res
   		enddo
   		
   		if(.not. this%normalize) mag = 1.0
   		
   		intens(start_dist:end_dist) = intens(start_dist:end_dist) + temp_humlicek(:)*abscoef/mag
   		
  -		
  +		!print *,mag
   		deallocate(temp_humlicek)
     	endif
     	if(.not. this%normalize) mag = 1.0
  
positive.code: 	
  	mag = this%m_mag
  	
  	start_dist = max(center_point - dist, ib)
  	end_dist = min(center_point + dist, ie)
  		
  	left_start = ib_rel
  
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)
    	
  