ENH: surfaceNoise - warn if no data available to process

This commit is contained in:
Andrew Heather
2017-03-06 12:09:05 +00:00
parent f3abf405a8
commit a2f54aef3c

View File

@ -617,53 +617,62 @@ void surfaceNoise::calculate()
scalarField PSDfAve(nFreq, 0); scalarField PSDfAve(nFreq, 0);
scalarField fOut(nFreq, 0); scalarField fOut(nFreq, 0);
for (label i = f0; i <= f1; ++i) if (nFreq == 0)
{ {
label freqI = i*fftWriteInterval_; WarningInFunction
fOut[i] = freq1[freqI]; << "No surface data available using a fftWriteInterval of "
const word gName = "fft"; << fftWriteInterval_ << endl;
}
else
{
for (label i = f0; i <= f1; ++i)
{
label freqI = i*fftWriteInterval_;
fOut[i] = freq1[freqI];
const word gName = "fft";
PrmsfAve[i] = writeSurfaceData PrmsfAve[i] = writeSurfaceData
( (
fNameBase, fNameBase,
gName, gName,
"Prmsf", "Prmsf",
freq1[freqI], freq1[freqI],
surfPrmsf[i], surfPrmsf[i],
procFaceOffset, procFaceOffset,
writePrmsf_ writePrmsf_
); );
PSDfAve[i] = writeSurfaceData PSDfAve[i] = writeSurfaceData
( (
fNameBase, fNameBase,
gName, gName,
"PSDf", "PSDf",
freq1[freqI], freq1[freqI],
surfPSDf[i], surfPSDf[i],
procFaceOffset, procFaceOffset,
writePSDf_ writePSDf_
); );
writeSurfaceData writeSurfaceData
( (
fNameBase, fNameBase,
gName, gName,
"PSD", "PSD",
freq1[freqI], freq1[freqI],
noiseFFT::PSD(surfPSDf[i]), noiseFFT::PSD(surfPSDf[i]),
procFaceOffset, procFaceOffset,
writePSD_ writePSD_
); );
writeSurfaceData writeSurfaceData
( (
fNameBase, fNameBase,
gName, gName,
"SPL", "SPL",
freq1[freqI], freq1[freqI],
noiseFFT::SPL(surfPSDf[i]*deltaf), noiseFFT::SPL(surfPSDf[i]*deltaf),
procFaceOffset, procFaceOffset,
writeSPL_ writeSPL_
); );
}
} }
graph Prmsfg graph Prmsfg