[kdarc]

Hi.
I have been able to create cellZones in my mesh, which as best as I can tell are simply sets of cell indices.

My question is about how-to:

How do I make use of cell zones??
For example, I want to give a value to a field only in a certain cellZone.

I imagine this operation to be akin to defining boundary conditions on only specific patches, but so far I have been unable to find any example code for cellZones.

I found something relevant looking in a file called mixerFvMesh.C, shown below:

const labelList& cellAddr =
cellZones()[cellZones().findZoneID("carboxen")];

forAll (cellAddr, cellI)
{
k_ads[cellI] = k0_ads*Foam::exp(-E_ads/(R_gas*T[cellI])) ; // rate of adsoption of HCHO onto carboxen [1/s]

}

Are there some .H files or other preparations I would need to get this code to compile?? A compile library??

my error output is:

cEqn.H: In function âint main(int, char**)â:
cEqn.H:5: error: no matching function for call to Foam::cellZone::cellZone()â

any help greatly appreciated.

[kdarc]

An in-house colleague kindly sorted this out for me.

The syntax is:

const label cellZoneID = mesh.cellZones().findZoneID("carboxen");
const labelList& cellAddr = mesh.cellZones()[cellZoneID];

etc etc...