[cig-commits] r22316 - seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/src/specfem3D

[xie.zhinan at geodynamics.org]

Author: xie.zhinan
Date: 2013-06-15 12:13:54 -0700 (Sat, 15 Jun 2013)
New Revision: 22316

Modified:
   seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/src/specfem3D/compute_element.F90
Log:
fix one error in computation of strain inside one element


Modified: seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/src/specfem3D/compute_element.F90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/src/specfem3D/compute_element.F90	2013-06-15 19:07:25 UTC (rev 22315)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/undo_att/src/specfem3D/compute_element.F90	2013-06-15 19:13:54 UTC (rev 22316)
@@ -1616,16 +1616,17 @@
         duzdxl_plus_duxdzl = duzdxl + duxdzl
         duzdyl_plus_duydzl = duzdyl + duydzl
 
+        eps_trace_over_3_loc(i,j,k) = ONE_THIRD * (duxdxl + duydyl + duzdzl) 
+        epsilondev_loc(1,i,j,k) = duxdxl - eps_trace_over_3_loc(i,j,k)
+        epsilondev_loc(2,i,j,k) = duydyl - eps_trace_over_3_loc(i,j,k)
+        epsilondev_loc(3,i,j,k) = 0.5 * duxdyl_plus_duydxl
+        epsilondev_loc(4,i,j,k) = 0.5 * duzdxl_plus_duxdzl
+        epsilondev_loc(5,i,j,k) = 0.5 * duzdyl_plus_duydzl
+
       enddo
     enddo
   enddo
 
-  eps_trace_over_3_loc(i,j,k) = ONE_THIRD * (duxdxl + duydyl + duzdzl) 
-  epsilondev_loc(1,i,j,k) = duxdxl - eps_trace_over_3_loc(i,j,k)
-  epsilondev_loc(2,i,j,k) = duydyl - eps_trace_over_3_loc(i,j,k)
-  epsilondev_loc(3,i,j,k) = 0.5 * duxdyl_plus_duydxl
-  epsilondev_loc(4,i,j,k) = 0.5 * duzdxl_plus_duxdzl
-  epsilondev_loc(5,i,j,k) = 0.5 * duzdyl_plus_duydzl
 
  end subroutine compute_element_strain_undo_att_Dev
 
@@ -1798,17 +1799,19 @@
         duzdxl_plus_duxdzl = duzdxl + duxdzl
         duzdyl_plus_duydzl = duzdyl + duydzl
 
+        templ = ONE_THIRD * (duxdxl + duydyl + duzdzl) 
+        epsilondev_loc_nplus1(1,i,j,k) = duxdxl - templ
+        epsilondev_loc_nplus1(2,i,j,k) = duydyl - templ
+        epsilondev_loc_nplus1(3,i,j,k) = 0.5 * duxdyl_plus_duydxl
+        epsilondev_loc_nplus1(4,i,j,k) = 0.5 * duzdxl_plus_duxdzl
+        epsilondev_loc_nplus1(5,i,j,k) = 0.5 * duzdyl_plus_duydzl
+
       enddo
     enddo
   enddo
 
-  templ = ONE_THIRD * (duxdxl + duydyl + duzdzl) 
-  epsilondev_loc_nplus1(1,i,j,k) = duxdxl - templ
-  epsilondev_loc_nplus1(2,i,j,k) = duydyl - templ
-  epsilondev_loc_nplus1(3,i,j,k) = 0.5 * duxdyl_plus_duydxl
-  epsilondev_loc_nplus1(4,i,j,k) = 0.5 * duzdxl_plus_duxdzl
-  epsilondev_loc_nplus1(5,i,j,k) = 0.5 * duzdyl_plus_duydzl
 
+
  end subroutine compute_element_strain_att_Dev
 !=====================================================================