diff --git a/src/TNL/Experimental/Hamilton-Jacobi/Solvers/hamilton-jacobi/HamiltonJacobiProblem.h b/src/TNL/Experimental/Hamilton-Jacobi/Solvers/hamilton-jacobi/HamiltonJacobiProblem.h
index a4144200028fccf124aad1730b2415fa1cd9e5ff..7f1bd4193960ba3e6bbc8f2cb1fb12f7c4e84cdf 100644
--- a/src/TNL/Experimental/Hamilton-Jacobi/Solvers/hamilton-jacobi/HamiltonJacobiProblem.h
+++ b/src/TNL/Experimental/Hamilton-Jacobi/Solvers/hamilton-jacobi/HamiltonJacobiProblem.h
@@ -25,6 +25,8 @@
 #include <solvers/pde/tnlLinearSystemAssembler.h>
 #include <functions/tnlMeshFunction.h>
 
+#include <TNL/Meshes/DistributedMeshes/DistributedGridIO.h>
+
 template< typename Mesh,
 		    typename DifferentialOperator,
 		    typename BoundaryCondition,
diff --git a/src/TNL/Experimental/Hamilton-Jacobi/Solvers/hamilton-jacobi/HamiltonJacobiProblem_impl.h b/src/TNL/Experimental/Hamilton-Jacobi/Solvers/hamilton-jacobi/HamiltonJacobiProblem_impl.h
index 3cc6388497c062865a812cca629c57d45d715d09..9244b18332e40887c131fbb68cabbffe417bfef2 100644
--- a/src/TNL/Experimental/Hamilton-Jacobi/Solvers/hamilton-jacobi/HamiltonJacobiProblem_impl.h
+++ b/src/TNL/Experimental/Hamilton-Jacobi/Solvers/hamilton-jacobi/HamiltonJacobiProblem_impl.h
@@ -123,12 +123,24 @@ setInitialCondition( const Config::ParameterContainer& parameters,
 {
   this->bindDofs( mesh, dofs );
   const String& initialConditionFile = parameters.getParameter< String >( "initial-condition" );
-  if( ! this->solution.boundLoad( initialConditionFile ) )
+  if(CommunicatorType::isDistributed())
   {
-    std::cerr << "I am not able to load the initial condition from the file " << initialConditionFile << "." <<std::endl;
-    return false;
+    std::cout<<"Nodes Distribution: " << uPointer->getMesh().getDistributedMesh()->printProcessDistr() << std::endl;
+    if(distributedIOType==Meshes::DistributedMeshes::MpiIO)
+      Meshes::DistributedMeshes::DistributedGridIO<MeshFunctionType,Meshes::DistributedMeshes::MpiIO> ::load(initialConditionFile, *uPointer );
+    if(distributedIOType==Meshes::DistributedMeshes::LocalCopy)
+      Meshes::DistributedMeshes::DistributedGridIO<MeshFunctionType,Meshes::DistributedMeshes::LocalCopy> ::load(initialConditionFile, *uPointer );
+    uPointer->template synchronize<CommunicatorType>();
   }
-  return true;
+  else
+  {
+    if( ! this->solution.boundLoad( initialConditionFile ) )
+    {
+      std::cerr << "I am not able to load the initial condition from the file " << initialConditionFile << "." <<std::endl;
+      return false;
+    }
+  }
+   return true;
 }
 
 template< typename Mesh,