Implementation of UMDFLUX subroutine

Is there any example fortran script for umdflux subroutine for the non uniform heat flux (goldak)?