I am using functional data with this. I processed it through proc.py, and the defaults on that appear to be smoothwm for surfA and pial for surfB, so it sounds like the slow_surf_clustsim.py defaults are fine, I just need to change my SurfClust surf_A option to smoothwm.