[cig-commits] r19633 - in seismo/3D/FAULT_SOURCE/branches/new_fault_db: . DATA
ampuero at geodynamics.org
ampuero at geodynamics.org
Tue Feb 14 17:26:11 PST 2012
Author: ampuero
Date: 2012-02-14 17:26:11 -0800 (Tue, 14 Feb 2012)
New Revision: 19633
Added:
seismo/3D/FAULT_SOURCE/branches/new_fault_db/DATA/STATIONS.in
Modified:
seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT
Log:
improved README after workshop with Surendra, Javier and Yingdi
Added: seismo/3D/FAULT_SOURCE/branches/new_fault_db/DATA/STATIONS.in
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/DATA/STATIONS.in (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/DATA/STATIONS.in 2012-02-15 01:26:11 UTC (rev 19633)
@@ -0,0 +1,4 @@
+str12dp00 SC 3000.0 12000.00 0.0 0.0
+str-12dp00 SC 3000.0 -12000.00 0.0 0.0
+str12dp75 SC 3000.0 12000.00 0.0 7500.0
+str-12dp75 SC 3000.0 -12000.00 0.0 7500.0
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT 2012-02-14 23:46:25 UTC (rev 19632)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT 2012-02-15 01:26:11 UTC (rev 19633)
@@ -16,7 +16,7 @@
It is not a replacement for the original manual, which you should read first.
Contents:
- Download
+ Download and updates
Installation
Running a simulation
Examples
@@ -126,9 +126,6 @@
for instance:
EXAMPLES/splay_faults/splay_faults.py
- To run CUBIT from different working directories it is convenient to create the
- environment variable PYTHONPATH (or add to an existing one) with the path SPECFEM3D/CUBIT
-
This creates several mesh files in directory CUBIT/MESH/:
absorbing_surface_file_* (5 files)
free_surface_file
@@ -139,7 +136,10 @@
The CUBIT graphics window should show the mesh (e.g. EXAMPLES/splay_faults/splay_faults.jpg)
+ To run CUBIT from different working directories it is convenient to
+ include the path to SPECFEM3D/CUBIT in an environment variable PYTHONPATH
+
3. Partition the mesh with the domain decomposition software SCOTCH.
>> cd ~/SPECFEM3D/decompose_mesh_SCOTH
@@ -165,6 +165,9 @@
>> mpirun -np nproc ./xgenerate_databases
or
>> ./run/run.xdatabases
+ or
+ >> qsub -l nodes=... -l walltime=... go_mesher
+ (the script go_mesher is in utils/)
This creates several binary mesh files for each processor (proc000***.bin)
in directory DATABASES_MPI.
@@ -173,7 +176,8 @@
6. If this is the first time you run a simulation with this mesh, set the parameter
DT (the time step) DATA/Par_file to a value smaller than the
"Maximum suggested time step" in OUTPUT_FILES/output_mesher.txt.
- Also edit the variables NSTEP and NTSTEP_* in DATA/Par_file accordingly.
+ Consider the note in section IV about the effect of Kelvin-Voigt damping on the critical timestep.
+ Edit the variables NSTEP and NTSTEP_* in DATA/Par_file accordingly.
7. Run the solver:
@@ -181,15 +185,18 @@
>> mpirun -np nproc ./xspecfem3D
or
>>./run/run.xspecfem3d
+ or
+ >> qsub -l nodes=... -l walltime=... go_solver
+ (the script go_solver is in utils/)
III. EXAMPLES
--------------
-The package includes a few examples, the SCEC benchmark problems:
+The package includes examples, the SCEC benchmark problems:
+ TPV5, a planar vertical strike-slip fault
- (+ TPV14 and TPV15, vertical strike-slip fault system with a fault branch)
+ + TPV14 and TPV15, vertical strike-slip fault system with a fault branch
+ Splay fault models from Wendt et al. (2009)
To run the examples:
@@ -327,6 +334,39 @@
The fault-id identifies the fault that contains the station.
It is the index of appearance in the faults list after line 2 of Par_file_faults
+
+DATA/FAULT/input_file.txt Heterogeneous stresses and friction parameters
+ Documented in page 10 of EXAMPLES/tpv16.crack/description/TPV16_17_Description_v03.pdf
+ To activate this feature, in fault_solver.f90 set
+ TPV16 = .true.
+ then re-compile the code:
+ cd SPECFEM3D
+ make
+
+
+Heterogeneous velocity models can be given in a regular grid.
+ In the mesh generation Python file (.py) or in the CUBIT journal file (.jou),
+ set the material 1st attribute to -1:
+ block * attribute index 1 -1
+ In module model_tomography.f90:
+ Set the parameter TOMO_FILENAME to the name of the file containing the velocity model.
+ Re-compile the code:
+ cd SPECFEM3D
+ make
+ The format of the velocity file is:
+ Line 1: ORIG_X, ORIG_Y, ORIG_Z, END_X, END_Y, END_Z
+ (coordinates of the two extreme corners of the box)
+ Line 2: SPACING_X, SPACING_Y, SPACING_Z
+ (regular grid spacing in each direction)
+ Line 3: NX, NY, NZ
+ (number of grid point in each direction)
+ Line 4: VP_MIN, VP_MAX, VS_MIN, VS_MAX, RHO_MIN, RHO_MAX
+ (min and max values of P wave speed, S wave speed and density)
+ Line 5 to 4+NX*NY*NZ:
+ x_tomo,y_tomo,z_tomo,vp_tomo,vs_tomo,rho_tomo
+ (position, P wave speed, S wave speed, density)
+
+
Several files are generated automatically by xgenerate_databases in directory DATABASES_MPI
and do not need to be modified by the user.
@@ -361,14 +401,15 @@
3. Rupture time files are named Rupture_time_FAULT-id. One file is generated for each fault.
The files are ascii and start with a header (12 lines long) followed by a data block
with the following format, one line per fault node:
- # Column #1 = X position (m)
- # Column #2 = Y position (m)
- # Column #3 = Z position (m)
+ # Column #1 = horizontal coordinate, distance along strike (m)
+ # Column #2 = vertical coordinate, distance down-dip (m)
# Column #3 = rupture time (s)
+4. Fault quantities (slip, slip rate, stresses, etc) at regular times are stored
+ in binary data files called Snapshot#it#.bin, where #it# is the timestep number.
+ These can be read in Matlab with the function Post-processing/FSEM3D_snapshot.m
-
V. POST-PROCESSING AND VISUALIZATION
-------------------------------------
@@ -398,9 +439,11 @@
FC = ifort
MPIFC = mpiifort
-Compiling SCOTCH:
+If you can't compile SCOTCH:
+ cp -r /home/surendra/decompose_mesh_SCOTCH/scotch_5.1 .
+or in SPECFEM3D/decompose_mesh_SCOTCH/Makefile set
+ SCOTCH_LIBS = -L/home/surendra/decompose_mesh_SCOTCH/scotch_5.1/lib/ -lscotch -lscotcherr
Submitting a job:
- See example script
- qsub
+ qsub -l nodes=48 -l walltime=30:00 -N simulation_name -m bae -M your at email.address go_mesher
More information about the CIG-COMMITS
mailing list