anchor.diff:      rnum = (7.21d0+0.01*dlog10(dble(hecMAT%N)))*10.0d0/dble(hecMAT%N)
       cond = 1.0d0
       !rnum = (7.25d0)*10.0d0/dble(hecMAT%N)
  -    
  +
       open(fio,file='nonzero.dat',status='replace')
  -      write(fio,"(a,f12.5,i)")"##magic number 10 : 7.2, ",rnum,hecMAT%N
  +      write(fio,"(a,f12.5,i0)")"##magic number 10 : 7.2, ",rnum,hecMAT%N
         do i= 1, n
           jS= hecMAT%indexL(i-1) + 1
           jE= hecMAT%indexL(i  )
  -        write(fio,"(i,a,i)")i,"    ",i
  +        write(fio,"(i0,a,i0)")i,"    ",i
           do j= jS, jE
             in = hecMAT%itemL(j)
  -          write(fio,"(i,a,i)")i,"    ",in
  -          write(fio,"(i,a,i)")in,"    ",i
  +          write(fio,"(i0,a,i0)")i,"    ",in
  +          write(fio,"(i0,a,i0)")in,"    ",i
           enddo
         enddo
       close(fio)
  -    
  +
       open(fio,file='nonzero.plt',status='replace')
         if(ftype == 4)then
           write(fio,"(a)")'set terminal png'
  
positive.code:     implicit none
      type (hecmwST_local_mesh) :: hecMESH
      type (hecmwST_matrix)     :: hecMAT
      
      integer(kind=kint) :: i, j, in, jS, jE, ftype, n, ndof, nnz, fio
      real(kind=kreal) :: rnum, dens, cond
  
  
negative.code:     !ftype = 2: eps
      !ftype = 4: png
      ftype = 2
  
      n = hecMAT%N
      ndof = 3*n
      nnz = 9*n + 9*2*hecMAT%indexL(n)
  