anchor.diff:     temp_L2(4:6,p%ngp*p%elem_total+2) = p%uuN0(temp_int-2:temp_int,p%elem_total)
      DEALLOCATE(temp_GLL)
      DEALLOCATE(SP_Coef)
  +!WRITE(*,*) 'uuN0'
  +!WRITE(*,*) p%uuN0(:,1)
   
      ! Compute sectional propertities ( 6 by 6 stiffness and mass matrices)
      ! at Gauss points
  
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:    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
  