gea-program 0.11.0__tar.gz → 0.12.0__tar.gz

{gea_program-0.11.0/gea_program.egg-info → gea_program-0.12.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: gea-program
-Version: 0.11.0
+Version: 0.12.0
 Summary: GEA-Program: downhole gauge monitoring - canonical measurement records, quality rules, drift evaluation with a re-fit change log, well-test validation, alarm management, model cards, store-and-forward, monthly SLA measurement, FAT/SAT protocols and a dashboard, generated from live objects in the client's reporting outline.
 Author: Daniel T. Murphy
 License: MPL-2.0
@@ -125,7 +125,7 @@ gea dashboard --catalog-well volve_f12_f14_production_excerpt:15/9-F-12:10000 --
 gea client-report --report accuracy --out client_report
 gea model-cards --out model_cards
 gea sbom --out sbom
-gea accept                     # the product gate (320 checks)
+gea accept                     # the product gate (324 checks)
 gea help drift                 # the help library, by the job (16 pages; the same text is on every dashboard page)
 gea guide                      # the click-by-click tester guide (docs/TESTER_GUIDE.md)
 gea gui                        # the desktop window (pip install "gea-program[desktop]")
@@ -566,6 +566,24 @@ operation between them written out; a point on an unknown datum gets no WGS
 84 coordinates and is named. Without the operator's partition, ACL and legal
 tag the manifest is NOT LOADABLE and says so. `gea help osdu`.
 
+### The QuakeML catalogue export
+
+The site's own events in the format the regulator's catalogue tools read.
+After `associate`, `gea workspace --action quakeml-export --site <id>
+[--agency "<name>"]` writes `reports/sites/<site>/quakeml/<site>_events.xml`
+as QuakeML 1.2 to the Basic Event Description schema: one event per
+associated event with its origin (the misfit region as the uncertainties,
+depth in metres and "operator assigned" because it was declared, the method
+and the earth model named, the quality block with the counts, the RMS, the
+azimuthal gap and the distances), one pick per station with the record's own
+FDSN codes, one arrival per pick with its residual. Every origin and pick is
+automatic and preliminary in the schema's own words; the catalogue's AGREES
+or DIFFERS is a comment with the catalogue's event id; no magnitude is
+written and a comment says why. The shape is checked the way a reader checks
+it first - identifiers, allowed children, required children, enumerations,
+references - and the file validates against the published XSD. `gea help
+quakeml`; `gea quakeml` runs the labelled scene.
+
 ### Audit / Update
 
 The Audit / Update tab holds three things. The whole audit log

{gea_program-0.11.0 → gea_program-0.12.0}/README.md

@@ -86,7 +86,7 @@ gea dashboard --catalog-well volve_f12_f14_production_excerpt:15/9-F-12:10000 --
 gea client-report --report accuracy --out client_report
 gea model-cards --out model_cards
 gea sbom --out sbom
-gea accept                     # the product gate (320 checks)
+gea accept                     # the product gate (324 checks)
 gea help drift                 # the help library, by the job (16 pages; the same text is on every dashboard page)
 gea guide                      # the click-by-click tester guide (docs/TESTER_GUIDE.md)
 gea gui                        # the desktop window (pip install "gea-program[desktop]")
@@ -527,6 +527,24 @@ operation between them written out; a point on an unknown datum gets no WGS
 84 coordinates and is named. Without the operator's partition, ACL and legal
 tag the manifest is NOT LOADABLE and says so. `gea help osdu`.
 
+### The QuakeML catalogue export
+
+The site's own events in the format the regulator's catalogue tools read.
+After `associate`, `gea workspace --action quakeml-export --site <id>
+[--agency "<name>"]` writes `reports/sites/<site>/quakeml/<site>_events.xml`
+as QuakeML 1.2 to the Basic Event Description schema: one event per
+associated event with its origin (the misfit region as the uncertainties,
+depth in metres and "operator assigned" because it was declared, the method
+and the earth model named, the quality block with the counts, the RMS, the
+azimuthal gap and the distances), one pick per station with the record's own
+FDSN codes, one arrival per pick with its residual. Every origin and pick is
+automatic and preliminary in the schema's own words; the catalogue's AGREES
+or DIFFERS is a comment with the catalogue's event id; no magnitude is
+written and a comment says why. The shape is checked the way a reader checks
+it first - identifiers, allowed children, required children, enumerations,
+references - and the file validates against the published XSD. `gea help
+quakeml`; `gea quakeml` runs the labelled scene.
+
 ### Audit / Update
 
 The Audit / Update tab holds three things. The whole audit log

{gea_program-0.11.0 → gea_program-0.12.0}/gea/__init__.py

@@ -122,7 +122,7 @@ from .profile_catalog import (
     read_pangaea_txt,
 )
 
-__version__ = "0.11.0"
+__version__ = "0.12.0"
 __all__ = [
     "aging_rate", "rate_psi_yr", "accuracy_psi",
     "Sensor", "SimulatorConfig", "DownholeEngine",

{gea_program-0.11.0 → gea_program-0.12.0}/gea/__main__.py

@@ -282,7 +282,7 @@ def main(argv=None) -> int:
                       choices=["init", "add-file", "add-catalog", "add-live", "remove", "list", "migrate", "refresh", "audit",
                                "add-seismic", "refresh-seismic", "remove-seismic", "sar-film", "refresh-all", "add-track", "refresh-track", "remove-track",
                                "seismic-dataset", "seismic-field", "add-site", "sites", "site-report", "remove-site", "rename",
-                               "set-well", "sra-define", "sra-report", "associate", "osdu-export"])
+                               "set-well", "sra-define", "sra-report", "associate", "osdu-export", "quakeml-export"])
     p_ws.add_argument("--stations", type=str, nargs="+", default=None, help="add-track: two or more array station ids of this workspace")
     p_ws.add_argument("--site", type=str, default=None, help="site-report/remove-site: the site id")
     p_ws.add_argument("--api", type=str, default=None, help="set-well: the API number")
