[CIG-SEISMO] SPECFEM3D

[Gula Aleksey]

Hello,



I have a several questions about SPECFEM3D:

I try to use specfem in local problems, more precisely, I compute models in
scale of seismic.

So, I have next questions:



1) How I can choose source time function, I want to use  “ricker”, for
example?

2) What about coordinates of receivers and source, if I use CUBIT for
meshing with next parameters:



cubit.cmd('brick x 2000 y 2000 z 2215')

cubit.cmd('volume 1 move x 10000 y 10000 z -1107,5')

cubit.cmd('brick x 2000 y 2000 z 4600')

cubit.cmd('volume 2 move x 1000 y 1000 z -2445')

cubit.cmd('brick x 2000 y 2000 z 5000')

cubit.cmd('volume 3 move x 1000 y 1000 z -2925')



SPECFEM3D use UTM zones, but in my case it doesn’t matter in which it
should to be located.



I made next CMTSOLUTION file



PDE  0000 01 01 00 00 00.00  1000 1000 0.200 0.0 0.0 hom_explosion
event name:       hom_explosion
time shift:       0.00
half duration:    0.5
latitude:      1000.0
longitude:     1000.0
depth:          0.2
Mrr:       1.000000e+23
Mtt:       1.000000e+23
Mpp:       1.000000e+23
Mrt:       0.000000
Mrp:       0.000000
Mtp:       0.000000



In result of computation was of the fact that coordinates of source change
to (300;300;200)



So I made next STATIONS file



R1 KG 500.00000 1000.00000 0.0 200.0
R2 KG 550.00000 1000.00000 0.0 200.0
R3 KG 600.00000 1000.00000 0.0 200.0
R4 KG 650.00000 1000.00000 0.0 200.0
R5 KG 700.00000 1000.00000 0.0 200.0
R6 KG 750.00000 1000.00000 0.0 200.0
R7 KG 800.00000 1000.00000 0.0 200.0
R8 KG 850.00000 1000.00000 0.0 200.0
R9 KG 900.00000 1000.00000 0.0 200.0
R10 KG 950.00000 1000.00000 0.0 200.0
R11 KG 1000.00000 1000.00000 0.0 200.0
R12 KG 1050.00000 1000.00000 0.0 200.0
R13 KG 1100.00000 1000.00000 0.0 200.0
R14 KG 1150.00000 1000.00000 0.0 200.0
R15 KG 1200.00000 1000.00000 0.0 200.0
R16 KG 1250.00000 1000.00000 0.0 200.0
R17 KG 1300.00000 1000.00000 0.0 200.0
R18 KG 1350.00000 1000.00000 0.0 200.0
R19 KG 1400.00000 1000.00000 0.0 200.0
R20 KG 1450.00000 1000.00000 0.0 200.0
R21 KG 1500.00000 1000.00000 0.0 200.0



But after computation I had next error, that no receivers didn’t  get into
the model



Best regards,

Aleksey
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/cig-seismo/attachments/20120224/93f3a5f9/attachment.htm