anchor.diff:  
   END SUBROUTINE BD_ComputeBladeMassNew
   
  -SUBROUTINE BD_ComputeElementMass(Nuu0,NGPpos,EMass0_GL,&
  -                                 ngp,gw,hhx,hpx,Jaco,&
  -                                 node_elem,dof_node,&
  +SUBROUTINE BD_ComputeElementMass(NGPpos,EMass0_GL,&
  +                                 ngp,gw,Jaco,&
                                    elem_mass,elem_CG,elem_IN,ErrStat,ErrMsg)
   
   !-------------------------------------------------------------------------------
   ! This subroutine total element forces and mass matrices
   !-------------------------------------------------------------------------------
   
  -   REAL(ReKi),INTENT(IN   )    :: Nuu0(:) ! Nodal initial position for each element !bjj: NOT USED
      REAL(ReKi),INTENT(IN   )    :: NGPpos(:,:)
      REAL(ReKi),INTENT(IN   )    :: EMass0_GL(:,:,:) ! Nodal material properties for each element
      REAL(ReKi),INTENT(  OUT)    :: elem_mass  ! Total element force (Fd, Fc, Fb)
  
positive.code: ! nodes along beam axis
  ! Nodal forces = K u
  !***************************************************************************************
     REAL(ReKi),        INTENT(IN   ):: uuN0(:,:) ! Initial position vector
     REAL(ReKi),        INTENT(IN   ):: rrN0(:,:) ! Initial position vector
     REAL(ReKi),        INTENT(IN   ):: Stif0(:,:,:) ! Element stiffness matrix
     REAL(ReKi),        INTENT(IN   ):: Mass0(:,:,:) ! Element stiffness matrix
     REAL(ReKi),        INTENT(IN   ):: gravity(:) !
     TYPE(BD_InputType),INTENT(IN   ):: u           ! Inputs at t
     REAL(ReKi),        INTENT(IN   ):: uuN(:) ! Displacement of Mass 1: m
     REAL(ReKi),        INTENT(IN   ):: vvN(:) ! Displacement of Mass 1: m
     INTEGER(IntKi),    INTENT(IN   ):: node_elem ! Node per element
     INTEGER(IntKi),    INTENT(IN   ):: dof_node ! Degrees of freedom per element
     INTEGER(IntKi),    INTENT(IN   ):: elem_total ! Total number of elements
     INTEGER(IntKi),    INTENT(IN   ):: dof_total ! Total number of degrees of freedom ! bjj: NOT USED
     INTEGER(IntKi),    INTENT(IN   ):: node_total ! Total number of nodes
     INTEGER(IntKi),    INTENT(IN   ):: ngp
     INTEGER(IntKi),    INTENT(IN   ):: quadrature
     REAL(ReKi),        INTENT(IN   ):: gw(:)
  
negative.code:        return
     end if
  
     CALL BD_GenerateDynamicElementAcc(uuN,vvN,Stif0,Mass0,gravity,u,&
                                       damp_flag,beta,&
                                       quadrature,gw,hhx,hpx,Jacobian,uu0,E10,                &
                                       elem_total,node_elem,dof_node,ngp,&
                                       RHS,MassM,ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
     DO j=1,node_total
  