anchor.diff:      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
   
  
positive.code:     !ftype = 2: eps
      !ftype = 4: png
      ftype = 2
      
      n = hecMAT%N
      ndof = 3*n
      nnz = 9*n + 9*2*hecMAT%indexL(n)
  
negative.code:         write(fio,"(a)")'set out "image.png"'
        else
          write(fio,"(a)")'set out "image.eps"'
  
          write(fio,"(a)"     )'set label 1 "Name" at graph 1.1,0.9'
          write(fio,"(a)")'set label 2 "N" at graph 1.1,0.85'
          write(fio,"(a)")'set label 3 "Non-Zero Elem." at graph 1.1,0.8'
          write(fio,"(a)")'set label 4 "Density [%]" at graph 1.1,0.75'
          write(fio,"(a)")'set label 9 "Condition Num." at graph 1.1,0.7'
  
          write(fio,"(a)"     )'set label 5 ":  matrix" at graph 1.4,0.9'
          write(fio,"(a,i0,a)")'set label 6 ":  ',ndof,'" at graph 1.4,0.85'
          write(fio,"(a,i0,a)")'set label 7 ":  ',nnz,'" at graph 1.4,0.8'
          write(fio,"(a,1pe9.2,a)")'set label 8 ": ',dens,'" at graph 1.4,0.75'
          write(fio,"(a,1pe9.2,a)")'set label 10 ": ',cond,'" at graph 1.4,0.7'
        endif
  
        write(fio,"(a,f12.5,a)")'plot "nonzero.dat" pointtype 5 pointsize ',rnum,' linecolor rgb "#F96566"'
      close(fio)
  
      write(*,*)''
      write(*,*)' ### Command recommendation'
      write(*,*)' gnuplot -persist "nonzero.plt"'
  
      !call system('gnuplot -persist "nonzero.plt"')
  
      !open(fio,file='nonzero.dat',status='old')
      !close(fio,status="delete")
      !open(fio,file='nonzero.plt',status='old')
  