( '=R', 'x2', [(), (), 1, ()], 'y_coarse'),
]
fieldio += [
( '=w', 'x1', [(), (), 1, ()], 'xx' ),
( '=w', 'x2', [(), (), 1, ()], 'yy' ),
]
# nucleation
delts = 1.2
rnucl = 2000
# output
_dnt = int( 0.04 / dt )
ihypo = 20000./dx[0], (nn[1] + 1) / 2, 1.5
k = ihypo[1]
fieldio += [
( '=w', 'sl', [(),k,1,(1,-1,_dnt)], 'sl' ),
( '=w', 'tsm', [(),k,1,(1,-1,_dnt)], 'tsm' ),
( '=w', 'svm', [(),k,1,(1,-1,_dnt)], 'svm' ),
( '=w', 'tnm', [(),k,1,(1,-1,_dnt)], 'tnm' ),
( '=w', 'psi', [(),k,1,(1,-1,_dnt)], 'psi' ),
( '=w', 'a1', [(),(),1,(1,-1,_dnt)], 'ax' ),
( '=w', 'a2', [(),(),1,(1,-1,_dnt)], 'ay' ),
( '=w', 'epm', [(),(),1,(1,-1,_dnt)], 'epm' ),
]
sord.run( locals() )
#!/usr/bin/env python
import sord # import the sord module
rundir = 'tmp' # directory location for output
dx = 100.0, 100.0, 100.0 # spatial step length in x, y, and z
dt = 0.0075 # time step length
nn = 61, 61, 61 # number of mesh nodes in x, y, and z
nt = 60 # number of time steps
fieldio = [ # field variable input and output
('=', 'rho', [], 2670.0), # material density
('=', 'vp', [], 6000.0), # material P-wave velocity
('=', 'vs', [], 3464.0), # material S-wave velocity
('=', 'gam', [], 0.3), # material viscosity
('=w', 'v1', [0, 0, 31, -1], 'vx'), # write X velocity slice output
('=w', 'v2', [0, 0, 31, -1], 'vy'), # write Y velocity slice output
]
ihypo = 31.0, 31.0, 31.0 # source location
source = 'potency' # source type
source1 = 1e6, 1e6, 1e6 # source normal components
source2 = 0.0, 0.0, 0.0 # source shear components
timefunction = 'brune' # source time function
period = 6 * dt # source dominant period
sord.run(locals()) # launch SORD job
