anchor.diff:  
      allocate(gamma_idx(0:JmaxAll,-1:1))
      call ArrayStart('gamma_idx',info,size(gamma_idx),kind(gamma_idx))
  +   allocate(gamma_comb(0:JmaxAll,-1:1))
  +   call ArrayStart('gamma_comb',info,size(gamma_comb),kind(gamma_comb))
  +   do i=0,JmaxAll
  +	  do j= max(0,i-1),min(JmaxAll,i+1)    
  +	  	gamma_comb(i,i-j) = get_Voigt_gamma_val(Nspecies,real(i,rk),real(j,rk))
  +	  enddo
  +   enddo
  +   
  +   
  +   call ArrayStart('gamma_idx',info,size(gamma_idx),kind(gamma_idx))
      if(proftype(1:5) == 'VOI-F') then
      	
      	gamma_idx = 1
  
positive.code:      
      ! f = gamma_comb(Jpp,Jp-Jpp)
       !return
       
       
       !
  
negative.code:      !
       f = halfwidth
       !
    end function get_Voigt_gamma_val
    !
    function voigt_faddeeva(nu,nu0,alpha,gamma) result(f)
     !
  