Skip to content
This repository was archived by the owner on Jun 7, 2024. It is now read-only.

Commit 46f0d22

Browse files
committed
Debug assistance in JDAS when excluding reactions from the ODE set.
This helped us figure out when PDep reactions are or are not excluded from the ODE solver because they duplicate existing reactions.
1 parent 0f824ba commit 46f0d22

File tree

1 file changed

+19
-5
lines changed

1 file changed

+19
-5
lines changed

source/RMG/jing/rxnSys/JDAS.java

Lines changed: 19 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -279,11 +279,25 @@ public void generatePDepReactionList(ReactionModel p_reactionModel,
279279
// check if this reaction is already in the list and also
280280
// check if this reaction has a reverse reaction which is already present in the list.
281281
if (rxn.reactantEqualsProduct()) continue;
282-
if (troeList.contains(rxn) || troeList.contains(reverse)) continue;
283-
if (thirdBodyList.contains(rxn) || thirdBodyList.contains(reverse)) continue;
284-
if (lindemannList.contains(rxn) || lindemannList.contains(reverse)) continue;
285-
if (seedList.contains(rxn) || seedList.contains(reverse)) continue; // exclude rxns already in seed mechanism
286-
282+
if (troeList.contains(rxn) || troeList.contains(reverse)) {
283+
//Logger.debug(String.format("Excluding FAME-estimated PDep rate for %s from ODEs because its Troe rate is in a reaction library or seed mechanism.",rxn));
284+
continue; // exclude rxns already in seed mechanism
285+
}
286+
else if (thirdBodyList.contains(rxn) || thirdBodyList.contains(reverse)) {
287+
//Logger.debug(String.format("Excluding FAME-estimated PDep rate for %s from ODEs because its 3-body rate is in a reaction library or seed mechanism.",rxn));
288+
continue; // exclude rxns already in seed mechanism
289+
}
290+
else if (lindemannList.contains(rxn) || lindemannList.contains(reverse)) {
291+
//Logger.debug(String.format("Excluding FAME-estimated PDep rate for %s from ODEs because its Lindemann rate is in a reaction library or seed mechanism.",rxn));
292+
continue; // exclude rxns already in seed mechanism
293+
}
294+
else if (seedList.contains(rxn) || seedList.contains(reverse)) {
295+
//Logger.debug(String.format("Excluding FAME-estimated PDep rate for %s from ODEs because it's in the seed mechanism",rxn));
296+
continue; // exclude rxns already in seed mechanism
297+
}
298+
else {
299+
//Logger.debug(String.format("Including FAME-estimated PDep rate for %s in ODEs because it's not in the seed mechanism, nor does it have a P-dep rate from a reaction library.",rxn));
300+
}
287301
if (!pDepList.contains(rxn) && !pDepList.contains(reverse)) {
288302
pDepList.add(rxn);
289303
} else if (pDepList.contains(rxn) && !pDepList.contains(reverse)) {

0 commit comments

Comments
 (0)