diff --git a/src/closure_mod.F90 b/src/closure_mod.F90
index 06eaf27633e371e355a1a1d09d40c0e857dc4fe1..8d57639cf5ca56ee86ea547a7f3bc873dc0a7236 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) .GT. dmax)
+                        .AND. (2*parray(ip)+jarray(ij) .LE. dmax)
       ENDDO
     ENDDO  
   CASE DEFAULT