[CIG-SEISMO] SPECFEM3D

[daniel peter]

Hi Aleksey,

1): if you like fortran coding, then depending if you use an acoustic  
or elastic domain & sources, look into the  
compute_add_sources_acoustic.f90 / compute_add_sources_elastic.f90  
files and change the comp_source_time_function***(..) routines to  
comp_source_time_function_rickr(..). Or you could just redefine the  
functions in comp_source_time_function.f90.

if you don't like fortran, set the half-duration to 0.0 in the  
CMTSOLUTION file and convolve the traces with your source time  
function afterwards.

2): set flag SUPPRESS_UTM_PROJECTION = .true. in Par_file,then it  
won't use, that is suppresses UTM conversions.

you find examples in the examples/ folder to get the idea, the  
documentation is described in manual_SPECFEM3D.pdf and there are  
tutorials on the website:
http://www.geodynamics.org/cig/software/specfem3d/tutorials

hopefully this helps out a bit,
best wishes,
daniel



On Feb 24, 2012, at 5:49 AM, Gula Aleksey wrote:

> 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
>
> _______________________________________________
> CIG-SEISMO mailing list
> CIG-SEISMO at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo

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