@@ -315,7 +315,8 @@ def main(argv=None) -> int:
     p_ws.add_argument("--legal-tag", type=str, nargs="*", default=None, help="osdu-export: legal tag(s)")
     p_ws.add_argument("--country", type=str, nargs="*", default=None, help="osdu-export: otherRelevantDataCountries (default US)")
     p_ws.add_argument("--operator-org", type=str, default="", help="osdu-export: the operator's Organisation record id on the platform")
-    p_ws.add_argument("--out", type=str, default=None, help="osdu-export: the folder to write the manifest and files into")
+    p_ws.add_argument("--out", type=str, default=None, help="osdu-export: the folder to write the manifest and files into; quakeml-export: the file to write")
+    p_ws.add_argument("--agency", type=str, default="", help="quakeml-export: the agencyID written in creationInfo (default: the site's name)")
     p_ws.add_argument("--depth-km", type=float, default=6.0, help="associate: the hypocentral depth the location is made under, km (declared)")
     p_ws.add_argument("--rms-tol", type=float, default=0.15, help="associate: the RMS residual an event must fit within, s")
     p_ws.add_argument("--wells", type=str, nargs="+", default=None, help="add-site: well ids of this workspace")
@@ -399,6 +400,9 @@ def main(argv=None) -> int:
 
     p_osdu = sub.add_parser("osdu", help="the OSDU-shaped export: the self-test on a labelled site")
     p_osdu.add_argument("--action", choices=["selftest"], default="selftest")
+    p_qml = sub.add_parser("quakeml", help="the QuakeML 1.2 catalogue export: the self-test on the association leg's labelled scene")
+    p_qml.add_argument("--action", choices=["selftest"], default="selftest")
+    p_qml.add_argument("--out", type=str, default=None, help="where to leave the self-test's file (default: a temporary folder)")
     p_sra = sub.add_parser("sra", help="the Seismicity Response Area packet: the self-test on a labelled scene")
     p_sra.add_argument("--action", choices=["selftest"], default="selftest")
     p_sra.add_argument("--json", action="store_true")
