[cig-commits] r15657 - in seismo/2D/SPECFEM2D/trunk: . DATA
cmorency at geodynamics.org
cmorency at geodynamics.org
Fri Sep 4 17:59:36 PDT 2009
Author: cmorency
Date: 2009-09-04 17:59:36 -0700 (Fri, 04 Sep 2009)
New Revision: 15657
Modified:
seismo/2D/SPECFEM2D/trunk/DATA/Par_file
seismo/2D/SPECFEM2D/trunk/specfem2D.F90
Log:
Fixed a bug in specfem2D.F90 on plotting diplacement ps files during forward calculations. The adjoint displacement filed was wrongly entered.
Thanks to Ali Can BAKIR for pointing that out.
Modified: seismo/2D/SPECFEM2D/trunk/DATA/Par_file
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file 2009-09-04 23:41:52 UTC (rev 15656)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file 2009-09-05 00:59:36 UTC (rev 15657)
@@ -73,7 +73,7 @@
# display parameters
NTSTEP_BETWEEN_OUTPUT_INFO = 2000 # display frequency in time steps
-output_postscript_snapshot = .false. # output Postscript snapshot of the results
+output_postscript_snapshot = .true. # output Postscript snapshot of the results
output_color_image = .true. # output color image of the results
imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
cutsnaps = 1. # minimum amplitude in % for snapshots
Modified: seismo/2D/SPECFEM2D/trunk/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/specfem2D.F90 2009-09-04 23:41:52 UTC (rev 15656)
+++ seismo/2D/SPECFEM2D/trunk/specfem2D.F90 2009-09-05 00:59:36 UTC (rev 15657)
@@ -6719,7 +6719,7 @@
if (myrank == 0) write(IOUT,*) 'drawing displacement vector as small arrows...'
- call compute_vector_whole_medium(potential_acoustic,b_displ_elastic,displs_poroelastic,&
+ call compute_vector_whole_medium(potential_acoustic,displ_elastic,displs_poroelastic,&
elastic,poroelastic,vector_field_display, &
xix,xiz,gammax,gammaz,ibool,hprime_xx,hprime_zz,nspec,npoin,numat,kmato,density,rhoext,assign_external_model)
More information about the CIG-COMMITS
mailing list