anchor.diff:         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(1)      = 0.0_BDKi
          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)
  
positive.code:    elk(:,:) = 0.0_BDKi
     elf(:)   = 0.0_BDKi
  
     fact = .TRUE.
  
  
  
negative.code:    REAL(BDKi),         INTENT(IN   ):: uu0(:,:)
     REAL(BDKi),         INTENT(IN   ):: E10(:,:)
     REAL(BDKi),         INTENT(  OUT):: Force(:)
     INTEGER(IntKi),     INTENT(  OUT):: ErrStat      ! Error status of the operation
     CHARACTER(*),       INTENT(  OUT):: ErrMsg       ! Error message if ErrStat /= ErrID_None
  
     INTEGER(IntKi)                   :: ErrStat2     ! Temporary Error status
     CHARACTER(ErrMsgLen)             :: ErrMsg2      ! Temporary Error message
     CHARACTER(*), PARAMETER          :: RoutineName = 'BD_DynamicSolutionForce'
  
     ErrStat = ErrID_None
     ErrMsg  = ""
  
  
     CALL BD_GenerateDynamicElementForce(uuN,vvN,aaN,                       &
                                         Stif0,Mass0,u,                     &
                                         damp_flag,beta,                    &
                                         elem_total,node_elem,dof_node,ngp, &
                                         quadrature,gw,hhx,hpx,Jaco,uu0,E10,&
                                         Force,ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  