[blaise]

Dear All,

i am facing a problem with topoSet when i try to use the field cellLevel to select cells. topoSet cannot find the cellLevel file. I need this because I want to do mesh refinement in dependence of the previous refinement level. This worked in OpenFOAM 2.2.x. Currently i am using build 2.3.x-2d875bf8e478.

Here is my topoSetDict:

actions
(
    {
        name    FS;
        type    cellSet;
        action  new;
        source boxToCell;
        sourceInfo
        {
            box (-1000 -1000 -1.7) (1000 1000 1.7);
        }
    }
    // values of field within certain range
    {
        name    FS;
        type    cellSet;
        action  subset;    
        source  fieldToCell;
        sourceInfo
        {
            fieldName   cellLevel;   
            min         4;
            max         4;
        }
    }    
);

This is the error message:

...
Create polyMesh for time = 66

Reading topoSetDict

Time = 66
    mesh not changed.
Created cellSet FS
    Applying source boxToCell
    Adding cells with center within boxes 1((-1000 -1000 -1.7) (1000 1000 1.7))
    cellSet FS now size 406636
Read set cellSet FS with size 406636
    Applying source fieldToCell
--> FOAM Warning : 
    From function fieldToCell::applyToSet(const topoSetSource::setAction, topoSet& set)
    in file sets/cellSources/fieldToCell/fieldToCell.C at line 203
    Cannot read field cellLevel from time 66
    cellSet FS now size 0

End

There is a cellLevel file in 66/polyMesh/, but i suppose this is not the place it should be, at least in comparison to version 2.2.x. In 2.2.x this file would sit at the top level of the respective time directory, and could even be read by paraFoam for display of cellLevels.

I thought this issue could depend on the debug setting in snappyHexMeshDict. Setting this to 2 is meant to give:

2 : write volScalarField with cellLevel for postprocessing

but gives strangely numbered time directories (like 66, probably intermediate mesh states).

So, how can topoSet get access to cellLevel information again?

Thank you!