[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