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: 
  END SUBROUTINE BD_UpdateDynamicGA2
  !-----------------------------------------------------------------------------------------------------------------------------------
  SUBROUTINE BD_GenerateDynamicElementAcc(uuN0,rrN0,uuN,vvN,Stif0,Mass0,gravity,u,    &
                                          damp_flag,beta,                             &
                                          quadrature,gw,hhx,hpx,Jacobian,uu0,E10,     &
                                          elem_total,node_elem,dof_total,dof_node,ngp,&
                                          RHS,MassM,ErrStat,ErrMsg)
  !----------------------------------------------------------------------------------------
  ! This subroutine computes Global mass matrix and force vector for the beam.
  !----------------------------------------------------------------------------------------
     REAL(ReKi),        INTENT(IN   ):: uuN0(:,:) ! Initial position vector
     REAL(ReKi),        INTENT(IN   ):: rrN0(:,:) ! Initial position vector
     REAL(ReKi),        INTENT(IN   ):: uuN(:) ! Displacement of Mass 1: m
     REAL(ReKi),        INTENT(IN   ):: vvN(:) ! Velocity of Mass 1: m/s
     REAL(ReKi),        INTENT(IN   ):: Stif0(:,:,:) ! Element stiffness matrix
  
negative.code: 
  END SUBROUTINE BD_BoundaryGA2
  !-----------------------------------------------------------------------------------------------------------------------------------
  SUBROUTINE BD_DynamicSolutionGA2( uuNf,vvNf,aaNf,xxNf,          &
                                    Stif0,Mass0,gravity,u,damp_flag,beta,   &
                                    node_elem,dof_node,elem_total,dof_total,&
                                    quadrature,gw,hhx,hpx,Jaco,uu0,E10,            &
  