Error writing output when scalar_moment is none

[mrosskopf]

Hello,

I am currently testing SourceSpec for several events in the BedrettoLab together with Fred and we encountered the Error shown below.
Is there anything we need to change or is an exemption needed in the main code?

Thanks for your help!
Best
Martina

---------------------------------------------------
Output written to file: sspec_out/2026aawgrlxfz/2026aawgrlxfz.ssp.out
Output written to file: sspec_out/2026aawgrlxfz/2026aawgrlxfz.ssp.yaml
Traceback (most recent call last):
  File "/srv/sspec/shared_env/sharedSspec/bin/source_spec", line 6, in <module>
    sys.exit(main())
             ^^^^^^
  File "/srv/sspec/shared_env/sharedSspec/lib/python3.11/site-packages/sourcespec/source_spec.py", line 78, in main
    write_output(config, sspec_output)
  File "/srv/sspec/shared_env/sharedSspec/lib/python3.11/site-packages/sourcespec/ssp_output.py", line 440, in write_output
    write_qml(config, sspec_output)
  File "/srv/sspec/shared_env/sharedSspec/lib/python3.11/site-packages/sourcespec/ssp_qml_output.py", line 136, in write_qml
    mt.scalar_moment = Mo_summary.value
    ^^^^^^^^^^^^^^^^
  File "/srv/sspec/shared_env/sharedSspec/lib/python3.11/site-packages/obspy/core/event/base.py", line 384, in __setattr__
    raise ValueError(msg)
ValueError: On MomentTensor object: Value 'nan' for 'scalar_moment' is not a finite floating point value.
ls -ltrh sspec_out/2026aawgrlxfz/2026aawgrlxfz.ssp.yaml

[claudiodsf]

Hi Martina,
thanks for the feedback.

This is a strange error, since the scalar moment should not be nan at this stage: the code should have thrown some error earlier...

Could you please share the log file for this event (probably called 2026aawgrlxfz.ssp.log)?

[mrosskopf]

Hi Claudio,

sure I attached the log file.

Best
Martina

2026aawgrlxfz.ssp.log

[claudiodsf]

Thanks,

there is something strange going on in the inversion: the retrieved Mw does not respect the given bounds.

Exemple:

2026-04-10 07:07:53,405 ssp_inversion        INFO     8R.V0401..FPH broadb: initial values: Mw_0: -2.083; fc_0: 249.98; t_star_0: 0.045
2026-04-10 07:07:53,407 ssp_inversion        INFO     8R.V0401..FPH broadb: bounds: Mw: -2.5715, -1.6749; fc: 100.0, 900.0; t_star: 0.001, 0.25
2026-04-10 07:08:06,836 ssp_inversion        INFO     8R.V0401..FPH broadb: optimal values: Mw: 2.2283; fc: 0.0678; t_star: 0.0004
2026-04-10 07:08:06,836 ssp_inversion        INFO     8R.V0401..FPH broadb: Normalized RMS misfit: 0.306
2026-04-10 07:08:06,836 ssp_inversion        INFO     8R.V0401..FPH broadb: Quality of fit: 73.615%
2026-04-10 07:08:06,837 ssp_inversion        INFO     8R.V0401..FPH broadb: spectral weight close to fc: 1.000

Bounds: Mw: -2.5715, -1.6749
Estimated value: Mw: 2.2283

Could you please zip the whole folder sspec_out/2026aawgrlxfz and share it here?

[mrosskopf]

Sure, I attached it. :)

2026aawgrlxfz.zip