anchor.diff:          !
           dtemp=maxtemp/real(npoints,rk)
           !
  -        if (specttype(1:4)=='PART') then 
  +        if (proftype(1:4)=='PART') then 
             !
             if (ipartf<0.or.ipartf>3) stop "illegal partfunc component, can be only 0,1,2,3"
             !
  
positive.code:               enddo
                !              
              enddo
              !$omp parallel end do
              !
              do iomp = 1,N_omp_procs
                 intens(:) = intens(:) + intens_omp(:,iomp)
  
negative.code:                !
               case ('LIFETIME')
                 !
                 !omp critical
                 !
                 if (Asum(ilevelf)<0) Asum(ilevelf) = 0 
                 !
                 Asum(ilevelf) = Asum(ilevelf) + acoef
                 !omp end critical
                 !
                 !print*,ilevelf,indexf,indexi,acoef
                 !
  