[cig-commits] r20852 - in seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103: . DATA DATA/FAULT
surendra at geodynamics.org
surendra at geodynamics.org
Thu Oct 18 12:42:01 PDT 2012
Author: surendra
Date: 2012-10-18 12:42:01 -0700 (Thu, 18 Oct 2012)
New Revision: 20852
Added:
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/CMTSOLUTION
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/FAULT/FAULT_STATIONS
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/FAULT/Par_file_faults
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/Par_file
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/STATIONS
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/TPV103.jou
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/TPV103.py
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/description.pdf
Log:
Added missing tpv103 files
Added: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/CMTSOLUTION
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/CMTSOLUTION (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/CMTSOLUTION 2012-10-18 19:42:01 UTC (rev 20852)
@@ -0,0 +1,13 @@
+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
Added: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/FAULT/FAULT_STATIONS
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/FAULT/FAULT_STATIONS (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/FAULT/FAULT_STATIONS 2012-10-18 19:42:01 UTC (rev 20852)
@@ -0,0 +1,11 @@
+9
+ 12000.0 0.0 -3000.0 faultst120dp030 1
+ 9000.0 0.0 -7500.0 faultst090dp075 1
+ 12000.0 0.0 -12000.0 faultst120dp120 1
+ -12000.0 0.0 -3000.0 faultst-120dp030 1
+ -9000.0 0.0 -7500.0 faultst-090dp075 1
+ -12000.0 0.0 -12000.0 faultst-120dp120 1
+ 0.0 0.0 -3000.0 faultst000dp030 1
+ 0.0 0.0 -7500.0 faultst000dp075 1
+ 0.0 0.0 -12000.0 faultst000dp120 1
+
Added: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/FAULT/Par_file_faults
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/FAULT/Par_file_faults (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/FAULT/Par_file_faults 2012-10-18 19:42:01 UTC (rev 20852)
@@ -0,0 +1,12 @@
+1
+0.00134 ! #tag #eta(damping) FAULT A
+1 ! 1 = dyn 2=kin
+1200 ! NTOUT : Number of time steps
+200 ! NTSNAP: time interation of snapshots
+-1e0 ! V_HEALING (-1 : Healing off)
+1e-3 ! V_RUPT
+&BEGIN_FAULT /
+&INIT_STRESS S1=40.0e6,S2=0.0e0,S3=-120.0e6,n1=0,n2=0,n3=0 /
+&RSF V0 =1.e-6, f0 =0.6, a =0.01, b =0.014, L =0.4, V_init =1.e-16, theta_init =0.563591842632738, C=0.e0, nV0 =0, nf0 =0, na=0, nb=0, nL=0, nV_init=0, ntheta_init=0 /
+&ASP Fload =0.e0, nFload =1 /
+&DIST2D shape='circle-exp', val =1.e0, xc = 0.e3, yc =0.e3, zc = -7.5e3, r =3.e3 /
Added: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/Par_file
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/Par_file (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/Par_file 2012-10-18 19:42:01 UTC (rev 20852)
@@ -0,0 +1,53 @@
+
+# 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 = 192
+
+# time step parameters
+NSTEP = 4800
+DT = 0.0025d0
+
+# 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 = .true.
+# save AVS or OpenDX movies
+MOVIE_SURFACE = .false.
+MOVIE_VOLUME = .false.
+NTSTEP_BETWEEN_FRAMES = 100
+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 = /home/surendra/SCEC_benchmarks/TPV101_DATABASES_MPI_serial_Fault_EXT2ABS_larger
+
+# interval at which we output time step info and max of norm of displacement
+NTSTEP_BETWEEN_OUTPUT_INFO = 200
+
+# interval in time steps for writing of seismograms
+NTSTEP_BETWEEN_OUTPUT_SEISMOS = 20000
+
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION = .false.
+
+## Here defined tags that will be taken for faults.
+## which friction laws .
+
+
Added: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/STATIONS
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/STATIONS (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/DATA/STATIONS 2012-10-18 19:42:01 UTC (rev 20852)
@@ -0,0 +1,6 @@
+body-060s-120dp000 SF -12000.0 -6000.0 0.0 0.0
+body060s-120dp000 SF -12000.0 +6000.0 0.0 0.0
+body-090st000dp000 SF 0.0 -9000.0 0.0 0.0
+body090st000dp000 SF 0.0 +9000.0 0.0 0.0
+body-060st120dp000 SF +12000.0 -6000.0 0.0 0.0
+body060st120dp000 SF +12000.0 +6000.0 0.0 0.0
Added: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/TPV103.jou
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/TPV103.jou (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/TPV103.jou 2012-10-18 19:42:01 UTC (rev 20852)
@@ -0,0 +1,45 @@
+reset
+
+set echo off
+set journal off
+
+brick x 60000 y 60000 z 46000
+move volume 1 z -7000
+webcut volume 1 with plane yplane
+
+surface 12 7 merge off
+merge tolerance 1e-3
+
+imprint all
+merge all
+
+
+
+surface all size 400
+mesh volume all
+
+
+
+
+group "upp" add node in surface 12
+group "upp" remove node with x_coord = 30000
+group "upp" remove node with x_coord = -30000
+group "upp" remove node with z_coord = -30000
+group "upp" remove node with z_coord = 16000
+nodeset 200 group upp
+nodeset 200 move 0.01 0 0
+
+group "lowr" add node in surface 7
+group "lowr" remove node with x_coord = 30000
+group "lowr" remove node with x_coord = -30000
+group "lowr" remove node with z_coord = -30000
+group "lowr" remove node with z_coord = 16000
+nodeset 201 group lowr
+nodeset 201 move -0.01 0 0
+
+
+
+
+
+
+
Added: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/TPV103.py
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/TPV103.py (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/TPV103.py 2012-10-18 19:42:01 UTC (rev 20852)
@@ -0,0 +1,74 @@
+#!python
+#!python
+#!/usr/bin/env python
+
+import cubit
+import cubit2specfem3d
+
+import os
+import sys
+from save_fault_nodes_elements import *
+from absorbing_boundary import *
+
+cubit.cmd('playback "TPV103.jou" ')
+
+os.system('mkdir -p MESH')
+
+
+xmin = [9,16]
+xmax = [11,13]
+ymin = [3]
+ymax = [5]
+zmax = [8,15]
+zmin = [10,14]
+entities=['face']
+define_bc(entities,xmin,xmax,ymin,ymax,zmin,zmax)
+
+cubit.cmd('block 1 name "elastic 1" ')
+cubit.cmd('block 1 attribute count 5')
+cubit.cmd('block 1 attribute index 1 1')
+cubit.cmd('block 1 attribute index 2 6000')
+cubit.cmd('block 1 attribute index 3 3464')
+cubit.cmd('block 1 attribute index 4 2670')
+cubit.cmd('block 1 attribute index 5 13')
+
+cubit.cmd('block 2 name "elastic 2" ')
+cubit.cmd('block 2 attribute count 5')
+cubit.cmd('block 2 attribute index 1 1')
+cubit.cmd('block 2 attribute index 2 6000')
+cubit.cmd('block 2 attribute index 3 3464')
+cubit.cmd('block 2 attribute index 4 2670')
+cubit.cmd('block 2 attribute index 5 13')
+
+#### Export to SESAME format using cubit2specfem3d.py of GEOCUBIT
+
+cubit2specfem3d.export2SESAME('MESH')
+
+
+
+
+
+Au = [8] # A_up
+Ad = [3] # A_down
+
+faultA = fault_input(1,Au,Ad)
+quads_Aup,quads_Adp = save_cracks(faultA.name,faultA.surface_u,faultA.surface_d)
+#Unpacking list.
+quads_Au=unpack_list(quads_Aup)
+quads_Ad=unpack_list(quads_Adp)
+
+print 'len(Au):',len(quads_Au)
+print 'len(Ad):',len(quads_Ad)
+
+if not (len(quads_Au)==len(quads_Ad)):
+ print 'Number of elements for each fauld side up and down do not concide'
+
+save_elements_nodes(faultA.name,quads_Au,quads_Ad)
+
+
+
+
+
+
+
+
Added: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/description.pdf
===================================================================
(Binary files differ)
Property changes on: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/description.pdf
___________________________________________________________________
Name: svn:mime-type
+ application/octet-stream
More information about the CIG-COMMITS
mailing list