anchor.diff:         case ("ENERMAX")
             !
             call readf(enermax)
  -          !
  -          ! energy threshold based on the uncertainty 
  -          !
  -       case('UNCERTAINTY_MAX','UNC_MAX','UNC')
  -          !
  -          uncertainty_filter = .true.
  -          !
  -          call readf(unc_threshold)
             !
          case('GAUSSIAN','GAUSS','DOPPL','DOPPLER','RECT','BOX','BIN','STICKS','STICK','GAUS0','DOPP0',&
               'LOREN','LORENTZIAN','LORENTZ','MAX','VOIGT','PSEUDO','PSE-ROCCO','PSE-LIU','VOI-QUAD','PHOENIX',&
  
positive.code:          call readi(gtot(i))
           call readf(jrot(i))
           !
           !if (uncertainty_filter) then 
           !  !
           !  call readf(unc)
           !  !
           !endif
           !
           if (mod(nint(jrot(i)*2),2)/=0) then
             write(b_fmt,"(f5.1)") jrot(i)             ; quantum_numbers(0,i) = adjustl(b_fmt)
           else
  
negative.code:    if (allocated(gtot)) call ArrayStop('gtot')
     if (specttype(1:4)=='PART') call ArrayStop('pf')
     !
     return
     !
     116 continue
     !
     write(out,"('Error read of unc: missing or wrongly specified UNC in QN')") 
     stop 'Error read of unc: missing or wrongly specified UNC in QN'
     !
    end subroutine intensity
    !
    !
  