[cig-commits] r18114 - in seismo/3D/SPECFEM3D/trunk: examples/noise_tomography src/meshfem3D

[yangl at geodynamics.org]

Author: yangl
Date: 2011-03-15 09:37:40 -0700 (Tue, 15 Mar 2011)
New Revision: 18114

Modified:
   seismo/3D/SPECFEM3D/trunk/examples/noise_tomography/README
   seismo/3D/SPECFEM3D/trunk/src/meshfem3D/meshfem3D.f90
Log:
fix a bug in meshfem3D.f90: line 650--spacing_x_interface_top changed to spacing_y_interface_top

Modified: seismo/3D/SPECFEM3D/trunk/examples/noise_tomography/README
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/noise_tomography/README	2011-03-15 15:26:06 UTC (rev 18113)
+++ seismo/3D/SPECFEM3D/trunk/examples/noise_tomography/README	2011-03-15 16:37:40 UTC (rev 18114)
@@ -10,7 +10,7 @@
 
 
 **************
-0.README
+0. README
 **************
 
    This example shows how one can construct noise cross-correlation sensitivity kernels.
@@ -76,6 +76,6 @@
 4. Results
 **************
 
-   At the end of those simulations, kernel files will be saved in SPECFEM3D/in_out_files/NOISE_TOMOGRAPHY/1st/ & SPECFEM3D/in_out_files/NOISE_TOMOGRAPHY/1st/
+   At the end of those simulations, kernel files will be saved in SPECFEM3D/in_out_files/NOISE_TOMOGRAPHY/1st/ & SPECFEM3D/in_out_files/NOISE_TOMOGRAPHY/2nd/
    If mesh2vtu.pl works, you will have 1st_beta_kernel.vtu & 2nd_beta_kernel.vtu in SPECFEM3D/in_out_files/NOISE_TOMOGRAPHY/, which can be directly viewed in Paraview
    Once you visualize them in Paraview, you should see pictures similar to those in SPECFEM3D/examples/noise_tomography/REF_KERNELS/

Modified: seismo/3D/SPECFEM3D/trunk/src/meshfem3D/meshfem3D.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/meshfem3D/meshfem3D.f90	2011-03-15 15:26:06 UTC (rev 18113)
+++ seismo/3D/SPECFEM3D/trunk/src/meshfem3D/meshfem3D.f90	2011-03-15 16:37:40 UTC (rev 18114)
@@ -647,7 +647,7 @@
 
 ! get coordinate of corner in bathy/topo model
         icornerlong = int((long - orig_x_interface_top) / spacing_x_interface_top) + 1
-        icornerlat = int((lat - orig_y_interface_top) / spacing_x_interface_top) + 1
+        icornerlat = int((lat - orig_y_interface_top) / spacing_y_interface_top) + 1
 
 ! avoid edge effects and extend with identical point if outside model
         if(icornerlong < 1) icornerlong = 1