[cig-commits] [commit] devel: added a stop statement if poroelastic materials are used for adjoint runs because Christina Morency mentioned that her routines for poroelastic materials are fully OK for forward runs but not debugged yet for adjoint runs (5482435)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/specfem3d

On branch  : devel
Link       : https://github.com/geodynamics/specfem3d/compare/5cd4f87b76191f953d3356cba3bcc7c082410955...5482435a44c766976b16ef4139e5b545834a7300

>---------------------------------------------------------------

commit 5482435a44c766976b16ef4139e5b545834a7300
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date:   Sun Dec 21 01:03:33 2014 +0100

    added a stop statement if poroelastic materials are used for adjoint runs because Christina Morency mentioned that her routines for poroelastic materials are fully OK for forward runs but not debugged yet for adjoint runs


>---------------------------------------------------------------

5482435a44c766976b16ef4139e5b545834a7300
 src/specfem3D/initialize_simulation.f90 | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/src/specfem3D/initialize_simulation.f90 b/src/specfem3D/initialize_simulation.f90
index 58b33c0..105de19 100644
--- a/src/specfem3D/initialize_simulation.f90
+++ b/src/specfem3D/initialize_simulation.f90
@@ -373,6 +373,7 @@
   use specfem_par_elastic
   use specfem_par_acoustic
   use specfem_par_poroelastic
+
   implicit none
 
   ! check simulation parameters
@@ -386,6 +387,8 @@
   if (SIMULATION_TYPE /= 1 .and. NSOURCES > 999999) &
     call exit_MPI(myrank,'for adjoint simulations, NSOURCES <= 999999, if you need more change i6.6 in write_seismograms.f90')
 
+  if (SIMULATION_TYPE /= 1 .and. POROELASTIC_SIMULATION) stop 'poroelastic simulations for adjoint runs not supported yet'
+
   ! snapshot file names: ADJOINT attenuation
   if (ATTENUATION .and. ((SIMULATION_TYPE == 1 .and. SAVE_FORWARD) .or. SIMULATION_TYPE == 3)) &
     call create_name_database(prname_Q,myrank,OUTPUT_FILES_PATH)
@@ -446,7 +449,7 @@
       stop 'GPU mode does not support N_SLS /= 3 yet'
   endif
 
-  if (POROELASTIC_SIMULATION) stop 'poroelastic simulations on GPU not supported yet'
+  if (POROELASTIC_SIMULATION) stop 'poroelastic simulations on GPUs not supported yet'
 
   ! initializes GPU and outputs info to files for all processes
   call initialize_cuda_device(myrank,ncuda_devices)