[CIG-SEISMO] PS

[Dimitri Komatitsch]

PS: but you are right, in the line that is commented out
"ibool(ispec_source,i_source)" should be "ibool(i_source,ispec_source)", 
and iglob_source should be declared.

Let me fix the problem and release a new version.

Thanks,
Dimitri.

On 04/30/2012 05:22 PM, Dimitri Komatitsch wrote:
>
> 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