Permeability Subroutine in Abaqus

If I gonna make spatially varying permeability value into abaqus using subroutine, what kind of subroutine I can use?. I have a .txt file containing x, y, z, and the value of permeabilty. How can I call this file into my fortran programming code? Do anyone have an example to define such this subroutine?.Thank you