diff --git a/mlprec/impl/solver/mld_c_mumps_solver_apply.F90 b/mlprec/impl/solver/mld_c_mumps_solver_apply.F90
index 7424e6aa..f489dbea 100644
--- a/mlprec/impl/solver/mld_c_mumps_solver_apply.F90
+++ b/mlprec/impl/solver/mld_c_mumps_solver_apply.F90
@@ -136,7 +136,7 @@ subroutine c_mumps_solver_apply(alpha,sv,x,beta,y,desc_data,&
   if (sv%ipar(1) == mld_local_solver_ ) then
     call psb_geaxpby(alpha,gx,beta,y,desc_data,info)
   else
-    call psb_scatter(gx, ww, desc_data, info, root=0)   
+    call psb_scatter(gx, ww, desc_data, info, root=izero)   
     if (info == psb_success_) then
       call psb_geaxpby(alpha,ww,beta,y,desc_data,info)
     end if
diff --git a/mlprec/impl/solver/mld_d_mumps_solver_apply.F90 b/mlprec/impl/solver/mld_d_mumps_solver_apply.F90
index 8b77ae6b..b6b704fe 100644
--- a/mlprec/impl/solver/mld_d_mumps_solver_apply.F90
+++ b/mlprec/impl/solver/mld_d_mumps_solver_apply.F90
@@ -136,7 +136,7 @@ subroutine d_mumps_solver_apply(alpha,sv,x,beta,y,desc_data,&
   if (sv%ipar(1) == mld_local_solver_ ) then
     call psb_geaxpby(alpha,gx,beta,y,desc_data,info)
   else
-    call psb_scatter(gx, ww, desc_data, info, root=0)   
+    call psb_scatter(gx, ww, desc_data, info, root=izero)   
     if (info == psb_success_) then
       call psb_geaxpby(alpha,ww,beta,y,desc_data,info)
     end if
diff --git a/mlprec/impl/solver/mld_s_mumps_solver_apply.F90 b/mlprec/impl/solver/mld_s_mumps_solver_apply.F90
index 01ba296d..e252c41f 100644
--- a/mlprec/impl/solver/mld_s_mumps_solver_apply.F90
+++ b/mlprec/impl/solver/mld_s_mumps_solver_apply.F90
@@ -136,7 +136,7 @@ subroutine s_mumps_solver_apply(alpha,sv,x,beta,y,desc_data,&
   if (sv%ipar(1) == mld_local_solver_ ) then
     call psb_geaxpby(alpha,gx,beta,y,desc_data,info)
   else
-    call psb_scatter(gx, ww, desc_data, info, root=0)   
+    call psb_scatter(gx, ww, desc_data, info, root=izero)   
     if (info == psb_success_) then
       call psb_geaxpby(alpha,ww,beta,y,desc_data,info)
     end if
diff --git a/mlprec/impl/solver/mld_z_mumps_solver_apply.F90 b/mlprec/impl/solver/mld_z_mumps_solver_apply.F90
index d39517f7..970be263 100644
--- a/mlprec/impl/solver/mld_z_mumps_solver_apply.F90
+++ b/mlprec/impl/solver/mld_z_mumps_solver_apply.F90
@@ -136,7 +136,7 @@ subroutine z_mumps_solver_apply(alpha,sv,x,beta,y,desc_data,&
   if (sv%ipar(1) == mld_local_solver_ ) then
     call psb_geaxpby(alpha,gx,beta,y,desc_data,info)
   else
-    call psb_scatter(gx, ww, desc_data, info, root=0)   
+    call psb_scatter(gx, ww, desc_data, info, root=izero)   
     if (info == psb_success_) then
       call psb_geaxpby(alpha,ww,beta,y,desc_data,info)
     end if