diff --git a/src/closure_mod.F90 b/src/closure_mod.F90
index 9e18d7040b3851afdf07e1b9c881a4cbb1697cb1..06eaf27633e371e355a1a1d09d40c0e857dc4fe1 100644
--- a/src/closure_mod.F90
+++ b/src/closure_mod.F90
@@ -45,7 +45,7 @@ SUBROUTINE set_closure_model
   CASE('max_degree')
     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)) &
+          evolve_mom(ip,ij) = (parray(ip).GE.0) .AND. (jarray(ij).GE.0) &
                         .AND. (2*parray(ip)+jarray(ij) .GT. dmax)
       ENDDO
     ENDDO