anchor.diff:  
      CALL BD_CopyInput(u, u_tmp, MESH_NEWCOPY, ErrStat2, ErrMsg2)
         CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  +!WRITE(*,*) 'u_Inic'
  +!DO k=1,3
  +!WRITE(*,*) u%RootMotion%Orientation(k,:,1)
  +!ENDDO
  +!WRITE(*,*) 'END u_Inic'
      CALL BD_InputGlobalLocal(p,u_tmp,ErrStat2,ErrMsg2)
         CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
      CALL BD_CalcIC(u_tmp,p,x,OtherState,ErrStat2,ErrMsg2)
         CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  -   x%q(:) = 0.0D0
   
      ! Define initial guess for the system outputs here:
   
  
positive.code:       CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
     CALL BD_CrvExtractCrv(p%GlbRot,temp_glb,ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
     DO i=1,p%elem_total
         IF( i .EQ. 1) THEN
             k = 1
  
negative.code:    ! Incorporate boundary conditions
     CALL BD_BoundaryGA2(x,p,u_interp,t+p%dt,OtherState,ErrStat2,ErrMsg2)
           call SetErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,RoutineName)
  !WRITE(*,*) 'x%q'
  !WRITE(*,*) x%q(:)
  !WRITE(*,*) 'x%dqdt'
  !WRITE(*,*) x%dqdt(:)
  !WRITE(*,*) 'OS%acc'
  !WRITE(*,*) OS_tmp%acc
  !WRITE(*,*) 'OS%xcc'
  !WRITE(*,*) OS_tmp%xcc
  
     ! find x, acc, and xcc at t+dt
     CALL BD_DynamicSolutionGA2( p%uuN0,x%q,x%dqdt,OtherState%acc,OtherState%xcc,&
  