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

[walter at geodynamics.org]

Author: walter
Date: 2007-12-07 15:47:36 -0800 (Fri, 07 Dec 2007)
New Revision: 8457

Modified:
   seismo/2D/SPECFEM2D/trunk/plotpost.f90
   seismo/2D/SPECFEM2D/trunk/specfem2D.f90
Log:
fixed iglobsource bug detected by Paul Cristini in plotpost.f90


Modified: seismo/2D/SPECFEM2D/trunk/plotpost.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/plotpost.f90	2005-01-26 22:17:23 UTC (rev 8456)
+++ seismo/2D/SPECFEM2D/trunk/plotpost.f90	2007-12-07 23:47:36 UTC (rev 8457)
@@ -11,7 +11,7 @@
 !
 !========================================================================
 
-  subroutine plotpost(displ,coord,vpext,iglob_source,st_xval,st_zval,it,dt,coorg, &
+  subroutine plotpost(displ,coord,vpext,x_source,z_source,st_xval,st_zval,it,dt,coorg, &
           xinterp,zinterp,shapeint, &
           Uxinterp,Uzinterp,flagrange,density,elastcoef,knods,kmato,ibool, &
           numabs,codeabs,anyabs,stitle,npoin,npgeo,vpmin,vpmax,nrec, &
@@ -31,7 +31,7 @@
   double precision, dimension(MAXCOLORS) :: red,green,blue
 
   integer it,nrec,nelemabs,numat,iptsdisp,nspec
-  integer i,iglob_source,npoin,npgeo,ngnod
+  integer i,npoin,npgeo,ngnod
 
   integer kmato(nspec),knods(ngnod,nspec)
   integer ibool(NGLLX,NGLLZ,nspec)
@@ -43,7 +43,7 @@
   double precision flagrange(NGLLX,iptsdisp)
   double precision density(numat),elastcoef(4,numat)
 
-  double precision dt,timeval
+  double precision dt,timeval,x_source,z_source
   double precision displ(NDIM,npoin),coord(NDIM,npoin)
   double precision vpext(npoin)
 
@@ -849,8 +849,8 @@
 !
 !----  write position of the source
 !
-  xw = coord(1,iglob_source)
-  zw = coord(2,iglob_source)
+  xw = x_source
+  zw = z_source
   xw = (xw-xmin)*rapp_page + orig_x
   zw = (zw-zmin)*rapp_page + orig_z
   xw = xw * centim

Modified: seismo/2D/SPECFEM2D/trunk/specfem2D.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/specfem2D.f90	2005-01-26 22:17:23 UTC (rev 8456)
+++ seismo/2D/SPECFEM2D/trunk/specfem2D.f90	2007-12-07 23:47:36 UTC (rev 8457)
@@ -447,7 +447,7 @@
 !
 
 ! version "propre mais lente" ou version "sale mais rapide"
-  if(fast_numbering) then
+  if(FAST_NUMBERING) then
     call createnum_fast(knods,ibool,shape2D,coorg,npoin,npgeo,nspec,ngnod)
   else
     call createnum_slow(knods,ibool,npoin,nspec,ngnod)
@@ -1635,7 +1635,7 @@
 ! for elastic medium
   if(ELASTIC .and. vecttype == 1) then
     write(IOUT,*) 'drawing displacement field...'
-    call plotpost(displ,coord,vpext,iglob_source,st_xval,st_zval, &
+    call plotpost(displ,coord,vpext,x_source,z_source,st_xval,st_zval, &
           it,deltat,coorg,xinterp,zinterp,shape2D_display, &
           Uxinterp,Uzinterp,flagrange,density,elastcoef,knods,kmato,ibool, &
           numabs,codeabs,anyabs,stitle,npoin,npgeo,vpmin,vpmax,nrec, &
@@ -1644,7 +1644,7 @@
 
   else if(ELASTIC .and. vecttype == 2) then
     write(IOUT,*) 'drawing velocity field...'
-    call plotpost(veloc,coord,vpext,iglob_source,st_xval,st_zval, &
+    call plotpost(veloc,coord,vpext,x_source,z_source,st_xval,st_zval, &
           it,deltat,coorg,xinterp,zinterp,shape2D_display, &
           Uxinterp,Uzinterp,flagrange,density,elastcoef,knods,kmato,ibool, &
           numabs,codeabs,anyabs,stitle,npoin,npgeo,vpmin,vpmax,nrec, &
@@ -1653,7 +1653,7 @@
 
   else if(ELASTIC .and. vecttype == 3) then
     write(IOUT,*) 'drawing acceleration field...'
-    call plotpost(accel,coord,vpext,iglob_source,st_xval,st_zval, &
+    call plotpost(accel,coord,vpext,x_source,z_source,st_xval,st_zval, &
           it,deltat,coorg,xinterp,zinterp,shape2D_display, &
           Uxinterp,Uzinterp,flagrange,density,elastcoef,knods,kmato,ibool, &
           numabs,codeabs,anyabs,stitle,npoin,npgeo,vpmin,vpmax,nrec, &
@@ -1669,7 +1669,7 @@
 ! for acoustic medium, compute gradient for display, displ represents the potential
     call compute_gradient_fluid(displ,vector_field_postscript, &
           xix,xiz,gammax,gammaz,ibool,hprime_xx,hprime_zz,NSPEC,npoin)
-    call plotpost(vector_field_postscript,coord,vpext,iglob_source,st_xval,st_zval, &
+    call plotpost(vector_field_postscript,coord,vpext,x_source,z_source,st_xval,st_zval, &
           it,deltat,coorg,xinterp,zinterp,shape2D_display, &
           Uxinterp,Uzinterp,flagrange,density,elastcoef,knods,kmato,ibool, &
           numabs,codeabs,anyabs,stitle,npoin,npgeo,vpmin,vpmax,nrec, &
@@ -1681,7 +1681,7 @@
 ! for acoustic medium, compute gradient for display, veloc represents the first derivative of the potential
     call compute_gradient_fluid(veloc,vector_field_postscript, &
           xix,xiz,gammax,gammaz,ibool,hprime_xx,hprime_zz,NSPEC,npoin)
-    call plotpost(vector_field_postscript,coord,vpext,iglob_source,st_xval,st_zval, &
+    call plotpost(vector_field_postscript,coord,vpext,x_source,z_source,st_xval,st_zval, &
           it,deltat,coorg,xinterp,zinterp,shape2D_display, &
           Uxinterp,Uzinterp,flagrange,density,elastcoef,knods,kmato,ibool, &
           numabs,codeabs,anyabs,stitle,npoin,npgeo,vpmin,vpmax,nrec, &