following the "Converting a problem description file to direct problem construction" thread, I have removed the ambiguity in numbers of components of field vs. variable. What was before
u = FieldVariable('u', 'unknown', field, mesh.dim)
u = FieldVariable('u', 'unknown', field)
The variable's dimensions are given by the field. Nothing changes in the problem description files.