adjtomo / seisflows

An automated workflow tool for full waveform inversion and adjoint tomography
http://seisflows.readthedocs.org
BSD 2-Clause "Simplified" License
172 stars 122 forks source link

Error at 'postprocess_event_kernels' stage #185

Closed raulleoncz closed 3 months ago

raulleoncz commented 5 months ago

Hello @bch0w, this is a new issue related to the problem #182. I ran the simulation step by step and I found out that I got the problem at the step 3 of the workflow. The problem is the same as in #182. I checked the file 'adj_solver.log' and I getting this: Captura de pantalla 2024-01-24 203933

In case of 'OUTPUT_FILES/output_solver.txt', there was not a file with that name. I also see if I have the same number of traces and there was the 10 receivers that I was expecting: Captura de pantalla 2024-01-24 204133

I add you the error I got from 'combine_vp.log' just in case: Captura de pantalla 2024-01-24 204025

Hope you can help me. Thank you.

bch0w commented 5 months ago

Hi @raulleoncz, sorry for a slow reply here. I think I understand where this is coming from now, the kernels output by SPECFEM are in the format proc??????_<parameter>_acoustic_kernel.bin, whereas SeisFlows is expecting files named proc??????_<parameter>_kernel.bin. So the file manipulation functions are failing to catch the right kernel files and the workflow is failing.

That additional acoustic tag was something I was not aware of as I have not explored the 2D acoustic workflow much. Thanks for bringing this to my attention, it will likely require me to run a 2D acoustic example on my own to ensure that the file manipulation parts of the workflow are able to successfully recognize and catch these new kernel files.

As for how you can get past this in the interim, this is a bit hacky and not guaranteed to work, but you might be able to move these kernel files to the correct location on your own. Here are some pseudo commands that you can run:

cd ~/seisflows/Tiber/scratch/solver/001/OUTPUT_FILES
mv *_kernel.bin ~/seisflows/Tiber/scratch/eval_grad/kernels/001
cd  ~/seisflows/Tiber/scratch/eval_grad/kernels/001
rename _acoustic '' *
cd ~/seisflows
seisflows resume

I am not 100% sure this will work as the parameter names are also different. Let me test out an acoustic example on my own and get back to you! I'll also try to work on the SU preprocessing at the same time.

bch0w commented 3 months ago

Hi @raulleoncz, I will close this issue and continue work on this in #202, please continue discussion in that issue. Thanks!