anchor.diff:         RHS(temp_id1+3) = kp_coord(kp_member,i)
          RHS(temp_id1+4) = 0.0D0
          CALL LAPACK_getrf( 4*(kp_member-1), 4*(kp_member-1), K,indx,&
  -                          ErrStat2, ErrMsg) 
  +                          ErrStat2, ErrMsg2) 
             CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
          CALL LAPACK_getrs( 'N',4*(kp_member-1), K,indx,RHS,&
  -                          ErrStat2, ErrMsg) 
  +                          ErrStat2, ErrMsg2) 
             CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
          DO j=1,kp_member-1
              DO m=1,4
  
positive.code:        ENDDO
         RHS(:) = RHS(:) + F_PointLoad(:)
  
  
         errf = 0.0D0
         feqv = 0.0D0
         DO j=1,dof_total-6
             feqv(j) = RHS(j+6)
             RHS_LU(j) = RHS(j+6)
             DO k=1,dof_total-6
                 StifK_LU(j,k) = StifK(j+6,k+6)
             ENDDO
         ENDDO
  
         CALL ludcmp(StifK_LU,dof_total-6,indx,d)
         CALL lubksb(StifK_LU,dof_total-6,indx,RHS_LU,ai_temp)
  
         ai = 0.0D0
         DO j=1,dof_total-6
             ai(j+6) = ai_temp(j)
         ENDDO
         CALL BD_UpdateDynamicGA2(ai,uuNf,vvNf,aaNf,xxNf,coef,node_total,dof_node)
         IF(i==1) THEN
             Eref = SQRT(abs(DOT_PRODUCT(ai_temp,feqv)))*TOLF
             IF(Eref .LE. TOLF) RETURN
         ENDIF
         IF(i .GT. 1) THEN
             Enorm = SQRT(abs(DOT_PRODUCT(ai_temp,feqv)))
             IF(Enorm .LE. Eref) RETURN
         ENDIF
  
         IF(i==niter) THEN
            CALL setErrStat( ErrID_Fatal, "Solution does not converge after the maximum number of iterations", ErrStat, ErrMsg, RoutineName)
            RETURN
         ENDIF
     ENDDO
  
     END SUBROUTINE BD_DynamicSolutionGA2
  
     SUBROUTINE BD_GenerateDynamicElementGA2(uuN0,uuNf,vvNf,aaNf,                 &
  
negative.code:    CHARACTER(ErrMsgLen)             :: ErrMsg2     ! Temporary Error message
     CHARACTER(*), PARAMETER          :: RoutineName = 'BD_DynamicSolutionGA2'         
     REAL(ReKi)                       :: errf
     REAL(ReKi),           ALLOCATABLE:: StifK(:,:)
     REAL(ReKi),           ALLOCATABLE:: StifK_LU(:,:)
     REAL(ReKi),           ALLOCATABLE:: RHS(:)
     REAL(ReKi),           ALLOCATABLE:: RHS_LU(:)
     REAL(ReKi),           ALLOCATABLE:: MassM(:,:)
     REAL(ReKi),           ALLOCATABLE:: DampG(:,:)
     REAL(ReKi),           ALLOCATABLE:: F_PointLoad(:)
     REAL(ReKi),           ALLOCATABLE:: ai(:)
     REAL(ReKi)                       :: Eref
     REAL(ReKi)                       :: Enorm
     REAL(ReKi)                       :: d
     INTEGER(IntKi),       ALLOCATABLE:: indx(:)
     INTEGER(IntKi)                   :: i
     INTEGER(IntKi)                   :: j
     INTEGER(IntKi)                   :: k
  