anchor.diff:        do l = 1,np+1
            
          if(abs(spts(j)-npts(l)).LE.eps) then
  -         Ps(l,j) = 1.
  +         Ps(l,j) = 1.0_BDKi
          else
  -         dnum = 1.
  -         den = 1.
  +         dnum = 1.0_BDKi
  +         den  = 1.0_BDKi
            do k = 1,np+1
              if (k.NE.l) then
  -             den = den*(npts(l) - npts(k))
  +             den  = den *(npts(l) - npts(k))
                dnum = dnum*(spts(j) - npts(k))
              endif
            enddo
  
positive.code:              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 )
         mmm  = 0.0_BDKi
         mEta = 0.0_BDKi
         rho  = 0.0_BDKi
         mmm          = EMass0_GL(1,1,igp)
         mEta(2)      = -EMass0_GL(1,6,igp)
         mEta(3)      =  EMass0_GL(1,5,igp)
         rho(1:3,1:3) = EMass0_GL(4:6,4:6,igp)
  
negative.code: ! This subroutine calculates elemetal internal forces
  !------------------------------------------------------------
  
     REAL(BDKi),     INTENT(IN   ):: Nuuu(:)          ! Nodal displacement of Mass 1 for each element
     REAL(BDKi),     INTENT(IN   ):: Nrrr(:)          ! Nodal rotation parameters for displacement of Mass 1
     REAL(BDKi),     INTENT(IN   ):: Nvvv(:)          ! Nodal velocity of Mass 1: m/s for each element
     REAL(BDKi),     INTENT(IN   ):: EStif0_GL(:,:,:) ! Nodal material properties for each element
     REAL(BDKi),     INTENT(IN   ):: EMass0_GL(:,:,:) ! Nodal material properties for each element
     INTEGER(IntKi), INTENT(IN   ):: damp_flag        ! Number of Gauss points
     REAL(BDKi),     INTENT(IN   ):: beta(:)
     INTEGER(IntKi), INTENT(IN   ):: ngp              ! Number of Gauss points
     REAL(BDKi),     INTENT(IN   ):: gw(:)
     REAL(BDKi),     INTENT(IN   ):: hhx(:,:)
     REAL(BDKi),     INTENT(IN   ):: hpx(:,:)
     REAL(BDKi),     INTENT(IN   ):: Jaco(:)
     REAL(BDKi),     INTENT(IN   ):: uu0(:)
     REAL(BDKi),     INTENT(IN   ):: E10(:)
     INTEGER(IntKi), INTENT(IN   ):: node_elem        ! Node per element
     INTEGER(IntKi), INTENT(IN   ):: dof_node         ! Degrees of freedom per node
     REAL(BDKi),     INTENT(  OUT):: elf(:)           ! Total element force (Fd, Fc, Fb)
     INTEGER(IntKi), INTENT(  OUT):: ErrStat          ! Error status of the operation
     CHARACTER(*),   INTENT(  OUT):: ErrMsg           ! Error message if ErrStat /= ErrID_None
  
     REAL(BDKi)                   :: temp_Naaa(dof_node*node_elem)
     REAL(BDKi)                   :: RR0(3,3)
  