anchor.diff:  !                                                                      !
   ! Software Name : FrontISTR Ver. 3.6                                   !
   !                                                                      !
  -!      Module Name : Library                                            !
  +!      Module Name : Library                                           !
   !                                                                      !
   !            Written by Yasuji Fukahori (Univ. of Tokyo)               !
   !                                                                      !
  
positive.code: !C Mesh Summary
  !C
        write(ILOG,"(a)") '###  Mesh Summary  ###'
        write(ILOG,"(a,i12)") 'Num of node:',hecMESH%n_node
        write(ILOG,"(a,i12)") 'Num of DOF :',hecMESH%n_node*hecMESH%n_dof
        ndof2  = hecMESH%n_dof**2
        ntdof2 = (hecMESH%n_node*hecMESH%n_dof)**2
        write(ILOG,"(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  )
          ic_type = hecMESH%elem_type_item(itype)
          write(ILOG,"(a,i4,a,i12)") 'Num of ',ic_type,':',jE-jS
        enddo
        nonzero = ndof2*(hecMAT%NP + hecMAT%NPU + hecMAT%NPL)
        write(ILOG,"(a,i12)") 'Num of NZ  :',nonzero
        write(ILOG,"(a,i12)") 'Num of DOF2:',ntdof2
        write(ILOG,"(a,1pe12.5,a)") 'Sparsity   :',100.0d0*dble(nonzero)/dble(ntdof2),"[%]"
  
  !C
  !C Output sparsity pattern
  
negative.code:       write(fio,"(a,f12.5,a)")'plot "nonzero.dat" pointtype 5 pointsize ',rnum,' linecolor rgb "#F96566"'
      close(fio)
      
      write(*,*)''
      write(*,*)'  gnuplot -persist "nonzero.plt"'
      write(*,*)''
      !call system('gnuplot -persist "nonzero.plt"')
      
      !open(fio,file='nonzero.dat',status='old')
  