anchor.diff:     DO i=1,p%elem_total
          DO j=1,p%node_elem
              temp_int = (j-1)*p%dof_node
  -           CALL BD_CrvMatrixR(p%uuN0(temp_int+4:temp_int+6,i),y%BldMotion%Orientation(:,:,(i-1)*p%node_elem+j))
  +           CALL BD_CrvMatrixR(p%uuN0(temp_int+4:temp_int+6,i),y%BldMotion%Orientation(:,:,(i-1)*p%node_elem+j),&
  +                                     ErrStat2,ErrMsg2)
          ENDDO
      ENDDO
      y%BldMotion%TranslationVel(:,:)  = 0.0D0
  
positive.code:    CHARACTER(ErrMsgLen)        :: ErrMsg2                      ! Temporary Error message
     CHARACTER(*), PARAMETER     :: RoutineName = 'BD_GaussPointDataAt0'
  
     uu0 = 0.0D0
     E10 = 0.0D0
     DO inode=1,node_elem
         hhi = hhx(inode)
         hpi = hpx(inode)
  
negative.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,node_elem,dof_node,elk,elf,elm,elg,  &
                                 ErrStat,ErrMsg)
  
     REAL(ReKi),     INTENT(IN   ):: Nuu0(:)
     REAL(ReKi),     INTENT(IN   ):: Nuuu(:)
  