anchor.diff:         Qe = 0.0D0
          Stif = 0.0D0
          cet = 0.0D0
  +       
  +       Fi = 0.0D0
  +       Mi = 0.0D0
  +       Gi = 0.0D0
  +       Ki = 0.0D0
  +       mEta = 0.0D0
  +       rho = 0.0D0
          CALL NodalDataAt0(node_elem,nelem,norder,dof_node,i,hhp,Nuu0,Jac,E10)
          CALL NodalData(Nuuu,Nrrr,Nuu0,Nrr0,E10,hhp,Stif0,Jac,&
                         &node_elem,nelem,i,norder,dof_node,&
                         &E1,RR0,kapa,Stif,cet)
          CALL ElasticForce(E1,RR0,kapa,Stif,cet,Fc,Fd,Oe,Pe,Qe)
  +       
  +       CALL NodalDataMass(RR0, mEta0, rho0, mEta, rho)
  +       CALL InertialForce(m00, mEta, rho, Nvvv, Naaa, Fi, Mi, Gi, Ki)
  +       Fc_elem(1:6,i) = Fc(1:6)
  +       Fd_elem(1:6,i) = Fd(1:6)
  +       Oe_elem(1:6,1:6,i) = Oe(1:6,1:6)
  +       Pe_elem(1:6,1:6,i) = Pe(1:6,1:6)
  +       Qe_elem(1:6,1:6,i) = Qe(1:6,1:6)
  +       Se_elem(1:6,1:6,i) = Stif(1:6,1:6)
  +       
          Fi_elem(1:6,i) = Fi(1:6)
          Ki_elem(1:6,1:6,i) = Ki(1:6,1:6)
          Gi_elem(1:6,1:6,i) = Gi(1:6,1:6)
  
positive.code:             IF(ALLOCATED(Pe_elem)) DEALLOCATE(Pe_elem)
              IF(ALLOCATED(Qe_elem)) DEALLOCATE(Qe_elem)
              IF(ALLOCATED(Se_elem)) DEALLOCATE(Se_elem)
          ENDIF
  
     END SUBROUTINE ElementMatrix
  
negative.code:            ENDDO
         ENDDO
     ENDDO
     
     DO i=1,node_elem
         DO j=1,6
             temp_id1 = (i-1)*dof_node+j
             DO k=1,6
                 temp_id2 = (i-1)*dof_node+k
                 elm(temp_id1,temp_id2) = elm(temp_id1, temp_id2) + w(i) * Mi_elem(j,k,i) * Jac
                 elg(temp_id1,temp_id2) = elg(temp_id1, temp_id2) + w(i) * Gi_elem(j,k,i) * Jac
                 elk(temp_id1,temp_id2) = elk(temp_id1, temp_id2) + w(i) * Ki_elem(j,k,i) * Jac
             ENDDO
         ENDDO
     ENDDO
  
     DO i=1,node_elem
         DO j=1,6
             temp_id1 = (i-1)*dof_node+j
             elf(temp_id1) = elf(temp_id1) - w(i) * Fi_elem(j,i) * Jac
         ENDDO
     ENDDO
  
  
  
  