#!/bin/tcsh -f

set subj = 620
set hdir = $PWD

@demo_prompt 'Make sure all SUMA and AFNI sessions are closed'
if ($status) goto END



suma  -spec SUMA/s${subj}_both.mini.spec     \
      -sv s${subj}_SurfVol_Alnd_Exp+orig. \
      -niml -dev  >& runlog   &
      
sleep 10
      
@demo_prompt 'Wait until SUMA is up'
if ($status) goto END

echo "Loading time series data on both sides" 
foreach hemi (lh rh)
   DriveSuma -com surf_cont   \
             -label     ${hemi}.smoothwm.asc \
             -I_range 500  \
             -T_val 0.25   \
             -load_dset ${hemi}.s${subj}_rest_r1.clean.wm.sm4.niml.dset  \
                  >& /dev/null
   sleep 5          
end

@demo_prompt 'Wait until Dsets are loaded'
if ($status) goto END
echo "Opening a graph window on the left hemisphere for the loaded dset"
DriveSuma   -com viewer_cont  \
               -key:v100 j \
               -key g                           


@demo_prompt 'When ready, go for InstaCorr'
if ($status) goto END
echo "Starting InstaCorr"
DriveSuma   -com viewer_cont  \
               -key:v100 j \
               -key D                           

@demo_prompt 'When ready, go setting the threshold'
if ($status) goto END
DriveSuma   -com surf_cont \
               -label     lh.smoothwm.asc \
               -T_val 0.25 \
            -com surf_cont \
               -label     rh.smoothwm.asc \
               -T_val 0.25 \
            
               
@demo_prompt 'When ready, start AFNI'
if ($status) goto END
afni -niml -yesplugouts &
sleep 4
DriveSuma   -com viewer_cont -key 't'
plugout_drive -com 'OPEN_WINDOW A.axialimage ' \
               -com "SWITCH_UNDERLAY s${subj}_rest_r1.clean"   \
               -com 'OPEN_WINDOW B.axialimage ' \
               -com "SWITCH_UNDERLAY B.s${subj}_SurfVol_Alnd_Exp"   \
               -quit
echo ""
echo "clicking in AFNI sends the time series to SUMA and performs correlation"
echo "with time series on surface"
echo ""