[cig-commits] r15487 - in seismo/3D/SPECFEM3D_GLOBE/tags: . v4.0.4

[leif at geodynamics.org]

Author: leif
Date: 2009-07-30 09:49:54 -0700 (Thu, 30 Jul 2009)
New Revision: 15487

Added:
   seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.4/
Modified:
   seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.4/write_seismograms.f90
Log:
Created SPECFEM3D GLOBE v4.0.4 tag, which incorporates a
SAC/rotated-output event location bug from Brian Savage.


Copied: seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.4 (from rev 15486, seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.3)

Modified: seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.4/write_seismograms.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.3/write_seismograms.f90	2009-07-30 14:23:31 UTC (rev 15486)
+++ seismo/3D/SPECFEM3D_GLOBE/tags/v4.0.4/write_seismograms.f90	2009-07-30 16:49:54 UTC (rev 15487)
@@ -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