Hi, Mikey-
Yes, please don't use 3dcalc for this... yikes!
If you are asking about taking one input 4D dataset called DSET_4D and calculate the sum across the 4th dimension (=time), then you can use:
3dTstat -sum -prefix dset_sum.nii.gz DSET_4D
For future reference, this program has lots of options for calculating voxelwise quantities across time.
Then you can use 3dcalc to divide the result by whatever you want, such as with this unlikely candidate:
3dcalc -a dset_sum.nii.gz -expr 'a/3.14159' -prefix dset_sum_div_pi.nii.gz
And indeed, if your question was to perform voxelwise addition over a set of 3D datasets, you could use the following:
3dmean -sum -prefix all_sum.nii.gz DSET1 DSET2 DSET3 ....
--pt