anchor.diff:  !** 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, iS, iE, itype
  +      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
  
positive.code: !                                                                      !
  ! Software Name : FrontISTR Ver. 3.6                                   !
  !                                                                      !
  !      Module Name : Library                                            !
  !                                                                      !
  !            Written by Yasuji Fukahori (Univ. of Tokyo)               !
  !                                                                      !
  
negative.code: !C Mesh Summary
  !C
        write(ILOG,"(a)") '###  Mesh Summary  ###'
        write(ILOG,"(a,i12)") '  Num of node:',hecMESH%n_node
        write(*   ,"(a,i12)") '  Num of node:',hecMESH%n_node
        write(ILOG,"(a,i12)") '  Num of DOF :',hecMESH%n_node*hecMESH%n_dof
        write(*   ,"(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
        write(*   ,"(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
          write(*   ,"(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(*   ,"(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),"[%]"
  
  !C
  !C Output sparsity pattern
  