[CIG-SEISMO] Potential bug in SPECFEM3D_Cartisen

[Junwei Huang]

Hi Dimitri,
I am afraid I don't know enough about the software to modify it without
introducing extra issues. In my case, I may work around it by hardwiring
the loop to
          imin = 1
          imax = NGLLX
          jmin = 1
          jmax = NGLLY
          kmin = 1
          kmax = NGLLZ
for both source and receivers. I wonder what was the reason to exclude the
edge nodes.

Junwei




On Mon, Jul 24, 2017 at 12:54 PM, Dimitri Komatitsch <
komatitsch at lma.cnrs-mrs.fr> wrote:

>
> Hi Junwei,
>
> Thank you very much. Very useful. Do you know how this should be fixed? If
> so, could you please send me the fix, or submit it directly as a pull
> request, following https://github.com/geodynamics
> /specfem3d/wiki/Using-Git-for-SPECFEM ?
>
> Thank you,
> Best regards,
>
> Dimitri.
>
> On 07/24/2017 06:11 PM, Junwei Huang wrote:
>
>> Hi developers,
>> Would like to report that this line of code may cause issues:
>>
>> Inline image 1
>>
>>
>> As iglob_is_surface_external_mesh is true for the external meshes only
>> and the loop is iterating within internal nodes,
>> iglob_is_surface_external_mesh is always false.
>>
>> It can cause issues when "SOURCES_CAN_BE_BURIED" or
>> "RECEIVERS_CAN_BE_BURIED" are false. Hope it makes sense.
>>
>> Best regards,
>> Junwei
>>
>>
>> _______________________________________________
>> CIG-SEISMO mailing list
>> CIG-SEISMO at geodynamics.org
>> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
>>
>>
> --
> Dimitri Komatitsch, CNRS Research Director (DR CNRS)
> Laboratory of Mechanics and Acoustics, Marseille, France
> http://komatitsch.free.fr
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20170724/0f9cc7fc/attachment.html>