anchor.diff:           call readi(itemp)
            call readf(energy)
            !
  -         if (energy>enermax) then
  -           i = i - 1
  -           cycle
  -         endif
  +         !if (energy>enermax) then
  +         !  i = i - 1
  +         !  cycle
  +         !endif
            !
            energies(i) = energy
            !
  
positive.code:      real(rk),intent(in) :: Jf,Ji
       real(rk) :: halfwidth,gamma_,n_,f
       integer(ik) :: ispecies,Jpp,Jp
       !
       halfwidth = 0
       !
  
negative.code:        case ("HWHM","HALFWIDTH")
            !
            call readf(halfwidth)
            species(1)%gamma = halfwidth
            !
            if_halfwidth_defined = .true.
            !
         case ("IOFFSET")
            !
  