anchor.diff:     INTEGER(IntKi)                                     :: ErrStat2   ! Temporary Error status
      CHARACTER(ErrMsgLen)                               :: ErrMsg2    ! Temporary Error message
      CHARACTER(*), PARAMETER                            :: RoutineName = 'BD_GA2'
  +   REAL(ReKi):: temp_3(3)
   !   INTEGER(IntKi)                                     :: i
   
      ! Initialize ErrStat
  
positive.code:    CALL BD_CalcIC(u_tmp,p,x,OtherState,ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  
     ! Define initial guess for the system outputs here:
  
     y%BldForce%Force(:,:)    = 0.0D0
  
negative.code:    ! enter known endpoints  [-1.0, 1.0]
     x(1) = -1.0D+00
     x(N1) = 1.0D+00
  
     DO i = 1, N1
        x_it = -COS(pi * FLOAT(i-1) / N) ! initial guess - chebyshev points
  