[cig-commits] r18684 - in seismo/3D/FAULT_SOURCE/branches/new_fault_db: EXAMPLES/splay_faults/DATA/FAULT EXAMPLES/splay_faults/post decompose_mesh_SCOTCH
percygalvez at geodynamics.org
percygalvez at geodynamics.org
Tue Jul 5 06:21:36 PDT 2011
Author: percygalvez
Date: 2011-07-05 06:21:35 -0700 (Tue, 05 Jul 2011)
New Revision: 18684
Added:
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/post/Plot_uz_surface.m
Modified:
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/decompose_mesh_SCOTCH/fault_scotch.f90
Log:
fault shilding at one node and input files for went et al w=50km, and scripts for mat
Modified: 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 2011-07-01 15:39:19 UTC (rev 18683)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults 2011-07-05 13:21:35 UTC (rev 18684)
@@ -2,37 +2,36 @@
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
--1e0 ! V_HEALING (-1 : Healing off)
-1e-3 ! V_RUPT: slip velocity threshold to define rupture front
+100 ! NTOUT : Number of time steps
+1000 ! NTSNAP: time interation of snapshots
+-1e0 ! V_HEALING (-1 : Healing off)
+1e-3 ! V_RUPT
&BEGIN_FAULT /
-&INIT_STRESS S1=0e6,S2=0e0,S3=-4.78e6,n1=0,n2=3,n3=3 /
-&DIST2D shape='rectangle', val = -2.16e6, xc = 113.5e3, yc =120e3, zc= -19.5e3, lx=227e3, ly=80e3, lz=30e3 /
-&DIST2D shape='rectangle', val = -2.88e6, xc = 166.521863e3, yc =148e3, zc= -7.96209773e3, lx=23.79468e3, ly=24e3, lz=3.132628e3 /
-&DIST2D shape='rectangle-taper', val = -2.16e6, valh =-0.216e6, xc = 227.872623e3, yc =120e3, zc= -4.5e3, lx=45.5745247e3, ly=80e3, lz=3e3 /
-&DIST2D shape='rectangle', val =-0.478e6, xc = 113.5e3, yc =120e3, zc= -15e3, lx=227e3, ly=80e3, lz=30e3 /
-&DIST2D shape='rectangle', val = -4.78e6, xc = 113.5e3, yc =120e3, zc= -19.5e3, lx=227e3, ly=80e3, lz=30e3 /
-&DIST2D shape='rectangle-taper', val = -4.78e6, valh=-0.478e6, xc = 226.407564e3, yc =120e3, zc= -4.5e3, lx=45.5745247e3, ly=80e3, lz=3e3 /
+&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=22787.2623e0, 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.50e3, lx=22787.2623e0, ly=80.0e3, lz=3.0e3 /
-&SWF mus=1000e0,mud=1000e0,dc=0.5e0,nmus=1,nmud=1 /
-&DIST2D shape='rectangle', val = 0.6e0, xc = 113.5e3, yc =140e3, zc= -15e3, lx=227e3, ly=40e3, lz=30e3 /
-&DIST2D shape='rectangle', val = 0.2e0, xc = 113.5e3, yc =140e3, zc= -15e3, lx=227e3, ly=40e3, lz=30e3 /
+&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 =135.0e3, zc= -15.0e3, lx=227.0e3, ly=50.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = 0.2e0, xc = 113.5e3, yc =135.0e3, zc= -15.0e3, lx=227.0e3, ly=50.0e3, lz=30000.0e0 /
-
&BEGIN_FAULT /
-&INIT_STRESS S1=0e6,S2=0e0,S3=-100e6,n1=0,n2=3,n3=4 /
-&DIST2D shape='rectangle', val = -5.69e6, xc = 70.5144284e3, yc =120e3, zc=-11.25e3, lx=23.3826859e3, ly=80e3, lz=13.5e3 /
-&DIST2D shape='rectangle', val = -3.54e6, xc = 43.8948822e3, yc =120e3, zc=-22e3, lx=29.8564065e3, ly=80e3, lz= 8e3 /
-&DIST2D shape='rectangle-taper', val = -5.69e6, valh=-0.569e6, xc = 55.8191065e3, yc =120e3, zc= -4.5e3, lx=45.5745247e3, ly=80e3, lz=3e3 /
-&DIST2D shape='rectangle', val =-1.238e6, xc = 74.4115427e3, yc =120e3, zc= -9e3, lx=31.1769145e3, ly=80e3, lz=18e3 /
-&DIST2D shape='rectangle', val =-12.38e6, xc = 70.5144284e3, yc =120e3, zc=-11.25e3, lx=23.3826859e3, ly=80e3, lz=13.5e3 /
-&DIST2D shape='rectangle-taper', val =-12.38e6, valh=-1.238e6, xc = 55.8191065e3, yc =120e3, zc= -4.5e3, lx=45.5745247e3, ly=80e3, lz=3e3 /
-&DIST2D shape='rectangle', val = -9.04e6, xc = 43.8948822e3, yc =120e3, zc=-22e3, lx=29.8564065e3, ly=80e3, lz= 8e3 /
+&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 = 82205.7714e0, yc =120.0e3, zc= -4.5e3, lx=5196.15242e0, 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 = 82205.7714e0, yc =120.0e3, zc= -4.5e3, lx=5196.15242e0, 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=1000e0,dc=0.5e0,nmus=1,nmud=1 /
-&DIST2D shape='rectangle', val = 0.6e0, xc = 60.1915e3, yc =120e3, zc=-13e3, lx=59.61698e3, ly=80e3, lz=26e3 /
-&DIST2D shape='rectangle', val = 0.2e0, xc = 60.1915e3, yc =120e3, zc=-13e3, lx=59.61698e3, ly=80e3, lz=26e3 /
+&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 /
Modified: 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.barrier 2011-07-01 15:39:19 UTC (rev 18683)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.barrier 2011-07-05 13:21:35 UTC (rev 18684)
@@ -2,37 +2,36 @@
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
--1e0 ! V_HEALING (-1 : Healing off)
-1e-3 ! V_RUPT: slip velocity threshold to define rupture front
+100 ! NTOUT : Number of time steps
+1000 ! NTSNAP: time interation of snapshots
+-1e0 ! V_HEALING (-1 : Healing off)
+1e-3 ! V_RUPT
&BEGIN_FAULT /
-&INIT_STRESS S1=0e6,S2=0e0,S3=-4.78e6,n1=0,n2=3,n3=3 /
-&DIST2D shape='rectangle', val = -2.16e6, xc = 113.5e3, yc =120e3, zc= -19.5e3, lx=227e3, ly=80e3, lz=30e3 /
-&DIST2D shape='rectangle', val = -2.88e6, xc = 166.521863e3, yc =148e3, zc= -7.96209773e3, lx=23.79468e3, ly=24e3, lz=3.132628e3 /
-&DIST2D shape='rectangle-taper', val = -2.16e6, valh =-0.216e6, xc = 226.407564e3, yc =120e3, zc= -4.5e3, lx=45.5745247e3, ly=80e3, lz=3e3 /
-&DIST2D shape='rectangle', val =-0.478e6, xc = 113.5e3, yc =120e3, zc= -15e3, lx=227e3, ly=80e3, lz=30e3 /
-&DIST2D shape='rectangle', val = -4.78e6, xc = 113.5e3, yc =120e3, zc= -19.5e3, lx=227e3, ly=80e3, lz=30e3 /
-&DIST2D shape='rectangle-taper', val = -4.78e6, valh=-0.478e6, xc = 226.407564e3, yc =120e3, zc= -4.5e3, lx=45.5745247e3, ly=80e3, lz=3e3 /
+&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=22787.2623e0, 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.50e3, lx=22787.2623e0, ly=80.0e3, lz=3.0e3 /
-&SWF mus=1000e0,mud=1000e0,dc=0.5e0,nmus=1,nmud=1 /
-&DIST2D shape='rectangle', val = 0.6e0, xc = 113.5e3, yc =140e3, zc= -15e3, lx=227e3, ly=40e3, lz=30e3 /
-&DIST2D shape='rectangle', val = 0.2e0, xc = 113.5e3, yc =140e3, zc= -15e3, lx=227e3, ly=40e3, lz=30e3 /
+&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 =135.0e3, zc= -15.0e3, lx=227.0e3, ly=50.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = 0.2e0, xc = 113.5e3, yc =135.0e3, zc= -15.0e3, lx=227.0e3, ly=50.0e3, lz=30000.0e0 /
-
&BEGIN_FAULT /
-&INIT_STRESS S1=0e6,S2=0e0,S3=-100e6,n1=0,n2=3,n3=4 /
-&DIST2D shape='rectangle', val = -5.69e6, xc = 70.5144284e3, yc =120e3, zc=-11.25e3, lx=23.3826859e3, ly=80e3, lz=13.5e3 /
-&DIST2D shape='rectangle', val = -3.54e6, xc = 43.8948822e3, yc =120e3, zc=-22e3, lx=29.8564065e3, ly=80e3, lz= 8e3 /
-&DIST2D shape='rectangle-taper', val = -5.69e6, valh=-0.569e6, xc = 55.8191065e3, yc =120e3, zc= -4.5e3, lx=45.5745247e3, ly=80e3, lz=3e3 /
-&DIST2D shape='rectangle', val =-1.238e6, xc = 74.4115427e3, yc =120e3, zc= -9e3, lx=31.1769145e3, ly=80e3, lz=18e3 /
-&DIST2D shape='rectangle', val =-12.38e6, xc = 70.5144284e3, yc =120e3, zc=-11.25e3, lx=23.3826859e3, ly=80e3, lz=13.5e3 /
-&DIST2D shape='rectangle-taper', val =-12.38e6, valh=-1.238e6, xc = 55.8191065e3, yc =120e3, zc= -4.5e3, lx=45.5745247e3, ly=80e3, lz=3e3 /
-&DIST2D shape='rectangle', val = -9.04e6, xc = 43.8948822e3, yc =120e3, zc=-22e3, lx=29.8564065e3, ly=80e3, lz= 8e3 /
+&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 = 82205.7714e0, yc =120.0e3, zc= -4.5e3, lx=5196.15242e0, 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 = 82205.7714e0, yc =120.0e3, zc= -4.5e3, lx=5196.15242e0, 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=1000e0,dc=0.5e0,nmus=1,nmud=1 /
-&DIST2D shape='rectangle', val = 0.6e0, xc = 60.1915e3, yc =120e3, zc=-13e3, lx=59.61698e3, ly=80e3, lz=26e3 /
-&DIST2D shape='rectangle', val = 0.2e0, xc = 60.1915e3, yc =120e3, zc=-13e3, lx=59.61698e3, ly=80e3, lz=26e3 /
+&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 /
Modified: 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/DATA/FAULT/Par_file_faults.homogeneous 2011-07-01 15:39:19 UTC (rev 18683)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.homogeneous 2011-07-05 13:21:35 UTC (rev 18684)
@@ -2,37 +2,36 @@
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
--1e0 ! V_HEALING (-1 : Healing off)
-1e-3 ! V_RUPT: slip velocity threshold to define rupture front
+100 ! NTOUT : Number of time steps
+1000 ! NTSNAP: time interation of snapshots
+-1e0 ! V_HEALING (-1 : Healing off)
+1e-3 ! V_RUPT
&BEGIN_FAULT /
-&INIT_STRESS S1=0e6,S2=0e0,S3=-4.78e6,n1=0,n2=3,n3=3 /
-&DIST2D shape='rectangle', val = -2.16e6, xc = 113.5e3, yc =120e3, zc= -19.5e3, lx=227e3, ly=80e3, lz=30e3 /
-&DIST2D shape='rectangle', val = -2.88e6, xc = 166.521863e3, yc =148e3, zc= -7.96209773e3, lx=23.79468e3, ly=24e3, lz=3.132628e3 /
-&DIST2D shape='rectangle-taper', val = -2.16e6, valh =-0.216e6, xc = 226.407564e3, yc =120e3, zc= -4.5e3, lx=45.5745247e3, ly=80e3, lz=3e3 /
-&DIST2D shape='rectangle', val =-0.478e6, xc = 113.5e3, yc =120e3, zc= -15e3, lx=227e3, ly=80e3, lz=30e3 /
-&DIST2D shape='rectangle', val = -4.78e6, xc = 113.5e3, yc =120e3, zc= -19.5e3, lx=227e3, ly=80e3, lz=30e3 /
-&DIST2D shape='rectangle-taper', val = -4.78e6, valh=-0.478e6, xc = 226.407564e3, yc =120e3, zc= -4.5e3, lx=45.5745247e3, ly=80e3, lz=3e3 /
+&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=22787.2623e0, 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.50e3, lx=22787.2623e0, ly=80.0e3, lz=3.0e3 /
-&SWF mus=1000e0,mud=1000e0,dc=0.5e0,nmus=1,nmud=1 /
-&DIST2D shape='rectangle', val = 0.6e0, xc = 113.5e3, yc =140e3, zc= -15e3, lx=227e3, ly=80e3, lz=30e3 /
-&DIST2D shape='rectangle', val = 0.2e0, xc = 113.5e3, yc =140e3, zc= -15e3, lx=227e3, ly=80e3, lz=30e3 /
+&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=0e6,S2=0e0,S3=-100e6,n1=0,n2=3,n3=4 /
-&DIST2D shape='rectangle', val = -5.69e6, xc = 70.5144284e3, yc =120e3, zc=-11.25e3, lx=23.3826859e3, ly=80e3, lz=13.5e3 /
-&DIST2D shape='rectangle', val = -3.54e6, xc = 43.8948822e3, yc =120e3, zc=-22e3, lx=29.8564065e3, ly=80e3, lz= 8e3 /
-&DIST2D shape='rectangle-taper', val = -5.69e6, valh=-0.569e6, xc = 55.8191065e3, yc =120e3, zc= -4.5e3, lx=45.5745247e3, ly=80e3, lz=3e3 /
-&DIST2D shape='rectangle', val =-1.238e6, xc = 74.4115427e3, yc =120e3, zc= -9e3, lx=31.1769145e3, ly=80e3, lz=18e3 /
-&DIST2D shape='rectangle', val =-12.38e6, xc = 70.5144284e3, yc =120e3, zc=-11.25e3, lx=23.3826859e3, ly=80e3, lz=13.5e3 /
-&DIST2D shape='rectangle-taper', val =-12.38e6, valh=-1.238e6, xc = 55.8191065e3, yc =120e3, zc= -4.5e3, lx=45.5745247e3, ly=80e3, lz=3e3 /
-&DIST2D shape='rectangle', val = -9.04e6, xc = 43.8948822e3, yc =120e3, zc=-22e3, lx=29.8564065e3, ly=80e3, lz= 8e3 /
+&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 = 82205.7714e0, yc =120.0e3, zc= -4.5e3, lx=5196.15242e0, 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 = 82205.7714e0, yc =120.0e3, zc= -4.5e3, lx=5196.15242e0, 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=1000e0,dc=0.5e0,nmus=1,nmud=1 /
-&DIST2D shape='rectangle', val = 0.6e0, xc = 60.1915e3, yc =120e3, zc=-13e3, lx=59.61698e3, ly=80e3, lz=26e3 /
-&DIST2D shape='rectangle', val = 0.2e0, xc = 60.1915e3, yc =120e3, zc=-13e3, lx=59.61698e3, ly=80e3, lz=26e3 /
+&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 /
Added: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/post/Plot_uz_surface.m
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/post/Plot_uz_surface.m (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/post/Plot_uz_surface.m 2011-07-05 13:21:35 UTC (rev 18684)
@@ -0,0 +1,65 @@
+
+clear all;
+close all;
+
+%Output file located in ~/SPECFEM3D/OUTPUT_FILES
+BinFile = 'moviedata033000';
+par = 3;
+scale = 3;
+SLIM = [-scale scale];
+
+NDAT = 6 ;
+fid = fopen(BinFile);
+BinRead = fread(fid,[1,inf],'single')';
+fclose(fid);
+
+BinRead = reshape( BinRead(:),length(BinRead)/(NDAT),NDAT);
+BinRead = BinRead(2:end-1,:);
+
+ d.X = BinRead(:,1)/1e3; % in km
+ d.Y = BinRead(:,2)/1e3; % in km
+ d.Z = BinRead(:,3)/1e3; % in km
+ d.dX = BinRead(:,4);
+ d.dY = BinRead(:,5);
+ d.dZ = BinRead(:,6);
+
+switch 1
+ case(par == 1)
+ d.dD = d.dX ;
+ case(par == 2)
+ d.dD = d.dY ;
+ case(par == 3)
+ d.dD = d.dZ ;
+end
+
+XLIM = [0 120];
+YLIM = [0 300];
+% Shift to Wendt's coordinate system
+d.Y = 180-d.Y;
+d.X = d.X;
+
+ind = find(d.Z == 0);
+
+dxsurf = d.Y(ind);
+dysurf = d.X(ind);
+ddDsurf = d.dD(ind);
+
+indrange =find((dxsurf > 0 & dxsurf < 120) & (dysurf > -50 & dysurf < 300));
+
+x = dxsurf(indrange);
+y = dysurf(indrange);
+dz = ddDsurf(indrange);
+
+scatter3(x,y,dz,[],dz,'.');
+axis('tight')
+xlabel('Along strike (km)');
+ylabel('Down dip (km)');
+
+zlim([-1 1]);
+caxis([-scale scale]);
+
+az= 55 ; el=70;
+view([az el]);
+
+colorbar
+
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 2011-07-01 15:39:19 UTC (rev 18683)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH/fault_scotch.f90 2011-07-05 13:21:35 UTC (rev 18684)
@@ -287,6 +287,7 @@
!pgb: leaving more space for fault elements.
!jpa: But if the number of fault elements is much smaller than nproc_null
! we will end up with a very UNbalanced proc 0 !
+ !pgb: Solution , fault parellelization .
ipart=1
if (nproc > 1) then
do i = 1, nproc_null
@@ -297,9 +298,9 @@
end if
deallocate(elem_proc_null)
endif
- call mesh2dual_ncommonnodes_fault(nelmnts, nnodes, nsize, sup_neighbour, &
- elmnts, xadj, adjncy, nnodes_elmnts, &
- nodes_elmnts, max_neighbour, 4, esize)
+! call mesh2dual_ncommonnodes_fault(nelmnts, nnodes, nsize, sup_neighbour, &
+! elmnts, xadj, adjncy, nnodes_elmnts, &
+! nodes_elmnts, max_neighbour, 4, esize)
! coupled elements
! ---------------
@@ -315,52 +316,58 @@
! 1 2
! Allocating elements with double shield layer
- print *, "Fault shield double-layer :"
- do el = 0, nelmnts-1
- if ( is_on_fault(el+1) ) then
- part(el) = 0
- do k1 = xadj(el), xadj(el+1) - 1
- el_1 = adjncy(k1)
- part(el_1) = 0
- do k2 = xadj(el_1), xadj(el_1+1) - 1
- el_2 = adjncy(k2)
- part(el_2) = 0
- enddo
- enddo
- endif
- enddo
+!
+! =========== FAULT SHIELD double-layer =========== !
+! print *, "Fault shield double-layer :"
+! do el = 0, nelmnts-1
+! if ( is_on_fault(el+1) ) then
+! part(el) = 0
+! do k1 = xadj(el), xadj(el+1) - 1
+! el_1 = adjncy(k1)
+! part(el_1) = 0
+! do k2 = xadj(el_1), xadj(el_1+1) - 1
+! el_2 = adjncy(k2)
+! part(el_2) = 0
+! enddo
+! enddo
+! endif
+! enddo
+! ===================================================== !
-! list of elements per node
-! nnodes_elmnts(i) = number of elements containing node #i (i=0:nnodes-1)
-! nodes_elmnts(nsize*i:nsize*i+nnodes_elmnts(i)-1) = index of elements (starting at 0) containing node #i
-! nnodes_elmnts(:) = 0
-! nodes_elmnts(:) = 0
-! do i = 0, esize*nelmnts-1
-! nodes_elmnts(elmnts(i)*nsize+nnodes_elmnts(elmnts(i))) = i/esize
-! nnodes_elmnts(elmnts(i)) = nnodes_elmnts(elmnts(i)) + 1
-! end do
+! Fault zone layer = the set of elements that contain at least one fault node
+ print *, "Fault zone layer :"
- ! fault zone layer = the set of elements that contain at least one fault node
-! print *, "Fault zone layer :"
-! do iflt=1,size(faults)
-! do e=1,faults(iflt)%nspec
-! do k=1,4
+! List of elements per node
+! nnodes_elmnts(i) = number of elements containing node #i (i=0:nnodes-1)
+! nodes_elmnts(nsize*i:nsize*i+nnodes_elmnts(i)-1) = index of elements (starting at 0) containing node #i
+! nsize = maximun number of elements in a node.
+! esize = nodes per element.
-! inode = faults(iflt)%inodes1(k,e)-1 ! node index, starting at 0
-! k1 = nsize*inode
-! k2 = k1 + nnodes_elmts(inode) -1
-! part( nodes_elmnts(k1:k2) ) = 0
+ nnodes_elmnts(:) = 0
+ nodes_elmnts(:) = 0
-! inode = faults(iflt)%inodes2(k,e)-1
-! k1 = nsize*inode
-! k2 = k1 + nnodes_elmts(inode) -1
-! part( nodes_elmnts(k1:k2) ) = 0
+ do i = 0, esize*nelmnts-1
+ nodes_elmnts(elmnts(i)*nsize+nnodes_elmnts(elmnts(i))) = i/esize
+ nnodes_elmnts(elmnts(i)) = nnodes_elmnts(elmnts(i)) + 1
+ end do
-! end do
-! end do
-! end do
-! jpa: if the previous section works then clean up obsolete stuff
+ do iflt=1,size(faults)
+ do e=1,faults(iflt)%nspec
+ do k=1,4
+ inode = faults(iflt)%inodes1(k,e)-1 ! node index, starting at 0
+ k1 = nsize*inode
+ k2 = k1 + nnodes_elmnts(inode) -1
+ part( nodes_elmnts(k1:k2) ) = 0
+ inode = faults(iflt)%inodes2(k,e)-1
+ k1 = nsize*inode
+ k2 = k1 + nnodes_elmnts(inode) -1
+ part( nodes_elmnts(k1:k2) ) = 0
+
+ end do
+ end do
+ end do
+
nproc_null_final = count( part == 0 )
print *, nproc_null_final
More information about the CIG-COMMITS
mailing list