[CIG-SEISMO] need help with define_external_model.f90 in specfem2D

[Dimitri Komatitsch]

Hi Ignazio, Hi all,

I think Alexis will know the answer because he actively uses SPECFEM2D, 
including for external models :-)

Best wishes,
Dimitri.

On 03/09/2017 05:21 PM, ignazio damiani wrote:
> Dear all,
> I'm using specfem2D and I'm trying to assign an external velocity model. The aim is to create a model with lateral variation of seismic wave velocities and density.
>  I'm trying to use define_external_model.f90 file but even though I suppress the routine for the AK135F, the function doesn't work.
> I always get the same error:
>
>> External model: ak135-f
>>  STOP wrong flag number
>
> This is what i'm doing:
> firstly, i set in Par_file MODEL=external.
> Secondly, in read_external_model.f90 i set: case ('external') call define_external_model_dummy().
>> And finally in define_external_model.f90 i leave the subroutine as define_external_model_dummy(), i suppress the AK135F routine and i modify the number of material_element(ispec) based on my coordinates and region.
>> What's wrong? because it remains the error.
>
> Thanks for your help.
> Thanks Dimitri for the help and for giving me your email.
>
> Best regards
>
> Ignazio Damiani
> --
> Msc in Exploration and Applied Geophysics. Pisa. Italy
>

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