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: 
  !** Local variables
        integer(kind=kint) :: nelem, mid, j, isect, nline, tline, icel, iiS
        integer(kind=kint) :: ndof2
        integer(kind=kint) :: ie, ia, jelem, ic_type, nn, jS, jE, itype
        integer(kind=kint) :: nodLOCAL(20),NTOTsum(1)
        integer(kind=8)    :: ntdof2,nonzero
        real(kind=kreal)   :: al, almin, almax, AA, thick, vol, avvol
        real(kind=kreal)   :: tvol, tvmax, tvmin, tlmax, tlmin, asp, aspmax
        real(kind=kreal)   :: xx(20),yy(20),zz(20)
  
negative.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,1pe12.5,a)") '  Sparsity   :',100.0d0*dble(nonzero)/dble(ntdof2),"[%]"
        write(*   ,"(a,1pe12.5,a)") '  Sparsity   :',100.0d0*dble(nonzero)/dble(ntdof2),"[%]"
  
  