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:    INTEGER(IntKi)                                     :: ErrStat2   ! Temporary Error status
     CHARACTER(ErrMsgLen)                               :: ErrMsg2    ! Temporary Error message
     CHARACTER(*), PARAMETER                            :: RoutineName = 'BD_GA2'
  !   INTEGER(IntKi)                                     :: i
  
     ! Initialize ErrStat
  
negative.code: !DO i=1,3
  !WRITE(*,*) u_interp%RootMotion%Orientation(i,:,1)
  !ENDDO
  !WRITE(*,*) u_interp%RootMotion%TranslationDisp(:,1)
  !WRITE(*,*) 'END u_interp'
  !CALL BD_CrvExtractCrv(TRANSPOSE(u_interp%RootMotion%Orientation(:,:,1)),temp_3,ErrStat2,ErrMsg2)
  !WRITE(*,*) temp_3(:)
  !WRITE(*,*) u_interp%RootMotion%RotationAcc(:,1)
                   
     ! GA2: prediction        
     CALL BD_TiSchmPredictorStep( x_tmp%q,x_tmp%dqdt,OS_tmp%acc,OS_tmp%xcc,             &
  