diff --git a/tutorials/sphereGranFlow/toteblender/settings/geometryDict b/tutorials/sphereGranFlow/toteblender/settings/geometryDict
new file mode 100644
index 00000000..49d486bd
--- /dev/null
+++ b/tutorials/sphereGranFlow/toteblender/settings/geometryDict
@@ -0,0 +1,151 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+objectName 	geometryDict;
+objectType 	dictionary;
+/* ------------------------------------------------------------------------- */
+// motion model: rotating object around an axis
+motionModel rotatingAxisMotion; 
+
+surfaces
+{
+	enterGate
+	{
+	// type of wall
+		type planeWall;
+	// coords of wall
+		p1 (-0.05    -0.05    0.3);
+		p2 (-0.05    0.05     0.3);
+		p3 ( 0.05    0.05     0.3);
+		p4 (0.05     -0.05    0.3);
+	// material of wall
+		material prop1;
+	// motion component name
+		motion rotAxis;	
+	}
+	
+    cylinderinlet
+	{
+	// type of the wall
+		type 		cylinderWall;
+	// begin point of cylinder axis 	
+		p1 			(0.0 0.0 0.28);
+	// end point of cylinder axis  
+		p2 			(0.0 0.0 0.3);
+	// radius at p1  
+		radius1 	0.03;
+	// radius at p2		
+		radius2 	0.03;
+	// number of divisions		
+		resolution 	36;
+	// material name of this wall
+		material 	prop1;
+	// motion component name   	
+		motion rotAxis;		
+	}
+
+	coneShelltop
+	{
+	// type of the wall	
+		type 		cylinderWall;
+	// begin point of cylinder axis  	
+		p1 			(0.0 0.0 0.2);
+	// end point of cylinder axis  
+		p2 			(0.0 0.0 0.28);
+	// radius at p1  
+		radius1 	0.1;
+	// radius at p2		
+		radius2 	0.03;
+	// number of divisions		
+		resolution 	36;
+	// material name of this wall      	
+		material 	prop1;
+	// motion component name   	
+		motion rotAxis;		
+	}
+
+	cylinderShell
+	{
+	// type of the wall
+		type 		cylinderWall;
+	// begin point of cylinder axis  	
+		p1 			(0.0 0.0 0.1);
+	// end point of cylinder axis
+		p2 			(0.0 0.0 0.2);
+	// radius at p1	  
+		radius1 	0.1;
+	// radius at p2			
+		radius2 	0.1;
+	// number of divisions		
+		resolution 	36;
+	// material name of this wall	      	
+		material 	prop1; 
+	// motion component name  	
+		motion rotAxis;		
+	}
+
+	coneShelldown
+	{
+	// type of the wall
+		type 		cylinderWall;
+	// begin point of cylinder axis  	
+		p1 			(0.0 0.0 0.02);
+	// end point of cylinder axis  
+		p2 			(0.0 0.0 0.1);
+	// radius at p1  
+		radius1 	0.03;
+	// radius at p2		
+		radius2 	0.1;
+	// number of divisions		
+		resolution 	36;
+	// material name of this wall	      	
+		material 	prop1;
+	// motion component name   	
+		motion rotAxis;		
+	}
+	/*
+	This is a plane wall at the exit of silo
+	*/
+
+	    cylinderoutlet
+	{
+	// type of the wall
+		type 		cylinderWall;  	
+	// begin point of cylinder axis	
+		p1 			(0.0 0.0 0.0);
+	// end point of cylinder axis	  
+		p2 			(0.0 0.0 0.02);
+	// radius at p1  
+		radius1 	0.03;
+	// radius at p2			
+		radius2 	0.03;
+	// number of divisions			
+		resolution 	36;
+	// material name of this wall	      	
+		material 	prop1;
+	// motion component name	   	
+		motion rotAxis;		
+	}
+	exitGate
+	{
+		type planeWall;
+		p1 (-0.05    -0.05    0);
+		p2 (-0.05    0.05     0);
+		p3 ( 0.05    0.05     0);
+		p4 (0.05     -0.05    0);
+		material prop1;
+		motion rotAxis;		
+	}
+		
+}
+// information for rotatingAxisMotion motion model 
+rotatingAxisMotionInfo
+{
+	rotAxis 
+	{
+		p1 (-0.1 0.0 0.15);	// first point for the axis of rotation 
+		p2 (0.1 0.0 0.15);	// second point for the axis of rotation
+		omega 3; 		// rotation speed (rad/s)
+	}
+}
diff --git a/tutorials/sphereGranFlow/toteblender/settings/particlesDict b/tutorials/sphereGranFlow/toteblender/settings/particlesDict
new file mode 100644
index 00000000..49a26e76
--- /dev/null
+++ b/tutorials/sphereGranFlow/toteblender/settings/particlesDict
@@ -0,0 +1,59 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+objectName 	particlesDict;
+objectType 	dictionary;
+/* ------------------------------------------------------------------------- */
+setFields
+{
+	/*
+	    Default value for fields defined for particles
+	These fields should always be defined for simulations with 
+	spherical particles.
+	*/
+
+	defaultValue 
+	{
+// linear velocity (m/s)	
+		velocity 		realx3 	(0 0 0); 
+// linear acceleration (m/s2)		
+		acceleration 	realx3 	(0 0 0);
+// rotational velocity (rad/s)		 
+		rotVelocity 	realx3 	(0 0 0);
+// name of the particle shape 		 
+		shapeName 		word	sphere1; 
+	}
+
+	selectors
+	{}
+}
+
+// positions particles 
+positionParticles
+{
+// ordered positioning
+	method positionOrdered;     
+// maximum number of particles in the simulation
+	maxNumberOfParticles 40000;
+// perform initial sorting based on morton code? 
+	mortonSorting Yes;             
+// box for positioning particles 
+	box  
+	{
+// lower corner point of the box 	
+		min (-0.06 -0.06 0.08);
+// upper corner point of the box 
+		max (0.06 0.06 0.18);
+	}
+
+	positionOrderedInfo
+	{
+// minimum space between centers of particles
+		diameter 0.004;
+// number of particles in the simulation 	 	
+		numPoints 20000;
+// axis order for filling the space with particles		 	
+		axisOrder (z y x);  
+	}
+}
diff --git a/tutorials/sphereGranFlow/toteblender/settings/settingsDict b/tutorials/sphereGranFlow/toteblender/settings/settingsDict
new file mode 100644
index 00000000..28b9a552
--- /dev/null
+++ b/tutorials/sphereGranFlow/toteblender/settings/settingsDict
@@ -0,0 +1,34 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+objectName 	settingsDict;
+objectType 	dictionary;;
+/*---------------------------------------------------------------------------*/
+run 			toteBlender;
+// time step for integration (s)
+dt 				0.00001; 	
+// start time for simulation 
+startTime 		0; 			
+// end time for simulation 
+endTime 		10; 		
+// time interval for saving the simulation
+saveInterval 	0.1; 		
+// maximum number of digits for time folder
+timePrecision   6;			 
+// gravity vector (m/s2)
+g 				(0 0 -9.8);  
+
+/* Simulation domain */
+/*	every particles that goes outside this domain is deleted. */
+domain 
+{
+	min (-0.5 -0.5 -0.5);
+	max (0.5 0.5 0.5);
+}
+// integration method
+integrationMethod 		AdamsBashforth2; 	 
+// report timers?
+timersReport 			Yes;  	
+// time interval for reporting timers
+timersReportInterval   	0.01;