anchor.diff:     rrp = 0.0D0
      DO inode=1,node_elem
          hhi = hhx(inode)
  -       hpi = hpx(inode)
  +       hpi = hpx(inode)/Jaco
          temp_id = (inode-1)*dof_node
          temp_id2 = (inode-1)*dof_node/2
          DO i=1,3
  
positive.code:        temp_id2 = (igp - 1) * (dof_node/2)
         Stif(:,:) = 0.0D0
         Stif(1:6,1:6) = EStif0_GL(1:6,1:6,igp)
         CALL BD_GaussPointData(hhx(:,igp),hpx(:,igp),Nuuu,Nrrr,&
               uu0(temp_id1+1:temp_id1+6),E10(temp_id2+1:temp_id2+3),node_elem,dof_node,&
               uuu,uup,E1,RR0,kapa,Stif,cet,ErrStat2,ErrMsg2)
            CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  
negative.code: !--------------------------------------------------------------------------
     REAL(ReKi),    INTENT(IN   ):: hhx(:)      ! Shape function
     REAL(ReKi),    INTENT(IN   ):: hpx(:)      ! Derivative of shape function
     REAL(ReKi),    INTENT(IN   ):: Jaco        ! Jacobian
     REAL(ReKi),    INTENT(IN   ):: Nuuu(:)     ! Element nodal displacement array
     REAL(ReKi),    INTENT(IN   ):: Nrrr(:)     ! Element nodal relative rotation array
     REAL(ReKi),    INTENT(IN   ):: uu0(:)      ! Initial position array at Gauss point
  