[cig-commits] [commit] devel: Fixes bcasting starting address for 'addressing' array (d8730de)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/specfem3d_globe

On branch  : devel
Link       : https://github.com/geodynamics/specfem3d_globe/compare/b9fb1aa33196d161098710455fadbb4ed91c5e47...897de40783bd1a4630c2aacd3fa5f8b016d4c189

>---------------------------------------------------------------

commit d8730de16a48e0cc4e6dbe5cbc1fb160eb38debd
Author: Matthieu Lefebvre <ml15 at princeton.edu>
Date:   Fri Dec 5 09:40:54 2014 -0500

    Fixes bcasting starting address for 'addressing' array


>---------------------------------------------------------------

d8730de16a48e0cc4e6dbe5cbc1fb160eb38debd
 src/specfem3D/read_mesh_databases.F90 | 9 ++++++---
 1 file changed, 6 insertions(+), 3 deletions(-)

diff --git a/src/specfem3D/read_mesh_databases.F90 b/src/specfem3D/read_mesh_databases.F90
index 7929c17..f88e731 100644
--- a/src/specfem3D/read_mesh_databases.F90
+++ b/src/specfem3D/read_mesh_databases.F90
@@ -888,7 +888,7 @@
   endif ! I_should_read_the_database
 
   call bcast_all_i_for_database(ichunk, 1)
-  call bcast_all_i_for_database(addressing(1,1,1), size(addressing))
+  call bcast_all_i_for_database(addressing(1,0,0), size(addressing))
 
   end subroutine read_mesh_databases_addressing
 
@@ -918,6 +918,9 @@
     else  
       call read_mesh_databases_MPI_CM()
     endif
+    !call synchronize_all()
+    !print *, "read_mesh_databases_MPI_CM done"
+    !call flush()
   endif
   call bcast_mesh_databases_MPI_CM()
 
@@ -989,7 +992,6 @@
     if (ier /= 0 ) call exit_mpi(myrank,'Error allocating array b_buffer_send_vector_inner_core etc.')
   endif
 
-
   ! user output
   if (myrank == 0) then
     write(IMAIN,*) '  for overlapping of communications with calculations:'
@@ -1687,7 +1689,6 @@
     if (ier /= 0 ) &
       call exit_mpi(myrank,'Error allocating array my_neighbours_crust_mantle etc.')
   endif
-
   if (num_interfaces_crust_mantle > 0) then
     call bcast_all_i_for_database(max_nibool_interfaces_cm, 1)
     if (.not. I_should_read_the_database) then 
@@ -1784,7 +1785,9 @@
               stat=ier)
       if (ier /= 0 ) call exit_mpi(myrank,'Error allocating array ibool_interfaces_outer_core')
     endif
+  endif
 
+  if (num_interfaces_outer_core > 0) then
     call bcast_all_i_for_database(my_neighbours_outer_core(1), size(my_neighbours_outer_core))
     call bcast_all_i_for_database(nibool_interfaces_outer_core(1), size(nibool_interfaces_outer_core))
     call bcast_all_i_for_database(ibool_interfaces_outer_core(1,1), size(ibool_interfaces_outer_core))