[cig-commits] r20952 - in seismo/3D/FAULT_SOURCE/branches/new_fault_db: . decompose_mesh_SCOTCH src
surendra at geodynamics.org
surendra at geodynamics.org
Fri Oct 26 20:29:54 PDT 2012
Author: surendra
Date: 2012-10-26 20:29:54 -0700 (Fri, 26 Oct 2012)
New Revision: 20952
Modified:
seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT
seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH/fault_scotch.f90
seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_dynamic.f90
seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_kinematic.f90
Log:
Removed the sub-directory FAULT within DATA
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT 2012-10-27 03:19:45 UTC (rev 20951)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT 2012-10-27 03:29:54 UTC (rev 20952)
@@ -228,14 +228,14 @@
DATA/STATIONS List of stations outside the fault (see manual page 23).
-DATA/FAULT/Par_file_faults contains parameters of the fault. The first part of this file
+DATA/Par_file_faults contains parameters of the fault. The first part of this file
has a strict format:
Line 1: Number of faults (NF)
Lines 2 to NF+1: Kelvin Voigt damping (in seconds) for each fault. (See below how to set this parameter)
Line NF+2: Type of simulation (1=dynamic , 2 = kinematic)
Line NF+3: Number of time steps between updates of the time series outputs at selected
- fault points (see DATA/FAULT/FAULT_STATIONS), usually a large number (100s or 1000s).
+ fault points (see DATA/FAULT_STATIONS), usually a large number (100s or 1000s).
Note that the sampling rate of the time series is usually much higher.
Line NF+4: Number of time steps between fault snapshot outputs (quantities at every fault
point exported at regular times), usually a large number (100s or 1000s).
@@ -332,14 +332,14 @@
and look for the "maximum suggested time step" in OUTPUT_FILES/output_mesher.txt.
-DATA/FAULT/FAULT_STATIONS Stations in the fault plane.
+DATA/FAULT_STATIONS Stations in the fault plane.
Line 1: number of stations.
Line 2 to end: 5 columns: X, Y, Z (-depth), station name, fault-id
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
+DATA/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.
@@ -381,7 +381,7 @@
Several output files are saved in ~/SPECFEM3D/OUTPUT_FILES:
-1. Seismograms for each station on the fault plane given in DATA/FAULT/FAULT_STATIONS.
+1. Seismograms for each station on the fault plane given in DATA/FAULT_STATIONS.
One output file is generated for each station, named after the station. The files
are ascii and start with a header (22 lines long) followed by a data block with the
following format, one line per time sample:
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH/fault_scotch.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH/fault_scotch.f90 2012-10-27 03:19:45 UTC (rev 20951)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH/fault_scotch.f90 2012-10-27 03:29:54 UTC (rev 20952)
@@ -33,7 +33,7 @@
character(len=256),intent(in) :: localpath_name
integer :: nbfaults, iflt, ier
- open(101,file='../DATA/FAULT/Par_file_faults',status='old',action='read',iostat=ier)
+ open(101,file='../DATA/Par_file_faults',status='old',action='read',iostat=ier)
if (ier==0) then
read(101,*) nbfaults
else
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_dynamic.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_dynamic.f90 2012-10-27 03:19:45 UTC (rev 20951)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_dynamic.f90 2012-10-27 03:29:54 UTC (rev 20952)
@@ -149,7 +149,7 @@
dummy_idfault = 0
- open(unit=IIN_PAR,file='DATA/FAULT/Par_file_faults',status='old',iostat=ier)
+ open(unit=IIN_PAR,file='DATA/Par_file_faults',status='old',iostat=ier)
if( ier /= 0 ) then
write(6,*) 'File Par_file_faults not found: assume no faults'
close(IIN_PAR)
@@ -289,7 +289,7 @@
integer :: relz_num,sub_relz_num, num_cell_str,num_cell_dip, hypo_cell_str,hypo_cell_dip
integer :: i
- open(unit=IIN_NUC,file='DATA/FAULT/input_file.txt',status='old',iostat=ier)
+ open(unit=IIN_NUC,file='DATA/input_file.txt',status='old',iostat=ier)
read(IIN_NUC,*) relz_num,sub_relz_num
read(IIN_NUC,*) num_cell_str,num_cell_dip,siz_str,siz_dip
read(IIN_NUC,*) hypo_cell_str,hypo_cell_dip,hypo_loc_str,hypo_loc_dip,rad_T_str,rad_T_dip
@@ -1016,7 +1016,7 @@
! requested coordinate
IIN = 251
- open(IIN,file='DATA/FAULT/FAULT_STATIONS',status='old',action='read',iostat=ier)
+ open(IIN,file='DATA/FAULT_STATIONS',status='old',action='read',iostat=ier)
read(IIN,*) np
DataT%npoin =0
do i=1,np
@@ -1031,7 +1031,7 @@
allocate(DataT%name(DataT%npoin))
allocate(DataT%dist(DataT%npoin)) !Surendra : for parallel fault
- open(IIN,file='DATA/FAULT/FAULT_STATIONS',status='old',action='read',iostat=ier)
+ open(IIN,file='DATA/FAULT_STATIONS',status='old',action='read',iostat=ier)
if( ier /= 0 ) stop 'error opening FAULT_STATIONS file'
read(IIN,*) np
k = 0
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_kinematic.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_kinematic.f90 2012-10-27 03:19:45 UTC (rev 20951)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_kinematic.f90 2012-10-27 03:29:54 UTC (rev 20952)
@@ -114,7 +114,7 @@
dummy_idfault = 0
- open(unit=IIN_PAR,file='DATA/FAULT/Par_file_faults',status='old',iostat=ier)
+ open(unit=IIN_PAR,file='DATA/Par_file_faults',status='old',iostat=ier)
if( ier /= 0 ) then
write(6,*) 'Have not found Par_file_faults: assume no faults'
return
@@ -413,7 +413,7 @@
! requested coordinate
IIN = 251
- open(IIN,file='DATA/FAULT/FAULT_STATIONS',status='old',action='read',iostat=ier)
+ open(IIN,file='DATA/FAULT_STATIONS',status='old',action='read',iostat=ier)
read(IIN,*) np
DataT%npoin =0
do i=1,np
@@ -428,7 +428,7 @@
allocate(DataT%name(DataT%npoin))
allocate(DataT%dist(DataT%npoin)) !Surendra : for parallel fault
- open(IIN,file='DATA/FAULT/FAULT_STATIONS',status='old',action='read',iostat=ier)
+ open(IIN,file='DATA/FAULT_STATIONS',status='old',action='read',iostat=ier)
if( ier /= 0 ) stop 'error opening FAULT_STATIONS file'
read(IIN,*) np
k = 0
More information about the CIG-COMMITS
mailing list