<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi Aleksey,<div><br></div><div>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 <span class="Apple-style-span" style="font-family: Monaco; font-size: 12px; ">comp_source_time_function***(..) </span>routines to <span class="Apple-style-span" style="font-family: Monaco; font-size: 12px; ">comp_source_time_function_rickr(..). </span>Or you could just redefine the functions in comp_source_time_function.f90.</div><div><br></div><div>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.</div><div><br></div><div>2): set flag <span class="Apple-style-span" style="font-family: Monaco; font-size: 12px; ">SUPPRESS_UTM_PROJECTION = .true. in Par_file,</span>then it won't use, that is suppresses UTM conversions. </div><div><br></div><div>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:</div><div><a href="http://www.geodynamics.org/cig/software/specfem3d/tutorials">http://www.geodynamics.org/cig/software/specfem3d/tutorials</a></div><div><br></div><div>hopefully this helps out a bit,</div><div>best wishes,</div><div>daniel</div><div><br></div><div><br></div><div><br><div><div>On Feb 24, 2012, at 5:49 AM, Gula Aleksey wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><meta http-equiv="Content-Type" content="text/html; charset=Windows-1252"><p class="MsoNormal"><span lang="EN-US">Hello,</span><span lang="EN-US"></span></p><div><span lang="EN-US"> </span><br class="webkit-block-placeholder"></div><p class="MsoNormal"><span lang="EN-US">I have a several questions about SPECFEM3D:</span></p><p class="MsoNormal"><span lang="EN-US">I try to use specfem in local problems, more precisely, I compute models in scale of seismic.</span></p><p class="MsoNormal"><span lang="EN-US">So, I have next questions:</span></p><div><span lang="EN-US"> </span><br class="webkit-block-placeholder"></div><p class="MsoNormal"><span lang="EN-US">1) How I can choose source time function, I want to use “ricker”, for example?</span></p><p class="MsoNormal"><span lang="EN-US">2) What about coordinates of receivers and source, if I use CUBIT for meshing with next parameters:</span></p><div><span lang="EN-US"> </span><br class="webkit-block-placeholder"></div><p class="MsoNormal"><span lang="EN-US">cubit.cmd('brick x 2000 y 2000 z 2215')</span></p><p class="MsoNormal"><span lang="EN-US">cubit.cmd('volume 1 move x 10000 y 10000 z -1107,5')</span></p><p class="MsoNormal"><span lang="EN-US">cubit.cmd('brick x 2000 y 2000 z 4600')</span></p><p class="MsoNormal"><span lang="EN-US">cubit.cmd('volume 2 move x 1000 y 1000 z -2445')</span></p><p class="MsoNormal"><span lang="EN-US">cubit.cmd('brick x 2000 y 2000 z 5000')</span></p><p class="MsoNormal"><span lang="EN-US">cubit.cmd('volume 3 move x 1000 y 1000 z -2925')</span></p><div><span lang="EN-US"> </span><br class="webkit-block-placeholder"></div><p class="MsoNormal"><span lang="EN-US">SPECFEM3D use UTM zones, but in my case it doesn’t matter in which it should to be located.</span></p><div><span lang="EN-US"> </span><br class="webkit-block-placeholder"></div><p class="MsoNormal"><span lang="EN-US">I made next CMTSOLUTION file </span></p><div><span lang="EN-US"> </span><br class="webkit-block-placeholder"></div><p class="MsoPlainText"><span lang="EN-US" style="font-family:"Courier New"">PDE 0000 01 01 00 00 00.00 1000 1000 0.200 0.0 0.0 hom_explosion<br> event name: hom_explosion<br> time shift: 0.00<br> half duration: 0.5<br> latitude: 1000.0<br> longitude: 1000.0<br> depth: 0.2<br> Mrr: 1.000000e+23<br> Mtt: 1.000000e+23<br> Mpp: 1.000000e+23<br> Mrt: 0.000000<br> Mrp: 0.000000<br> Mtp: 0.000000</span><span lang="EN-US"></span></p><div><span lang="EN-US"> </span><br class="webkit-block-placeholder"></div><p class="MsoNormal"><span lang="EN-US">In result of computation was of the fact that coordinates of source change to (300;300;200) </span></p><div><span lang="EN-US"> </span><br class="webkit-block-placeholder"></div><p class="MsoNormal"><span lang="EN-US">So I made next STATIONS file</span></p><div><span lang="EN-US"> </span><br class="webkit-block-placeholder"></div><p class="MsoPlainText"><span lang="EN-US" style="font-family:"Courier New"">R1 KG 500.00000 1000.00000 0.0 200.0<br> R2 KG 550.00000 1000.00000 0.0 200.0<br> R3 KG 600.00000 1000.00000 0.0 200.0<br> R4 KG 650.00000 1000.00000 0.0 200.0<br> R5 KG 700.00000 1000.00000 0.0 200.0<br> R6 KG 750.00000 1000.00000 0.0 200.0<br> R7 KG 800.00000 1000.00000 0.0 200.0<br> R8 KG 850.00000 1000.00000 0.0 200.0<br> R9 KG 900.00000 1000.00000 0.0 200.0<br> R10 KG 950.00000 1000.00000 0.0 200.0<br> R11 KG 1000.00000 1000.00000 0.0 200.0<br> R12 KG 1050.00000 1000.00000 0.0 200.0<br> R13 KG 1100.00000 1000.00000 0.0 200.0<br> R14 KG 1150.00000 1000.00000 0.0 200.0<br> R15 KG 1200.00000 1000.00000 0.0 200.0<br> R16 KG 1250.00000 1000.00000 0.0 200.0<br> R17 KG 1300.00000 1000.00000 0.0 200.0<br> R18 KG 1350.00000 1000.00000 0.0 200.0<br> R19 KG 1400.00000 1000.00000 0.0 200.0<br> R20 KG 1450.00000 1000.00000 0.0 200.0<br> R21 KG 1500.00000 1000.00000 0.0 200.0</span><span lang="EN-US"></span></p><div><span lang="EN-US"> </span><br class="webkit-block-placeholder"></div><p class="MsoNormal"><span lang="EN-US">But after computation I had next error, that no receivers didn’t get into the model</span></p><div><span lang="EN-US"> </span><br class="webkit-block-placeholder"></div><p class="MsoNormal"><span lang="EN-US">Best regards,</span></p><p class="MsoNormal"><span lang="EN-US">Aleksey</span></p> _______________________________________________<br>CIG-SEISMO mailing list<br><a href="mailto:CIG-SEISMO@geodynamics.org">CIG-SEISMO@geodynamics.org</a><br>http://geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo<br></blockquote></div><br></div></body></html>