anchor.diff:  
   END SUBROUTINE BD_GyroForce
   !-----------------------------------------------------------------------------------------------------------------------------------
  -SUBROUTINE BD_ElementMatrixForce(Nuu0,Nuuu,Nrr0,Nrrr,Nvvv,&
  +SUBROUTINE BD_ElementMatrixForce(Nuuu,Nrrr,Nvvv,&
                                    EStif0_GL,EMass0_GL,     &
                                    damp_flag,beta,          &
                                    ngp,gw,hhx,hpx,Jaco,uu0,E10,&
  
positive.code:        mEta(3)      =  EMass0_GL(1,5,igp)
         rho(1:3,1:3) = EMass0_GL(4:6,4:6,igp)
         CALL BD_GaussPointDataMass(hhx(:,igp),hpx(:,igp),Jaco(igp),Nvvv,temp_Naaa,RR0,&
               node_elem,dof_node, vvv,temp_aaa,vvp,mmm,mEta,rho,ErrStat2,ErrMsg2)
            CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
         CALL BD_ElasticForce(E1,RR0,kapa,Stif,cet,fact,Fc,Fd,Oe,Pe,Qe,ErrStat2,ErrMsg2)
            CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  
negative.code:    REAL(ReKi),     INTENT(  OUT):: vvv(:)
     REAL(ReKi),     INTENT(  OUT):: vvp(:)
     REAL(ReKi),     INTENT(  OUT):: aaa(:)
     REAL(ReKi),     INTENT(INOUT):: mEta(:)
     REAL(ReKi),     INTENT(INOUT):: rho(:,:)
     INTEGER(IntKi), INTENT(  OUT):: ErrStat       ! Error status of the operation
  