@@ -1496,6 +1500,14 @@ def main(argv=None) -> int:
             return 0
         print(f"notifications: {'configured' if n.cfg else 'not configured'}" + (f" - {len(n.cfg['channels'])} channel(s), {len(n.cfg['rules'])} rule(s), quiet {n.cfg['quiet_s']} s" if n.cfg else " (gea notify --example)"))
         return 0
+    elif a.cmd == "quakeml":
+        from . import quakeml as _Q
+        r = _Q.selftest(a.out)
+        print(f"quakeml selftest [{r['label']}]: {r['status']}")
+        for k, v in r["checks"].items():
+            print(f"  {'ok ' if v else 'BAD'} {k}")
+        print(_Q.report_text({**r, 'status': 'WRITTEN' if r['status'] == 'OK' else 'FAILED', 'agency': 'SIMULATION_SELF_TEST site', 'refused': 2}))
+        return 0 if r["status"] == "OK" else 1
     elif a.cmd == "osdu":
         from . import osdu as _O
         r = _O.selftest()
@@ -1753,6 +1765,12 @@ def main(argv=None) -> int:
                                    countries=a.country, operator_org_id=a.operator_org, out_dir=a.out, actor=a.actor)
                 print(_osdu_text(r)); v = r["validate"]
                 print(f"  structure: {'ok' if v['ok'] else 'PROBLEMS'} ({v['records']} record(s))" + ('' if v['ok'] else '; ' + '; '.join(v['problems'][:5]))); return 0
+            elif a.action == "quakeml-export":
+                from .quakeml import report_text as _qml_text
+                if not a.site:
+                    raise SystemExit("quakeml-export needs --site <id>")
+                r = ws.quakeml_export(a.site, agency=a.agency, out=a.out, actor=a.actor)
+                print(_qml_text(r)); return 0 if r["status"] == "WRITTEN" else 1
             elif a.action == "seismic-dataset":
                 r = ws.write_seismic_dataset_report(actor=a.actor)
                 inv = r["inventory"]

{gea_program-0.11.0 → gea_program-0.12.0}/gea/acceptance_tests.py

@@ -4473,6 +4473,7 @@ def section_bb_supervisor(tmp: str) -> None:
        "BB4 a code that is neither 0 nor 86 is a death, not a decision: the launcher reads the operator's answer from records/auto_restart.flag, "
        "tries at most five times in a row so a panel that crashes while starting hands over the error instead of flickering all night, and reads "
        "that file outside an if-block because a redirect inside parentheses is parsed before the block runs")
+    cmd_real = cmd                                   # the launcher as shipped (template or installed file), before the stand-in below
     if cmd is None:                                  # a Linux kit: BB6 needs a Windows launcher text to build its good and bad kits from
         cmd = ("@echo off\ntitle GEA - Operations Control Panel\n:gea_run\npython\\python.exe -m gea serve\nset GEA_RC=%ERRORLEVEL%\n"
                "if \"%GEA_RC%\"==\"86\" goto gea_run\n:gea_read_flag\nset /p GEA_AUTO=<\"%GEA_WORKSPACE%\\records\\auto_restart.flag\"\n"
@@ -4636,12 +4637,13 @@ def section_bb_supervisor(tmp: str) -> None:
     ok('"--no-browser", action="store_true"' in main_src
        and "svc.open_browser = not a.no_browser" in main_src and "svc.console = True" in main_src
        and hasattr(Service, "_open_browser_later") and Service(wsp, host="127.0.0.1", port=0, scheduler=False).open_browser is False
-       and (cmd is None or ("%GEA_BROWSER%" in cmd and "timeout /t 4" not in cmd and 'if "%1"=="--no-browser" set GEA_BROWSER=--no-browser' in cmd
-                            and "set GEA_BROWSER=--no-browser" in cmd))
+       and (cmd_real is None or ("%GEA_BROWSER%" in cmd_real and "timeout /t 4" not in cmd_real
+                                 and 'if "%1"=="--no-browser" set GEA_BROWSER=--no-browser' in cmd_real and "set GEA_BROWSER=--no-browser" in cmd_real))
        and (sh is None or ("$browser" in sh and 'browser="--no-browser"' in sh)),
        "BB13 `gea serve` opens the control panel in the browser itself once its port is listening, and the launchers no longer race it with a "
        "browser of their own - opened first, it showed a failure page from the second before the server was up; a scheduled start passes "
-       "--no-browser through, a restart from the panel does not open a second page, and the service run in-process by this gate opens nothing")
+       "--no-browser through, a restart from the panel does not open a second page, and the service run in-process by this gate opens nothing; a Linux "
+       "kit, which carries no Windows launcher, is judged on its shell launcher alone")
 
 
 def section_bc_sra(tmp: str) -> None:
