[CIG-SEISMO] [SW4] difficulties to implement material model

Sun Feb 18 11:21:10 PST 2018

Dear all,

I’m trying implement material model using command block, but SW4 only take in consideration the first line of the block. In consequence the error is the same for all points of computation grid. In the following example I used the block model of example file named artie-coarse.in and I got the same error with my material model:

Point (i,j,k)=(9, 152, -1) in grid g=0
with (x,y,z)=(5.333333e+02,1.006667e+04,-1.333333e+02) and depth=-1.333333e+02
is outside the block domain: -1.000000e+05<= x <= 2.000000e+05, -1.000000e+05 <= y <= 2.000000e+05, 2.100000e+04 <= depth <= 8.000000e+04

The block bounds are:
********Parsing source command*********
Cartesian coordinates of source at (lon, lat)=(-1.552100e+02, 1.933000e+01) is (x,y)=(31959.8, 72480.5)
Moment source at x=3.195983e+04, y=7.248048e+04, z=8.000000e+03 is centered at grid point i=480, j=1088, k=121, in grid=0
********Done parsing source command*********
*** Ignoring command: 'absdepth=0'
Block 1 has bounds -1.00000000e+05 2.00000000e+05 -1.00000000e+05 2.00000000e+05 2.10000000e+04 8.00000000e+04
Block 2 has bounds -1.00000000e+05 2.00000000e+05 -1.00000000e+05 2.00000000e+05 0.00000000e+00 2.00000000e+00
Block 3 has bounds -1.00000000e+05 2.00000000e+05 -1.00000000e+05 2.00000000e+05 2.00000000e+00 6.00000000e+00
Block 4 has bounds -1.00000000e+05 2.00000000e+05 -1.00000000e+05 2.00000000e+05 6.00000000e+00 1.20000000e+01
Block 5 has bounds -1.00000000e+05 2.00000000e+05 -1.00000000e+05 2.00000000e+05 1.20000000e+01 2.00000000e+01
Block 6 has bounds -1.00000000e+05 2.00000000e+05 -1.00000000e+05 2.00000000e+05 2.00000000e+01 3.00000000e+01
Block 7 has bounds -1.00000000e+05 2.00000000e+05 -1.00000000e+05 2.00000000e+05 3.00000000e+01 1.00000000e+02
Block 8 has bounds -1.00000000e+05 2.00000000e+05 -1.00000000e+05 2.00000000e+05 1.00000000e+02 3.00000000e+02
Block 9 has bounds -1.00000000e+05 2.00000000e+05 -1.00000000e+05 2.00000000e+05 3.00000000e+02 5.00000000e+02
Block 10 has bounds -1.00000000e+05 2.00000000e+05 -1.00000000e+05 2.00000000e+05 5.00000000e+02 7.00000000e+02
Block 11 has bounds -1.00000000e+05 2.00000000e+05 -1.00000000e+05 2.00000000e+05 7.00000000e+02 1.00000000e+03
Block 12 has bounds -1.00000000e+05 2.00000000e+05 -1.00000000e+05 2.00000000e+05 1.00000000e+03 3.00000000e+03
Block 13 has bounds -1.00000000e+05 2.00000000e+05 -1.00000000e+05 2.00000000e+05 3.00000000e+03 5.00000000e+03
Block 14 has bounds -1.00000000e+05 2.00000000e+05 -1.00000000e+05 2.00000000e+05 5.00000000e+03 6.00000000e+03
Block 15 has bounds -1.00000000e+05 2.00000000e+05 -1.00000000e+05 2.00000000e+05 6.00000000e+03 1.10000000e+04
Block 16 has bounds -1.00000000e+05 2.00000000e+05 -1.00000000e+05 2.00000000e+05 1.10000000e+04 1.60000000e+04
Block 17 has bounds -1.00000000e+05 2.00000000e+05 -1.00000000e+05 2.00000000e+05 1.60000000e+04 2.10000000e+04

********Done reading the input file*********
Receiver INFO for station hawaii:
     initial location (x,y,z) = 7.98272135e+04 3.87629303e+04 0.00000000e+00 zTopo= 0.00000000e+00
     nearest grid point (x,y,z) = 7.98000000e+04 3.87333333e+04 0.00000000e+00 h= 6.66666667e+01 with indices (i,j,k)= 1198 582 1 in grid 0
   Execution time, reading input file 2 minutes 1.58231549e+01 seconds 
  Using Bjorn's fast (parallel) IO library
