[cig-commits] r13805 - seismo/3D/SPECFEM3D_SESAME/trunk

[nlegoff at geodynamics.org]

Author: nlegoff
Date: 2009-01-06 06:54:21 -0800 (Tue, 06 Jan 2009)
New Revision: 13805

Modified:
   seismo/3D/SPECFEM3D_SESAME/trunk/create_movie_AVS_DX.f90
Log:
fixed a bug in conversion of shakemaps to DX format (create_movie_AVS_DX.f90) for lowres shakemaps. Shakemap generation seems to be buggy, investigating.

Modified: seismo/3D/SPECFEM3D_SESAME/trunk/create_movie_AVS_DX.f90
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/trunk/create_movie_AVS_DX.f90	2009-01-06 08:43:23 UTC (rev 13804)
+++ seismo/3D/SPECFEM3D_SESAME/trunk/create_movie_AVS_DX.f90	2009-01-06 14:54:21 UTC (rev 13805)
@@ -466,12 +466,12 @@
 ! for shaking map, norm of U stored in ux, V in uy and A in uz
           if(plot_shaking_map) then
 !!!! NL NL mute value near source
-              if ( (sqrt(((x(i,j) - (X_SOURCE_EXT_MESH))**2 + &
-                   (y(i,j) - (Y_SOURCE_EXT_MESH))**2 + &
-                   (z(i,j) - (Z_SOURCE_EXT_MESH))**2)) < RADIUS_TO_MUTE) &
+              if ( (sqrt(((dble(xcoord) - (X_SOURCE_EXT_MESH))**2 + &
+                   (dble(ycoord) - (Y_SOURCE_EXT_MESH))**2 + &
+                   (dble(zcoord) - (Z_SOURCE_EXT_MESH))**2)) < RADIUS_TO_MUTE) &
                    .and. MUTE_SOURCE) then
 
-                display(i,j) = 0.
+                field_display(ilocnum+ieoff) = 0. 
               else