anchor.diff:  	   do i=0,JmaxAll
   	   	do j= max(0,i-1),min(JmaxAll,i+1) 
   	   		
  -	   		call fast_voigt%generate_indices(get_Voigt_gamma_n(Nspecies,real(i,rk),real(j,rk)),gamma_idx(i,i-j),JmaxAll)
  +	   		call fast_voigt%generate_indices(gamma_comb(i,i-j),gamma_idx(i,i-j),JmaxAll)
   	   	enddo
   	   enddo  
      	 
  
positive.code:      
      ! f = gamma_comb(Jpp,Jp-Jpp)
       !return
       
       
       !
  
negative.code: 
     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
  