diff --git a/src/optimisation/adjointOptimisation/adjoint/objectives/incompressible/objectiveUniformityCellZone/objectiveUniformityCellZone.C b/src/optimisation/adjointOptimisation/adjoint/objectives/incompressible/objectiveUniformityCellZone/objectiveUniformityCellZone.C
index 7314896024b6cb98268bfab4486067c863f553d0..0e57c1ff45076a32ce43fc6e82e82efc377d0bd6 100644
--- a/src/optimisation/adjointOptimisation/adjoint/objectives/incompressible/objectiveUniformityCellZone/objectiveUniformityCellZone.C
+++ b/src/optimisation/adjointOptimisation/adjoint/objectives/incompressible/objectiveUniformityCellZone/objectiveUniformityCellZone.C
@@ -5,8 +5,8 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
- Copyright (C) 2007-2023 PCOpt/NTUA
- Copyright (C) 2013-2023 FOSS GP
+ Copyright (C) 2007-2025 PCOpt/NTUA
+ Copyright (C) 2013-2025 FOSS GP
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
@@ -118,9 +118,9 @@ scalar objectiveUniformityCellZone::J()
const volVectorField& U = vars_.UInst();
const scalarField& V = mesh_.V().field();
- for (const label zI : zones_)
+ forAll(zones_, zI)
{
- const cellZone& zoneI = mesh_.cellZones()[zI];
+ const cellZone& zoneI = mesh_.cellZones()[zones_[zI]];
scalarField VZone(V, zoneI);
vectorField UZone(U.primitiveField(), zoneI);
volZone_[zI] = gSum(VZone);
@@ -138,9 +138,9 @@ void objectiveUniformityCellZone::update_dJdv()
{
const volVectorField& U = vars_.U();
- for (const label zI : zones_)
+ forAll(zones_, zI)
{
- const cellZone& zoneI = mesh_.cellZones()[zI];
+ const cellZone& zoneI = mesh_.cellZones()[zones_[zI]];
for (const label cellI : zoneI)
{
dJdvPtr_()[cellI] = (U[cellI] - UMean_[zI])/volZone_[zI];
@@ -154,9 +154,9 @@ void objectiveUniformityCellZone::update_divDxDbMultiplier()
volScalarField& divDxDbMult = divDxDbMultPtr_();
const volVectorField& U = vars_.U();
- for (const label zI : zones_)
+ forAll(zones_, zI)
{
- const cellZone& zoneI = mesh_.cellZones()[zI];
+ const cellZone& zoneI = mesh_.cellZones()[zones_[zI]];
for (const label cellI : zoneI)
{
divDxDbMult[cellI] =