anchor.diff:  
   ! 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
  +  bl(1)  = 180.     ! Lower bound is 180 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
  +  bl(2)  = 180.     ! Lower bound is 180 cm^-1
     bu(2)  = 4000.    ! Upper bound is 4000 cm^-1
   
   ! x(3) corresponds to the degeneracy of the first hindered rotor
  
positive.code: !  Define the initial guess for the solution
  !  Note that the two pseudo-frequencies should be different so that the
  !  Jacobian of x(1) is not automatically zero
    x(1) = 800.0      !First harm osc. frequency: 800 cm^-1
    x(2) = 1200.0      ! Second harm osc. frequency: 1200 cm^-1
    x(3) = floor(REAL(N_rot)/2.)! Hind. rot. degen.:  1/2 the # of hind. rot.
    x(4) = 1000.0      ! First hind. rot. pseudo-barrier: 1000 cm^-1
    x(5) = 100.0     ! First hind. rot.osc. freq: 100 cm^-1
    x(6) = 2000.0      ! Second hind. rot. pseudo-barrier: 2000 cm^-1
  
  !  Tell how much storage we gave the solver.
    iwork(1) = lwork
  
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 first hindered rotor pseudo-barrier
  