diff --git a/src/postProcessing/functionObjects/forces/forceCoeffs/forceCoeffs.C b/src/postProcessing/functionObjects/forces/forceCoeffs/forceCoeffs.C index 7494702f5f8ae49e19d1cc133abe7f38b76cec36..7f094c2e3fd06268a72301b5876ad6fde15920c7 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 9830719c86cfe7e970fd6217c36010d1c41da5e5..d20434eb52ca9705d2aa7c53a53375e33338e440 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<Field<vector> > lf(3); - List<Field<vector> > 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<Field<vector> > lf(3); + List<Field<vector> > 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; }