#!/bin/tcsh -f

if ("$1" == "-clean") then
   \rm -f v2s.lh.DEL.1D.dset
   \rm -f v2s.lh.DEL.1D.col
   goto END
endif 

if ( ! -f v2s.lh.DEL.niml.dset ) then
   echo "Missing data. Need to run first:"
   echo "   tcsh run_3dVol2Surf"
   echo ""
   goto BEND
endif

#ScaleToMap is not all that useful anymore.
#You typically colorize the data interactively in SUMA.
#

#Create a .1D version of the .niml dset
ConvertDset    -input v2s.lh.DEL.niml.dset \
               -prefix v2s.lh.DEL.1D.dset \
               -prepend_node_index_1D

#delete output file (from previous run)
if ( -f v2s.lh.DEL.1D.col) rm -f v2s.lh.DEL.1D.col
#run the command for mapping data values to color map
ScaleToMap     -input   v2s.lh.DEL.1D.dset 0  1         \
               -cmap    afni_p20                \
               -apr     30                      \
               -msk     -3.0 -0.1               \
               -nomsk_col                       \
               > v2s.lh.DEL.1D.col

END: