anchor.diff:            !
             if (mod(npoints,2)==0) npoints = npoints + 1
             !
  -        case ("ABSORPTION","LIFETIME","OSCILLATOR")
  +        case ("ABSORPTION","LIFETIME","EMISSION")
             !
             specttype = trim(w)
             !
  
positive.code:               enddo
                !              
              enddo
              !$omp parallel end do
              !
              do iomp = 1,N_omp_procs
                 intens(:) = intens(:) + intens_omp(:,iomp)
  
negative.code:            iswap_ = 1
             abscoef_ram = 0
             !
             !omp parallel do private(iswap,indexf,indexi,acoef,ilevelf,ileveli,energyf,energyi,ifilter,jf,ji,tranfreq,tranfreq0,offset,abscoef,cutoff)&
             !omp&  schedule(static) reduction(+:iswap_,intband) shared(ilevelf_ram,ileveli_ram,abscoef_ram,acoef_ram,nu_ram,Asum)
             loop_swap : do iswap = 1,nswap_
               !
  