diff --git a/src/closure_mod.F90 b/src/closure_mod.F90
index 8d57639cf5ca56ee86ea547a7f3bc873dc0a7236..9a80a8bdb57cd8bc8c2cfef5e0da1510efd3095f 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