anchor.diff:  
   END SUBROUTINE BD_GenerateDynamicElementForce
   !-----------------------------------------------------------------------------------------------------------------------------------
  -SUBROUTINE BD_DynamicSolutionForce(uuN,vvN,aaN,                                      &
  +SUBROUTINE BD_DynamicSolutionForce(uuN,vvN,aaN,                                   &
                                      Stif0,Mass0,u,                                 &
  -                                   damp_flag,beta,                                        &
  +                                   damp_flag,beta,                                &
                                      node_elem,dof_node,elem_total,ngp,&
                                      quadrature,gw,hhx,hpx,Jaco,uu0,E10,&
                                      Force,ErrStat,ErrMsg)
  
positive.code:        CALL BD_ElasticForce(E1,RR0,kapa,Stif,cet,fact,Fc,Fd,Oe,Pe,Qe,ErrStat2,ErrMsg2)
            CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
         mmm          = EMass0_GL(1,1,igp)
         mEta(2)      =-EMass0_GL(1,6,igp)
         mEta(3)      = EMass0_GL(1,5,igp)
         rho(1:3,1:3) = EMass0_GL(4:6,4:6,igp)
  
negative.code:    CHARACTER(ErrMsgLen)         :: ErrMsg2                      ! Temporary Error message
     CHARACTER(*), PARAMETER      :: RoutineName = 'BD_ElememntMatrixForce'
  
     ErrStat   = ErrID_None
     ErrMsg    = ""
     elf       = 0.0_BDKi
     temp_Naaa = 0.0_BDKi
     
     fact = .FALSE.
  
  