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

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2009-07-21 16:50:13 -0700 (Tue, 21 Jul 2009)
New Revision: 15469

Modified:
   seismo/2D/SPECFEM2D/trunk/define_external_model.f90
Log:
fixed missing "myrank" argument in define_external_model


Modified: seismo/2D/SPECFEM2D/trunk/define_external_model.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/define_external_model.f90	2009-07-21 06:01:21 UTC (rev 15468)
+++ seismo/2D/SPECFEM2D/trunk/define_external_model.f90	2009-07-21 23:50:13 UTC (rev 15469)
@@ -40,7 +40,7 @@
 !
 !========================================================================
 
-  subroutine define_external_model(x,y,iflag_element,rho,vp,vs)
+  subroutine define_external_model(x,y,iflag_element,rho,vp,vs,myrank)
 
   implicit none
 
@@ -49,14 +49,14 @@
 ! user can modify this routine to assign any different external Earth model (rho, vp, vs)
 ! based on the x and y coordinates of that grid point and the flag of the region it belongs to
 
-  integer, intent(in) :: iflag_element
+  integer, intent(in) :: iflag_element,myrank
 
   double precision, intent(in) :: x,y
 
   double precision, intent(out) :: rho,vp,vs
 
 ! dummy routine here, just to demonstrate how the model can be assigned
-  if(iflag_element == 1 .or. x < 1700.d0 .or. y >= 2300.d0) then
+  if(myrank == 0 .and. iflag_element == 1 .or. x < 1700.d0 .or. y >= 2300.d0) then
     rho = 2000.d0
     vp = 3000.d0
     vs = vp / sqrt(3.d0)