I wonder whether there is an h-adaptive multi-dimensional integration package, containing the algorithm like hcubarture package of Julia language, for modern Fortran?