From 4a37b52dca2ba09c927bf46a95d013a7087a140b Mon Sep 17 00:00:00 2001 From: andy Date: Wed, 25 Sep 2013 16:12:06 +0100 Subject: [PATCH] Revert "ENH: Updated forces and forceCoeffs bin writing" This reverts commit 1e711927a76bb17b8e4dad6c8c7cc60ab3a027f5. --- .../forces/forceCoeffs/forceCoeffs.C | 25 ++++---- .../functionObjects/forces/forces/forces.C | 60 +++++++++++-------- 2 files changed, 49 insertions(+), 36 deletions(-) diff --git a/src/postProcessing/functionObjects/forces/forceCoeffs/forceCoeffs.C b/src/postProcessing/functionObjects/forces/forceCoeffs/forceCoeffs.C index 7494702f5f..7f094c2e3f 100644 --- a/src/postProcessing/functionObjects/forces/forceCoeffs/forceCoeffs.C +++ b/src/postProcessing/functionObjects/forces/forceCoeffs/forceCoeffs.C @@ -63,14 +63,16 @@ void Foam::forceCoeffs::writeFileHeader(const label i) << "# bins : " << nBin_ << nl << "# start : " << binMin_ << nl << "# delta : " << binDx_ << nl - << "# direction : " << binDir_ << nl; + << "# direction : " << binDir_ << nl + << "# Time"; - file(i) - << "# Time" - << tab << "bin" - << tab << "Cm" - << tab << "Cd" - << tab << "Cl"; + for (label j = 0; j < nBin_; j++) + { + const word jn = Foam::name(j); + + file(i) + << tab << "Cm" << jn << tab << "Cd" << jn << tab << "Cl" << jn; + } } else { @@ -214,16 +216,17 @@ void Foam::forceCoeffs::write() } } + file(1)<< obr_.time().value(); + forAll(coeffs[0], i) { file(1) - << obr_.time().value() - << tab << i << tab << coeffs[2][i] << tab << coeffs[1][i] - << tab << coeffs[0][i] - << endl; + << tab << coeffs[0][i]; } + + file(1) << endl; } if (log_) diff --git a/src/postProcessing/functionObjects/forces/forces/forces.C b/src/postProcessing/functionObjects/forces/forces/forces.C index 9830719c86..d20434eb52 100644 --- a/src/postProcessing/functionObjects/forces/forces/forces.C +++ b/src/postProcessing/functionObjects/forces/forces/forces.C @@ -94,19 +94,29 @@ void Foam::forces::writeFileHeader(const label i) << "# bins : " << nBin_ << nl << "# start : " << binMin_ << nl << "# delta : " << binDx_ << nl - << "# direction : " << binDir_ << nl; + << "# direction : " << binDir_ << nl + << "# Time"; - file(i) - << "# Time" - << tab << "bin" - << tab << "forces(pressure,viscous,porous)" - << tab << "moment(pressure,viscous,porous)"; + for (label j = 0; j < nBin_; j++) + { + const word jn = Foam::name(j); + file(i) + << tab + << "forces" << jn << "(pressure,viscous,porous) " + << "moment" << jn << "(pressure,viscous,porous)"; + } if (localSystem_) { - file(i) - << tab << "localForces(pressure,viscous,porous)" - << tab << "localMoment(pressure,viscous,porous)"; + for (label j = 0; j < nBin_; j++) + { + const word jn = Foam::name(j); + + file(i) + << tab + << "localForces" << jn << "(pressure,viscous,porous) " + << "localMoments" << jn << "(pressure,viscous,porous)"; + } } } else @@ -425,17 +435,27 @@ void Foam::forces::writeBins() f[0][i] += f[0][i-1]; f[1][i] += f[1][i-1]; f[2][i] += f[2][i-1]; + m[0][i] += m[0][i-1]; m[1][i] += m[1][i-1]; m[2][i] += m[2][i-1]; } } - List > lf(3); - List > lm(3); + file(1) << obr_.time().value(); + + forAll(f[0], i) + { + file(1) + << tab + << "(" << f[0][i] << "," << f[1][i] << "," << f[2][i] << ") " + << "(" << m[0][i] << "," << m[1][i] << "," << m[2][i] << ")"; + } if (localSystem_) { + List > lf(3); + List > lm(3); lf[0] = coordSys_.localVector(force_[0]); lf[1] = coordSys_.localVector(force_[1]); lf[2] = coordSys_.localVector(force_[2]); @@ -455,27 +475,17 @@ void Foam::forces::writeBins() lm[2][i] += lm[2][i-1]; } } - } - forAll(f[0], i) - { - file(1) - << obr_.time().value() - << tab << i - << tab << "(" << f[0][i] << "," << f[1][i] << "," << f[2][i] << ")" - << tab << "(" << m[0][i] << "," << m[1][i] << "," << m[2][i] << ")"; - - if (localSystem_) + forAll(lf[0], i) { file(1) << tab - << "(" << lf[0][i] << "," << lf[1][i] << "," << lf[2][i] << ")" - << tab + << "(" << lf[0][i] << "," << lf[1][i] << "," << lf[2][i] << ") " << "(" << lm[0][i] << "," << lm[1][i] << "," << lm[2][i] << ")"; } - - file(1) << endl; } + + file(1) << endl; }