Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Pass eventID to NLLoc Hypocenter-Phase file #59

Open
loicbachelot opened this issue Sep 26, 2024 · 2 comments
Open

Pass eventID to NLLoc Hypocenter-Phase file #59

loicbachelot opened this issue Sep 26, 2024 · 2 comments

Comments

@loicbachelot
Copy link

Hello,

I have multiple events in the same Phase file. I separate the events with a blank line and I added a comment with the event id before the list of arrivals.
What is the best approach to have this event id in the hypocenter-Phase file, either as the filename or as a field in the file itself?
Having this eventID will be very valuable to keep track of each event separately.
For this use case, it it better to keep each event in a different Phase file?

Thanks

@alomax
Copy link
Collaborator

alomax commented Sep 27, 2024

Hello @loicbachelot

If you are using NLL-phase format NLLOC_OBS or HYPO71or HYPOELLIPSE for input, you can add one or two NLL-hypocenter-phase lines to add event id, magnitude and other information which will be propagated into the output. The lines for each event are of the form:

PUBLIC_ID smi:org.gfz-potsdam.de/geofon/mynet2024suwa
QUALITY  Pmax -1 MFmin -1 MFmax -1 RMS -1 Nphs -1 Gap -1 Dist -1 Mamp 1.75 39 Mdur -9.90 0
PHASE ID Ins Cmp On Pha  FM Date     HrMn   Sec     Err  ErrMag    Coda      Amp       Per       PriorWt
NN_SS15 ?    ?    ? P      ? 20240924 1327 12.4000 GAU  3.00e-01 -1.00e+00 -1.00e+00 -1.00e+00  1.00e+00
NN_SS15 ?    ?    ? S      ? 20240924 1327 14.9323 GAU  3.00e-01 -1.00e+00 -1.00e+00 -1.00e+00  1.00e+00
NN_SS09 ?    N    ? S      ? 20240924 1327 15.8282 GAU  2.00e-01 -1.00e+00 -1.00e+00 -1.00e+00  1.00e+00
...

With NLLOC_OBS format, the NLL-hypocenter-phase FOCALMECH line is also read if present.

If you add the FILENAME_PUBLIC_ID identifier in the LOCHYPOUT control statement, the public_id will be appended to the data-time in the NLL event filename.

I hope these (not really documented) options help!

Best regards,
Anthony

@loicbachelot
Copy link
Author

Hello @alomax,

The combination of FILENAME_PUBLIC_ID identifier in the LOCHYPOUT control statement and the PUBLIC_ID in the NLLOC_OBS is exactly the behavior I was looking for, thanks for the tip!

Best,
Loïc

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants