[Lits]

Has anyone ever output the coordinates of all the grids in a celset in openfoam, or any other information about the grids?

#include "fvCFD.H"
//#include "dynamicFvMesh.H"
#include "CMULES.H"
#include "EulerDdtScheme.H"
#include "localEulerDdtScheme.H"
#include "CrankNicolsonDdtScheme.H"
#include "subCycle.H"
//#include "immiscibleIncompressibleTwoPhaseMixture.H"
//#include "turbulentTransportModel.H"
#include "pimpleControl.H"
#include "fvOptions.H"
#include "CorrectPhi.H"
#include "fvcSmooth.H"
#include <iomanip>
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

// -------------INPUT PARAMETERS ---------------------


int main(int argc, char *argv[])
{
timeSelector::addOptions();

#include "setRootCase.H"
#include "createTime.H"
instantList timeDirs = timeSelector::select0(runTime, args);
#include "createMesh.H"
//#include "cellZoneSet.H"
#include "mapPolyMesh.H"
#include "polyMesh.H"

//=========finding grids ========

labelList cellZoneID = mesh.cellZones().findIndices("anodeGDL");
Info << "cellZoneID = " << cellZoneID << endl;
label cellZoneID2 = mesh.cellZones().findZoneID("flow");
Info << "cellZoneID2 = " << cellZoneID2 << endl;
const cellZone& cz = mesh.cellZones()[cellZoneID];
const labelList& cells = mesh.cellZones()[cellZoneID[0]];;
forAll(cells, i)
{

label celli = cells[i];
Info << "1" << endl;
}



return 0;
//}
}

I tried the above code and it compiles successfully, but it runs with an error:


#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 ? in "/home/lts/OpenFOAM/lts-6/platforms/linux64GccDPInt32Opt/bin/setBINDER2_invert"
#3 ? in "/home/lts/OpenFOAM/lts-6/platforms/linux64GccDPInt32Opt/bin/setBINDER2_invert"
#4 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#5 ? in "/home/lts/OpenFOAM/lts-6/platforms/linux64GccDPInt32Opt/bin/setBINDER2_invert"
Aborted (core dumped)

I will be appreciate for any advice.