[cig-commits] r20534 - in seismo/3D/SPECFEM3D/trunk: examples/homogeneous_poroelastic/in_data_files src/generate_databases src/specfem3D

[cmorency at geodynamics.org]

Author: cmorency
Date: 2012-07-22 08:46:42 -0700 (Sun, 22 Jul 2012)
New Revision: 20534

Modified:
   seismo/3D/SPECFEM3D/trunk/examples/homogeneous_poroelastic/in_data_files/Par_file
   seismo/3D/SPECFEM3D/trunk/src/generate_databases/save_arrays_solver.f90
   seismo/3D/SPECFEM3D/trunk/src/specfem3D/read_mesh_databases.f90
Log:
Updated Par_file in example/homogeneous_poroelastic/
Corrected arguments in check_mesh_resolution_poro() in /specfem3D/read_mesh_database.f90
Suppressed stop flag for forward poro calculation in generate_databases/save_arrays_solver.f90


Modified: seismo/3D/SPECFEM3D/trunk/examples/homogeneous_poroelastic/in_data_files/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/homogeneous_poroelastic/in_data_files/Par_file	2012-07-21 16:48:13 UTC (rev 20533)
+++ seismo/3D/SPECFEM3D/trunk/examples/homogeneous_poroelastic/in_data_files/Par_file	2012-07-22 15:46:42 UTC (rev 20534)
@@ -5,18 +5,18 @@
 SIMULATION_TYPE                 = 1   
 # 0 = earthquake simulation,  1/2/3 = three steps in noise simulation
 NOISE_TOMOGRAPHY                = 0   
-SAVE_FORWARD                    = .false.
+SAVE_FORWARD                    = .true.
 
 # UTM projection parameters
 UTM_PROJECTION_ZONE             = 11
 SUPPRESS_UTM_PROJECTION         = .true.
 
 # number of MPI processors
-NPROC                           = 16
+NPROC                           = 14
 
 # time step parameters
-NSTEP                           = 4000 #4000
-DT                              = 3.d-5 #0.5d-3
+NSTEP                           = 4000
+DT                              = 5.d-5 #0.5d-3
 
 # models:
 # available options are:
@@ -33,6 +33,7 @@
 ATTENUATION                     = .false.
 USE_OLSEN_ATTENUATION           = .false.
 ANISOTROPY                      = .false.
+GRAVITY                         = .false.
 
 # absorbing boundary conditions for a regional simulation
 ABSORBING_CONDITIONS            = .false.
@@ -53,12 +54,17 @@
 LOCAL_PATH                      = ../in_out_files/DATABASES_MPI
 
 # interval at which we output time step info and max of norm of displacement
-NTSTEP_BETWEEN_OUTPUT_INFO      = 500
+NTSTEP_BETWEEN_OUTPUT_INFO      = 250
 
 # interval in time steps for writing of seismograms
 NTSTEP_BETWEEN_OUTPUT_SEISMOS   = 10000
 
+# interval in time steps for reading adjoint traces
+# # 0 = read the whole adjoint sources at the same time
+NTSTEP_BETWEEN_READ_ADJSRC      = 0
+
 # print source time function
 PRINT_SOURCE_TIME_FUNCTION      = .true.
 
-
+# set to true to use GPUs 
+GPU_MODE                        = .false.

Modified: seismo/3D/SPECFEM3D/trunk/src/generate_databases/save_arrays_solver.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/generate_databases/save_arrays_solver.f90	2012-07-21 16:48:13 UTC (rev 20533)
+++ seismo/3D/SPECFEM3D/trunk/src/generate_databases/save_arrays_solver.f90	2012-07-22 15:46:42 UTC (rev 20534)
@@ -124,8 +124,6 @@
 ! poroelastic
   call any_all_l( ANY(ispec_is_poroelastic), POROELASTIC_SIMULATION )
   if( POROELASTIC_SIMULATION ) then
-    stop 'not implemented yet: write rmass_solid_poroelastic .. '
-
     write(IOUT) rmass_solid_poroelastic
     write(IOUT) rmass_fluid_poroelastic
     write(IOUT) rhoarraystore

Modified: seismo/3D/SPECFEM3D/trunk/src/specfem3D/read_mesh_databases.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/specfem3D/read_mesh_databases.f90	2012-07-21 16:48:13 UTC (rev 20533)
+++ seismo/3D/SPECFEM3D/trunk/src/specfem3D/read_mesh_databases.f90	2012-07-22 15:46:42 UTC (rev 20534)
@@ -554,7 +554,6 @@
     rho_vp = 0.0_CUSTOM_REAL
     rho_vs = 0.0_CUSTOM_REAL
     call check_mesh_resolution_poro(myrank,NSPEC_AB,NGLOB_AB,ibool,xstore,ystore,zstore, &
-                                    rho_vp,rho_vs, &
                                     DT,model_speed_max,min_resolved_period, &
                                     phistore,tortstore,rhoarraystore,rho_vpI,rho_vpII,rho_vsI, &
                                     LOCAL_PATH,SAVE_MESH_FILES)