! begin file my_e_matrix_inc
! ...................................................................
! ** GET_APPLICATION_E_MATRIX PROBLEM DEPENDENT STATEMENTS FOLLOW **
!                 For REAL (DP) :: E (N_R_B, N_R_B) 
!                 Given INTEGER,   INTENT(IN)  :: IE
!                       REAL (DP), INTENT(IN)  :: XYZ (N_SPACE)
!   Given modules: Elem_Type_Data, 
!     System_Constants, Sys_Properties_Data
! ...................................................................
! Define any new array or variable types, then give statements
! INTEGER,  SAVE ::                           ! Needed on next call
! REAL(DP), SAVE ::                           ! Needed on next call
! LOGICAL,  SAVE ::                           ! Needed on next call
!
! INTEGER  ::                                 ! Local use
! REAL(DP) ::                                 ! Local use
! LOGICAL  ::                                 ! Local use
  INTEGER, SAVE :: NOTE = 0
    IF ( NOTE == 0 ) THEN ; NOTE = 1;
      IF ( DEBUG_E .OR. DEBUG_INCLUDE ) PRINT *, &
        'NOTE: USED my_e_matrix_inc'
    END IF

!       APPLICATION DEPENDENT CONSTANT CONSTITUTIVE MATRIX 
    PRINT *,'NO APPLICATION my_e_matrix_inc PROVIDED'
    CALL REAL_IDENTITY (N_R_B, E) ! default
!                           or
!   CALL REAL_IDENTITY (N_R_B, E) !   DEFAULT TO IDENTITY MATRIX
!                           or
!--->  DEFINE PROPERTIES (PLANE STRESS ELASTICITY)
!     Y  = GET_FLT_LP (IE, 1)  ! ELASTIC MODULUS
!     PR = GET_FLT_LP (IE, 2)  ! POISSON'S RATIO
!
!      FORM THE CONSTITUTIVE MATRIX
! CALL PLANE_STRESS_E_MATRIX (GET_FLT_LP (IE, 1),  &
!                             GET_FLT_LP (IE, 2), E)
!                           or
!--->  DEFINE PROPERTIES (POISSON EQUATION)
!     K  = GET_FLT_LP (IE, 1)  ! CONDUCTIVITY
!
!      FORM THE CONSTITUTIVE MATRIX
! CALL POISSON_ISOTROPIC_E_MATRIX (GET_FLT_LP (IE, 1),  E)
!
!                           etc
! end file my_e_matrix_inc