sfunif3 now accepts an optional layers= parameter to input a file with lateral distribution of properties inside each layer. Note that the number of layers is one more than the number of interfaces.
See rsf/rsf/unif3 for examples.