anchor.diff:            call read_line(eof) ; if (eof) exit
             !call readu(w)
             !
  -          ifilter = 0 
  +          ifilter = 0
             !
             call readu(w)
             !
  
positive.code:       Jrot = 0
        gtot = 0
        !
        ! In case of specttype = partfunc the partition function will be computed for a series of temperatures. 
        ! npoints stands for the number of temperature steps, and the frequency limits as temperature limits
        !
        if (proftype(1:4)=='PART'.or.proftype(1:4)=='COOL') then
          !
          if (trim(enrfilename)=="none") then 
             !
             write(out,"('PARTITION FUNCTION ERROR: no energy file is specified')") 
             stop "no energy file is specified"
             !
          endif
          !
          dtemp=maxtemp/real(npoints,rk)
          !
          if (proftype(1:4)=='PART') then 
            !
            if (ipartf<0.or.ipartf>3) stop "illegal partfunc component, can be only 0,1,2,3"
            !
  
negative.code:        !
       enddo
       !
       close(tunit,status='keep')
       !
     elseif (any( trim(proftype(1:3))==(/'DOP','GAU','REC','BIN','BOX','LOR','VOI','MAX','PSE'/)) ) then
       !
       write(ioname, '(a)') 'Cross sections or intensities'
       call IOstart(trim(ioname),tunit)
  