The func2hcp.py script uses the -volume-parcel-resampling-generic command that we took out of our pipelines, in order to instead do cubic spline resampling, resolving an issue of it creating a grid pattern in the noise variance (since this command uses a form of linear resampling). We now use -cifti-resample after making a subject-specific subcortical cifti file at the acquired resolution. See here:
From Tim Coalson