afni.nimh.nih.gov

That -glueto approach could be simplified up a bit...

If one of your 3d+time datasets is 710 TRs long, then
it would be enough to do:

3dcalc -a dset.710+tlrc -expr 0 -prefix zero.710

If not, presumably 4blocks_zeros_ones_for_convolution.1D
is still 710 TRs. long. In that case it would be enough to
take a single volume and that 710 TR 1D time series as
inputs, and still output zeros:

3dcalc -a mean_iso_HRF+orig"[0]"                \
       -b 4blocks_zeros_ones_for_convolution.1D \
       -expr 0 -prefix zero.710

- rick