Hi all, 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) is now u = FieldVariable('u', 'unknown', field) The variable's dimensions are given by the field. Nothing changes in the problem description files. r.