ecs.Q().rmap(Region::CONDUCTOR);
            ecs.Q().mapInternalField(Region::THERMAL);
            //fvc::extrapolate(ecs.Q()[Region::THERMAL]);

I have this mapping of region fields The fields overlap, but when doing extrapolate I get floating point errors, I am also not so sure if I need to extrapolate as mentioned they overlap, and the values are cellCentered for instance this is volumetric heating W/m^3