1) This works well using modulo arithmetic. Try
something like:

1deval -expr 'bool(mod(t+1,4))' -num 240


2) One method for adding the extra 60 zero volumes
to the 180 volume time series would be to first
create a 181 volume time series by pre-pending
(or appending, whichever) the extra zero volume,
making it time index 0, leaving 1..180 as the real
data.

Then you just need a time series that counts from
1 to 180, but inserts a 0 after every 3 values,
and this can be use for sub-brick selection.

Note that you also need to pad a 1D time series
with zeros in the case of motion parameters and
possibly the enorm time series. Use 1d_tool.py
for that, as is done in this class data example,
AFNI_data6/FT_analysis/s17.proc.FT.rest.11.

Search for the string "zero pad censored TRs",
where 1d_tool.py -censor_fill_parent is applied.

- rick