[cig-commits] [commit] master: these have the same realkind, so why typecast? (30efcd9)

cig_noreply at geodynamics.org cig_noreply at geodynamics.org
Fri Oct 17 05:29:12 PDT 2014


Repository : https://github.com/geodynamics/axisem

On branch  : master
Link       : https://github.com/geodynamics/axisem/compare/607f803cf074063627513d235f9ed0837fc1dd44...b6457db24acdde4a4e1c08935ae1b22adf87f5bf

>---------------------------------------------------------------

commit 30efcd9bdc1a38784800203653ccb9ad31d1e48c
Author: martinvandriel <martin at vandriel.de>
Date:   Wed Oct 15 14:35:33 2014 +0200

    these have the same realkind, so why typecast?


>---------------------------------------------------------------

30efcd9bdc1a38784800203653ccb9ad31d1e48c
 SOLVER/seismograms.f90 | 20 ++++++++++----------
 1 file changed, 10 insertions(+), 10 deletions(-)

diff --git a/SOLVER/seismograms.f90 b/SOLVER/seismograms.f90
index 13157fb..02d56f6 100644
--- a/SOLVER/seismograms.f90
+++ b/SOLVER/seismograms.f90
@@ -909,23 +909,23 @@ subroutine nc_compute_recfile_seis_bare(disp, iseismo)
 
   if (src_type(1) == 'monopole') then
      do i=1, num_rec
-          disp_rec(1,i) = real(disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),1))
+          disp_rec(1,i) = disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),1)
           disp_rec(2,i) = 0
-          disp_rec(3,i) = real(disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),3))  
+          disp_rec(3,i) = disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),3)
      enddo
   elseif (src_type(1) == 'dipole') then
      do i=1, num_rec
-          disp_rec(1,i) = real(disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),1) &
-                             + disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),2))
-          disp_rec(2,i) = real(disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),1) &
-                             - disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),2))
-          disp_rec(3,i) = real(disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),3))
+          disp_rec(1,i) = disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),1) &
+                        + disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),2)
+          disp_rec(2,i) = disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),1) &
+                        - disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),2)
+          disp_rec(3,i) = disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),3)
      enddo
   elseif (src_type(1) == 'quadpole') then
      do i=1, num_rec
-          disp_rec(1,i) = real(disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),1))
-          disp_rec(2,i) = real(disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),2))
-          disp_rec(3,i) = real(disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),3))  
+          disp_rec(1,i) = disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),1)
+          disp_rec(2,i) = disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),2)
+          disp_rec(3,i) = disp(recfile_el(i,2), recfile_el(i,3), recfile_el(i,1),3)
      enddo
   end if !src_type(1)
 



More information about the CIG-COMMITS mailing list