[cig-commits] r16264 - seismo/2D/SPECFEM2D/trunk

[cmorency at geodynamics.org]

Author: cmorency
Date: 2010-02-15 17:16:43 -0800 (Mon, 15 Feb 2010)
New Revision: 16264

Modified:
   seismo/2D/SPECFEM2D/trunk/specfem2D.F90
Log:
Corrected the unit of "anglerec" when using rec_normal_to_surface = .false., that is, replacing cosrot_irec(:)=cos(anglerec) by cosrot_irec(:)=cos(anglerec_irec(:)), where anglerec_irec(:) = anglerec *pi/180 is in radian.


Modified: seismo/2D/SPECFEM2D/trunk/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/specfem2D.F90	2010-02-15 18:46:19 UTC (rev 16263)
+++ seismo/2D/SPECFEM2D/trunk/specfem2D.F90	2010-02-16 01:16:43 UTC (rev 16264)
@@ -1876,8 +1876,8 @@
     stop 'anglerec should be zero when receivers are normal to the topography'
 
   anglerec_irec(:) = anglerec * pi / 180.d0
-  cosrot_irec(:) = cos(anglerec)
-  sinrot_irec(:) = sin(anglerec)
+  cosrot_irec(:) = cos(anglerec_irec(:))
+  sinrot_irec(:) = sin(anglerec_irec(:))
 endif
 
 !