[cig-commits] r13013 - seismo/3D/SPECFEM3D/branches/update_temporary

[nlegoff at geodynamics.org]

Author: nlegoff
Date: 2008-10-09 08:57:54 -0700 (Thu, 09 Oct 2008)
New Revision: 13013

Modified:
   seismo/3D/SPECFEM3D/branches/update_temporary/constants.h.in
   seismo/3D/SPECFEM3D/branches/update_temporary/locate_receivers.f90
Log:
added choice if seismograms are recorded along normal to surface.

Modified: seismo/3D/SPECFEM3D/branches/update_temporary/constants.h.in
===================================================================
--- seismo/3D/SPECFEM3D/branches/update_temporary/constants.h.in	2008-10-09 15:33:09 UTC (rev 13012)
+++ seismo/3D/SPECFEM3D/branches/update_temporary/constants.h.in	2008-10-09 15:57:54 UTC (rev 13013)
@@ -130,6 +130,9 @@
   logical, parameter :: RECVS_CAN_BE_BURIED_EXT_MESH = .true.
   logical, parameter :: SOURCES_CAN_BE_BURIED_EXT_MESH = .true.
 
+! the seismograms are normal to surface
+  logical, parameter :: EXT_MESH_RECV_NORMAL = .false.
+
 !
   logical, parameter :: EXTERNAL_MESH_MOVIE_SURFACE = .false.
   logical, parameter :: EXTERNAL_MESH_CREATE_SHAKEMAP = .false.

Modified: seismo/3D/SPECFEM3D/branches/update_temporary/locate_receivers.f90
===================================================================
--- seismo/3D/SPECFEM3D/branches/update_temporary/locate_receivers.f90	2008-10-09 15:33:09 UTC (rev 13012)
+++ seismo/3D/SPECFEM3D/branches/update_temporary/locate_receivers.f90	2008-10-09 15:57:54 UTC (rev 13013)
@@ -478,6 +478,7 @@
     v_vector(3) = w_vector(1)*u_vector(2) - w_vector(2)*u_vector(1)
 
 ! build rotation matrice nu for seismograms
+    if (EXT_MESH_RECV_NORMAL) then
 !     East (u)
       nu(1,1,irec) = u_vector(1)
       nu(1,2,irec) = v_vector(1)
@@ -492,7 +493,23 @@
       nu(3,1,irec) = u_vector(3)
       nu(3,2,irec) = v_vector(3)
       nu(3,3,irec) = w_vector(3)
+      else
+!     East
+      nu(1,1,irec) = 1.d0
+      nu(1,2,irec) = 0.d0
+      nu(1,3,irec) = 0.d0
 
+!     North
+      nu(2,1,irec) = 0.d0
+      nu(2,2,irec) = 1.d0
+      nu(2,3,irec) = 0.d0
+
+!     Vertical
+      nu(3,1,irec) = 0.d0
+      nu(3,2,irec) = 0.d0
+      nu(3,3,irec) = 1.d0
+      endif
+
   endif ! of if (USE_EXTERNAL_MESH .and. (.not. RECVS_CAN_BE_BURIED_EXT_MESH))
 
 ! end of loop on all the stations