I'm having difficulty figuring out what I'm doing wrong. I'm getting the following output/error:

++ '-stim_times_AM1 1 stimuli/happy.1D' has 1 auxiliary values per time point
++ '-stim_times_AM1 1': basis function model 'dmBLOCK' uses 1 parameters,
    out of the 1 found in timing file 'stimuli/happy.1D'
++ '-stim_times_AM1 2 stimuli/angry.1D' has 1 auxiliary values per time point
++ '-stim_times_AM1 2': basis function model 'dmBLOCK' uses 1 parameters,
    out of the 1 found in timing file 'stimuli/angry.1D'
++ '-stim_times_AM1 3 stimuli/sad.1D' has 1 auxiliary values per time point
++ '-stim_times_AM1 3': basis function model 'dmBLOCK' uses 1 parameters,
    out of the 1 found in timing file 'stimuli/sad.1D'
++ '-stim_times_AM1 4 stimuli/neutral.1D' has 1 auxiliary values per time point
++ '-stim_times_AM1 4': basis function model 'dmBLOCK' uses 1 parameters,
    out of the 1 found in timing file 'stimuli/neutral.1D'
** FATAL ERROR: 3dDeconvolve dies: Use -num_glt option to specify number of GLTs when more than 10

I've tried it with and without the -num_glt option and I get the same error each time, so I'm wondering if the problem might be related to my 1D files instead?

Here is what a sample 1D file contains:

2.994:0.605 17.987:0.6 29.999:0.531 39.001:0.531 41.995:0.604 65.989:0.509 92.993:0.603 95.989:0.631 101.995:0.681 113.991:0.583 125.988:0.584 146.986:0.577 158.999:0.615 168.0:0.591 179.997:0.594 188.998:0.651 194.989:0.636 198.0:0.584 203.989:0.57 215.987:0.578 218.998:0.545 239.996:0.561 248.997:0.605 251.992:0.663 281.992:0.619 302.989:0.627 323.988:0.624 326.999:0.657 338.996:0.652 368.995:0.632 383.987:0.613 386.998:0.588 407.996:0.61 410.991:0.615 419.993:0.58 422.988:0.641 435.001:0.596 449.992:0.684 470.99:0.545 473.986:0.613

I would appreciate any advice!

And you're right about selecting HA[1..4]. I was modifying an older script that used a tent function. Thanks for pointing that out.