From 0404bd0cf06397560ba95ed225404a513b26c27b Mon Sep 17 00:00:00 2001
From: Ian Jauslin <ian.jauslin@rutgers.edu>
Date: Tue, 15 Oct 2024 11:41:31 -0400
Subject: [PATCH] Print negative values of alpha

---
 src/navier-stokes.c | 6 ++++++
 1 file changed, 6 insertions(+)

diff --git a/src/navier-stokes.c b/src/navier-stokes.c
index 3ae4ece..7b74a5f 100644
--- a/src/navier-stokes.c
+++ b/src/navier-stokes.c
@@ -235,6 +235,12 @@ int enstrophy(
       prevtime=time;
     }
 
+    // print alpha when it gets negative
+    if(alpha<0){
+      fprintf(stderr,"## negative alpha: % .8e at time % .8e\n",alpha, time);
+      printf("## negative alpha: % .15e at time % .15e\n",alpha, time);
+    }
+
     // get ready for next step
     step=next_step;