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