[cig-commits] r12604 - seismo/3D/SPECFEM3D_GLOBE/trunk

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2008-08-11 13:36:00 -0700 (Mon, 11 Aug 2008)
New Revision: 12604

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/create_regions_mesh.f90
Log:
fixed a small issue (incorrect flag) when computing the local height of oceans in the case of a 3D Earth


Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/create_regions_mesh.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/create_regions_mesh.f90	2008-08-11 01:47:21 UTC (rev 12603)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/create_regions_mesh.f90	2008-08-11 20:36:00 UTC (rev 12604)
@@ -1604,8 +1604,8 @@
 
         iglobnum=ibool(ix_oceans,iy_oceans,iz_oceans,ispec_oceans)
 
-! compute local height of oceans
-        if(ISOTROPIC_3D_MANTLE) then
+! if 3D Earth, compute local height of oceans
+        if(CASE_3D) then
 
 ! get coordinates of current point
           xval = xstore(ix_oceans,iy_oceans,iz_oceans,ispec_oceans)
@@ -1636,6 +1636,7 @@
           endif
 
         else
+! if 1D Earth, use oceans of constant thickness everywhere
           height_oceans = THICKNESS_OCEANS_PREM
         endif