Using clusters with 3dRemlfit

Good morning
I would like to perform an analysis with 3dREMLfit and 3dMEMA based on permutations. I mean, I would like to perform the estimation many times on surrogate data, in order to build the distribution of the statistic voxel by voxel.
Is there any advantage I can get from cluster computing? I have the possibility to work on clusters, but I’d like to know if this Is supported by AFNI functions