Issue using external tomo file for plane wave simulation

[kenenig]

I was having issue running the tomographic_ocean_model example. It give back an immediate backtrace error when I run it either serially or with multiple processors. I am using an Macbook with a M1 processor so my guess is that there is something wrong when first compiling Specfem but I am not sure. When I tried creating my own tomographic model and running a similar plane wave simulation, the tomo file was read successfully but returned the same backtrace error when running the internal mesher. Has anyone else had a similar issue using an external tomography file?

[homnath]

Which version are you using? It works fine with the recent devel branch. But I tried on Linux and WINDOWS.

[kenenig]

I am running 8.1.0 version and I am using the master branch. I just looked into the devel branch too but it said that it had issues on Mac.

As an update, I am able to run this example. However, for the simulation which has my external tomography file but a different issue occurs. The tomography file appears to be successfully read (see screenshot below) but the last velocity is the only one used and the material is made homogenous. I was trying to test if I can use a tomography file for a non-plane wave simulation and it was giving me the same issue. I included my Par_file, SOURCE, and external tomography file used for the run in the attached zip file for context.

SF_tomo_issue.zip

[homnath]

How did you find, it is using the homogeneous model? Can you also upload all mesh files?

…

--------------------------------------- Hom Nath Gharti, PhD Assistant Professor | Digital Earth Scientist https://www.digitalearthscience.com/ Department of Geological Sciences and Geological Engineering Miller Hall 314, 36 Union St Queen’s University Kingston, ON K7L 3N6, Canada Queen’s University is situated on traditional Anishinaabe and Haudenosaunee Territory.<https://www.queensu.ca/encyclopedia/t/traditional-territories>

________________________________ From: kenenig ***@***.***> Sent: Wednesday, May 8, 2024 11:55 PM To: SPECFEM/specfem3d ***@***.***> Cc: Hom Nath Gharti ***@***.***>; Comment ***@***.***> Subject: Re: [SPECFEM/specfem3d] Issue using external tomo file for plane wave simulation (Discussion #1694) I am running 8.1.0 version and I am using the master branch. I just looked into the devel branch too but it said that it had issues on Mac. As an update, I am able to run this example. However, for the simulation which has my external tomography file but a different issue occurs. The tomography file appears to be successfully read (see screenshot below) but the last velocity is the only one used and the material is made homogenous. I was trying to test if I can use a tomography file for a non-plane wave simulation and it was giving me the same issue. I included my Par_file, SOURCE, and external tomography file used for the run in the attached zip file for context. [SF_tomo_issue.zip](https://github.c<https://github.c/> SF_tomo_terminal_output.png (view on web)<https://github.com/SPECFEM/specfem3d/assets/140196629/9c22140c-3ff5-4b90-8fb4-f7ac8e8bd121> om/SPECFEM/specfem3d/files/15257309/SF_tomo_issue.zip) — Reply to this email directly, view it on GitHub<#1694 (reply in thread)>, or unsubscribe<https://github.com/notifications/unsubscribe-auth/ABMCQ4TWLM3Z3QUM2CRFXTDZBLXSBAVCNFSM6AAAAABGZQAAS2VHI2DSMVQWIX3LMV43SRDJONRXK43TNFXW4Q3PNVWWK3TUHM4TGNRTGAYTO>. You are receiving this because you commented.Message ID: ***@***.***>

[kenenig]

Yes here they are. I used the plot_model_from_gll_or_binary_file.py script to plot the vs, vp, and rho and they were all one number.

mesh_0000.zip
mesh_0001.zip
mesh_0002.zip
mesh_0003_rho_vp.zip
mesh_0003vs_0004rho.zip
mesh_0004_vs_vp.zip

[homnath]

I see. Can you send MOR1_interfaces_outsideonly.dat file?

[kenenig]

Yes, here it is. I only defined the outer portions of the model space while troubleshooting. MOR1_interfaces_outsideonly.dat.zip

[homnath]

I think it is correctly implemeting the tomographic model. Does your model resemble the one shown in the snapshot?

[homnath]