anchor.diff:  
   END SUBROUTINE BD_GenerateStaticElement
   !-----------------------------------------------------------------------------------------------------------------------------------
  -SUBROUTINE BD_StaticElementMatrix(Nuu0,Nuuu,Nrr0,Nrrr,Distr_GL,gravity,&
  +SUBROUTINE BD_StaticElementMatrix(Nuuu,Nrrr,Distr_GL,gravity,&
                                     EMass0_GL,EStif0_GL,&
                                     ngp,gw,hhx,hpx,Jaco,uu0,E10,&
                                     node_elem,dof_node,elk,elf,&
                                     ErrStat,ErrMsg)
   
  -   REAL(ReKi),    INTENT(IN   ):: Nuu0(:) !bjj: NOT USED
      REAL(ReKi),    INTENT(IN   ):: Nuuu(:)
  -   REAL(ReKi),    INTENT(IN   ):: Nrr0(:) !bjj: NOT USED
      REAL(ReKi),    INTENT(IN   ):: Nrrr(:)
      REAL(ReKi),    INTENT(IN   ):: Distr_GL(:,:)
      REAL(ReKi),    INTENT(IN   ):: gravity(:)
  
positive.code: 
  END SUBROUTINE BD_NodalRelRot
  !-----------------------------------------------------------------------------------------------------------------------------------
  SUBROUTINE BD_ElementMatrixGA2(Nuu0,Nuuu,Nrr0,Nrrr,Nvvv,Naaa,           &
                                 EStif0_GL,EMass0_GL,gravity,DistrLoad_GL,&
                                 damp_flag,beta,                          &
                                 ngp,gw,hhx,hpx,Jaco,uu0,E10,             &
                                 node_elem,dof_node,fact,elk,elf,elm,elg, &
                                 ErrStat,ErrMsg)
  
     REAL(ReKi),     INTENT(IN   ):: Nuu0(:) !bjj: NOT USED
     REAL(ReKi),     INTENT(IN   ):: Nuuu(:)
     REAL(ReKi),     INTENT(IN   ):: Nrr0(:) !bjj: NOT USED
     REAL(ReKi),     INTENT(IN   ):: Nrrr(:)
     REAL(ReKi),     INTENT(IN   ):: Nvvv(:)
     REAL(ReKi),     INTENT(IN   ):: Naaa(:)
  
negative.code:        mEta(3)      = EMass0_GL(1,5,igp)
         rho(1:3,1:3) = EMass0_GL(4:6,4:6,igp)
         CALL BD_GaussPointDataMass(hhx(:,igp),hpx(:,igp),Jaco(igp),temp_Nvvv,temp_Naaa,RR0,&
               node_elem,dof_node,vvv,aaa,vvp,mEta,rho,ErrStat2,ErrMsg2)
            CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
         CALL BD_GravityForce(mmm,mEta,gravity,Fg,ErrStat2,ErrMsg2)
            CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  