[CIG-SEISMO] assign external model specfem2D issue

[Eve Tsang-Hin-Sun]

Dear all,

I would like to run some simulations with solid and fluid layers. I saw that there is an example provided in the EXAMPLES folder and was trying to work on it to understand the details. I actually encountered some problems when trying to assign an external velocity model. My first goal is very simple as I just try to adapt the "define_external_model_dummy » routine to the «  fluid_solid »  case.

As I understand it the variable « material_element » is a vector containing the model number for each spectral element (so 1 or 2 in my case). In the routine, each element is tested to read its model number and then something happens (eg assigning some velocity values).

In the Par_file_fluid_solid file, the fluid layer is numbered 2 and the solid layer is numbered 1.


Any way, each time I try to run the code, I always get the same error :
STOP error: non acoustic material in DATA/Par_file or external mesh redefined as acoustic in define_external_model()

meaning that something is wrong in my Par_file or in my routine. Could anyone explain to me why such an error would occur or what did I miss? I also join the two files.

My second problem concerns the «  read_external_file routine »  : in the function, the « define_external_routine » is called. However the inputs of the routine do not coincide with the ones given in the dummy example (for example, gravity or Nsq). But the code compiles anyway…How is it possible ? I admit that I am not fully familiar with f90…

Thanks for your help.

Regards,
Eve
--
Eve Tsang-Hin-Sun
Doctorante
Institut Universitaire Européen de la Mer (IUEM)
Laboratoires Domaines Océaniques (LDO)
Place Nicolas Copernic - 29280 Plouzané

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20150109/edc39e6d/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: define_external_model.f90
Type: application/octet-stream
Size: 2818 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20150109/edc39e6d/attachment.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20150109/edc39e6d/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Par_file_fluid_solid
Type: application/octet-stream
Size: 13839 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20150109/edc39e6d/attachment-0001.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20150109/edc39e6d/attachment-0002.html>