anchor.diff:              use_width_offset = .true.
             endif
             !
  +          if (offset<=0) stop 'offset cannot be negative or zero'
  +          !
          case ("MASS")
             !
             call readf(meanmass)
  
positive.code:                 !
                  !if (abscoef<cutoff) cycle
                  !
                  ipoint =  max(nint( ( tranfreq-freql)/dfreq )+1,1)
                  !
                  intens_omp(ipoint,iomp) = intens_omp(ipoint,iomp)+abscoef
                  !
  
negative.code:     real(rk)     :: freql=-small_
      real(rk)     :: freqr= 200000.0
      real(rk)     :: dfreq
      real(rk)     :: offset = -100000.0_rk
    end type gridT
    !
    type(selectT),save :: upper(filtermax),lower(filtermax)
  