diff --git a/src/overset/cellCellStencil/cellVolumeWeight/cellVolumeWeightCellCellStencil.C b/src/overset/cellCellStencil/cellVolumeWeight/cellVolumeWeightCellCellStencil.C index 4f02c9d466c48c7dffa7d2fcf728ef939b51c584..60464e71eb592426112c5fceb8340c3ca9063668 100644 --- a/src/overset/cellCellStencil/cellVolumeWeight/cellVolumeWeightCellCellStencil.C +++ b/src/overset/cellCellStencil/cellVolumeWeight/cellVolumeWeightCellCellStencil.C @@ -49,7 +49,7 @@ namespace cellCellStencils } Foam::scalar -Foam::cellCellStencils::cellVolumeWeight::defaultOverlapTolerance_ = 1e-9; +Foam::cellCellStencils::cellVolumeWeight::defaultOverlapTolerance_ = 1e-6; // * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * // @@ -953,7 +953,7 @@ bool Foam::cellCellStencils::cellVolumeWeight::update() // Interpolate. Check if enough overlap scalar v = mesh_.V()[cellI]; scalar overlapVol = sum(allWeights[cellI]); - if (overlapVol > (1.0-overlapTolerance_)*v) + if (overlapVol > overlapTolerance_*v) { allCellTypes[cellI] = INTERPOLATED; } @@ -971,7 +971,7 @@ bool Foam::cellCellStencils::cellVolumeWeight::update() } } - +/* // Knock out cell with insufficient interpolation weights forAll(allCellTypes, cellI) { @@ -989,7 +989,7 @@ bool Foam::cellCellStencils::cellVolumeWeight::update() } } } - +*/ if (debug) { tmp<volScalarField> tfld