[cig-commits] r15497 - in seismo/3D/SPECFEM3D_GLOBE/trunk: . version41_beta/src

[savage13 at geodynamics.org]

Author: savage13
Date: 2009-07-31 18:31:47 -0700 (Fri, 31 Jul 2009)
New Revision: 15497

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/write_seismograms.F90
   seismo/3D/SPECFEM3D_GLOBE/trunk/write_seismograms.f90
Log:
Fixed bug when outputting seismograms to set the source location to the CMT location, not the PDE location and, if requested, to rotate to radial and tangential using the CMT location

Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/write_seismograms.F90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/write_seismograms.F90	2009-07-31 22:39:46 UTC (rev 15496)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/version41_beta/src/write_seismograms.F90	2009-08-01 01:31:47 UTC (rev 15497)
@@ -594,15 +594,19 @@
  STEL = stele(irec)
  STDP = undef    !stdep(irec)
 !event values (hypocenter):
- EVLA   = elat
- EVLO   = elon
+! EVLA   = elat
+! EVLO   = elon
+ EVLA   = cmt_lat
+ EVLO   = cmt_lon
  EVEL   = undef  !not defined
- EVDP   = depth
+ EVDP   = cmt_depth
 
 !cmt location values (different from hypocenter location, usually):
- USER0  = cmt_lat
- USER1  = cmt_lon
- USER2  = cmt_depth
+! USER0  = cmt_lat
+! USER1  = cmt_lon
+ USER0  = elat
+ USER1  = elon
+ USER2  = depth
 
  USER3  = cmt_hdur !half duration from CMT if not changed to hdur=0.d0 (point source)
 
@@ -678,9 +682,9 @@
  KCMPNM = chn(3:3)           ! 3A8
  KNETWK = network_name(irec) !  A6
 
- KUSER0 = 'CMT_LAT_'          !  A8
- KUSER1 = 'CMT_LON_'          !  A8
- KUSER2 = 'CMTDEPTH'          !  A8
+ KUSER0 = 'PDE_LAT_'          !  A8
+ KUSER1 = 'PDE_LON_'          !  A8
+ KUSER2 = 'PDEDEPTH'          !  A8
 !----------------------------------
 
   if (OUTPUT_SEISMOS_SAC_ALPHANUM) then

Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/write_seismograms.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/write_seismograms.f90	2009-07-31 22:39:46 UTC (rev 15496)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/write_seismograms.f90	2009-08-01 01:31:47 UTC (rev 15497)
@@ -418,7 +418,8 @@
 
           ! BS BS calculate backazimuth needed to rotate East and North
           ! components to Radial and Transverse components
-          call get_backazimuth(elat,elon,stlat(irec),stlon(irec),backaz)
+!          call get_backazimuth(elat,elon,stlat(irec),stlon(irec),backaz)
+          call get_backazimuth(cmt_lat,cmt_lon,stlat(irec),stlon(irec),backaz)
 
           if (backaz>180.) then
              phi=backaz-180.
@@ -529,15 +530,19 @@
  STEL = stele(irec)
  STDP = undef    !stdep(irec)
 !event values (hypocenter):
- EVLA   = elat
- EVLO   = elon
+! EVLA   = elat
+! EVLO   = elon
+ EVLA   = cmt_lat
+ EVLO   = cmt_lon
  EVEL   = undef  !not defined
- EVDP   = depth
+ EVDP   = cmt_depth
 
 !cmt location values (different from hypocenter location, usually):
- USER0  = cmt_lat
- USER1  = cmt_lon
- USER2  = cmt_depth
+! USER0  = cmt_lat
+! USER1  = cmt_lon
+ USER0  = elat
+ USER1  = elon
+ USER2  = depth
 
  USER3  = cmt_hdur !half duration from CMT if not changed to hdur=0.d0 (point source)
 
@@ -613,9 +618,9 @@
  KCMPNM = chn(3:3)           ! 3A8
  KNETWK = network_name(irec) !  A6
 
- KUSER0 = 'CMT_LAT_'          !  A8
- KUSER1 = 'CMT_LON_'          !  A8
- KUSER2 = 'CMTDEPTH'          !  A8
+ KUSER0 = 'PDE_LAT_'          !  A8
+ KUSER1 = 'PDE_LON_'          !  A8
+ KUSER2 = 'PDEDEPTH'          !  A8
 !----------------------------------
 
   if (OUTPUT_SEISMOS_SAC_ALPHANUM) then