anchor.diff:             cycle
           end if
           !
  -        ib =  max(nint( ( tranfreq-halfwidth*ioffset-freql)/dfreq )+1,1)
  -        ie =  min(nint( ( tranfreq+halfwidth*ioffset-freql)/dfreq )+1,npoints)
  +        if (abscoef<abscoef_thresh) cycle
  +        !
  +        ib =  max(nint( ( tranfreq-offset-freql)/dfreq )+1,1)
  +        ie =  min(nint( ( tranfreq+offset-freql)/dfreq )+1,npoints)
           !
           select case (trim(proftype(1:5)))
               !
  
positive.code:         !
        endif 
        !
        indices(itemp) = i
        !
        j = 2*int(jrot(i))+1
        energy = energies(i)
  
negative.code:      !
       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)
  