diff --git a/openfpm_numerics b/openfpm_numerics
index 515ec9bf8f1b482ad0b3d2c45ed9c36673f15a8e..d0f3b7819273efc599f5667de5c7c73fb9f5be3c 160000
--- a/openfpm_numerics
+++ b/openfpm_numerics
@@ -1 +1 @@
-Subproject commit 515ec9bf8f1b482ad0b3d2c45ed9c36673f15a8e
+Subproject commit d0f3b7819273efc599f5667de5c7c73fb9f5be3c
diff --git a/src/Grid/staggered_dist_grid.hpp b/src/Grid/staggered_dist_grid.hpp
index 7c0a7e875d752ce0f7b833512245bd29d032f2bd..c336da0b7fa6e40babaf2827c8d1eeb7377b0255 100644
--- a/src/Grid/staggered_dist_grid.hpp
+++ b/src/Grid/staggered_dist_grid.hpp
@@ -104,6 +104,23 @@ public:
 	:grid_dist_id<dim,St,T,Decomposition,Memory,device_grid>(g_sz,domain,ghost)
 	{}
 
+
+    /*! It construct a grid of a specified size, defined on a specified Box space, having a specified ghost size and periodicity
+     *
+     * \param g_sz grid size on each dimension
+     * \param domain Box that contain the grid
+     * \param g Ghost part of the domain (given in grid units)
+     * \param p periodicity
+     *
+     * \warning In very rare case the ghost part can be one point bigger than the one specified
+     *
+     */
+	staggered_grid_dist(const size_t (& g_sz)[dim],const Box<dim,St> & domain,
+			     const Ghost<dim,long int> & g, const periodicity<dim> & p)
+	:grid_dist_id<dim,St,T,Decomposition,Memory,device_grid>(g_sz,domain,g,p)
+	{
+	}
+
 	/*! \brief set the staggered positions of the properties
 	 *
 	 * \tparam property p