3dLME documentation states, "3dLME only allows for one voxelwise covariate per subject". If I have multiple visits per subject, is it possible to use voxelwise covariates?
Currently, when I attempt this by specifying a separate voxelwise covariate file per subject visit the design matrix appears to process correctly but 3dLME halts after loading inputs with this error:

***** End of data structure information *****
++++++++++++++++++++++++++++++++++++++++++++++++++++
Reading input files now...
Reading input files: Done!
Range of input data: [-0.460, 0.891]

If the program hangs here for more than, for example, half an hour,
kill the process because the model specification or the GLT coding
is likely inappropriate.

[1] "Great, test run passed at voxel (54614, 1, 1)!"
[1] "Start to compute 1 slices along Z axis. You can monitor the progress"
[1] "and estimate the total run time as shown below."
[1] "05/03/21 19:12:31.496"
Error in rbind(inData, array(0, dim = c(fill, NoFile))) : 
  number of columns of matrices must match (see arg 2)
Execution halted

In a separate test, I instead specified the subject mean voxelwise covariate at each visit, ie the same mean file for each visit. This works for qVars but does not appear to work for vVars.

Here is an abbreviated command:

3dLME
                -model "DX*VISIT+COVARFILE" \
                -vVars "COVARFILE" \
                -ranEff '~1' \
                -SS_type 3 \

If i remove the voxelwise covariate from the model but not the datatable, the model runs fine.

thanks