anchor.diff:  END SUBROUTINE BD_GaussPointDataAt0
   !-----------------------------------------------------------------------------------------------------------------------------------
      SUBROUTINE BD_GaussPointData(hhx,hpx,Nuuu,Nrrr,uu0,E10,node_elem,dof_node,&
  -                                uuu,uup,E1,RR0,kapa,Stif,cet)
  +                                uuu,uup,E1,RR0,kapa,Stif,cet,ErrStat,ErrMsg)
      !--------------------------------------------------------------------------
      ! This subroutine computes Gauss point values: 1) uuu, 2) uup, 3) E1
      ! 4) RR0, 5) kapa, 6) Stif, and 7) cet
  
positive.code:            DistrLoad_GL(4:6,j)  = u%DistrLoad%Moment(1:3,temp_id+j+1)
         ENDDO
  
         CALL BD_NodalRelRot(Nuu0,node_elem,dof_node,Nrr0)
         CALL BD_NodalRelRot(Nuuu,node_elem,dof_node,Nrrr)
  
         CALL BD_ElemNodalDisp(vvN,node_elem,dof_node,nelem,Nvvv,ErrStat2,ErrMsg2)
  
  
negative.code:    REAL(ReKi),    INTENT(  OUT):: kapa(:)     ! Curvature starin vector at Gauss point
     REAL(ReKi),    INTENT(  OUT):: cet         ! Extension-torsion coefficient at Gauss point
     REAL(ReKi),    INTENT(  OUT):: Stif(:,:)   ! C/S stiffness matrix resolved in inertial frame at Gauss point
     INTEGER(IntKi),INTENT(  OUT):: ErrStat       ! Error status of the operation
     CHARACTER(*),  INTENT(  OUT):: ErrMsg        ! Error message if ErrStat /= ErrID_None
  
     REAL(ReKi)                  :: rrr(3)
     REAL(ReKi)                  :: rrp(3)
  