anchor.diff:  
   ! x(5) corresponds to the second 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
   
   ! x(6) corresponds to the third harmonic oscillator pseudo-frequency
     ind(6) = 3        ! Has both an upper and lower bound
  -  bl(6)  = 400.     ! Lower bound is 400 cm^-1
  +  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) = 500.0 !Harm. osc. frequency
  +  x(1) = 300.0 !Harm. osc. frequency
     x(2) = FLOOR(REAL(N_rot)/2.) !Degeneracy of first pseudo-frequency: 1/2 N_rot
     x(3) = 500.0 ! First pseudo-barrier height: 500 cm^-1
     x(4) = 3000. ! Second pseudo-barrier height: 3000 cm^-1
  -  x(5) = 1000.0 !Harm. osc. frequency
  -  x(6) = 2000.0 !Harm. osc. frequency
  +  x(5) = 1500.0 !Harm. osc. frequency
  +  x(6) = 3000.0 !Harm. osc. frequency
   
   !  Tell how much storage we gave the solver.
     iwork(1) = lwork
  
positive.code: 
  ! x(1) corresponds to the first harmonic oscillator pseudo-frequency
    ind(1) = 3        ! Has both an upper and lower bound
    bl(1)  = 400.     ! Lower bound is 400 cm^-1
    bu(1)  = 4000.    ! Upper bound is 4000 cm^-1
  
  ! x(2) corresponds to the second harmonic oscillator pseudo-frequency
    ind(2) = 3        ! Has both an upper and lower bound
    bl(2)  = 400.     ! Lower bound is 400 cm^-1
    bu(2)  = 4000.    ! Upper bound is 4000 cm^-1
  
  ! x(3) corresponds to the degeneracy of the first hindered rotor
  
negative.code: 
  ! x(2) corresponds to the first harmonic oscillator pseudo-frequency
    ind(2) = 3        ! Has both an upper and lower bound
    bl(2)  = 180.     ! Lower bound is 180 cm^-1
    bu(2)  = 4000.    ! Upper bound is 4000 cm^-1
  
  ! x(3) corresponds to the second harmonic oscillator pseudo-frequency
    ind(3) = 3        ! Has both an upper and lower bound
    bl(3)  = 180.     ! Lower bound is 180 cm^-1
    bu(3)  = 4000.    ! Upper bound is 4000 cm^-1
  
  ! x(4) corresponds to the degeneracy of the first hindered rotor
  