Hello, <div>I am trying to use specfem3D for a sort of unusual purpose compared to its traditional use. I am replicating small scale laboratory experiments about wave-propagation on metal-plate with scatterers and 3D features embedded or attached to it. I have been working for a while and the results seems overall good. One thing still remain not clear to me: how to set the right force magnitude to the source?</div>
<div><div> </div><div>The source I want to emulate acts in the vertical direction (perpendicular to the plate, Z-axis) and it applies a load of the order of 0.1N. How can i set this value? The code seems to ignore the moment tensor parameter when i am using a point force and rather using the &quot;<span class="Apple-style-span" style="font-family:Menlo;font-size:11px">FACTOR_FORCE_SOURCE&quot;</span>.  I have been following a recent discussion (september), which confirmed that by tweaking the code myself i went to the right direction, I disabled the original Ricker for the Heaviside. Those are the settings I have been using in the constant.h file so far:</div>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo"><br></p><p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo">logical, parameter :: USE_FORCE_POINT_SOURCE = .<span style="color:#bc319c">true</span>.</p>
<p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo"><span style="color:#bc319c">double</span> precision, parameter :: FACTOR_FORCE_SOURCE = <span style="color:#2834cf">1.</span>d4 </p>
<div><span class="Apple-style-span" style="font-family:Menlo;font-size:11px">integer, parameter :: COMPONENT_FORCE_SOURCE = <span style="color:#2834cf">3</span>  ! takes direction <span style="color:#bc319c">in</span> comp E/N/Z = <span style="color:#2834cf">1</span>/<span style="color:#2834cf">2</span>/<span style="color:#2834cf">3</span></span> </div>
<div><br></div><div>and regarding the receivers:</div><div><p style="margin:0.0px 0.0px 0.0px 0.0px;font:11.0px Menlo">logical, parameter :: EXT_MESH_RECV_NORMAL = .<span style="color:#bc319c">true</span>.</p></div><div><br>
</div><div>I have no UTM projection active and the source is a Dirac pulse (or better an Heaviside) as it is originally implemented in the code, i.e. reading the netlib error function. Frequency are high (10Khz).</div><div>
Thanks for the help and best regards</div><div> </div><div>Andrea Colombi</div><div><br></div><div><br></div><div><br></div><div>    </div></div>