[CIG-SEISMO] SPECFEM1D: Bug in Local to Global numbering of the source

[Dimitri Komatitsch]

Dear Martin,

Thanks for your feedback. I checked the source code, I think it is fine, 
both i_source and ispec_source are defined above, in the beginning of 
the code:

! set the source
   ispec_source = (NSPEC+1)/2
   i_source = (NGLL+1)/2


Best regards,
Dimitri.

On 04/24/2012 12:15 PM, Martin van Driel wrote:
> Dear all,
>
> I think the source implementation in SPECFEM1D is buggy (it is commented
> out anyway, but it might save others some time to have it correct
> nonetheless):
>
> the variables iglob_source and stf are not defined and the indices of
> the local to global numbering are swapped. This way it seems to work:
>
> ! add source at global level
>      accel(ibool(i_source,ispec_source)) =&
> 	accel(ibool(i_source,ispec_source)) +&
>           source_time_function(dble(it - 1) * DT - hdur,hdur) * source_amp
>
>
> All the best,
> Marti
> _______________________________________________
> CIG-SEISMO mailing list
> CIG-SEISMO at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo

-- 
Dimitri Komatitsch - komatitsch aT lma.cnrs-mrs.fr
CNRS Research Director (DR CNRS), Laboratory of Mechanics and Acoustics,
UPR 7051, Marseille, France    http://komatitsch.free.fr