[cig-commits] r20200 - seismo/2D/SPECFEM2D/trunk/src/specfem2D

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2012-05-20 14:37:15 -0700 (Sun, 20 May 2012)
New Revision: 20200

Modified:
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/define_external_model.f90
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_external_model.f90
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
Log:
no need for argument myrank in define_external_model(), as suggested by Paul Cristini


Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/define_external_model.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/define_external_model.f90	2012-05-20 17:48:51 UTC (rev 20199)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/define_external_model.f90	2012-05-20 21:37:15 UTC (rev 20200)
@@ -43,7 +43,7 @@
 !========================================================================
 
 
-  subroutine define_external_model(x,y,iflag_element,myrank,rho,vp,vs,QKappa_attenuation,&
+  subroutine define_external_model(x,y,iflag_element,rho,vp,vs,QKappa_attenuation,&
        Qmu_attenuation,c11,c13,c15,c33,c35,c55 )
 
   implicit none
@@ -53,7 +53,7 @@
 ! 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,myrank
+  integer, intent(in) :: iflag_element
 
   double precision, intent(in) :: x,y
 
@@ -62,7 +62,7 @@
   double precision, intent(out) :: c11,c15,c13,c33,c35,c55
 
 ! dummy routine here, just to demonstrate how the model can be assigned
-   if(myrank == 0 .and. iflag_element == 1 .or. x < 1700.d0 .or. y >= 2300.d0) then
+   if(iflag_element == 1 .or. x < 1700.d0 .or. y >= 2300.d0) then
      rho = 2000.d0
      vp = 3000.d0
      vs = vp / sqrt(3.d0)
@@ -89,3 +89,4 @@
    endif
 
   end subroutine define_external_model
+

Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_external_model.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_external_model.f90	2012-05-20 17:48:51 UTC (rev 20199)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_external_model.f90	2012-05-20 21:37:15 UTC (rev 20200)
@@ -49,14 +49,14 @@
                 f0_attenuation,inv_tau_sigma_nu1_sent,phi_nu1_sent, &
                 inv_tau_sigma_nu2_sent,phi_nu2_sent,Mu_nu1_sent,Mu_nu2_sent, &
                 inv_tau_sigma_nu1,inv_tau_sigma_nu2,phi_nu1,phi_nu2,Mu_nu1,Mu_nu2,&
-                coord,kmato,myrank,rhoext,vpext,vsext, &
+                coord,kmato,rhoext,vpext,vsext, &
                 QKappa_attenuationext,Qmu_attenuationext, &
                 c11ext,c13ext,c15ext,c33ext,c35ext,c55ext,READ_EXTERNAL_SEP_FILE)
 
   implicit none
   include "constants.h"
 
-  integer :: nspec,myrank,nglob
+  integer :: nspec,nglob
   double precision  :: f0_attenuation
 
   ! Mesh
@@ -117,7 +117,7 @@
         do i = 1,NGLLX
 
           iglob = ibool(i,j,ispec)
-          call define_external_model(coord(1,iglob),coord(2,iglob),kmato(ispec),myrank,&
+          call define_external_model(coord(1,iglob),coord(2,iglob),kmato(ispec),&
                                     rhoext(i,j,ispec),vpext(i,j,ispec),vsext(i,j,ispec), &
                                     QKappa_attenuationext(i,j,ispec),Qmu_attenuationext(i,j,ispec),&
                                     c11ext(i,j,ispec),c13ext(i,j,ispec),c15ext(i,j,ispec), &

Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90	2012-05-20 17:48:51 UTC (rev 20199)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90	2012-05-20 21:37:15 UTC (rev 20200)
@@ -1990,7 +1990,7 @@
                 f0_attenuation,inv_tau_sigma_nu1_sent,phi_nu1_sent, &
                 inv_tau_sigma_nu2_sent,phi_nu2_sent,Mu_nu1_sent,Mu_nu2_sent, &
                 inv_tau_sigma_nu1,inv_tau_sigma_nu2,phi_nu1,phi_nu2,Mu_nu1,Mu_nu2,&
-                coord,kmato,myrank,rhoext,vpext,vsext, &
+                coord,kmato,rhoext,vpext,vsext, &
                 QKappa_attenuationext,Qmu_attenuationext, &
                 c11ext,c13ext,c15ext,c33ext,c35ext,c55ext,READ_EXTERNAL_SEP_FILE)
   end if