anchor.diff:  
   ! x(5) corresponds to the harmonic oscillator pseudo-frequency
     ind(5) = 3        ! Has both an upper and lower bound
  -  bl(5)  = 400.     ! Lower bound is 400 cm^-1
  +  bl(5)  = 180.     ! Lower bound is 180 cm^-1
     bu(5)  = 4000.    ! Upper bound is 4000 cm^-1
   
   !  Define the initial guess for the solution
  
positive.code: 
  ! x(3) corresponds to the first harmonic oscillator pseudo-frequency
    ind(3) = 3        ! Has both an upper and lower bound
    bl(3)  = 400.     ! Lower bound is 400 cm^-1
    bu(3)  = 4000.    ! Upper bound is 4000 cm^-1
  
  ! x(4) corresponds to the harmonic oscillator degeneracy
  
negative.code: 
  ! x(5) corresponds to the third harmonic oscillator pseudo-frequency
    ind(5) = 3        ! Has both an upper and lower bound
    bl(5)  = 180.     ! Lower bound is 180 cm^-1
    bu(5)  = 4000.    ! Upper bound is 4000 cm^-1
  
  ! x(6) corresponds to the fourth harmonic oscillator pseudo-frequency
    ind(6) = 3        ! Has both an upper and lower bound
    bl(6)  = 180.     ! Lower bound is 180 cm^-1
    bu(6)  = 4000.    ! Upper bound is 4000 cm^-1
  
  !  Define the initial guess for the solution
    x(1) = 250.0 !Harm. osc. frequency
    x(2) = FLOOR(REAL(N_vib)/3.) !Degeneracy
    x(3) = 800.0 !Harm. osc. frequency
    x(4) = FLOOR(REAL(N_vib)/3.) !Degeneracy
    x(5) = 1600.0 !Harm. osc. frequency
    x(5) = 4000.0 !Harm. osc. frequency
  
  !  Tell how much storage we gave the solver.
    iwork(1) = lwork
  