hi_constant.f90

module hi_const
    
	integer,parameter   :: ngr = -10
    integer             :: ngl = -10 ! how many gaussian point on 

    real(8),parameter   :: dlt(3,3) = RESHAPE((/1.,0.,0.,0.,1.,0.,0.,0.,1./),(/3,3/))
    ! d-matrix--unix matix
    
    real(8),parameter   :: hi_PI  = 4.D0*DATAN(1.D0)
    real(8),parameter   :: cnu = 0.3
    real(8),parameter   :: TOLGP=1.D-14


    real(8),parameter :: cnr_lcl_mtx(18) =(/-1.,-1.,1.,-1.,1.,1.,-1.,1.,0.,-1.,1.,0.,0.,1.,          &
        &         -1.,0.,0.,0./) ! 9-point element

    integer,parameter :: node_grp_by_edge(12) = (/1,2,5, 2,3,6, 3,4,7, 4,1,8/)

    real(8),parameter :: TOL=1.D-5 !tolerence?

    integer,parameter :: KSB(4)=(/-2,1,2,-1/) !some matrix,delcared
    ! IF(DABS(XIP(IABS(KS))-DBLE(KS)/DABS(DBLE(KS))).LT.TOL) GOTO 80  
    ! edge 1 , xip(2) != -1
    ! edge 2 , xip(1) != 1
    ! edge 3 , xip(2) != 1
    ! edge 4 , xip(1) != -1
    ! why?

end module hi_const