from smriprep.workflows.surfaces import init_morph_grayords_wf
wf = init_morph_grayords_wf(grayord_density="91k", omp_nthreads=1)