anchor.diff:           if (allocated(indx       )) deallocate(indx       )
   
         end subroutine Cleanup
  -   END SUBROUTINE BD_DynamicSolutionGA2
  +END SUBROUTINE BD_DynamicSolutionGA2
   
      SUBROUTINE BD_GenerateDynamicElementGA2(uuN0,uuNf,vvNf,aaNf,                 &
                                              Stif0,Mass0,gravity,u,damp_flag,beta,&
  
positive.code:    Twist_Angle = 0.0D0
     Twist_Angle = Coef(4,1) + Coef(4,2)*eta + Coef(4,3)*eta*eta + Coef(4,4)*eta*eta*eta
  
     END SUBROUTINE BD_ComputeIniNodalPosition
  
     SUBROUTINE BD_CrvExtractCrv(Rr,cc,ErrStat,ErrMsg)
     !--------------------------------------------------
     ! This subroutine computes the CRV parameters given
     ! the components of the rotation matrix
     !--------------------------------------------------
  
     REAL(ReKi),    INTENT(IN   ):: Rr(:,:)       ! Rotation Matrix
     REAL(ReKi),    INTENT(  OUT):: cc(:)         ! Crv paramteres
  
negative.code: 
     ENDDO
  
  END SUBROUTINE BD_TiSchmPredictorStep
  
  SUBROUTINE BD_TiSchmComputeCoefficients(rhoinf,deltat,coef,ErrStat,ErrMsg)
  !----------------------------------------------------------------------
  ! This subroutine calculates the coefficients used in generalized-alpha
  ! time integrator
  !----------------------------------------------------------------------
  
     REAL(DbKi),    INTENT(IN   ):: rhoinf
     REAL(DbKi),    INTENT(IN   ):: deltat
  