[cig-commits] r21282 - seismo/3D/SPECFEM3D/trunk/src/specfem3D

dkomati1 at geodynamics.org dkomati1 at geodynamics.org
Sat Jan 19 10:05:34 PST 2013


Author: dkomati1
Date: 2013-01-19 10:05:34 -0800 (Sat, 19 Jan 2013)
New Revision: 21282

Modified:
   seismo/3D/SPECFEM3D/trunk/src/specfem3D/pml_compute_memory_variables.f90
Log:
added "else if" instead of different independent "if"s to speed up the comparisons


Modified: seismo/3D/SPECFEM3D/trunk/src/specfem3D/pml_compute_memory_variables.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/specfem3D/pml_compute_memory_variables.f90	2013-01-19 18:03:21 UTC (rev 21281)
+++ seismo/3D/SPECFEM3D/trunk/src/specfem3D/pml_compute_memory_variables.f90	2013-01-19 18:05:34 UTC (rev 21282)
@@ -317,9 +317,8 @@
         enddo
      enddo
   enddo
-  endif
 
-  if( CPML_regions(ispec_CPML) == 2 ) then
+  else if( CPML_regions(ispec_CPML) == 2 ) then
   do k=1,NGLLZ
      do j=1,NGLLY
         do i=1,NGLLX
@@ -571,9 +570,8 @@
         enddo
      enddo
   enddo
-  endif
 
- if( CPML_regions(ispec_CPML) == 3 ) then
+ else if( CPML_regions(ispec_CPML) == 3 ) then
   do k=1,NGLLZ
      do j=1,NGLLY
         do i=1,NGLLX
@@ -825,9 +823,8 @@
         enddo
      enddo
   enddo
-  endif
 
-  if( CPML_regions(ispec_CPML) == 4 ) then
+  else if( CPML_regions(ispec_CPML) == 4 ) then
   do k=1,NGLLZ
      do j=1,NGLLY
         do i=1,NGLLX
@@ -1411,9 +1408,8 @@
         enddo
      enddo
   enddo
-  endif
 
-  if( CPML_regions(ispec_CPML) == 6 ) then
+  else if( CPML_regions(ispec_CPML) == 6 ) then
 
   do k=1,NGLLZ
      do j=1,NGLLY
@@ -1703,9 +1699,8 @@
         enddo
      enddo
   enddo
-  endif
 
-  if( CPML_regions(ispec_CPML) == 7 ) then
+  else if( CPML_regions(ispec_CPML) == 7 ) then
 
   do k=1,NGLLZ
      do j=1,NGLLY
@@ -2184,6 +2179,11 @@
         enddo
      enddo
   enddo
+
+  else
+
+    stop 'wrong PML flag in PML memory variable calculation routine'
+
   endif
 
 end subroutine pml_compute_memory_variables



More information about the CIG-COMMITS mailing list