I should also note that instead of calculating nonlinear alignment within afni_proc.py via:
-tlrc_base MNI152_2009_template.nii.gz \
... we typically nowadays recommend running @SSwarper to calculate both the skullstripping (SS) and nonlinear warp to standard space.
The command would be run ahead of time, e.g.:
@SSwarper \
-input DSET_ANAT \
-base MNI152_2009_template_SSW.nii.gz \
-subid ${subj} \
-odir OUT_DIR
... and then, as noted here in this sliiightly modified clip of the @SSwarper help file, you provide both the skullstripped anatomical (and turning off further skullstripping) and the warp dsets:
| afni_proc.py \
| [...other stuff here: processing blocks, options...] \
| -copy_anat OUT_DIR/anatSS.${subj}.nii \
| -anat_has_skull no \
| ....
| -volreg_tlrc_warp -tlrc_base MNI152_2009_template_SSW.nii.gz \
| -tlrc_NL_warp \
| -tlrc_NL_warped_dsets \
| OUT_DIR/anatQQ.${subj}.nii \
| OUT_DIR/anatQQ.${subj}.aff12.1D \
| OUT_DIR/anatQQ.${subj}_WARP.nii
--pt