anchor.diff:  
   END SUBROUTINE BD_GenerateDynamicElementForce
   !-----------------------------------------------------------------------------------------------------------------------------------
  -SUBROUTINE BD_DynamicSolutionForce(uuN0,rrN0,uuN,vvN,aaN,                                      &
  -                                   Stif0,Mass0,gravity,u,                                 &
  +SUBROUTINE BD_DynamicSolutionForce(uuN,vvN,aaN,                                      &
  +                                   Stif0,Mass0,u,                                 &
                                      damp_flag,beta,                                        &
  -                                   node_elem,dof_node,elem_total,dof_total,node_total,ngp,&
  +                                   node_elem,dof_node,elem_total,ngp,&
                                      quadrature,gw,hhx,hpx,Jaco,uu0,E10,&
                                      Force,ErrStat,ErrMsg)
   !***************************************************************************************
   ! This subroutine calculates the finite-element nodal forces along the beam
   ! Nodal forces = C \dot{u} + K u
   !***************************************************************************************
  -   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   ):: aaN(:) ! Velocity of Mass 1: m/s
      REAL(ReKi),         INTENT(IN   ):: Stif0(:,:,:) ! Element stiffness matrix
      REAL(ReKi),         INTENT(IN   ):: Mass0(:,:,:) ! Element stiffness matrix
  -   REAL(ReKi),         INTENT(IN   ):: gravity(:) !
      INTEGER(IntKi),     INTENT(IN   ):: damp_flag ! Number of Gauss points
      REAL(ReKi),         INTENT(IN   ):: beta(:)
      TYPE(BD_InputType), INTENT(IN   ):: u           ! Inputs at t
      INTEGER(IntKi),     INTENT(IN   ):: node_elem ! Node per element
  -   INTEGER(IntKi),     INTENT(IN   ):: dof_node ! Degrees of freedom per element ! bjj: NOT USED
  +   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
  -   INTEGER(IntKi),     INTENT(IN   ):: node_total ! Total number of nodes  ! bjj: NOT USED
      INTEGER(IntKi),     INTENT(IN   ):: ngp ! Number of Gauss points
      INTEGER(IntKi),     INTENT(IN   ):: quadrature
      REAL(ReKi),         INTENT(IN   ):: gw(:)
  
positive.code:    REAL(ReKi),     INTENT(  OUT):: vvv(:)
     REAL(ReKi),     INTENT(  OUT):: vvp(:)
     REAL(ReKi),     INTENT(  OUT):: aaa(:)
     REAL(ReKi),     INTENT(INOUT):: mmm  ! bjj: NOT USED
     REAL(ReKi),     INTENT(INOUT):: mEta(:)
     REAL(ReKi),     INTENT(INOUT):: rho(:,:)
     INTEGER(IntKi), INTENT(  OUT):: ErrStat       ! Error status of the operation
  
negative.code:    REAL(ReKi),        INTENT(IN   ):: E10(:,:)
     INTEGER(IntKi),    INTENT(IN   ):: elem_total ! Total number of elements
     INTEGER(IntKi),    INTENT(IN   ):: node_elem ! Node per element
     INTEGER(IntKi),    INTENT(IN   ):: dof_node ! Degrees of freedom per node
     INTEGER(IntKi),    INTENT(IN   ):: ngp ! Number of Gauss points
     REAL(ReKi),        INTENT(  OUT):: MassM(:,:) ! Mass matrix
  