anchor.diff:  
   END SUBROUTINE BD_GenerateStaticElementForce
   !-----------------------------------------------------------------------------------------------------------------------------------
  -SUBROUTINE BD_StaticElementMatrixForce(Nuu0,Nuuu,Nrr0,Nrrr,Nvvv,EStif0_GL,EMass0_GL,&
  -              gravity,DistrLoad_GL,ngp,gw,hhx,hpx,Jaco,uu0,E10,                     &
  +SUBROUTINE BD_StaticElementMatrixForce(Nuuu,Nrrr,EStif0_GL,&
  +              ngp,gw,hhx,hpx,Jaco,uu0,E10,                     &
                 node_elem,dof_node,elf,ErrStat,ErrMsg)
   !-------------------------------------------------------------------------------
   ! This subroutine calculates elemental internal node force for static analysis
   !-------------------------------------------------------------------------------
  -   REAL(ReKi),    INTENT(IN   ):: Nuu0(:) ! Nodal initial position for each element !bjj: NOT USED
      REAL(ReKi),    INTENT(IN   ):: Nuuu(:) ! Nodal displacement of Mass 1 for each element
  -   REAL(ReKi),    INTENT(IN   ):: Nrr0(:) ! Nodal rotation parameters for initial position !bjj: NOT USED
      REAL(ReKi),    INTENT(IN   ):: Nrrr(:) ! Nodal rotation parameters for displacement of Mass 1
  -   REAL(ReKi),    INTENT(IN   ):: Nvvv(:) ! Nodal velocity of Mass 1: m/s for each element ! bjj: NOT USED
      REAL(ReKi),    INTENT(IN   ):: EStif0_GL(:,:,:) ! Nodal material properties for each element
  -   REAL(ReKi),    INTENT(IN   ):: EMass0_GL(:,:,:) ! Nodal material properties for each element ! bjj: NOT USED
  -   REAL(ReKi),    INTENT(IN   ):: gravity(:) ! bjj: NOT USED
  -   REAL(ReKi),    INTENT(IN   ):: DistrLoad_GL(:,:) ! Nodal material properties for each element ! bjj: NOT USED
      INTEGER(IntKi),INTENT(IN   ):: ngp ! Number of Gauss points
      REAL(ReKi),    INTENT(IN   ):: gw(:)
      REAL(ReKi),    INTENT(IN   ):: hhx(:,:)
  
positive.code:            DistrLoad_GL(1:3,j) = u%DistrLoad%Force(1:3,temp_id+j)
             DistrLoad_GL(4:6,j) = u%DistrLoad%Moment(1:3,temp_id+j)
         ENDDO
         CALL BD_StaticElementMatrixForce(uuN0(:,nelem),Nuuu,rrN0(:,nelem),Nrrr,Nvvv,&
                 EStif0_GL,EMass0_GL,gravity,DistrLoad_GL,&
                 ngp,gw,hhx,hpx,Jaco(:,nelem),uu0(:,nelem),E10(:,nelem),&
                 node_elem,dof_node,elf,ErrStat2,ErrMsg2)
            CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  
  
negative.code: ! Given states (u,v) and accelerations (acc,xcc) at the initial of a time step (t_i),
  ! it returns the values of states and accelerations at the end of a time step (t_f)
  !------------------------------------------------------------------------------------
     REAL(ReKi),         INTENT(IN   ):: Stif0(:,:,:)
     REAL(ReKi),         INTENT(IN   ):: Mass0(:,:,:)
     REAL(ReKi),         INTENT(IN   ):: gravity(:)
  