anchor.diff:        DO J = 0,p%TipNode ! Loop through the blade nodes / elements
   
            LinAccES(:,J,K) = m%RtHS%LinAccESt(:,K,J)
  -!SP_start
            AngAccEK(:,J,K) = m%RtHs%AngAccEKt(:,J,K)
  -!SP_end
            DO I = 1,p%DOFs%NPSE(K)  ! Loop through all active (enabled) DOFs that contribute to the QD2T-related linear accelerations of blade K
               LinAccES(:,J,K) = LinAccES(:,J,K) + m%RtHS%PLinVelES(K,J,p%DOFs%PSE(K,I),0,:)*m%QD2T(p%DOFs%PSE(K,I))
  -!SP_start
               AngAccEK(:,J,K) = AngAccEK(:,J,K) + m%RtHS%PAngVelEM(K,J,p%DOFs%PSE(K,I),0,:)*m%QD2T(p%DOFs%PSE(K,I))
  -!SP_end
            ENDDO             ! I - All active (enabled) DOFs that contribute to the QD2T-related linear accelerations of blade K
   
         ENDDO             ! J - Blade nodes / elements
  
positive.code:    REAL(R8Ki)                   :: TmpVec2   (NDims)                               ! A temporary vector.
  
     REAL(ReKi)                   :: LinAccES (NDims,0:p%TipNode,p%NumBl)            ! Total linear acceleration of a point on a   blade (point S) in the inertia frame (body E for earth).
  !SP_start
     REAL(ReKi)                   :: AngAccEK (NDims,0:p%TipNode,p%NumBl)            ! Total linear acceleration of a point on a   blade (point S) in the inertia frame (body E for earth).
  !SP_end
     REAL(ReKi)                   :: LinAccET (NDims,0:p%TwrNodes)                   ! Total linear acceleration of a point on the tower (point T) in the inertia frame (body E for earth).
     REAL(ReKi)                   :: AngAccEF (NDims,0:p%TwrNodes)                   ! Total angular acceleration of tower element J (body F) in the inertia frame (body E for earth).
     REAL(ReKi)                   :: FrcS0B   (NDims,p%NumBl)                        ! Total force at the blade root (point S(0)) due to the blade.
  
negative.code:       ErrMsg = ' Error allocating memory for LinAccESt.'
        RETURN
     ENDIF
     ALLOCATE ( RtHS%AngAccEKt( Dims, 0:p%TipNode, p%NumBl ) , STAT=ErrStat )
     IF ( ErrStat /= 0_IntKi )  THEN
        ErrStat = ErrID_Fatal
        ErrMsg = ' Error allocating memory for AngAccEKt.'
        RETURN
     ENDIF
  
     ALLOCATE(RtHS%AngPosHM(Dims, p%NumBl, 0:p%TipNode), STAT=ErrStat )
     IF ( ErrStat /= 0_IntKi )  THEN
  