Assuming a PARALLEL file system
Writing images from (up to) 8 procs
Setting up SBP boundary stencils
Info: Grid #0 min z-coordinate: 0.00000000e+00
Info: mGridSize[0]=6.66666667e+01

The input parameters are the following ones:

# SW4 
# input RUN:teste
#
fileio path=teste  pfs=1 verbose=5
#
# Grid coords
grid x=100.0e3 y=100.0e3 z=40.0e3 nx=1501 lat=19.042900 lon=-155.900000 az=0
#
supergrid gp=60
# source event
source lat=19.33 lon=-155.21 depth=8.0e3 strike=175 dip=60 rake=-30 m0=3.13e16 mxx=4.17e12 myy=9.36e12 mzz=-1.35e13 mxy=-2.43e13 t0=1.0 freq=6.28 type=Gaussian
#
#time steps=0
time t=60 
#
#globalmaterial  vpmin=866 vsmin=500
#
supergrid gp=30
#
# Earth Model
# block model 1D
absdepth=0
block z1=21000.0 vs=3800.0 vp=6800.0 rho=2950.0 qs=190.0 qp=380.0
block z1=0.0 z2=2.0 vs=450.0 vp=1700.0 rho=2000.0 qs=22.5 qp=45.0
block z1=2.0 z2=6.0 vs=650.0 vp=1800.0 rho=2100.0 qs=32.5 qp=65.0
block z1=6.0 z2=12.0 vs=850.0 vp=1800.0 rho=2100.0 qs=42.5 qp=85.0
block z1=12.0 z2=20.0 vs=950.0 vp=1900.0 rho=2100.0 qs=47.5 qp=95.0
block z1=20.0 z2=30.0 vs=1150.0 vp=2000.0 rho=2200.0 qs=57.5 qp=115.0
block z1=30.0 z2=100.0 vs=1200.0 vp=2400.0 rho=2200.0 qs=60.0 qp=120.0
block z1=100.0 z2=300.0 vs=1400.0 vp=2800.0 rho=2300.0 qs=70.0 qp=140.0
block z1=300.0 z2=500.0 vs=1600.0 vp=3100.0 rho=2400.0 qs=80.0 qp=160.0
block z1=500.0 z2=700.0 vs=1800.0 vp=3400.0 rho=2450.0 qs=90.0 qp=180.0
block z1=700.0 z2=1000.0 vs=2100.0 vp=3700.0 rho=2500.0 qs=105.0 qp=210.0
block z1=1000.0 z2=3000.0 vs=2400.0 vp=4400.0 rho=2600.0 qs=120.0 qp=240.0
block z1=3000.0 z2=5000.0 vs=2800.0 vp=5100.0 rho=2700.0 qs=140.0 qp=280.0
block z1=5000.0 z2=6000.0 vs=3150.0 vp=5600.0 rho=2750.0 qs=157.5 qp=315.0
block z1=6000.0 z2=11000.0 vs=3600.0 vp=6150.0 rho=2825.0 qs=180.0 qp=360.0
block z1=11000.0 z2=16000.0 vs=3650.0 vp=6320.0 rho=2850.0 qs=182.5 qp=365.0
block z1=16000.0 z2=21000.0 vs=3700.0 vp=6550.0 rho=2900.0 qs=185.0 qp=370.0
#
# restricted output
rec lat=19.76 lon=-155.53 depth=0 sta=POHA file=hawaii usgsformat=0 sacformat=1 variables=velocity
#
# Image output z=0
#image z=0 mode=hmax file=image time=59.4
#image z=0 mode=hmax file=image timeInterval=0.2
#image z=0 mode=vmax file=image time=59.4

Sincerely yours,
João Fontiela

*********************
*********************
João Fontiela
PhD Student
skype: fontiela
email: fontiela at gmail.com
mobile phone +351 917 862 830

Profissional contacts:
email: jfontiela at uevora.pt
Tel: (+351) 266 745 354 ext: 5395 or 5413
Fax: (+351) 266 745 394 

Profissional address:
Universidade de Évora 
Instituto de Ciências da Terra
Colégio Luís Verney 
Rua Romão Ramalho, 59
7000-671 Évora - Portugal

ResearchGate: https://www.researchgate.net/profile/Joao_Fontiela

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20180218/06d40c2d/attachment-0001.html>