anchor.diff:         call Cleanup()
          return
      end if
  -   NGPpos(:,:)  = 0.0D0
  -   EMass0_GL(:,:,:)  = 0.0D0
  -   elem_mass= 0.0D0
  -   elem_CG(:)= 0.0D0
  -   elem_IN(:,:)= 0.0D0
  +   NGPpos(:,:)  = 0.0_BDKi
  +   EMass0_GL(:,:,:)  = 0.0_BDKi
  +   elem_mass= 0.0_BDKi
  +   elem_CG(:)= 0.0_BDKi
  +   elem_IN(:,:)= 0.0_BDKi
   
      DO nelem=1,elem_total
   
  
positive.code:    eed(1:3) = eed(1:3) + MATMUL(BD_Tilde(E1),ome)
  
     ! Compute damping matrix
     temp_b(:,:) = 0.0D0
     DO i=1,6
         temp_b(i,i) = beta(i)
     ENDDO
  
negative.code: ! It also calcuates the linearized matrices Mi, Gi, and Ki for N-R algorithm
  !---------------------------------------------------------------------------
  
     REAL(BDKi),    INTENT(IN   ):: m00
     REAL(BDKi),    INTENT(IN   ):: mEta(:)
     REAL(BDKi),    INTENT(IN   ):: rho(:,:)
     REAL(BDKi),    INTENT(IN   ):: vvv(:)
     REAL(BDKi),    INTENT(IN   ):: aaa(:)
     LOGICAL,       INTENT(IN   ):: fact
     REAL(BDKi),    INTENT(  OUT):: Fi(6)
     REAL(BDKi),    INTENT(  OUT):: Mi(6,6)
     REAL(BDKi),    INTENT(  OUT):: Gi(6,6)
     REAL(BDKi),    INTENT(  OUT):: Ki(6,6)
     INTEGER(IntKi),INTENT(  OUT):: ErrStat       ! Error status of the operation
     CHARACTER(*),  INTENT(  OUT):: ErrMsg        ! Error message if ErrStat /= ErrID_None
  
     REAL(BDKi)                  :: beta(3)
     REAL(BDKi)                  :: gama(3)
     REAL(BDKi)                  :: nu(3)
     REAL(BDKi)                  :: epsi(3,3)
     REAL(BDKi)                  :: mu(3,3)
     REAL(BDKi)                  :: ome(3)
     REAL(BDKi)                  :: omd(3)
     REAL(BDKi)                  :: tempV(3)
     REAL(BDKi)                  :: tempA(3)
     INTEGER(IntKi)              :: i
  
     ErrStat = ErrID_None
     ErrMsg  = ""
     Fi(:)   = 0.0_BDKi
     Mi(:,:) = 0.0_BDKi
     Gi(:,:) = 0.0_BDKi
     Ki(:,:) = 0.0_BDKi
  
     ome(:) = vvv(4:6)
     omd(:) = aaa(4:6)
  