#!/bin/tcsh
set suma_dir = mysumadir       # usually a directory called SUMA
set fslabel_dir = myFSlabeldir  # this is a parallel directory to your existing SUMA output directory called "label"
set hands = (lh rh)

cd $fslabel_dir
# make the annotation labels on the original FreeSurfer surfaces for each hemisphere
foreach hand ( $hands )
  FSread_annot -input ${hand}.aparc.DKTatlas40.annot -roi_1D $suma_dir/${hand}.aparc.DKTatlas40.annot.1D.roi \
  -dset  $suma_dir/${hand}.aparc.DKTatlas40.annot.niml.dset -cmap_1D $suma_dir/${hand}.aparc.DKTatlas40.annot.1D.cmap \
  -FScmap aparc.annot.DKTatlas40.ctab
end

cd $suma_dir
# transform labels to standard meshes
foreach hand ( $hands )
   MapIcosahedron -spec $suma_dir/FS_${hand}.spec -ld 60 -NN_dset_map $suma_dir/${hand}.aparc.DKTatlas40.annot.niml.dset \
    -morph $suma_dir/${hand}.sphere.reg.gii -prefix std.60.DKT.{$hand}.
end