anchor.diff:  
      y%BldMotion%TranslationDisp(:,:) = 0.0D0
      y%BldMotion%Orientation(:,:,:)   = 0.0D0
  -   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
      y%BldMotion%RotationVel(:,:)     = 0.0D0
      y%BldMotion%TranslationAcc(:,:)  = 0.0D0
  
positive.code:            temp_id2= (i-1)*p%node_elem+j
             IF(i .EQ. 1 .AND. j .EQ. 1) THEN
                 temp6(:) = 0.0D0
                 temp6(1:3) = u%RootMotion%TranslationAcc(1:3,1)
                 temp6(4:6) = u%RootMotion%RotationAcc(1:3,1)
             ELSE
                 temp6(:) = 0.0D0
                 temp6(1:3) = y%BldMotion%TranslationAcc(1:3,temp_id2)
                 temp6(4:6) = y%BldMotion%RotationAcc(1:3,temp_id2)
                 temp6(1:3) = MATMUL(TRANSPOSE(p%GlbRot),temp6(1:3))
                 temp6(4:6) = MATMUL(TRANSPOSE(p%GlbRot),temp6(4:6))
             ENDIF
             OtherState%Acc(temp_id+1:temp_id+6) = temp6(1:6)
         ENDDO
     ENDDO
  
negative.code:        u%DistrLoad%Moment(:,k) = 0.0D0
     ENDDO
  
     CALL BD_CopyInput(u, u_tmp, MESH_NEWCOPY, ErrStat2, ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
     CALL BD_InputGlobalLocal(p,u_tmp,ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
     CALL BD_CalcIC(u_tmp,p,x,OtherState,ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
     ! Define initial guess for the system outputs here:
  
  