From bb7730bfd873feb6cae8b3c8eb5bc1d37623188b Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Tue, 4 Apr 2023 10:55:19 +0200
Subject: [PATCH] typo, closure by max_degree works now

---
 src/closure_mod.F90 | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/closure_mod.F90 b/src/closure_mod.F90
index 8d57639c..9a80a8bd 100644
--- a/src/closure_mod.F90
+++ b/src/closure_mod.F90
@@ -46,7 +46,7 @@ SUBROUTINE set_closure_model
     DO ip = 1,local_np+ngp
       DO ij = 1, local_nj+ngj
           evolve_mom(ip,ij) = (parray(ip).GE.0) .AND. (jarray(ij).GE.0) &
-                        .AND. (2*parray(ip)+jarray(ij) .LE. dmax)
+                        .AND. (parray(ip)+2*jarray(ij) .LE. dmax)
       ENDDO
     ENDDO  
   CASE DEFAULT
-- 
GitLab