anchor.diff:        nonzero = ndof2*(hecMAT%NP + hecMAT%NPU + hecMAT%NPL)
         write(ILOG,"(a,i12)") '  Num of NZ  :',nonzero
         write(*   ,"(a,i12)") '  Num of NZ  :',nonzero
  -      write(ILOG,"(a,i12)") '  Num of DOF2:',ntdof2
  -      write(*   ,"(a,i12)") '  Num of DOF2:',ntdof2
         write(ILOG,"(a,1pe12.5,a)") '  Sparsity   :',100.0d0*dble(nonzero)/dble(ntdof2),"[%]"
         write(*   ,"(a,1pe12.5,a)") '  Sparsity   :',100.0d0*dble(nonzero)/dble(ntdof2),"[%]"
   
  
positive.code:       ntdof2 = (hecMESH%n_node*hecMESH%n_dof)**2
        write(ILOG,"(a,i12)") '  Num of elem:',hecMESH%n_elem
        write(*   ,"(a,i12)") '  Num of elem:',hecMESH%n_elem
        do itype = 1, hecMESH%n_elem_type
          jS = hecMESH%elem_type_index(itype-1)
          jE = hecMESH%elem_type_index(itype  )
  
negative.code:     integer(kind=kint) :: i, j, in, jS, jE, ftype, n, ndof, nnz, fio
      real(kind=kreal) :: rnum, dens, cond
      character :: fileid*3
  
      fio = 70 + hecMESH%my_rank
      write(fileid,"(i3.3)")hecMESH%my_rank
  
  