[cig-commits] r18680 - in seismo/3D/FAULT_SOURCE/branches/new_fault_db: EXAMPLES/splay_faults EXAMPLES/splay_faults/DATA EXAMPLES/splay_faults/DATA/FAULT EXAMPLES/splay_faults/post EXAMPLES/tpv5.crack/DATA EXAMPLES/tpv5.crack/DATA/FAULT Post-processing
ampuero at geodynamics.org
ampuero at geodynamics.org
Thu Jun 30 05:51:34 PDT 2011
Author: ampuero
Date: 2011-06-30 05:51:34 -0700 (Thu, 30 Jun 2011)
New Revision: 18680
Added:
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/FAULT_STATIONS
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.barrier
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.homogeneous
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/post/
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/post/Plot_snapshot_splay_faults.m
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/FAULT/FAULT_STATIONS
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/FAULT/Par_file_faults
Removed:
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/CMTSOLUTION.in
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/FAULT_STATIONS.in
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.barrier.in
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.homogeneous.in
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.in
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/Par_file.in
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/STATIONS.in
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/CMTSOLUTION.in
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/FAULT/FAULT_STATIONS.in
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/FAULT/Par_file_faults.in
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/Par_file.in
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/STATIONS.in
seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/Plot_snapshot_splay_faults.m
Log:
changed input file names, cleaned up example input files
Deleted: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/CMTSOLUTION.in
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/CMTSOLUTION.in 2011-06-30 07:50:49 UTC (rev 18679)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/CMTSOLUTION.in 2011-06-30 12:51:34 UTC (rev 18680)
@@ -1,13 +0,0 @@
-PDE 1999 01 01 00 00 00.00 1000000 1000000 -1000000 1 1 dynamic_rupture
-event name: dynamic_rupture
-time shift: 0.0000
-half duration: 1.0
-latitude: 1000000.0
-longitude: 1000000.0
-depth: -1000000.0
-Mrr: 1.258900e+0
-Mtt: 1.258900e+0
-Mpp: 1.258900e+0
-Mrt: 0.000000e0
-Mrp: 0.000000e0
-Mtp: 0.000000e0
Copied: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/FAULT_STATIONS (from rev 18677, seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/FAULT_STATIONS.in)
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/FAULT_STATIONS (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/FAULT_STATIONS 2011-06-30 12:51:34 UTC (rev 18680)
@@ -0,0 +1,4 @@
+3
+176521.863 128000.0 -7962.0977 splayfaultA 1
+ 43894.8822 120000.0 -22000.0 splayfaultB 2
+ 74411.5427 120000.0 -9000.0 splayfaultC 2
Deleted: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/FAULT_STATIONS.in
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/FAULT_STATIONS.in 2011-06-30 07:50:49 UTC (rev 18679)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/FAULT_STATIONS.in 2011-06-30 12:51:34 UTC (rev 18680)
@@ -1,4 +0,0 @@
-3
-176521.863 128000.0 -7962.0977 splayfaultA 1
- 43894.8822 120000.0 -22000.0 splayfaultB 2
- 74411.5427 120000.0 -9000.0 splayfaultC 2
Copied: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults (from rev 18679, seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.in)
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults 2011-06-30 12:51:34 UTC (rev 18680)
@@ -0,0 +1,38 @@
+2
+0.00134 ! #tag #eta(damping) FAULT A
+0.00134 ! #tag #eta(damping) FAULT BC
+1 ! 1 = dyn 2=kin
+100 ! NTOUT : number of time steps between seismogram updates
+1000 ! NTSNAP: number of time steps between snapshot outputs
+-1.0e0 ! V_HEALING (-1 : Healing off)
+1e-3 ! V_RUPT: slip velocity threshold to define rupture front
+
+&BEGIN_FAULT /
+
+&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-4.78e6,n1=0,n2=3,n3=3 /
+&DIST2D shape='rectangle', val = -2.16e6, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = -2.88e6, xc = 166521.863e0, yc =148.0e3, zc= -7962.09773e0, lx=23794.68e0, ly=24.0e3, lz=3132.628e0 /
+&DIST2D shape='rectangle-taper', val = -2.16e6, valh =-0.216e6, xc = 226407.564e0, yc =120.0e3, zc= -4.50e3, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val =-0.478e6, xc = 113.5e3, yc =120.0e3, zc= -15000.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = -4.78e6, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle-taper', val = -4.78e6, valh=-0.478e6, xc = 226407.564e0, yc =120.0e3, zc= -4.5e3, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+
+&SWF mus=1000.0e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1 /
+&DIST2D shape='rectangle', val = 0.6e0, xc = 113.5e3, yc =140.0e3, zc= -15.0e3, lx=227.0e3, ly=40.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = 0.2e0, xc = 113.5e3, yc =140.0e3, zc= -15.0e3, lx=227.0e3, ly=40.0e3, lz=30000.0e0 /
+
+
+&BEGIN_FAULT /
+
+&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-100.0e6,n1=0,n2=3,n3=4 /
+&DIST2D shape='rectangle', val = -5.69e6, valh=0.0e0, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
+&DIST2D shape='rectangle', val = -3.54e6, valh=0.0e0, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
+&DIST2D shape='rectangle-taper', val = -5.69e6, valh=-0.569e6, xc = 55819.1065e0, yc =120.0e3, zc= -4.5e3, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val =-1.238e6, xc = 74411.5427e0, yc =120.0e3, zc= -9.0e3, lx=31176.9145e0, ly=80.0e3, lz=18.0e3 /
+&DIST2D shape='rectangle', val =-12.38e6, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
+&DIST2D shape='rectangle-taper', val =-12.38e6, valh=-1.238e6, xc = 55819.1065e0, yc =120.0e3, zc= -4.5e3, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val = -9.04e6, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
+
+&SWF mus=1000.e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1 /
+&DIST2D shape='rectangle', val = 0.6e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3,r=0e0, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
+&DIST2D shape='rectangle', val = 0.2e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3,r=0e0, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
Copied: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.barrier (from rev 18679, seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.in)
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.barrier (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.barrier 2011-06-30 12:51:34 UTC (rev 18680)
@@ -0,0 +1,38 @@
+2
+0.00134 ! #tag #eta(damping) FAULT A
+0.00134 ! #tag #eta(damping) FAULT BC
+1 ! 1 = dyn 2=kin
+100 ! NTOUT : Number of time steps
+1000 ! NTSNAP: time interation of snapshots
+-1.0e0 ! V_HEALING (-1 : Healing off)
+1e-3 ! V_RUPT
+
+&BEGIN_FAULT /
+
+&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-4.78e6,n1=0,n2=3,n3=3 /
+&DIST2D shape='rectangle', val = -2.16e6, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = -2.88e6, xc = 166521.863e0, yc =148.0e3, zc= -7962.09773e0, lx=23794.68e0, ly=24.0e3, lz=3132.628e0 /
+&DIST2D shape='rectangle-taper', val = -2.16e6, valh =-0.216e6, xc = 226407.564e0, yc =120.0e3, zc= -4.50e3, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val =-0.478e6, xc = 113.5e3, yc =120.0e3, zc= -15000.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = -4.78e6, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle-taper', val = -4.78e6, valh=-0.478e6, xc = 226407.564e0, yc =120.0e3, zc= -4.5e3, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+
+&SWF mus=1000.0e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1 /
+&DIST2D shape='rectangle', val = 0.6e0, xc = 113.5e3, yc =140.0e3, zc= -15.0e3, lx=227.0e3, ly=40.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = 0.2e0, xc = 113.5e3, yc =140.0e3, zc= -15.0e3, lx=227.0e3, ly=40.0e3, lz=30000.0e0 /
+
+
+&BEGIN_FAULT /
+
+&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-100.0e6,n1=0,n2=3,n3=4 /
+&DIST2D shape='rectangle', val = -5.69e6, valh=0.0e0, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
+&DIST2D shape='rectangle', val = -3.54e6, valh=0.0e0, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
+&DIST2D shape='rectangle-taper', val = -5.69e6, valh=-0.569e6, xc = 55819.1065e0, yc =120.0e3, zc= -4.5e3, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val =-1.238e6, xc = 74411.5427e0, yc =120.0e3, zc= -9.0e3, lx=31176.9145e0, ly=80.0e3, lz=18.0e3 /
+&DIST2D shape='rectangle', val =-12.38e6, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
+&DIST2D shape='rectangle-taper', val =-12.38e6, valh=-1.238e6, xc = 55819.1065e0, yc =120.0e3, zc= -4.5e3, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val = -9.04e6, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
+
+&SWF mus=1000.e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1 /
+&DIST2D shape='rectangle', val = 0.6e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3,r=0e0, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
+&DIST2D shape='rectangle', val = 0.2e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3,r=0e0, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
Deleted: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.barrier.in
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.barrier.in 2011-06-30 07:50:49 UTC (rev 18679)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.barrier.in 2011-06-30 12:51:34 UTC (rev 18680)
@@ -1,32 +0,0 @@
-2
-0.00134 ! #tag #eta(damping) FAULT A
-0.00134 ! #tag #eta(damping) FAULT BC
-1 ! 1 = dyn 2=kin
-100 ! NTOUT : Number of time steps
-1000 ! NTSNAP: time interation of snapshots
--1.0e0 ! V_HEALING (-1 : Healing off)
-1e-3 ! V_RUPT
-&BEGIN_FAULT dummy_idfault=1 /
-&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-4.78e6,n1=0,n2=3,n3=3 /
-&DIST2D shape='rectangle', val = -2.16e6, valh =0.0e0, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, r=0e0, l=0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
-&DIST2D shape='rectangle', val = -2.88e6, valh =0.0e0, xc = 166521.863e0, yc =148.0e3, zc= -7962.09773e0, r=0e0, l=0.0e0, lx=23794.68e0, ly=24.0e3, lz=3132.628e0 /
-&DIST2D shape='rectangle-taper', val = -2.16e6, valh =-0.216e6, xc = 226407.564e0, yc =120.0e3, zc= -4.50e3, r=0e0, l=0e0, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
-&DIST2D shape='rectangle', val =-0.478e6, valh =0.0e0, xc = 113.5e3, yc =120.0e3, zc= -15000.0e0, r=0e0, l=0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
-&DIST2D shape='rectangle', val = -4.78e6, valh =0.0e0, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, r=0e0, l=0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
-&DIST2D shape='rectangle-taper', val = -4.78e6, valh=-0.478e6, xc = 226407.564e0, yc =120.0e3, zc= -4.5e3, r=0e0, l=0e0, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
-&SWF mus=1000.0e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1,ndc=0 /
-&DIST2D shape='rectangle', val = 0.6e0, valh = 0.0e0, xc = 113.5e3, yc =140.0e3, zc= -15.0e3, r=0e0, l=0e0, lx=227.0e3, ly=40.0e3, lz=30000.0e0 /
-&DIST2D shape='rectangle', val = 0.2e0, valh = 0.0e0, xc = 113.5e3, yc =140.0e3, zc= -15.0e3, r=0e0, l=0e0, lx=227.0e3, ly=40.0e3, lz=30000.0e0 /
-&BEGIN_FAULT dummy_idfault=2 /
-&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-100.0e6,n1=0,n2=3,n3=4 /
-&DIST2D shape='rectangle', val = -5.69e6, valh=0.0e0, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3,r=0e0, l=0e0, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
-&DIST2D shape='rectangle', val = -3.54e6, valh=0.0e0, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, r=0e0, l=0e0, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
-&DIST2D shape='rectangle-taper', val = -5.69e6, valh=-0.569e6, xc = 55819.1065e0, yc =120.0e3, zc= -4.5e3, r=0e0, l=0e0, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
-&DIST2D shape='rectangle', val =-1.238e6, valh=0.0e0, xc = 74411.5427e0, yc =120.0e3, zc= -9.0e3, r=0e0, l=0e0, lx=31176.9145e0, ly=80.0e3, lz=18.0e3 /
-&DIST2D shape='rectangle', val =-12.38e6, valh=0.0e0, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3,r=0e0, l=0e0, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
-&DIST2D shape='rectangle-taper', val =-12.38e6, valh=-1.238e6, xc = 55819.1065e0, yc =120.0e3, zc= -4.5e3, r=0e0, l=0e0, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
-&DIST2D shape='rectangle', val = -9.04e6, valh=0.0e0, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, r=0e0, l=0e0, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
-&SWF mus=1000.e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1,ndc=0 /
-&DIST2D shape='rectangle', val = 0.6e0, valh=0.0e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3,r=0e0, l=0e0, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
-&DIST2D shape='rectangle', val = 0.2e0, valh=0.0e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3,r=0e0, l=0e0, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
-&KINPAR kindt=10.0e-3
Copied: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.homogeneous (from rev 18677, seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.homogeneous.in)
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.homogeneous (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.homogeneous 2011-06-30 12:51:34 UTC (rev 18680)
@@ -0,0 +1,37 @@
+2
+0.00134 ! #tag #eta(damping) FAULT A
+0.00134 ! #tag #eta(damping) FAULT BC
+1 ! 1 = dyn 2=kin
+100 ! NTOUT : number of time steps between seismogram updates
+1000 ! NTSNAP: number of time steps between snapshot outputs
+-1.0e0 ! V_HEALING (-1 : Healing off)
+1e-3 ! V_RUPT: slip velocity threshold to define rupture front
+
+&BEGIN_FAULT /
+
+&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-4.78e6,n1=0,n2=3,n3=3 /
+&DIST2D shape='rectangle', val = -2.16e6, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = -2.88e6, xc = 166521.863e0, yc =148.0e3, zc= -7962.09773e0, lx=23794.68e0, ly=24.0e3, lz=3132.628e0 /
+&DIST2D shape='rectangle-taper', val = -2.16e6, valh =-0.216e6, xc = 226407.564e0, yc =120.0e3, zc= -4.50e3, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val =-0.478e6, xc = 113.5e3, yc =120.0e3, zc= -15000.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = -4.78e6, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle-taper', val = -4.78e6, valh=-0.478e6, xc = 226407.564e0, yc =120.0e3, zc= -4.5e3, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+
+&SWF mus=1000.0e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1,ndc=0 /
+&DIST2D shape='rectangle', val = 0.6e0, xc = 113.5e3, yc =120.0e3, zc= -15.0e3, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = 0.2e0, xc = 113.5e3, yc =120.0e3, zc= -15.0e3, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+
+&BEGIN_FAULT /
+
+&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-100.0e6,n1=0,n2=3,n3=4 /
+&DIST2D shape='rectangle', val = -5.69e6, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
+&DIST2D shape='rectangle', val = -3.54e6, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
+&DIST2D shape='rectangle-taper', val = -5.69e6, valh=-0.569e6, xc = 55819.1065e0, yc =120.0e3, zc= -4.5e3, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val =-1.238e6, xc = 74411.5427e0, yc =120.0e3, zc= -9.0e3, lx=31176.9145e0, ly=80.0e3, lz=18.0e3 /
+&DIST2D shape='rectangle', val =-12.38e6, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
+&DIST2D shape='rectangle-taper', val =-12.38e6, valh=-1.238e6, xc = 55819.1065e0, yc =120.0e3, zc= -4.5e3, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val = -9.04e6, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
+
+&SWF mus=1000.e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1,ndc=0 /
+&DIST2D shape='rectangle', val = 0.6e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
+&DIST2D shape='rectangle', val = 0.2e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
Deleted: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.homogeneous.in
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.homogeneous.in 2011-06-30 07:50:49 UTC (rev 18679)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.homogeneous.in 2011-06-30 12:51:34 UTC (rev 18680)
@@ -1,32 +0,0 @@
-2
-0.00134 ! #tag #eta(damping) FAULT A
-0.00134 ! #tag #eta(damping) FAULT BC
-1 ! 1 = dyn 2=kin
-100 ! NTOUT : Number of time steps
-1000 ! NTSNAP: time interation of snapshots
--1.0e0 ! V_HEALING (-1 : Healing off)
-1e-3 ! V_RUPT
-&BEGIN_FAULT dummy_idfault=1 /
-&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-4.78e6,n1=0,n2=3,n3=3 /
-&DIST2D shape='rectangle', val = -2.16e6, valh =0.0e0, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, r=0e0, l=0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
-&DIST2D shape='rectangle', val = -2.88e6, valh =0.0e0, xc = 166521.863e0, yc =148.0e3, zc= -7962.09773e0, r=0e0, l=0.0e0, lx=23794.68e0, ly=24.0e3, lz=3132.628e0 /
-&DIST2D shape='rectangle-taper', val = -2.16e6, valh =-0.216e6, xc = 226407.564e0, yc =120.0e3, zc= -4.50e3, r=0e0, l=0e0, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
-&DIST2D shape='rectangle', val =-0.478e6, valh =0.0e0, xc = 113.5e3, yc =120.0e3, zc= -15000.0e0, r=0e0, l=0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
-&DIST2D shape='rectangle', val = -4.78e6, valh =0.0e0, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, r=0e0, l=0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
-&DIST2D shape='rectangle-taper', val = -4.78e6, valh=-0.478e6, xc = 226407.564e0, yc =120.0e3, zc= -4.5e3, r=0e0, l=0e0, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
-&SWF mus=1000.0e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1,ndc=0 /
-&DIST2D shape='rectangle', val = 0.6e0, valh = 0.0e0, xc = 113.5e3, yc =120.0e3, zc= -15.0e3, r=0e0, l=0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
-&DIST2D shape='rectangle', val = 0.2e0, valh = 0.0e0, xc = 113.5e3, yc =120.0e3, zc= -15.0e3, r=0e0, l=0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
-&BEGIN_FAULT dummy_idfault=2 /
-&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-100.0e6,n1=0,n2=3,n3=4 /
-&DIST2D shape='rectangle', val = -5.69e6, valh=0.0e0, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3,r=0e0, l=0e0, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
-&DIST2D shape='rectangle', val = -3.54e6, valh=0.0e0, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, r=0e0, l=0e0, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
-&DIST2D shape='rectangle-taper', val = -5.69e6, valh=-0.569e6, xc = 55819.1065e0, yc =120.0e3, zc= -4.5e3, r=0e0, l=0e0, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
-&DIST2D shape='rectangle', val =-1.238e6, valh=0.0e0, xc = 74411.5427e0, yc =120.0e3, zc= -9.0e3, r=0e0, l=0e0, lx=31176.9145e0, ly=80.0e3, lz=18.0e3 /
-&DIST2D shape='rectangle', val =-12.38e6, valh=0.0e0, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3,r=0e0, l=0e0, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
-&DIST2D shape='rectangle-taper', val =-12.38e6, valh=-1.238e6, xc = 55819.1065e0, yc =120.0e3, zc= -4.5e3, r=0e0, l=0e0, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
-&DIST2D shape='rectangle', val = -9.04e6, valh=0.0e0, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, r=0e0, l=0e0, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
-&SWF mus=1000.e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1,ndc=0 /
-&DIST2D shape='rectangle', val = 0.6e0, valh=0.0e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3,r=0e0, l=0e0, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
-&DIST2D shape='rectangle', val = 0.2e0, valh=0.0e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3,r=0e0, l=0e0, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
-&KINPAR kindt=10.0e-3
Deleted: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.in
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.in 2011-06-30 07:50:49 UTC (rev 18679)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.in 2011-06-30 12:51:34 UTC (rev 18680)
@@ -1,31 +0,0 @@
-2
-0.00134 ! #tag #eta(damping) FAULT A
-0.00134 ! #tag #eta(damping) FAULT BC
-1 ! 1 = dyn 2=kin
-100 ! NTOUT : Number of time steps
-1000 ! NTSNAP: time interation of snapshots
--1.0e0 ! V_HEALING (-1 : Healing off)
-1e-3 ! V_RUPT
-&BEGIN_FAULT /
-&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-4.78e6,n1=0,n2=3,n3=3 /
-&DIST2D shape='rectangle', val = -2.16e6, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
-&DIST2D shape='rectangle', val = -2.88e6, xc = 166521.863e0, yc =148.0e3, zc= -7962.09773e0, lx=23794.68e0, ly=24.0e3, lz=3132.628e0 /
-&DIST2D shape='rectangle-taper', val = -2.16e6, valh =-0.216e6, xc = 226407.564e0, yc =120.0e3, zc= -4.50e3, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
-&DIST2D shape='rectangle', val =-0.478e6, xc = 113.5e3, yc =120.0e3, zc= -15000.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
-&DIST2D shape='rectangle', val = -4.78e6, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
-&DIST2D shape='rectangle-taper', val = -4.78e6, valh=-0.478e6, xc = 226407.564e0, yc =120.0e3, zc= -4.5e3, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
-&SWF mus=1000.0e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1,ndc=0 /
-&DIST2D shape='rectangle', val = 0.6e0, xc = 113.5e3, yc =140.0e3, zc= -15.0e3, lx=227.0e3, ly=40.0e3, lz=30000.0e0 /
-&DIST2D shape='rectangle', val = 0.2e0, xc = 113.5e3, yc =140.0e3, zc= -15.0e3, lx=227.0e3, ly=40.0e3, lz=30000.0e0 /
-&BEGIN_FAULT /
-&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-100.0e6,n1=0,n2=3,n3=4 /
-&DIST2D shape='rectangle', val = -5.69e6, valh=0.0e0, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
-&DIST2D shape='rectangle', val = -3.54e6, valh=0.0e0, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
-&DIST2D shape='rectangle-taper', val = -5.69e6, valh=-0.569e6, xc = 55819.1065e0, yc =120.0e3, zc= -4.5e3, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
-&DIST2D shape='rectangle', val =-1.238e6, xc = 74411.5427e0, yc =120.0e3, zc= -9.0e3, lx=31176.9145e0, ly=80.0e3, lz=18.0e3 /
-&DIST2D shape='rectangle', val =-12.38e6, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
-&DIST2D shape='rectangle-taper', val =-12.38e6, valh=-1.238e6, xc = 55819.1065e0, yc =120.0e3, zc= -4.5e3, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
-&DIST2D shape='rectangle', val = -9.04e6, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
-&SWF mus=1000.e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1,ndc=0 /
-&DIST2D shape='rectangle', val = 0.6e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3,r=0e0, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
-&DIST2D shape='rectangle', val = 0.2e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3,r=0e0, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
Deleted: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/Par_file.in
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/Par_file.in 2011-06-30 07:50:49 UTC (rev 18679)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/Par_file.in 2011-06-30 12:51:34 UTC (rev 18680)
@@ -1,53 +0,0 @@
-
-# forward or adjoint simulation
-SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint, 3 = both simultaneously
-SAVE_FORWARD = .false.
-
-# UTM projection parameters
-UTM_PROJECTION_ZONE = 11
-SUPPRESS_UTM_PROJECTION = .true.
-
-# number of MPI processors
-NPROC = 100
-
-# time step parameters
-NSTEP = 62000
-DT = 0.00325d0
-
-# parameters describing the model
-OCEANS = .false.
-ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
-ANISOTROPY = .false.
-
-# absorbing boundary conditions for a regional simulation
-ABSORBING_CONDITIONS = .true.
-ABSORB_FREE_SURFACE = .false.
-# save AVS or OpenDX movies
-MOVIE_SURFACE = .true.
-MOVIE_VOLUME = .false.
-NTSTEP_BETWEEN_FRAMES = 16500
-CREATE_SHAKEMAP = .false.
-SAVE_DISPLACEMENT = .true.
-USE_HIGHRES_FOR_MOVIES = .false.
-HDUR_MOVIE = 0.0
-
-# save AVS or OpenDX mesh files to check the mesh
-SAVE_MESH_FILES = .false.
-
-# path to store the local database file on each node
-LOCAL_PATH = /cluster/work/erdw/gpercy
-
-# interval at which we output time step info and max of norm of displacement
-NTSTEP_BETWEEN_OUTPUT_INFO = 500
-
-# interval in time steps for writing of seismograms
-NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
-
-# print source time function
-PRINT_SOURCE_TIME_FUNCTION = .false.
-
-## Here defined tags that will be taken for faults.
-## which friction laws .
-
-
Deleted: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/STATIONS.in
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/STATIONS.in 2011-06-30 07:50:49 UTC (rev 18679)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/STATIONS.in 2011-06-30 12:51:34 UTC (rev 18680)
@@ -1,3 +0,0 @@
-splayfaultA_surf SC 128000.0 176521.863 0.0 0.0
-splayfaultB_surf SC 120000.0 43894.8822 0.0 0.0
-splayfaultc_surf SC 120000.0 74411.5427 0.0 0.0
Copied: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/post/Plot_snapshot_splay_faults.m (from rev 18679, seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/Plot_snapshot_splay_faults.m)
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/post/Plot_snapshot_splay_faults.m (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/post/Plot_snapshot_splay_faults.m 2011-06-30 12:51:34 UTC (rev 18680)
@@ -0,0 +1,64 @@
+addpath ~/SPECFEM3D/Post-processing
+
+clear all
+close all
+ITSNAP = 33000;
+XLIM = [0 120]; % range along-strike
+SLIM = [-8 8]; % slip range
+
+dA = FSEM3D_snapshot(ITSNAP,1,'.',1);
+dBC = FSEM3D_snapshot(ITSNAP,1,'.',2);
+
+% along-dip slip is negative for thrust faulting
+dA.Dz = -dA.Dz;
+dBC.Dz = -dBC.Dz;
+
+% shift to Wendt's coordinate system
+dA.Y = 180-dA.Y;
+dA.X = 227-dA.X;
+dBC.Y = 180-dBC.Y;
+dBC.X = 90-dBC.X;
+
+figure(1)
+subplot(221)
+plot(dA.Y,dA.Dx,'.')
+set(gca,'XLim',XLIM)
+ylabel('Dx (m)')
+subplot(223)
+plot(dA.Y,dA.Dz,'.')
+set(gca,'XLim',XLIM)
+ylabel('Dz (m)')
+xlabel('Y along-strike (km)')
+subplot(222)
+plot(dBC.Y,dBC.Dx,'.')
+set(gca,'XLim',XLIM)
+subplot(224)
+plot(dBC.Y,dBC.Dz,'.')
+set(gca,'XLim',XLIM)
+
+figure(2)
+subplot(121)
+[hc,hf]=plotclr(dA.Y,dA.X,dA.Dz,'s',SLIM);
+set(hf,'YDir','reverse')
+xlabel('Along strike (km)')
+ylabel('Down dip (km)')
+axis equal
+axis tight
+set(hf,'XLim',XLIM)
+
+subplot(122)
+[hc,hf]=plotclr(dBC.Y,dBC.X,dBC.Dz,'s',SLIM);
+set(hf,'YDir','reverse')
+xlabel('Along strike (km)')
+ylabel('Down dip (km)')
+axis equal
+axis tight
+set(hf,'XLim',XLIM)
+
+% This is the colormap of Wendt et al, but it does not work well yet with plotclr:
+%colormap('jet');
+%map1 = colormap;
+%l1 = length(map1);
+%map2 = map1(l1/2+1:end,:);
+%colormap(map2);
+
Deleted: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/CMTSOLUTION.in
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/CMTSOLUTION.in 2011-06-30 07:50:49 UTC (rev 18679)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/CMTSOLUTION.in 2011-06-30 12:51:34 UTC (rev 18680)
@@ -1,13 +0,0 @@
-PDE 1999 01 01 00 00 00.00 1000000 1000000 -1000000 1 1 dynamic_rupture
-event name: dynamic_rupture
-time shift: 0.0000
-half duration: 0.0
-latitude: 1000000.0
-longitude: 1000000.0
-depth: -1000000.0
-Mrr: 1.258900e+0
-Mtt: 1.258900e+0
-Mpp: 1.258900e+0
-Mrt: 0.000000e0
-Mrp: 0.000000e0
-Mtp: 0.000000e0
Copied: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/FAULT/FAULT_STATIONS (from rev 18677, seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/FAULT/FAULT_STATIONS.in)
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/FAULT/FAULT_STATIONS (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/FAULT/FAULT_STATIONS 2011-06-30 12:51:34 UTC (rev 18680)
@@ -0,0 +1,17 @@
+16
+ -4500.0 0.0 0.0 faultst-045dp000 1
+ -4500.0 0.0 -7500.0 faultst-045dp075 1
+ -7500.0 0.0 0.0 faultst-075dp000 1
+ -7500.0 0.0 -7500.0 faultst-075dp075 1
+-12000.0 0.0 0.0 faultst-120dp000 1
+-12000.0 0.0 -7500.0 faultst-120dp075 1
+ 0.0 0.0 0.0 faultst000dp000 1
+ 0.0 0.0 -4500.0 faultst000dp045 1
+ 0.0 0.0 -7500.0 faultst000dp075 1
+ 0.0 0.0 -12000.0 faultst000dp120 1
+ 4500.0 0.0 0.0 faultst045dp000 1
+ 4500.0 0.0 -7500.0 faultst045dp075 1
+ 7500.0 0.0 0.0 faultst075dp000 1
+ 7500.0 0.0 -7500.0 faultst075dp075 1
+ 12000.0 0.0 0.0 faultst120dp000 1
+ 12000.0 0.0 -7500.0 faultst120dp075 1
Deleted: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/FAULT/FAULT_STATIONS.in
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/FAULT/FAULT_STATIONS.in 2011-06-30 07:50:49 UTC (rev 18679)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/FAULT/FAULT_STATIONS.in 2011-06-30 12:51:34 UTC (rev 18680)
@@ -1,17 +0,0 @@
-16
- -4500.0 0.0 0.0 faultst-045dp000 1
- -4500.0 0.0 -7500.0 faultst-045dp075 1
- -7500.0 0.0 0.0 faultst-075dp000 1
- -7500.0 0.0 -7500.0 faultst-075dp075 1
--12000.0 0.0 0.0 faultst-120dp000 1
--12000.0 0.0 -7500.0 faultst-120dp075 1
- 0.0 0.0 0.0 faultst000dp000 1
- 0.0 0.0 -4500.0 faultst000dp045 1
- 0.0 0.0 -7500.0 faultst000dp075 1
- 0.0 0.0 -12000.0 faultst000dp120 1
- 4500.0 0.0 0.0 faultst045dp000 1
- 4500.0 0.0 -7500.0 faultst045dp075 1
- 7500.0 0.0 0.0 faultst075dp000 1
- 7500.0 0.0 -7500.0 faultst075dp075 1
- 12000.0 0.0 0.0 faultst120dp000 1
- 12000.0 0.0 -7500.0 faultst120dp075 1
Copied: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/FAULT/Par_file_faults (from rev 18677, seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/FAULT/Par_file_faults.in)
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/FAULT/Par_file_faults (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/FAULT/Par_file_faults 2011-06-30 12:51:34 UTC (rev 18680)
@@ -0,0 +1,14 @@
+1
+0.00134 ! #tag #eta(damping)
+1 ! 1 = dyn 2=kin
+1000 ! NTOUT : Number of time steps
+100 ! NTSNAP: time interation of snapshots
+-1.0e0 ! V_HEALING (-1 : Healing off)
+1e-3 ! V_RUPT
+&BEGIN_FAULT /
+&INIT_STRESS S1=70.0e6, n1=3, S2=0.0e0,S3=-120.0e6 /
+&DIST2D shape='square', val = 78.0e6, xc = -7500.0e0, yc =0e0, zc= -7500.0e0, l=3000.0e0 /
+&DIST2D shape='square', val = 81.6e6, xc = 0e0, yc =0e0, zc= -7500.0e0, l=3000.0e0 /
+&DIST2D shape='square', val = 62.0e6, xc = 7500.0e0, yc =0e0, zc= -7500.0e0, l=3000.0e0 /
+&SWF mus=10000.0e0, nmus=1, mud=0.525e0,dc=0.4e0/
+&DIST2D shape='rectangle', val = 0.677e0, xc = 0e0, yc =0e0, zc= -7500.0e0, lx=30000.0e0, ly=1.0e0, lz=15000.0e0 /
Deleted: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/FAULT/Par_file_faults.in
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/FAULT/Par_file_faults.in 2011-06-30 07:50:49 UTC (rev 18679)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/FAULT/Par_file_faults.in 2011-06-30 12:51:34 UTC (rev 18680)
@@ -1,15 +0,0 @@
-1
-0.00134 ! #tag #eta(damping)
-1 ! 1 = dyn 2=kin
-1000 ! NTOUT : Number of time steps
-100 ! NTSNAP: time interation of snapshots
--1.0e0 ! V_HEALING (-1 : Healing off)
-1e-3 ! V_RUPT
-&BEGIN_FAULT dummy_idfault=1 /
-&INIT_STRESS S1=70.0e6,S2=0.0e0,S3=-120.0e6,n1=3,n2=0,n3=0 /
-&DIST2D shape='square', val = 78.0e6, xc = -7500.0e0, yc =0e0, zc= -7500.0e0, r=0e0, l=3000.0e0, lx=0e0, ly=0e0, lz=0e0 /
-&DIST2D shape='square', val = 81.6e6, xc = 0e0, yc =0e0, zc= -7500.0e0, r=0e0, l=3000.0e0, lx=0e0, ly=0e0, lz=0e0 /
-&DIST2D shape='square', val = 62.0e6, xc = 7500.0e0, yc =0e0, zc= -7500.0e0, r=0e0, l=3000.0e0, lx=0e0, ly=0e0, lz=0e0 /
-&SWF mus=10000.0e0,mud=0.525e0,dc=0.4e0,nmus=1,nmud=0,ndc=0 /
-&DIST2D shape='rectangle', val = 0.677e0, xc = 0e0, yc =0e0, zc= -7500.0e0, r=0e0, l=0e0, lx=30000.0e0, ly=0.0e0, lz=15000.0e0 /
-&KINPAR kindt=10.0e-3
Deleted: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/Par_file.in
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/Par_file.in 2011-06-30 07:50:49 UTC (rev 18679)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/Par_file.in 2011-06-30 12:51:34 UTC (rev 18680)
@@ -1,53 +0,0 @@
-
-# forward or adjoint simulation
-SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint, 3 = both simultaneously
-SAVE_FORWARD = .false.
-
-# UTM projection parameters
-UTM_PROJECTION_ZONE = 11
-SUPPRESS_UTM_PROJECTION = .true.
-
-# number of MPI processors
-NPROC = 100
-
-# time step parameters
-NSTEP = 12000
-DT = 0.001d0
-
-# parameters describing the model
-OCEANS = .false.
-ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
-ANISOTROPY = .false.
-
-# absorbing boundary conditions for a regional simulation
-ABSORBING_CONDITIONS = .true.
-ABSORB_FREE_SURFACE = .false.
-# save AVS or OpenDX movies
-MOVIE_SURFACE = .false.
-MOVIE_VOLUME = .false.
-NTSTEP_BETWEEN_FRAMES = 200
-CREATE_SHAKEMAP = .false.
-SAVE_DISPLACEMENT = .false.
-USE_HIGHRES_FOR_MOVIES = .false.
-HDUR_MOVIE = 0.0
-
-# save AVS or OpenDX mesh files to check the mesh
-SAVE_MESH_FILES = .false.
-
-# path to store the local database file on each node
-LOCAL_PATH = /home/galvez/FAULT_SOURCE/DATABASES_MPI
-
-# interval at which we output time step info and max of norm of displacement
-NTSTEP_BETWEEN_OUTPUT_INFO = 500
-
-# interval in time steps for writing of seismograms
-NTSTEP_BETWEEN_OUTPUT_SEISMOS = 10000
-
-# print source time function
-PRINT_SOURCE_TIME_FUNCTION = .false.
-
-## Here defined tags that will be taken for faults.
-## which friction laws .
-
-
Deleted: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/STATIONS.in
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/STATIONS.in 2011-06-30 07:50:49 UTC (rev 18679)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/DATA/STATIONS.in 2011-06-30 12:51:34 UTC (rev 18680)
@@ -1,4 +0,0 @@
-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
Deleted: seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/Plot_snapshot_splay_faults.m
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/Plot_snapshot_splay_faults.m 2011-06-30 07:50:49 UTC (rev 18679)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/Plot_snapshot_splay_faults.m 2011-06-30 12:51:34 UTC (rev 18680)
@@ -1,62 +0,0 @@
-clear all
-close all
-ITSNAP = 33000;
-XLIM = [0 120]; % range along-strike
-SLIM = [-8 8]; % slip range
-
-dA = FSEM3D_snapshot(ITSNAP,1,'.',1);
-dBC = FSEM3D_snapshot(ITSNAP,1,'.',2);
-
-% along-dip slip is negative for thrust faulting
-dA.Dz = -dA.Dz;
-dBC.Dz = -dBC.Dz;
-
-% shift to Wendt's coordinate system
-dA.Y = 180-dA.Y;
-dA.X = 227-dA.X;
-dBC.Y = 180-dBC.Y;
-dBC.X = 90-dBC.X;
-
-figure(1)
-subplot(221)
-plot(dA.Y,dA.Dx,'.')
-set(gca,'XLim',XLIM)
-ylabel('Dx (m)')
-subplot(223)
-plot(dA.Y,dA.Dz,'.')
-set(gca,'XLim',XLIM)
-ylabel('Dz (m)')
-xlabel('Y along-strike (km)')
-subplot(222)
-plot(dBC.Y,dBC.Dx,'.')
-set(gca,'XLim',XLIM)
-subplot(224)
-plot(dBC.Y,dBC.Dz,'.')
-set(gca,'XLim',XLIM)
-
-figure(2)
-subplot(121)
-[hc,hf]=plotclr(dA.Y,dA.X,dA.Dz,'s',SLIM);
-set(hf,'YDir','reverse')
-xlabel('Along strike (km)')
-ylabel('Down dip (km)')
-axis equal
-axis tight
-set(hf,'XLim',XLIM)
-
-subplot(122)
-[hc,hf]=plotclr(dBC.Y,dBC.X,dBC.Dz,'s',SLIM);
-set(hf,'YDir','reverse')
-xlabel('Along strike (km)')
-ylabel('Down dip (km)')
-axis equal
-axis tight
-set(hf,'XLim',XLIM)
-
-% This is the colormap of Wendt et al, but it does not work well yet with plotclr:
-%colormap('jet');
-%map1 = colormap;
-%l1 = length(map1);
-%map2 = map1(l1/2+1:end,:);
-%colormap(map2);
-
More information about the CIG-COMMITS
mailing list