anchor.diff:     !
      if (trim(enrfilename)/="none".and..not.hitran_do) then
   
  -      if (verbose>=2) print('(a,a,a/)'),'Reading energies from ',trim(enrfilename),'...'
  +      if (verbose>=2) print('(a/)'),'Reading energies from .states ..'
         !
         write(ioname, '(a)') 'Energy file'
         call IOstart(trim(ioname),enunit)
  
positive.code:      real(rk),intent(in) :: Jf,Ji
       real(rk) :: halfwidth,gamma_,n_,f
       integer(ik) :: ispecies,Jpp,Jp
       !
       ! if no species are specified in input the hwm value will be used for halfwidth instead 
       !if (Nspecies==0) then 
       !  f = hwhm
       !  return
       !endif
       !
       halfwidth = 0
       !
       do ispecies =1,Nspecies
         !
  
negative.code:             enddo
              !$omp end parallel do
              !
  
              !
          case ('DOPPL')
              !
  
              !
              !$omp parallel do private(iomp,iswap,abscoef,tranfreq,halfwidth0) shared(intens_omp) schedule(dynamic)
              do iomp = 1,N_omp_procs
  