anchor.diff:        !
       endif
       !
  +    if (proftype(1:3)=='BIN'.or.proftype(1:3)=='MAX') then
  +      offset = 0
  +    endif
  +    !
       if (proftype(1:3)/='BIN'.and.microns) then
         !
         write(out,"('Microns or um is currently only working with BIN, not ',a)") proftype(1:3)
         call report ("Illegal profile for Microns"//trim(w),.true.)
  +      !
  +    endif
  +    !
  +    if (proftype(1:3)=='BIN'.and.microns) then
  +      !
  +      proftype = 'BIN-MICRON'
  +      !
  +    endif
  +    !
  +    if (proftype(1:3)=='BIN'.and.use_resolving_power) then
  +      !
  +      proftype = 'BIN-R'
         !
       endif
       !
  
positive.code: 
  
  
    subroutine do_gf_oscillator_strength_France(Ntrans,iso,Jf,Ji,nu_ram,energyi_ram,abscoef_r)
      !
      implicit none 
      !
      integer(ik),intent(in) :: Ntrans,iso
      real(rk),intent(in) :: Jf,Ji,nu_ram(:),energyi_ram(:),abscoef_r(:)
      integer(ik) :: i,iloggf
      real(rk) :: lambda,energyi
  
     call get_Voigt_gamma_n(Nspecies,Jf,Ji)
  
     do i = 1,Ntrans
         !
         lambda = 1e7/nu_ram(i)
         !
         energyi = energyi_ram(i)
         !
         iloggf = nint(log(abscoef_r(i)))
         !
         write(out,"(i12,1x,f12.6,1x,i6,1x,f12.4,1x,f11.4,1x,f11.4,1x,f11.4)") &
                     iso,lambda,iloggf,energyi,&
                     species(1)%gamma,species(2)%gamma,species(1)%n
         !
  
     enddo
  
  !ielion, wl, gf-value, xi, igr, igs, igw
  !where ielion: internal code of the molecule for Phoenix, 
  !wl:wavelength, 
  !xi:lower state energy, 
  !gf-value: log gf coded in integer format,
  !igr: natural width constant, 
  !igs: Stark broadening constant, 
  !igw: van der Waals damping constants.
  !igr for gamma0 and igs for n (Voigt parameters) 
  
  
    end subroutine do_gf_oscillator_strength_France
  
  
  
negative.code: 
    end subroutine  do_Voigt
    !
  
  
    subroutine  do_rect(tranfreq,abscoef,dfreq,freq,offset,freql,intens)
       !
       implicit none
       !
       real(rk),intent(in) :: tranfreq,abscoef,dfreq,offset,freql
       real(rk),intent(in) :: freq(:)
       real(rk),intent(out) :: intens(:)
       real(rk) :: tranfreq_i
       integer(ik) :: ib,ie,ipoint
        ! 
        ib =  max(nint( ( tranfreq-offset-freql)/dfreq )+1,1)
        ie =  min(nint( ( tranfreq+offset-freql)/dfreq )+1,npoints)
        !
        !$omp parallel do private(ipoint,tranfreq_i) shared(intens) schedule(dynamic)
        do ipoint=ib,ie
           !
           tranfreq_i = freq(ipoint)
           !
           intens(ipoint)=intens(ipoint)+abscoef/dfreq
           !
        enddo
        !$omp end parallel do      
  
    end subroutine  do_rect
  
  
    !
    subroutine do_stick(sunit,nswap,nlines,maxitems,energies,Jrot,gtot,ilevelf_ram,ileveli_ram,abscoef_ram,acoef_ram,nchars_quanta,quantum_numbers)
  
  