readKMat Subroutine

    subroutine readKMat()
        implicit none
        character(*), parameter :: t_r = "readKMat"
        integer :: ierr
        integer(int64) :: kMatSize
        call kMatLin%readKMatFromFile("KDUMPLIN")
        call kMatSq%readKMatFromFile("KDUMPSQ")

        ! now, take care of the projected energy contribution if required
        if (tLogKMatProjE) then
            kMatSize = determineKMatSize()
            allocate(kMatProjEContrib(kMatSize), stat=ierr)
            kMatProjEContrib = 0.0_dp
            if (ierr /= 0) call stop_all(t_r, "Allocation failed")
        end if

    end subroutine readKMat