cccccccccccccccccccccccccccccccccccccccccc
c     com6.h : common var.
cccccccccccccccccccccccccccccccccccccccccc

c------------------------------------------------------
c     np: diameter of the disk
c-------------------------------------------------------
      integer msy,all,allw,bnd,bndw,nbndw
      parameter(msy=4)

      integer cntr,ncl,nsw
      integer a(4001),b(4001),aw(4001),bw(4001)
      integer ccl(4001)
      real*8 nl(4001)

c-------------------------------------------------------
c     n: triangular coordinate → particle number
c     nx,ny: particle number → triangular coordinate
c     ccl: particle num. of the lower edge of the disk
c-------------------------------------------------------

      common/iso1/ccl
      common/iso2/nl

c------------------------------------------------------
c     cntr: number of center particle of the disk
c     ncl: number of the edge particles of the disk
c------------------------------------------------------

      common/var/cntr,ncl,nsw,nbndw
      common/var2/all,allw,bnd,bndw

c------------------------------------------------------
c     a,b: particle number at the end of the bond
c------------------------------------------------------

      common/edg/a,b
      common/edg2/aw,bw

c-------------------------------------------------------
c     nwx: number of particles in one layer(must be even)
c     nwy: number of layer
c     nb: number of particles between surface particles
c-------------------------------------------------------

c      integer nwx,nwy
c      parameter (nwx=75,nwy=16)

c-------------------<NEIGHBOR LIST>----------------------

      integer lgt(4001),rgt(4001)
      common/block4/lgt,rgt

c-----------------------<ENERGY>---------------------------

      real*8 U(4001)
      real*8 uw(4001)
      common/block2/U
      common/block3/uw

c-----------------------<FORCE>----------------------------

      real*8 fxw(4001),fyw(4001),fx(4001),fy(4001)

      common /wforce/ fxw,fyw
      common/posi/fx,fy

c------------------<parameters for SI>--------------------
      real*8  c(msy),d(msy)
      common /symph/ c,d

c---------------<mean value of contact point>------------------
      real*8 xm,ym,vxm,vym
      common/mn/xm,ym,vxm,vym

c---------------<initial bond length>------------------
c      real*8 lh(nwy,nwx+2),lv(nwy,nwx+2),lrt(nwy,nwx+2),llt(nwy,nwx+2)
c      common/len/lh,lv,lrt,llt
      real*8 lw(4001)
      common/len2/lw

c----------------<SN2NN>---------------------
      integer sn2nn(4001)
      common/sn/sn2nn