[cig-commits] r17196 - seismo/3D/SPECFEM3D/trunk

[yangl at geodynamics.org]

Author: yangl
Date: 2010-09-15 13:54:02 -0700 (Wed, 15 Sep 2010)
New Revision: 17196

Modified:
   seismo/3D/SPECFEM3D/trunk/compute_forces_acoustic.f90
   seismo/3D/SPECFEM3D/trunk/save_arrays_solver.f90
Log:
fix a bug for reconstructed forward wavefield in acoustic adjoint simulations. (check deltat & b_deltat carefully)

Modified: seismo/3D/SPECFEM3D/trunk/compute_forces_acoustic.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/compute_forces_acoustic.f90	2010-09-15 12:43:49 UTC (rev 17195)
+++ seismo/3D/SPECFEM3D/trunk/compute_forces_acoustic.f90	2010-09-15 20:54:02 UTC (rev 17196)
@@ -286,7 +286,7 @@
 
   ! adjoint simulations  
   if (SIMULATION_TYPE == 3) &
-    b_potential_dot_acoustic(:) = b_potential_dot_acoustic(:) + deltatover2*b_potential_dot_dot_acoustic(:)
+    b_potential_dot_acoustic(:) = b_potential_dot_acoustic(:) + b_deltatover2*b_potential_dot_dot_acoustic(:)
 
   ! updates potential_dot_acoustic and potential_dot_dot_acoustic inside PML region for plotting seismograms/movies
   if(PML) call PML_acoustic_update_potentials(NGLOB_AB,NSPEC_AB, &

Modified: seismo/3D/SPECFEM3D/trunk/save_arrays_solver.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/save_arrays_solver.f90	2010-09-15 12:43:49 UTC (rev 17195)
+++ seismo/3D/SPECFEM3D/trunk/save_arrays_solver.f90	2010-09-15 20:54:02 UTC (rev 17196)
@@ -328,6 +328,15 @@
                         xstore_dummy,ystore_dummy,zstore_dummy,ibool, &
                         v_tmp,filename)
 
+!<YANGL
+!!! output density model for check
+    v_tmp = 0.0
+    where( rho_vp /= 0._CUSTOM_REAL ) v_tmp = rho_vp**2 / (FOUR_THIRDS * mustore + kappastore)
+    open(unit=27,file=prname(1:len_trim(prname))//'rho.bin',status='unknown',form='unformatted')
+    write(27) v_tmp
+    close(27)
+!>YANGL
+
     ! VTK file output
     ! saves attenuation flag assigned on each gll point into a vtk file 
     filename = prname(1:len_trim(prname))//'attenuation_flag'
@@ -853,4 +862,4 @@
 !
 !!=============================================================
     
-  
\ No newline at end of file
+