Does 3dLMEr -vVar option accept surface (.niml.dset)?

Based on your current setting, the model should be specified as follows:

-model 'TREAT*BL+AGE+TREAT*VIS_ORDER+(1|Subj)'

For a more detailed breakdown of model specification and syntax within 3dLMEr, I recommend checking out this blog post.

Gang Chen