Cluster correction in nilearn

Hi everybody,
I am doing a whole brian searchlight analysis using nilearn and I used two step permutation test to calculate the p vals in the second level.
I want to apply cluster correction but I see it is not implemeted in nilearn and I wonder is implemented in some other python package beside PyMVPA? Have someone already implemetend it using nilearn tools?
If not, than if I need to implement it myself what is must common straight forward way of doing cluster correction? Should I just use the correction from Stelzer et al. 2013 with initial pval of 0.001?
Would appriciate any help or suggestion.
Thanks,
Tamir