anchor.diff:                case('NPOINTS')
                   !
                   call readi(grid(igrid)%npoints)
  +                !
  +              case('OFFSET')
  +                !
  +                call readf(grid(igrid)%offset)
                   !
                 case default
                   !
  
positive.code:     real(rk)     :: freql=-small_
      real(rk)     :: freqr= 200000.0
      real(rk)     :: dfreq
    end type gridT
    !
    type(selectT),save :: upper(filtermax),lower(filtermax)
  
negative.code:             !
              select case(w)
              !
              case('TEMPMAX','MAXTEMP','MAX-TEMPERATURE',"TMAX")
                !
                call readf(maxtemp)
                !
                if (maxtemp<small_) call report ("Illegal Max Temperature maxtem=0 "//trim(w),.true.)
                !
              case ("NTEMPS","NPOINTS","NSTEPS","NUMBER-OF-TEMPERATURES","N")
                !
                call readi(npoints)
                !
  