diff --git a/src/postProcessing/functionObjects/forces/forceCoeffs/forceCoeffs.C b/src/postProcessing/functionObjects/forces/forceCoeffs/forceCoeffs.C index 7f094c2e3fd06268a72301b5876ad6fde15920c7..7494702f5f8ae49e19d1cc133abe7f38b76cec36 100644 --- a/src/postProcessing/functionObjects/forces/forceCoeffs/forceCoeffs.C +++ b/src/postProcessing/functionObjects/forces/forceCoeffs/forceCoeffs.C @@ -63,16 +63,14 @@ void Foam::forceCoeffs::writeFileHeader(const label i) << "# bins : " << nBin_ << nl << "# start : " << binMin_ << nl << "# delta : " << binDx_ << nl - << "# direction : " << binDir_ << nl - << "# Time"; + << "# direction : " << binDir_ << nl; - for (label j = 0; j < nBin_; j++) - { - const word jn = Foam::name(j); - - file(i) - << tab << "Cm" << jn << tab << "Cd" << jn << tab << "Cl" << jn; - } + file(i) + << "# Time" + << tab << "bin" + << tab << "Cm" + << tab << "Cd" + << tab << "Cl"; } else { @@ -216,17 +214,16 @@ 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]; + << tab << coeffs[0][i] + << endl; } - - file(1) << endl; } if (log_) diff --git a/src/postProcessing/functionObjects/forces/forces/forces.C b/src/postProcessing/functionObjects/forces/forces/forces.C index 38ab09ada640f9eee81f01ba1443b4180f626fef..ee739be63e92e8190f49a8a913449fc7c2c09c08 100644 --- a/src/postProcessing/functionObjects/forces/forces/forces.C +++ b/src/postProcessing/functionObjects/forces/forces/forces.C @@ -94,29 +94,19 @@ void Foam::forces::writeFileHeader(const label i) << "# bins : " << nBin_ << nl << "# start : " << binMin_ << nl << "# delta : " << binDx_ << nl - << "# direction : " << binDir_ << nl - << "# Time"; + << "# direction : " << binDir_ << nl; - for (label j = 0; j < nBin_; j++) - { - const word jn = Foam::name(j); + file(i) + << "# Time" + << tab << "bin" + << tab << "forces(pressure,viscous,porous)" + << tab << "moment(pressure,viscous,porous)"; - file(i) - << tab - << "forces" << jn << "(pressure,viscous,porous) " - << "moment" << jn << "(pressure,viscous,porous)"; - } if (localSystem_) { - 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)"; - } + file(i) + << tab << "localForces(pressure,viscous,porous)" + << tab << "localMoment(pressure,viscous,porous)"; } } else @@ -435,27 +425,17 @@ 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]; } } - 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] << ")"; - } + List<Field<vector> > lf(3); + List<Field<vector> > lm(3); 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]); @@ -475,17 +455,27 @@ 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] << ")"; - forAll(lf[0], i) + if (localSystem_) { file(1) << tab - << "(" << lf[0][i] << "," << lf[1][i] << "," << lf[2][i] << ") " + << "(" << lf[0][i] << "," << lf[1][i] << "," << lf[2][i] << ")" + << tab << "(" << lm[0][i] << "," << lm[1][i] << "," << lm[2][i] << ")"; } - } - file(1) << endl; + file(1) << endl; + } }