anchor.diff:             iswap_ = 1
              abscoef_ram = 0
              !
  -           !omp parallel do private(iswap,indexf,indexi,acoef,ilevelf,ileveli,ifilter,energyf,energyi,jf,ji,tranfreq,tranfreq0,offset,abscoef,cutoff)
  +           !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_
                !
  
positive.code:               enddo
                !              
              enddo
              !$omp parallel end do
              !
              do iomp = 1,N_omp_procs
                 intens(:) = intens(:) + intens_omp(:,iomp)
  
negative.code:           call readf(enermax)
            !
         case('GAUSSIAN','GAUSS','DOPPL','DOPPLER','RECT','BOX','BIN','STICKS','STICK','GAUS0','DOPP0',&
              'LOREN','LORENTZIAN','LORENTZ','MAX','VOIGT','PSEUDO','PSE-ROCCO','PSE-LIU','VOI-QUAD','PHOENIX')
            !
            proftype = trim(w)
            !
  