anchor.diff:      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
   
  
positive.code:       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),"[%]"
  
  
negative.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
        nelem_wo_mpc = 0
        do itype = 1, hecMESH%n_elem_type
          jS = hecMESH%elem_type_index(itype-1)
          jE = hecMESH%elem_type_index(itype  )
          ic_type = hecMESH%elem_type_item(itype)
          if(.not. hecmw_is_etype_link(ic_type)) nelem_wo_mpc = nelem_wo_mpc + jE-jS
        enddo
        write(ILOG,"(a,i12)") '   ** w/o MPC:',nelem_wo_mpc
        write(*   ,"(a,i12)") '   ** w/o MPC:',nelem_wo_mpc
        do itype = 1, hecMESH%n_elem_type
          jS = hecMESH%elem_type_index(itype-1)
          jE = hecMESH%elem_type_index(itype  )
  