[cig-commits] r18131 - seismo/3D/SPECFEM3D/trunk/src/specfem3D

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2011-03-24 07:44:25 -0700 (Thu, 24 Mar 2011)
New Revision: 18131

Modified:
   seismo/3D/SPECFEM3D/trunk/src/specfem3D/read_mesh_databases.f90
Log:
fixed the unallocated adjoint arrays bug found by Camille Mazoyer-Kervazo


Modified: seismo/3D/SPECFEM3D/trunk/src/specfem3D/read_mesh_databases.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/specfem3D/read_mesh_databases.f90	2011-03-24 00:57:39 UTC (rev 18130)
+++ seismo/3D/SPECFEM3D/trunk/src/specfem3D/read_mesh_databases.f90	2011-03-24 14:44:25 UTC (rev 18131)
@@ -467,6 +467,7 @@
 
 ! allocates adjoint arrays for acoustic simulations
   if( ACOUSTIC_SIMULATION .and. SIMULATION_TYPE == 3 ) then
+
     ! backward potentials
     allocate(b_potential_acoustic(NGLOB_ADJOINT))
     allocate(b_potential_dot_acoustic(NGLOB_ADJOINT))
@@ -484,6 +485,25 @@
     allocate(b_buffer_send_scalar_ext_mesh(max_nibool_interfaces_ext_mesh,num_interfaces_ext_mesh))
     allocate(b_buffer_recv_scalar_ext_mesh(max_nibool_interfaces_ext_mesh,num_interfaces_ext_mesh))
 
+  else
+
+    ! backward potentials
+    allocate(b_potential_acoustic(1))
+    allocate(b_potential_dot_acoustic(1))
+    allocate(b_potential_dot_dot_acoustic(1))
+
+    ! kernels
+    allocate(rho_ac_kl(1,1,1,1))
+    allocate(rhop_ac_kl(1,1,1,1))
+    allocate(kappa_ac_kl(1,1,1,1))
+    allocate(alpha_ac_kl(1,1,1,1))
+
+    ! MPI handling
+    allocate(b_request_send_scalar_ext_mesh(1))
+    allocate(b_request_recv_scalar_ext_mesh(1))
+    allocate(b_buffer_send_scalar_ext_mesh(1,1))
+    allocate(b_buffer_recv_scalar_ext_mesh(1,1))
+
   endif
 
 ! ADJOINT moho