[jans]

Hi,
For a particular application, I have written a code that creates the mesh object using the constructor

fvMesh (const IOobject &io, const Xfer< pointField > &points, const Xfer< faceList > &faces, const Xfer< labelList > &allOwner, const Xfer< labelList > &allNeighbour, const bool syncPar=true)

instead of the usual mesh constructor that reads in usual OpenFOAM files.

However, the fvSchemes information are absent from the resulting mesh object, i.e., mesh.schemesDict() will return empty braces {}. I expected this, and thought I can afterwards insert this missing information into the mesh object, but I am unable to find the right method to do this.

Data members like mesh.ddtSchemes_ etc are private and so cannot be edited later. E.g.,

error #308: member "Foam::fvSchemes::ddtSchemes_" (declared at line 57 of "/work/01255/siliu/OpenFOAM/OpenFOAM-2.3.0/src/finiteVolume/lnInclude/fvSchemes.H") is inaccessible
      mesh.ddtSchemes_.add(keyType("ref"),7.0);

mesh.schemesDict() returns a constant pointer, and hence cannot be edited later too. E.g.:

error: no instance of overloaded function "Foam::dictionary::add" matches the argument list and object (the object has cv-qualifiers that prevent a match)
In file included from icoFoam_iohdf5.C(44):
            argument types are: (Foam::keyType, double)
In file included from icoFoam_iohdf5.C(44):
            object type is: const Foam::dictionary
In file included from icoFoam_iohdf5.C(44):
      mesh.schemesDict().add(keyType("ref"),7.0);

Is there a way to do the above?
Thanks much in advance
Anirban

[marupio]

It looks like fvSchemes is hard-coded to read from a file, and if you construct the mesh with an IOobject whose read option is NO_READ, fvSchemes will be an empty dictionary. After it is constructed, you can add entries to its dictionary, but it never uses those... the schemes dictionaries it keeps are copies of the dictionary it had on construction. If you call read(), it will, again, try reading the file, and return false, indicating the read option is NO_READ. The quickest way I can see to bootstrap fvSchemes is to make fvSchemes::read(const dictionary& dict) a public function. Then you create a dictionary, add entries to it, and call mesh.fvSchemes::read(dict);

You may also experience problems with fvSolution.

Good luck!