@@ -4965,6 +4967,143 @@ def section_be_osdu(tmp: str) -> None:
        "file is a named gap; the structure validates; every export is audited; `gea osdu` runs the self-test from the command line")
 
 
+
+def section_bf_quakeml(tmp: str) -> None:
+    """Section BF - the QuakeML catalogue export: the site's own events in the exchange format a regulator's
+    catalogue tool reads, built to the Basic Event Description schema, automatic and preliminary in the
+    schema's own words, with no magnitude and with the catalogue's standing as a comment."""
+    import csv as _csv
+    import subprocess as _sp
+    import xml.etree.ElementTree as _ET
+    from . import quakeml as Q
+    from . import seismic as S
+    from . import seismic_assoc as SA
+    from .workspace import Workspace, WorkspaceError
+    from . import helplib as H
+    pkg = Path(__file__).parent
+    d = Path(tmp, "bf"); d.mkdir()
+    B = lambda t: "{%s}%s" % (Q.NS_BED, t)
+
+    # -- the self-test on the labelled scene, read back ------------------------------------------
+    r = Q.selftest(str(d / "self"))
+    root = _ET.parse(r["path"]).getroot()
+    origins = root.findall(".//" + B("origin"))
+    o = origins[0]
+    lat_u = float(o.find(B("latitude")).find(B("uncertainty")).text); lon_u = float(o.find(B("longitude")).find(B("uncertainty")).text)
+    hu = float(o.find(B("originUncertainty")).find(B("horizontalUncertainty")).text)
+    q = o.find(B("quality"))
+    ok(r["status"] == "OK" and all(r["checks"].values()) and r["label"] == "SIMULATION_SELF_TEST"
+       and root.tag == "{%s}quakeml" % Q.NS_Q and len(origins) == 2
+       and 0 < lat_u < 0.1 and 0 < lon_u < 0.1 and 0 < hu < 10000 and abs(hu - max(lat_u * Q.KM_PER_DEG, 0) * 1000) < 1.0 + hu * 0.5
+       and q.find(B("usedStationCount")).text == "5" and q.find(B("usedPhaseCount")).text == "5"
+       and float(q.find(B("standardError")).text) < 0.05 and 0 < float(q.find(B("azimuthalGap")).text) < 180
+       and 0 < float(q.find(B("minimumDistance")).text) < float(q.find(B("maximumDistance")).text) < 0.2
+       and o.find(B("methodID")).text.startswith("smi:local/gea/method/") and "vp-" in o.find(B("earthModelID")).text
+       and all(float(a.find(B("timeWeight")).text) == 1.0 and abs(float(a.find(B("timeResidual")).text)) < 0.15 and 0 < float(a.find(B("distance")).text) < 0.2
+               for a in o.findall(B("arrival"))),
+       "BF1 the self-test on the association leg's labelled scene: two events written and read back under the QuakeML 1.2 namespaces, five picks "
+       "and five arrivals each, depth in metres with depthType 'operator assigned' because it was declared, the misfit region as the latitude and "
+       "longitude uncertainties in degrees and the horizontal uncertainty in metres, the quality block with the counts, the RMS, the azimuthal gap "
+       "and the epicentral distances in degrees, every arrival with its residual and unit weight, the method and the earth model named, every "
+       "origin and pick automatic and preliminary, the catalogue's AGREES and DIFFERS as comments, the record's own FDSN codes in every waveformID, "
+       "and no magnitude anywhere")
+
+    # -- the shape check catches what a reader would refuse ------------------------------------------
+    text = Path(r["path"]).read_text(encoding="utf-8")
+    broken = []
+    bad1 = text.replace('publicID="smi:local/gea/SYNTHETIC/event/', 'publicID="bad id ', 1); broken.append(Q.validate(_ET.ElementTree(_ET.fromstring(bad1))))
+    bad2 = text.replace("<pickID>smi:local/gea/SYNTHETIC/pick/", "<pickID>smi:local/gea/SYNTHETIC/nopick/", 1); broken.append(Q.validate(_ET.ElementTree(_ET.fromstring(bad2))))
+    bad3 = text.replace("<evaluationStatus>preliminary</evaluationStatus>", "<evaluationStatus>reported</evaluationStatus>", 1); broken.append(Q.validate(_ET.ElementTree(_ET.fromstring(bad3))))
+    bad4 = text.replace("<depthType>", "<depthKind>").replace("</depthType>", "</depthKind>"); broken.append(Q.validate(_ET.ElementTree(_ET.fromstring(bad4))))
+    first_origin_end = text.index("</origin>")
+    bad5 = text[:first_origin_end] + '<magnitude publicID="smi:local/gea/x/mag/1"><mag><value>2.0</value></mag></magnitude>' + text[first_origin_end:]
+    broken.append(Q.validate(_ET.ElementTree(_ET.fromstring(bad5))))
+    bad6 = text.replace("<value>2026-03-05T11:03:", "<value>2026-03-05 11:03:", 1); broken.append(Q.validate(_ET.ElementTree(_ET.fromstring(bad6))))
+    ok(Q.validate(r["path"])["status"] == "VALID_SHAPE" and all(b["status"] == "INVALID" for b in broken)
+       and any("ResourceIdentifier" in p for p in broken[0]["problems"]) and any("pickID" in p for p in broken[1]["problems"])
+       and any("enumeration" in p for p in broken[2]["problems"]) and any("not in the schema" in p for p in broken[3]["problems"])
+       and any("magnitude" in p for p in broken[4]["problems"]) and any("dateTime" in p for p in broken[5]["problems"])
+       and all(Q.RID_PATTERN.match(x) for x in (Q.rid("Pad 3", "event", "assoc 1"), Q.rid("a", "x", authority="ab"), Q.rid("s", "comment", "a/b")))
+       and not Q.RID_PATTERN.match("smi:ab/x") and not Q.RID_PATTERN.match("urn:local/x"),
+       "BF2 the shape check is the first thing a QuakeML reader does, done here: an identifier off the schema's pattern, an arrival whose pick is "
+       "not in its event, an evaluation status outside the enumeration, a child element the schema does not know, a magnitude this program did "
+       "not estimate, and a time that is not an xs:dateTime are each named as a problem; and the identifier builder cleans a site name with a "
+       "space and a short authority to the pattern")
+
+    # -- through the workspace and the command line ---------------------------------------------------
+    sc = SA.synthetic_scene()
+    for sid, tr in sc["traces"].items():
+        S.write_mseed([tr], str(d / f"{sid}.mseed"))
+    cat = d / "cat.csv"
+    with open(cat, "w", newline="") as f:
+        w = _csv.writer(f); w.writerow(["EventID", "Origin Time", "Local Magnitude", "Latitude", "Longitude", "Depth"])
+        for c in sc["catalog"]:
+            w.writerow([c["event_id"], c["time"], c["magnitude"], c["lat"], c["lon"], 6.0])
+    ws = Workspace.create(str(d / "ws"), "Pad 3", actor="tester")
+    ids = [ws.add_seismic_station(sid, [str(d / f"{sid}.mseed")], st.lat, st.lon, "tester", datum=st.datum, band=[2.0, 20.0])["id"]
+           for sid, st in sc["stations"].items()]
+    ws.add_site("QML block", seismic=ids, client="Acme", actor="tester")
+    try:
+        ws.quakeml_export("QML-block"); early = False
+    except WorkspaceError as e:
+        early = "associate" in str(e)
+    run = lambda *a: _sp.run([sys.executable, "-m", "gea", "workspace", "--path", ws.path, "--actor", "tester", *a], capture_output=True, text=True)
+    r1 = run("--action", "associate", "--site", "QML-block", "--catalog", str(cat))
+    aj = json.loads(Path(ws.path, "reports", "sites", "QML-block", "association.json").read_text(encoding="utf-8"))
+    r2 = run("--action", "quakeml-export", "--site", "QML-block", "--agency", "Acme Operating")
+    xml_path = Path(ws.path, "reports", "sites", "QML-block", "quakeml", "QML-block_events.xml")
+    summ = json.loads(Path(ws.path, "reports", "sites", "QML-block", "quakeml", "export_summary.json").read_text(encoding="utf-8"))
+    root2 = _ET.parse(str(xml_path)).getroot()
+    agencies = {c.find(B("agencyID")).text for c in root2.iter(B("creationInfo"))}
+    wids = root2.findall(".//" + B("waveformID"))
+    audit = Path(ws.path, "records", "audit.jsonl").read_text(encoding="utf-8")
+    r3 = _sp.run([sys.executable, "-m", "gea", "quakeml", "--out", str(d / "cli")], capture_output=True, text=True)
+    help_md = (pkg / "help" / "quakeml.md").read_text(encoding="utf-8")
+    ok(early and r1.returncode == 0 and all((v.get("stream") or {}).get("station") == sid for sid, v in aj["stations"].items())
+       and r2.returncode == 0 and "WRITTEN - 2 event(s), 10 pick(s), 10 arrival(s)" in r2.stdout and "VALID_SHAPE" in r2.stdout and "gap:" not in r2.stdout
+       and xml_path.is_file() and summ["status"] == "WRITTEN" and summ["schema"] == "QuakeML-BED-1.2" and summ["n_events"] == 2 and not summ["gaps"]
+       and summ["validate"]["status"] == "VALID_SHAPE" and summ["association_generated_utc"] == aj["generated_utc"]
+       and agencies == {"Acme Operating"} and len(wids) == 10 and all(w.get("networkCode") == "XX" and w.get("channelCode") == "HHZ" for w in wids)
+       and "site.quakeml_export" in audit and r3.returncode == 0 and "SIMULATION_SELF_TEST" in r3.stdout and (d / "cli" / "SYNTHETIC_events.xml").is_file()
+       and H.VIEW_TOPIC.get("quakeml", "seismic") == "seismic" and any(t["topic"] == "quakeml" for t in H.topics()) and not H.check("quakeml")
+       and "operator assigned" in help_md and "preliminary" in help_md and "magnitude" in help_md.split("will not call a measurement")[1],
+       "BF3 `--action quakeml-export` refuses a site that has not been associated, and after `--action associate` (which now keeps each record's "
+       "FDSN codes) writes reports/sites/<site>/quakeml/<site>_events.xml with the agency the operator named in every creationInfo, the record's "
+       "codes in every waveformID, a summary beside it that carries the association's own timestamp, and audits it; `gea quakeml` runs the "
+       "labelled scene from the command line; the help page names the declared depth, the preliminary status and the magnitude it does not write")
+
+    # -- a station without codes or position is named as a gap, never invented ------------------------
+    doc = {"association": {"events": [{"id": "e1", "origin_utc": "2026-03-05T11:03:00.000Z", "lat": 31.95, "lon": -102.25, "depth_km": 6.0,
+                                        "depth_basis": "declared", "rms_s": 0.02, "n_stations": 3, "stations": ["A", "B", "C"], "dropped": [],
+                                        "picks": {"A": {"time": "2026-03-05T11:03:01.000Z", "snr": 10.0, "residual_s": 0.01},
+                                                  "B": {"time": "2026-03-05T11:03:01.200Z", "snr": 9.0, "residual_s": -0.01},
+                                                  "C": {"time": "2026-03-05T11:03:01.400Z", "snr": 8.0, "residual_s": 0.0}},
+                                        "region_km": {"extent_east_km": 1.0, "extent_north_km": 2.0, "nodes": 8, "reaches_grid_edge": False},
+                                        "magnitude": None, "magnitude_basis": "not estimated",
+                                        "model": {"vp_km_s": 5.8, "depth_km": 6.0, "pick_sigma_s": 0.05, "basis": "flat earth"}}],
+                           "refused": [], "n_unused_picks": 0, "model": {}, "basis": "b", "not_a_measurement": []}, "catalogue": None}
+    stations = {"A": {"id": "A", "lat": 31.99, "lon": -102.30, "stream": {"network": "TX", "station": "PB01", "location": "00", "channel": "HHZ"}},
+                "B": {"id": "B", "lat": 31.90, "lon": -102.20},                       # no codes
+                "C": {"id": "C"}}                                                       # no codes, no position
+    b = Q.build(doc, {"id": "s", "name": "Gap site"}, stations, version="0.0")
+    out = Path(d, "gaps.xml"); Q.write(b["tree"], str(out))
+    v = Q.validate(str(out))
+    rootg = _ET.parse(str(out)).getroot()
+    ws_ids = {w.get("stationCode"): w for w in rootg.findall(".//" + B("waveformID"))}
+    arrs = {a.find(B("pickID")).text.rsplit("/", 1)[1]: a for a in rootg.findall(".//" + B("arrival"))}
+    og = rootg.find(".//" + B("origin"))
+    ok(v["status"] == "VALID_SHAPE" and len(b["gaps"]) == 3 and sum("no network and station code" in g for g in b["gaps"]) == 2
+       and sum("no position" in g for g in b["gaps"]) == 1
+       and ws_ids["PB01"].get("networkCode") == "TX" and ws_ids["PB01"].get("locationCode") == "00"
+       and ws_ids["B"].get("networkCode") == "XX" and ws_ids["C"].get("networkCode") == "XX"
+       and arrs["A"].find(B("azimuth")) is not None and arrs["C"].find(B("azimuth")) is None and arrs["C"].find(B("distance")) is None
+       and abs(float(og.find(B("originUncertainty")).find(B("horizontalUncertainty")).text) - 1000.0) < 1e-6
+       and abs(float(og.find(B("latitude")).find(B("uncertainty")).text) - 1.0 / Q.KM_PER_DEG) < 1e-6
+       and og.find(B("quality")).find(B("azimuthalGap")) is not None,
+       "BF4 a station whose record carried no network and station code goes out under the test network XX with its own id and is named as a gap; "
+       "a station with no position gets a pick but an arrival without azimuth or distance and is named as a gap; neither is invented; the misfit "
+       "region of 1 km by 2 km becomes a 1000 m horizontal uncertainty and a latitude uncertainty of one kilometre in degrees")
+
 def main() -> int:
     print("GEA-Program - ACCEPTANCE SUITE (the product gate)")
     # a handle a job subprocess still holds at the end must not turn a finished gate into a traceback on Windows:
