[cig-commits] r13024 - seismo/2D/SPECFEM2D/branches/BIOT

[cmorency at geodynamics.org]

Author: cmorency
Date: 2008-10-14 10:47:17 -0700 (Tue, 14 Oct 2008)
New Revision: 13024

Modified:
   seismo/2D/SPECFEM2D/branches/BIOT/Makefile
   seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90
Log:
Save acoustic potential for adjoint source calculation


Modified: seismo/2D/SPECFEM2D/branches/BIOT/Makefile
===================================================================
--- seismo/2D/SPECFEM2D/branches/BIOT/Makefile	2008-10-14 10:14:38 UTC (rev 13023)
+++ seismo/2D/SPECFEM2D/branches/BIOT/Makefile	2008-10-14 17:47:17 UTC (rev 13024)
@@ -54,18 +54,18 @@
 #FLAGS_CHECK=-fast -Mbounds -Mneginfo -Mdclchk -Minform=warn
 
 # Intel
-F90 = ifort
-CC = gcc
-FLAGS_NOCHECK=-O0 -implicitnone -warn stderrors -warn truncated_source -warn argument_checking -warn unused -warn declarations -std95 -assume byterecl -check nobounds
-FLAGS_CHECK = $(FLAGS_NOCHECK)
+#F90 = ifort
+#CC = gcc
+#FLAGS_NOCHECK=-O0 -implicitnone -warn stderrors -warn truncated_source -warn argument_checking -warn unused -warn declarations -std95 -assume byterecl -check nobounds
+#FLAGS_CHECK = $(FLAGS_NOCHECK)
 
 # GNU gfortran
 #F90 = /opt/openmpi-1.2.1/gfortran64/bin/mpif90 -DUSE_MPI -DUSE_METIS -DUSE_SCOTCH
-#F90 = gfortran
-#CC = gcc
+F90 = gfortran
+CC = gcc
 #FLAGS_NOCHECK = -O3 -march=opteron -m64 -mfpmath=sse,387
-#FLAGS_NOCHECK = -std=gnu -fimplicit-none -frange-check -O2 -Wunused-labels -Waliasing -Wampersand -Wsurprising -Wline-truncation -Wunderflow
-#FLAGS_CHECK = $(FLAGS_NOCHECK) -fbounds-check
+FLAGS_NOCHECK = -std=gnu -fimplicit-none -frange-check -O2 -Wunused-labels -Waliasing -Wampersand -Wsurprising -Wline-truncation -Wunderflow
+FLAGS_CHECK = $(FLAGS_NOCHECK) -fbounds-check
 
 # IBM
 #F90 = xlf_r

Modified: seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90	2008-10-14 10:14:38 UTC (rev 13023)
+++ seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90	2008-10-14 17:47:17 UTC (rev 13024)
@@ -4813,6 +4813,12 @@
           dxd = vector_field_element(1,i,j)
           dzd = vector_field_element(2,i,j)
 
+        else if(.not. elastic(ispec) .and. .not. poroelastic(ispec) &
+                           .and. save_forward) then
+
+          dxd = potential_acoustic(iglob)
+          dzd = ZERO
+
         else if(seismotype == 1) then
 
              if(poroelastic(ispec)) then