diff --git a/src/opt_slip_plane.f90 b/src/opt_slip_plane.f90 index e607ae8..80da930 100644 --- a/src/opt_slip_plane.f90 +++ b/src/opt_slip_plane.f90 @@ -86,7 +86,7 @@ module opt_slip_plane end do !Check to make sure all lattice points exist for the current element - if(any(lat_points(m:m+esize-1, n:n+esize-1, o:o+esize-1) == 0)) cycle latloop + if(any(.not.lat_points(m:m+esize-1, n:n+esize-1, o:o+esize-1))) cycle latloop !Check to see if the plane intersects this element if not then add it maxp = maxval(rfill(sdim,1:basisnum(lat_ele(ie)),1:ng_node(lat_ele(ie)))) minp = minval(rfill(sdim,1:basisnum(lat_ele(ie)),1:ng_node(lat_ele(ie))))