Print negative values of alpha

2024-10-15 11:41:31 -04:00 · 2024-10-15 11:41:31 -04:00 · 0404bd0cf0
commit 0404bd0cf0
parent 7675148e7d
1 changed files with 6 additions and 0 deletions
--- 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;