From e534882f5aed12fa5e2592babd55c915e26d5ad5 Mon Sep 17 00:00:00 2001
From: andy <a.heather@opencfd.co.uk>
Date: Tue, 17 Mar 2009 18:33:27 +0000
Subject: [PATCH] adding warning messages
---
.../Templates/KinematicParcel/KinematicParcel.C | 13 +++++++++++++
.../Templates/ReactingParcel/ReactingParcel.C | 12 ++++++++++++
.../parcels/Templates/ThermoParcel/ThermoParcel.C | 13 +++++++++++++
3 files changed, 38 insertions(+)
diff --git a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C
index 637bf15fa4c..a6a3865ba18 100644
--- a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C
@@ -39,6 +39,19 @@ void Foam::KinematicParcel<ParcelType>::updateCellQuantities
)
{
rhoc_ = td.rhoInterp().interpolate(this->position(), cellI);
+ if (rhoc_ < SMALL)
+ {
+ WarningIn
+ (
+ "void Foam::KinematicParcel<ParcelType>::updateCellQuantities"
+ "("
+ "TrackData&, "
+ "const scalar, "
+ "const label"
+ ")"
+ ) << "Density < " << SMALL << " in cell " << cellI << nl << endl;
+ }
+
Uc_ = td.UInterp().interpolate(this->position(), cellI);
muc_ = td.muInterp().interpolate(this->position(), cellI);
diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
index 4ff61b6602a..8b59c5017a8 100644
--- a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
@@ -40,6 +40,18 @@ void Foam::ReactingParcel<ParcelType>::updateCellQuantities
ThermoParcel<ParcelType>::updateCellQuantities(td, dt, cellI);
pc_ = td.pInterp().interpolate(this->position(), cellI);
+ if (pc_ < SMALL)
+ {
+ WarningIn
+ (
+ "void Foam::ReactingParcel<ParcelType>::updateCellQuantities"
+ "("
+ "TrackData&, "
+ "const scalar, "
+ "const label"
+ ")"
+ ) << "Pressure < " << SMALL << " in cell " << cellI << nl << endl;
+ }
}
diff --git a/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C b/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C
index d8b999d0f58..b22a8e57e67 100644
--- a/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C
@@ -40,6 +40,19 @@ void Foam::ThermoParcel<ParcelType>::updateCellQuantities
KinematicParcel<ParcelType>::updateCellQuantities(td, dt, cellI);
Tc_ = td.TInterp().interpolate(this->position(), cellI);
+ if (Tc_ < SMALL)
+ {
+ WarningIn
+ (
+ "void Foam::ThermoParcel<ParcelType>::updateCellQuantities"
+ "("
+ "TrackData&, "
+ "const scalar, "
+ "const label"
+ ")"
+ ) << "Temperature < " << SMALL << " in cell " << cellI << nl << endl;
+ }
+
cpc_ = td.cpInterp().interpolate(this->position(), cellI);
}
--
GitLab