[cig-commits] r15665 - seismo/2D/SPECFEM2D/trunk

[cmorency at geodynamics.org]

Author: cmorency
Date: 2009-09-11 07:15:49 -0700 (Fri, 11 Sep 2009)
New Revision: 15665

Modified:
   seismo/2D/SPECFEM2D/trunk/adj_seismogram.f90
   seismo/2D/SPECFEM2D/trunk/specfem2D.F90
Log:
Fixed bug on acoustic kernels calculation. kappa & density called when building the kernels had wrong indices. Thanks to Yang Luo for pointing it out.


Modified: seismo/2D/SPECFEM2D/trunk/adj_seismogram.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/adj_seismogram.f90	2009-09-09 23:16:34 UTC (rev 15664)
+++ seismo/2D/SPECFEM2D/trunk/adj_seismogram.f90	2009-09-11 14:15:49 UTC (rev 15665)
@@ -85,8 +85,8 @@
 !      Nnorm = deltat * sum(time_window(:) * seism_veloc(:) * seism_veloc(:))
 ! cross-correlation traveltime adjoint source
       if(abs(Nnorm) > EPS) then
-!      ft_bar(:) = - seism_veloc(:) * time_window(:) / Nnorm
-      ft_bar(:) = seism_veloc(:) * time_window(:) / Nnorm
+      ft_bar(:) = - seism_veloc(:) * time_window(:) / Nnorm
+!      ft_bar(:) = seism_veloc(:) * time_window(:) / Nnorm
       print*,'Norm =', Nnorm
       else
       print *, 'norm < EPS for file '

Modified: seismo/2D/SPECFEM2D/trunk/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/specfem2D.F90	2009-09-09 23:16:34 UTC (rev 15664)
+++ seismo/2D/SPECFEM2D/trunk/specfem2D.F90	2009-09-11 14:15:49 UTC (rev 15665)
@@ -6452,8 +6452,8 @@
       do k = 1, NGLLZ
           do i = 1, NGLLX
             iglob = ibool(i,k,ispec)
-    kappal_ac_global(iglob) = poroelastcoef(1,2,kmato(ispec))
-    rhol_ac_global(iglob) = density(2,kmato(ispec))
+    kappal_ac_global(iglob) = poroelastcoef(3,1,kmato(ispec))
+    rhol_ac_global(iglob) = density(1,kmato(ispec))
           enddo
       enddo
      endif