[cig-commits] r16452 - seismo/3D/SPECFEM3D_GLOBE/tags/v5.0.1

[danielpeter at geodynamics.org]

Author: danielpeter
Date: 2010-03-24 12:22:49 -0700 (Wed, 24 Mar 2010)
New Revision: 16452

Modified:
   seismo/3D/SPECFEM3D_GLOBE/tags/v5.0.1/create_movie_GMT_global.f90
   seismo/3D/SPECFEM3D_GLOBE/tags/v5.0.1/write_seismograms.f90
Log:
copying changes for tagged version: updating version number in SAC headers and create_movie_GMT_global.f90

Modified: seismo/3D/SPECFEM3D_GLOBE/tags/v5.0.1/create_movie_GMT_global.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/tags/v5.0.1/create_movie_GMT_global.f90	2010-03-24 19:02:13 UTC (rev 16451)
+++ seismo/3D/SPECFEM3D_GLOBE/tags/v5.0.1/create_movie_GMT_global.f90	2010-03-24 19:22:49 UTC (rev 16452)
@@ -56,8 +56,8 @@
 
   ! muting source region
   logical, parameter :: MUTE_SOURCE = .true.
-  real(kind=CUSTOM_REAL) :: RADIUS_TO_MUTE = 10.0    ! start radius in degrees  
-  real(kind=CUSTOM_REAL) :: STARTTIME_TO_MUTE = 20.0 ! factor times hdur_movie
+  real(kind=CUSTOM_REAL) :: RADIUS_TO_MUTE = 1.0    ! start radius in degrees  
+  real(kind=CUSTOM_REAL) :: STARTTIME_TO_MUTE = 2.0 ! factor times hdur_movie
 
   ! normalizes output values
   logical, parameter :: NORMALIZE_VALUES = .true.
@@ -393,23 +393,23 @@
             ! approximate wavefront travel distance in degrees (~3.5 km/s wave speed for surface waves)
             mute_factor = 3.5 * (it-1)*DT / 6371. * 180./PI
             
-            ! approximate distance to source
+            ! approximate distance to source (in degrees)
             do while ( mute_factor > 360. )
               mute_factor = mute_factor - 360.
             enddo
             if( mute_factor > 180. ) mute_factor = 360. - mute_factor
 
-            ! limit size around source
-            if( mute_factor < 10. ) then
-              mute_factor = 0.0
-            endif
+            ! limit size around source (in degrees)
+            !if( mute_factor < 10. ) then
+            !  mute_factor = 0.0
+            !endif
             if( mute_factor > 80. ) then
               mute_factor = 80.0
             endif
             
-            print*,'muting radius: ',0.5 * mute_factor
+            print*,'muting radius: ',0.7 * mute_factor
                         
-            RADIUS_TO_MUTE = 0.5 * mute_factor * PI/180.
+            RADIUS_TO_MUTE = 0.7 * mute_factor * PI/180.
             
           else
             ! mute_factor used at the beginning for scaling displacement values

Modified: seismo/3D/SPECFEM3D_GLOBE/tags/v5.0.1/write_seismograms.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/tags/v5.0.1/write_seismograms.f90	2010-03-24 19:02:13 UTC (rev 16451)
+++ seismo/3D/SPECFEM3D_GLOBE/tags/v5.0.1/write_seismograms.f90	2010-03-24 19:22:49 UTC (rev 16452)
@@ -693,7 +693,7 @@
       ! indicates SEM synthetics
       ! by Ebru
       KUSER0 = 'SEM'          !  A8
-      KUSER1 = 'v5.0.0 '
+      KUSER1 = 'v5.0.1'
       KUSER2 = 'Tiger' ! aka. awesome (princeton) tiger version :)
 
       !KUSER0 = 'PDE_LAT_'          !  A8