diff --git a/blast/src/DBoundaryLayer.cpp b/blast/src/DBoundaryLayer.cpp
index 9886a798f26989f6efee41f6e8daa49668a7059e..804d7802f68f5e19e4805c1f297aab75e816c27e 100644
--- a/blast/src/DBoundaryLayer.cpp
+++ b/blast/src/DBoundaryLayer.cpp
@@ -157,6 +157,10 @@ void BoundaryLayer::computeBlwVel()
         blowingVelocity[iPoint - 1] = (blEdge->getRhoe(iPoint) * blEdge->getVt(iPoint) * deltaStar[iPoint] - blEdge->getRhoe(iPoint - 1) * blEdge->getVt(iPoint - 1) * deltaStar[iPoint - 1]) / (blEdge->getRhoe(iPoint) * (mesh->getLoc(iPoint) - mesh->getLoc(iPoint - 1)));
         if (std::isnan(blowingVelocity[iPoint - 1]))
         {
+            if (std::isnan(blEdge->getVt(iPoint)))
+                std::cout << "Velocity is nan at position " << iPoint << std::endl;
+            if (std::isnan(blEdge->getVt(iPoint-1)))
+                std::cout << "Velocity is nan at position " << iPoint-1 << std::endl;
             if (blEdge->getRhoe(iPoint) == 0.)
                 std::cout << "Density is zero at point " << iPoint << std::endl;
             if ((mesh->getLoc(iPoint) - mesh->getLoc(iPoint - 1)) == 0.)