diff --git a/src/thermophysicalModels/chemistryModel/chemistryModel/psiChemistryModel/newPsiChemistryModel.C b/src/thermophysicalModels/chemistryModel/chemistryModel/psiChemistryModel/newPsiChemistryModel.C
index 5195801f25..6a5606d0e7 100644
--- a/src/thermophysicalModels/chemistryModel/chemistryModel/psiChemistryModel/newPsiChemistryModel.C
+++ b/src/thermophysicalModels/chemistryModel/chemistryModel/psiChemistryModel/newPsiChemistryModel.C
@@ -35,7 +35,7 @@ Foam::autoPtr<Foam::psiChemistryModel> Foam::psiChemistryModel::New
 {
     word psiChemistryModelType;
     word thermoTypeName;
-    word userSel;
+    word userModel;
 
     // Enclose the creation of the chemistrtyProperties to ensure it is
     // deleted before the chemistrtyProperties is created otherwise the
@@ -53,37 +53,54 @@ Foam::autoPtr<Foam::psiChemistryModel> Foam::psiChemistryModel::New
             )
         );
 
-        chemistryPropertiesDict.lookup("psiChemistryModel") >> userSel;
+        chemistryPropertiesDict.lookup("psiChemistryModel") >> userModel;
 
         // construct chemistry model type name by inserting first template
         // argument
-        label tempOpen = userSel.find('<');
-        label tempClose = userSel.find('>');
+        label tempOpen = userModel.find('<');
+        label tempClose = userModel.find('>');
 
-        word className = userSel(0, tempOpen);
-        thermoTypeName = userSel(tempOpen + 1, tempClose - tempOpen - 1);
+        word className = userModel(0, tempOpen);
+        thermoTypeName = userModel(tempOpen + 1, tempClose - tempOpen - 1);
 
         psiChemistryModelType =
             className + '<' + typeName + ',' + thermoTypeName + '>';
     }
 
-    Info<< "Selecting psiChemistryModel " << userSel << endl;
+    if (debug)
+    {
+        Info<< "Selecting psiChemistryModel " << psiChemistryModelType << endl;
+    }
+    else
+    {
+        Info<< "Selecting psiChemistryModel " << userModel << endl;
+    }
 
     fvMeshConstructorTable::iterator cstrIter =
         fvMeshConstructorTablePtr_->find(psiChemistryModelType);
 
     if (cstrIter == fvMeshConstructorTablePtr_->end())
     {
-        wordList models = fvMeshConstructorTablePtr_->toc();
-        forAll(models, i)
+        if (debug)
         {
-            models[i] = models[i].replace(typeName + ',', "");
+            FatalErrorIn("psiChemistryModelBase::New(const mesh&)")
+                << "Unknown psiChemistryModel type " << psiChemistryModelType
+                << nl << nl << "Valid psiChemistryModel types are:" << nl
+                << fvMeshConstructorTablePtr_->toc() << nl << exit(FatalError);
         }
+        else
+        {
+            wordList models = fvMeshConstructorTablePtr_->toc();
+            forAll(models, i)
+            {
+                models[i] = models[i].replace(typeName + ',', "");
+            }
 
-        FatalErrorIn("psiChemistryModelBase::New(const mesh&)")
-            << "Unknown psiChemistryModel type " << userSel
-            << nl << nl << "Valid psiChemistryModel types are:" << nl
-            << models << nl << exit(FatalError);
+            FatalErrorIn("psiChemistryModelBase::New(const mesh&)")
+                << "Unknown psiChemistryModel type " << userModel
+                << nl << nl << "Valid psiChemistryModel types are:" << nl
+                << models << nl << exit(FatalError);
+        }
     }
 
     return autoPtr<psiChemistryModel>
diff --git a/src/thermophysicalModels/chemistryModel/chemistryModel/rhoChemistryModel/newRhoChemistryModel.C b/src/thermophysicalModels/chemistryModel/chemistryModel/rhoChemistryModel/newRhoChemistryModel.C
index 938b841221..15a7fbe9bf 100644
--- a/src/thermophysicalModels/chemistryModel/chemistryModel/rhoChemistryModel/newRhoChemistryModel.C
+++ b/src/thermophysicalModels/chemistryModel/chemistryModel/rhoChemistryModel/newRhoChemistryModel.C
@@ -35,7 +35,7 @@ Foam::autoPtr<Foam::rhoChemistryModel> Foam::rhoChemistryModel::New
 {
     word rhoChemistryModelType;
     word thermoTypeName;
-    word userSel;
+    word userModel;
 
     // Enclose the creation of the chemistrtyProperties to ensure it is
     // deleted before the chemistrtyProperties is created otherwise the
@@ -53,37 +53,54 @@ Foam::autoPtr<Foam::rhoChemistryModel> Foam::rhoChemistryModel::New
             )
         );
 
-        chemistryPropertiesDict.lookup("rhoChemistryModelType") >> userSel;
+        chemistryPropertiesDict.lookup("rhoChemistryModelType") >> userModel;
 
         // construct chemistry model type name by inserting first template
         // argument
-        label tempOpen = userSel.find('<');
-        label tempClose = userSel.find('>');
+        label tempOpen = userModel.find('<');
+        label tempClose = userModel.find('>');
 
-        word className = userSel(0, tempOpen);
-        thermoTypeName = userSel(tempOpen + 1, tempClose - tempOpen - 1);
+        word className = userModel(0, tempOpen);
+        thermoTypeName = userModel(tempOpen + 1, tempClose - tempOpen - 1);
 
         rhoChemistryModelType =
             className + '<' + typeName + ',' + thermoTypeName + '>';
     }
 
-    Info<< "Selecting rhoChemistryModel " << userSel << endl;
+    if (debug)
+    {
+        Info<< "Selecting rhoChemistryModel " << rhoChemistryModelType << endl;
+    }
+    else
+    {
+        Info<< "Selecting rhoChemistryModel " << userModel << endl;
+    }
 
     fvMeshConstructorTable::iterator cstrIter =
         fvMeshConstructorTablePtr_->find(rhoChemistryModelType);
 
     if (cstrIter == fvMeshConstructorTablePtr_->end())
     {
-        wordList models = fvMeshConstructorTablePtr_->toc();
-        forAll(models, i)
+        if (debug)
         {
-            models[i] = models[i].replace(typeName + ',', "");
+            FatalErrorIn("rhoChemistryModelBase::New(const mesh&)")
+                << "Unknown rhoChemistryModel type " << rhoChemistryModelType
+                << nl << nl << "Valid rhoChemistryModel types are:" << nl
+                << fvMeshConstructorTablePtr_->toc() << nl << exit(FatalError);
         }
+        else
+        {
+            wordList models = fvMeshConstructorTablePtr_->toc();
+            forAll(models, i)
+            {
+                models[i] = models[i].replace(typeName + ',', "");
+            }
 
-        FatalErrorIn("rhoChemistryModelBase::New(const mesh&)")
-            << "Unknown rhoChemistryModel type " << userSel
-            << nl << nl << "Valid rhoChemistryModel types are:" << nl
-            << models << nl << exit(FatalError);
+            FatalErrorIn("rhoChemistryModelBase::New(const mesh&)")
+                << "Unknown rhoChemistryModel type " << userModel
+                << nl << nl << "Valid rhoChemistryModel types are:" << nl
+                << models << nl << exit(FatalError);
+        }
     }
 
     return autoPtr<rhoChemistryModel>