@@ -5025,6 +5164,7 @@ def main() -> int:
         section_bc_sra(tmp)
         section_bd_association(tmp)
         section_be_osdu(tmp)
+        section_bf_quakeml(tmp)
     if _FAILS:
         print(f"[ACCEPTANCE] {len(_FAILS)} FAILURES ({_PASS} passed):")
         for f in _FAILS:

gea_program-0.12.0/gea/help/quakeml.md

@@ -0,0 +1,11 @@
+# The QuakeML catalogue export
+
+> The site's own events - what its stations located together - in QuakeML 1.2, the exchange format a regulator's catalogue tool, a university network's review desk and every seismological package read.
+
+**The command**: `gea workspace --path C:\site --action quakeml-export --site <id> [--agency "<name>"] [--out <file.xml>]`, after `--action associate` has run on the site (the export is built from `reports/sites/<site>/association.json`). `--agency` is the agencyID written in every creationInfo; the site's name by default. `gea quakeml` runs the self-test on the association leg's labelled scene and leaves the file where `--out` says.
+
+**What it writes**: `reports/sites/<site>/quakeml/<site>_events.xml` - a `q:quakeml` document in the QuakeML 1.2 namespaces, built to the Basic Event Description schema (QuakeML-BED-1.2.xsd): one `event` per associated event with its `origin` (time; latitude and longitude with the misfit region's half-extent as their uncertainty in degrees; depth in metres with `depthType` "operator assigned", because the depth was declared; `methodID` and `earthModelID` naming the grid search and the one-velocity flat-earth model; `quality` with the station and phase counts, the RMS as standardError, the azimuthal gap and the epicentral distances in degrees; `originUncertainty` with the misfit region in metres), one `pick` per station that made the event (the onset time with the pick uncertainty, the record's own FDSN codes in `waveformID`, phaseHint P), and one `arrival` per pick on the origin (azimuth, distance, the residual, weight 1). Every origin and pick is `evaluationMode` automatic and `evaluationStatus` preliminary. The catalogue's standing - AGREES, DIFFERS with both positions, NOT IN CATALOGUE - is a comment on the event with the catalogue's own event id. `export_summary.json` beside it carries the counts, the gaps and the shape check. Every export is in `records/audit.jsonl`.
+
+**The number to check**: `shape VALID_SHAPE` on the second line - the first checks a QuakeML reader makes, made here first: the namespaces; every publicID to the schema's ResourceIdentifier pattern and used once; every child element one the schema allows for its parent; time, latitude and longitude on every origin; time and waveformID on every pick; every arrival's pickID a pick of the same event; every event's preferredOriginID one of its origins; every enumerated value in its enumeration. A `gap:` line names a station whose record carried no network and station code (written under the test network XX with the station id) or no position (its arrival carries no azimuth or distance).
+
+**What this page will not call a measurement**: a magnitude - the association does not estimate one, the file carries no `magnitude` element, and a comment on every event says why. The catalogue's events as this site's - they are the catalogue's to publish; the file names them in comments and carries only what the site's own stations located. A reviewed or final event - a person has not looked at these, and the file says automatic and preliminary in the schema's own words. A located depth - it was declared, and `depthType` says so. A full schema validation - the shape check is what a reader reads first, not a validator; the file was validated against the published XSD when this export was written, and the check here holds it to that shape.

{gea_program-0.11.0 → gea_program-0.12.0}/gea/helplib.py

@@ -40,6 +40,7 @@ INDEX: List[tuple] = [
     ('notifications', 'Notifications', 'rules, channels, the delivery log'),
     ('sra', 'The Seismicity Response Area', 'the packet an operator inside an SRA puts in front of the Commission'),
     ('osdu', 'The OSDU-shaped export', 'the site as the manifest an operator\'s data platform loads'),
+    ('quakeml', 'The QuakeML catalogue export', 'the site\'s own events in the format the regulator\'s catalogue tools read'),
     ('doctor', 'Which code is running', 'gea doctor and the serve start-up gate'),
     ('audit-update', 'Audit / Update', 'the whole audit log with filters, the program against PyPI, every report against its source'),
     ('upkeep', 'Upkeep', 'housekeeping, sessions, the TLS proxy, the load test'),