anchor.diff:             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_NodalRelRot(Nuu0,node_elem,dof_node,Nrr0,ErrStat2,ErrMsg2)
  +       CALL BD_NodalRelRot(Nuuu,node_elem,dof_node,Nrrr,ErrStat2,ErrMsg2)
   
          CALL BD_ElemNodalDisp(vvN,node_elem,dof_node,nelem,Nvvv,ErrStat2,ErrMsg2)
   
  
positive.code:    fff = 0.0D0
     fff = MATMUL(Stif,eee)
  
     Wrk = 0.0D0
     Wrk = MATMUL(TRANSPOSE(RR0),tempS)
     e1s = Wrk(1)      !epsilon_{11} in material basis
  
     Wrk = 0.0D0
     Wrk = MATMUL(TRANSPOSE(RR0),tempK)
     k1s = Wrk(1)      !kapa_{1} in material basis
  
     DO i=1,3
         fff(i) = fff(i) + 0.5D0*cet*k1s*k1s*RR0(i,1)
         fff(i+3) = fff(i+3) + cet*e1s*k1s*RR0(i,1)
     ENDDO
  
     Fc = 0.0D0
     Fc = fff
     Wrk = 0.0D0
     Wrk(1:3) = fff(1:3)
     Fd = 0.0D0
     Fd(4:6) = MATMUL(TRANSPOSE(BD_Tilde(E1)),Wrk)
  
     C11 = 0.0D0
     C12 = 0.0D0
     C21 = 0.0D0
     C22 = 0.0D0
     C11(1:3,1:3) = Stif(1:3,1:3)
     C12(1:3,1:3) = Stif(1:3,4:6)
     C21(1:3,1:3) = Stif(4:6,1:3)
  
negative.code:    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),&
                                       ErrStat2,ErrMsg2)
         ENDDO
     ENDDO
     y%BldMotion%TranslationVel(:,:)  = 0.0D0
  