gea-program 0.34.0__tar.gz → 0.36.0__tar.gz

{gea_program-0.34.0/gea_program.egg-info → gea_program-0.36.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: gea-program
-Version: 0.34.0
+Version: 0.36.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
@@ -648,6 +648,40 @@ against it; `gea seedlink --config ... --hello` prints what a server says of
 itself; `gea help seedlink`. A TexNet station is served by EarthScope at
 `rtserve.iris.washington.edu:18000`, network `TX`.
 
+### The well in 3D
+
+The 3D tab draws each well's path from its directional survey by the
+minimum-curvature method, below sea level from its declared rig floor, among
+the mapped formations - interpreted horizon surfaces, or a formation-top
+map's contours at their elevations - with its picks on it. Colour the wells
+by formation, inclination, dogleg severity, or a drilling recorder's gamma,
+rate of penetration or weight on bit along the bit depth; replay the bit
+going down the hole from the recorder's bit depth or the survey stations'
+own times. Turn, pan and zoom by pointer, touch or keys. The Volve field
+(Equinor's open data: three original surveys, the field's formation picks,
+five mapped horizons) ships as the worked example, checked against itself.
+Until then the workspace's formation-top contours are drawn on their own.
+A workspace well is drawn once it has a position, a reference elevation and
+a survey: `gea trajectory` attaches one, verbatim, with its azimuth
+reference stated. `gea well3d`, `gea help well3d`.
+
+### The control room
+
+The Control room tab opens the wall: a screen read from across the room, one
+tile per well, worst first, each with its state in a shape and a word -
+Serious, Watch, Good, No data - the reasons in plain words, the alarms, the
+gauges, and for a drilling recorder the bottom, the bit, the rate of
+penetration and the hours on bottom at the end of the record. It refreshes
+itself every 30 seconds, says how old its last answer is, and puts a notice
+across the screen when the service stops answering; a recorded well is
+described at the end of its record and never called stale. F is full screen;
+the screen is kept awake where the browser allows it; when the wells do not
+fit, the pages turn, the worst first. The outgoing operator writes the shift
+handover - what is open, what is being watched, what was done - and it stays
+on the wall until the incoming operator acknowledges it; both go in the audit
+log. A wall pass lets a screen show the wall and nothing else with no one's
+account signed in. `gea wall`, `gea handover`, `gea help control-room`.
+
 ### The site map
 
 The Map tab puts everything in the workspace that has a position in one

{gea_program-0.34.0 → gea_program-0.36.0}/README.md

@@ -609,6 +609,40 @@ against it; `gea seedlink --config ... --hello` prints what a server says of
 itself; `gea help seedlink`. A TexNet station is served by EarthScope at
 `rtserve.iris.washington.edu:18000`, network `TX`.
 
+### The well in 3D
+
+The 3D tab draws each well's path from its directional survey by the
+minimum-curvature method, below sea level from its declared rig floor, among
+the mapped formations - interpreted horizon surfaces, or a formation-top
+map's contours at their elevations - with its picks on it. Colour the wells
+by formation, inclination, dogleg severity, or a drilling recorder's gamma,
+rate of penetration or weight on bit along the bit depth; replay the bit
+going down the hole from the recorder's bit depth or the survey stations'
+own times. Turn, pan and zoom by pointer, touch or keys. The Volve field
+(Equinor's open data: three original surveys, the field's formation picks,
+five mapped horizons) ships as the worked example, checked against itself.
+Until then the workspace's formation-top contours are drawn on their own.
+A workspace well is drawn once it has a position, a reference elevation and
+a survey: `gea trajectory` attaches one, verbatim, with its azimuth
+reference stated. `gea well3d`, `gea help well3d`.
+
+### The control room
+
+The Control room tab opens the wall: a screen read from across the room, one
+tile per well, worst first, each with its state in a shape and a word -
+Serious, Watch, Good, No data - the reasons in plain words, the alarms, the
+gauges, and for a drilling recorder the bottom, the bit, the rate of
+penetration and the hours on bottom at the end of the record. It refreshes
+itself every 30 seconds, says how old its last answer is, and puts a notice
+across the screen when the service stops answering; a recorded well is
+described at the end of its record and never called stale. F is full screen;
+the screen is kept awake where the browser allows it; when the wells do not
+fit, the pages turn, the worst first. The outgoing operator writes the shift
+handover - what is open, what is being watched, what was done - and it stays
+on the wall until the incoming operator acknowledges it; both go in the audit
+log. A wall pass lets a screen show the wall and nothing else with no one's
+account signed in. `gea wall`, `gea handover`, `gea help control-room`.
+
 ### The site map
 
 The Map tab puts everything in the workspace that has a position in one

{gea_program-0.34.0 → gea_program-0.36.0}/gea/__init__.py

@@ -123,7 +123,7 @@ from .profile_catalog import (
     read_pangaea_txt,
 )
 
-__version__ = "0.34.0"
+__version__ = "0.36.0"
 __all__ = [
     "aging_rate", "rate_psi_yr", "accuracy_psi",
     "Sensor", "SimulatorConfig", "DownholeEngine",

{gea_program-0.34.0 → gea_program-0.36.0}/gea/__main__.py

@@ -475,6 +475,38 @@ def main(argv=None) -> int:
     p_ck.add_argument("--note", type=str, default="", help="--confirm: what it was checked against (required)")
     p_ck.add_argument("--actor", type=str, default="cli")
 
+    p_tj = sub.add_parser("trajectory", help="attach a well's directional survey (WITSML trajectory or a table of MD, Inc, Azi), kept verbatim beside the well, so the 3D tab can draw its path")
+    p_tj.add_argument("--workspace", type=str, required=True)
+    p_tj.add_argument("--well", type=str, required=True)
+    p_tj.add_argument("--file", type=str, required=True)
+    p_tj.add_argument("--azimuth-reference", type=str, required=True, choices=("grid", "true", "magnetic"))
+    p_tj.add_argument("--grid", type=str, default=None, help="grid azimuths: the grid, e.g. TX_CENTRAL or UTM14N")
+    p_tj.add_argument("--declination", type=float, default=None, help="magnetic azimuths: the declination used, degrees east positive")
+    p_tj.add_argument("--md-unit", type=str, default="ft", choices=("ft", "m"))
+    p_tj.add_argument("--depth-reference", type=str, default=None, help="KB, DF or GL: must match the well's declared reference elevation")
+    p_tj.add_argument("--source", type=str, required=True, help='where the survey came from, e.g. "RRC imaged records, directional survey"')
+    p_tj.add_argument("--actor", type=str, default="cli")
+
+    p_3d = sub.add_parser("well3d", help="the 3D scene in words: each well's path, how it was checked, the surfaces and the replay (the 3D tab draws it)")
+    p_3d.add_argument("--scene", type=str, default="volve", choices=("volve", "workspace"))
+    p_3d.add_argument("--workspace", type=str, default=None)
+    p_3d.add_argument("--json", action="store_true")
+
+    p_wl = sub.add_parser("wall", help="the control-room board as text: a line per well, worst first, its state and why (the same board the wall screen shows)")
+    p_wl.add_argument("--workspace", type=str, default=None, help="the workspace folder, e.g. C:\\site")
+    p_wl.add_argument("--json", action="store_true")
+
+    p_ho = sub.add_parser("handover", help="the shift handover: list the notes, sign one (--write) or acknowledge one (--ack)")
+    p_ho.add_argument("--workspace", type=str, default=None, help="the workspace folder, e.g. C:\\site")
+    p_ho.add_argument("--write", action="store_true", help="sign a note: give at least one of --open, --watching, --done")
+    p_ho.add_argument("--open", type=str, default="", help="--write: what the next shift must pick up")
+    p_ho.add_argument("--watching", type=str, default="", help="--write: what is being watched, and why")
+    p_ho.add_argument("--done", type=str, default="", help="--write: what this shift did")
+    p_ho.add_argument("--ack", type=str, default=None, help="acknowledge the note with this id (not your own)")
+    p_ho.add_argument("--comment", type=str, default="", help="--ack: an optional comment")
+    p_ho.add_argument("--actor", type=str, default=None, help="who signs or acknowledges (required with --write and --ack)")
+    p_ho.add_argument("--json", action="store_true")
+
     p_sv = sub.add_parser("serve", help="serve the dashboard over a workspace: the page, the reports and the JSON API (standard library HTTP)")
     p_sv.add_argument("--workspace", type=str, required=True, help="the workspace folder (gea workspace --action init)")
     p_sv.add_argument("--host", type=str, default="127.0.0.1", help="bind address; 0.0.0.0 exposes it on the site network (put TLS in front)")
@@ -1922,6 +1954,79 @@ def main(argv=None) -> int:
         else:
             print(_SC.text(r))
         return 0 if r["ready"] else 2
+    elif a.cmd in ("trajectory", "well3d"):
+        import json as _json
+        from . import well3d as _W3
+        from .cli import _safe_stdout
+        from .workspace import Workspace, WorkspaceError
+        _safe_stdout()
+        ws = None
+        if a.cmd == "trajectory" or a.scene == "workspace":
+            if not a.workspace:
+                raise SystemExit(f"{a.cmd} needs --workspace, e.g. --workspace C:\\site")
+            try:
+                ws = Workspace(a.workspace)
+            except WorkspaceError as e:
+                raise SystemExit(f"{a.cmd}: {e}")
+        if a.cmd == "trajectory":
+            try:
+                with open(a.file, "rb") as f:
+                    content = f.read()
+                r = _W3.attach_trajectory(ws, a.well, a.file, content, a.actor, a.azimuth_reference, a.md_unit, a.depth_reference, a.grid,
+                                          a.declination, a.source)
+            except (ValueError, OSError, WorkspaceError) as e:
+                raise SystemExit(f"trajectory: {e}")
+            print(f"{a.well}: {r['stations']} stations, {r['md_range_m'][0]:,.1f} to {r['md_range_m'][1]:,.1f} m MD, azimuths {r['azimuth_reference']}"
+                  + (f" ({r['grid']})" if r.get('grid') else "") + f"; kept as wells/{a.well}/trajectory/{r['file']}")
+            return 0
+        sc = _W3.scene(ws, a.scene)
+        print(_json.dumps(sc, indent=1) if a.json else _W3.text(sc))
+        return 0
+    elif a.cmd in ("wall", "handover"):
+        import json as _json
+        from . import control_room as _CR
+        from .cli import _safe_stdout
+        from .workspace import Workspace, WorkspaceError
+        _safe_stdout()
+        if not a.workspace:
+            raise SystemExit(f"{a.cmd} needs --workspace, e.g. --workspace C:\\site")
+        try:
+            ws = Workspace(a.workspace)
+        except WorkspaceError as e:
+            raise SystemExit(f"{a.cmd}: {e}")
+        if a.cmd == "wall":
+            b = _CR.board(ws)
+            print(_json.dumps(b, indent=1) if a.json else _CR.text(b))
+            return 0
+        try:
+            if a.write:
+                if not a.actor:
+                    raise SystemExit("handover --write needs --actor, the name the note is signed with")
+                b = _CR.board(ws)
+                counts = {**b["counts"], "alarms_unacked": sum(((t.get("alarms") or {}).get("unacked") or 0) for t in b["tiles"])}
+                n = _CR.write_note(ws, a.actor, a.open, a.watching, a.done, counts)
+                print(f"handover {n['id']} signed by {n['by']} at {n['utc']}")
+                return 0
+            if a.ack:
+                if not a.actor:
+                    raise SystemExit("handover --ack needs --actor, the incoming operator's name")
+                r = _CR.ack_note(ws, a.ack, a.actor, a.comment)
+                print(f"handover {a.ack} acknowledged by {r['by']} at {r['utc']}")
+                return 0
+        except (ValueError, KeyError, PermissionError) as e:
+            raise SystemExit(f"handover: {e.args[0] if e.args else e}")
+        ns = _CR.notes(ws, 0)
+        if a.json:
+            print(_json.dumps(ns, indent=1))
+        elif not ns:
+            print("no handover notes yet")
+        else:
+            for n in ns:
+                print(f"{n['id']}  {n['utc']}  by {n['by']}  " + (f"acknowledged by {n['ack']['by']} at {n['ack']['utc']}" if n.get('ack') else "NOT YET ACKNOWLEDGED"))
+                for k in _CR.NOTE_FIELDS:
+                    if n.get(k):
+                        print(f"    {k}: {n[k]}")
+        return 0
     elif a.cmd == "serve":
         from .service import Service
         from .workspace import WorkspaceError

{gea_program-0.34.0 → gea_program-0.36.0}/gea/acceptance_tests.py

@@ -8448,6 +8448,358 @@ def section_cc_site_map(tmp: str) -> None:
        "bar, the not-placed wells with the way to place them, what is at a click; the help page"
        + (f" [failed: {', '.join(failed3)}]" if failed3 else ""))
 
+
+def section_cd_control_room(tmp: str) -> None:
+    """Section CD - the control room: the wall board judged by its stated rule (a shape and a word for each state, worst
+    first, a recorded well never called stale, a live feed's age measured), the shift handover signed by one operator and
+    acknowledged by another with both in the audit log, wall passes that show the wall and nothing else, the page."""
+    import http.cookiejar
+    import urllib.request
+    import urllib.error
+    from datetime import datetime as _dt, timedelta as _td, timezone as _tz
+    from . import control_room as CR
+    from . import edr as E
+    from . import helplib as H
+    from .workspace import Workspace
+    from .service import Service, Users
+    pkg = Path(__file__).parent
+    d = Path(tmp, "cd"); d.mkdir()
+
+    ev = E.selftest_events(str(d))
+    wsp = str(d / "ws"); ws = Workspace.create(wsp, "Pad 3", actor="tester")
+    wd = ws.add_well_file(ev["path"], display="Drilling well", actor="tester")
+    wa = ws.add_well_file(ev["path"], display="Alarm well", actor="tester")
+    wc = ws.add_well_catalog("volve_f12_f14_production_excerpt", "15/9-F-12", 10000.0, actor="tester")
+    ws.refresh_dashboard(actor="tester")
+    wn = ws.add_well_file(ev["path"], display="New well", actor="tester")          # added after the refresh: no reports yet
+    ra = Path(ws.reports_dir, "wells", wa["id"])
+    (ra / "alarm_event_report.json").write_text(json.dumps({"active": [
+        {"alarm_id": "SPP.HH_RANGE", "state": "ACTIVE_UNACKED", "priority": "P2", "tag_id": "SPP"},
+        {"alarm_id": "HKLD.QUALITY", "state": "ACTIVE_ACKED", "priority": "P4", "tag_id": "HKLD"}], "shelved": []}), encoding="utf-8")
+    now = _dt(2026, 10, 15, 12, 0, 0, tzinfo=_tz.utc)
+    iso = lambda s: (now - _td(seconds=s)).strftime("%Y-%m-%dT%H:%M:%SZ")
+    feed = lambda age, status="CONNECTED", wid=wn["id"]: [{"name": "p1", "well_id": wid, "status": status, "last_sample_utc": (iso(age) if age is not None else None), "patch": {"stale_after_s": 120.0}}]
+    b0 = CR.board(ws, [], now=now)
+    t0 = {t["well_id"]: t for t in b0["tiles"]}
+    dr = t0[wd["id"]]["drilling"]
+    live = {k: {t["well_id"]: t for t in CR.board(ws, feed(*a), now=now)["tiles"]}[wn["id"]] for k, a in
+            {"fresh": (30,), "late": (400,), "silent": (1000,), "down": (60, "DOWN"), "none": (None,)}.items()}
+    c1 = {
+        "states": t0[wa["id"]]["state"] == "SERIOUS" and t0[wa["id"]]["reasons"][:2] == ["1 unacknowledged P2 alarm(s)", "1 other active alarm(s)"]
+                  and t0[wc["id"]]["state"] == "SERIOUS" and "no explanation" in t0[wc["id"]]["reasons"][0]
+                  and t0[wd["id"]]["state"] == "GOOD" and t0[wd["id"]]["reasons"] == ["nothing needs attention"]
+                  and t0[wn["id"]]["state"] == "NO_DATA",
+        "shape_word": all(t["shape"] == CR.SHAPE[t["state"]] and t["word"] == CR.WORD[t["state"]] for t in b0["tiles"])
+                      and len(set(CR.SHAPE.values())) == 4 and len(set(CR.WORD.values())) == 4,
+        "order": [t["state"] for t in b0["tiles"]] == sorted((t["state"] for t in b0["tiles"]), key=CR.STATES.index)
+                 and b0["counts"] == {"SERIOUS": 2, "WATCH": 0, "GOOD": 1, "NO_DATA": 1},
+        "recorded": t0[wd["id"]]["live"] is None and t0[wd["id"]]["recorded"]["record_end_utc"] == dr["record_end_utc"]
+                    and "never called stale" in b0["rule"],
+        "read_only": not Path(wsp, "handover").exists() and not Path(wsp, "wall").exists(),
+        "drilling": dr["cache"] and dr["td_md_ft"] == 5750.0 and dr["bottom_ft"] == 5750.0 and dr["on_bottom_hours"] is not None
+                    and (dr["rop_ft_h"] is not None) == bool(dr["drilling_now"]),
+        "live": live["fresh"]["state"] == "GOOD" and live["late"]["state"] == "WATCH" and "past the 2 min stale limit" in live["late"]["reasons"][0]
+                and live["silent"]["state"] == "SERIOUS" and live["silent"]["reasons"][0].startswith("no live sample for 16 min")
+                and live["down"]["state"] == "SERIOUS" and live["none"]["state"] == "WATCH" and abs(live["fresh"]["live"]["age_s"] - 30) < 1e-6,
+    }
+    failed1 = [k for k, v in c1.items() if not v]
+    ok(not failed1,
+       "CD1 the wall board by its stated rule, read without writing anything: an unacknowledged P2 alarm and a gauge off with no explanation SERIOUS, a well with reports and "
+       "nothing else GOOD, a well with neither reports nor a live sample NO DATA, each state its own shape and word, worst first; a recorded well "
+       "described at the end of its record and never called stale (bottom 5,750 ft, the ROP only while drilling at the end); a live feed's "
+       "measured age against its stale limit - 30 s good, 400 s watch, 1,000 s serious, a patch down serious"
+       + (f" [failed: {', '.join(failed1)}]" if failed1 else ""))
+
+    users = Users(str(Path(wsp, "users.json")))
+    for n_, r_ in (("adm", "admin"), ("op1", "operator"), ("op2", "operator"), ("v", "viewer")):
+        users.add(n_, f"{n_}-pass-123", r_)
+    svc = Service(wsp, port=0, scheduler=False).start()
+    base_url = svc.url.rstrip("/")
+
+    def opener():
+        return urllib.request.build_opener(urllib.request.HTTPCookieProcessor(http.cookiejar.CookieJar()))
+
+    def call(op, method, path, body=None):
+        req = urllib.request.Request(base_url + path, method=method, data=(json.dumps(body).encode() if body is not None else None))
+        if method == "POST":
+            req.add_header("Content-Type", "application/json"); req.add_header("X-GEA-Action", "1")
+        try:
+            with op.open(req, timeout=120) as rr_:
+                return rr_.status, json.loads(rr_.read())
+        except urllib.error.HTTPError as e:
+            return e.code, json.loads(e.read())
+
+    def signin(name):
+        op = opener(); call(op, "POST", "/api/login", {"name": name, "password": f"{name}-pass-123"}); return op
+    try:
+        anon = opener()
+        st_anon, _ = call(anon, "GET", "/api/wall")
+        V, O1, O2, A = signin("v"), signin("op1"), signin("op2"), signin("adm")
+        st_v, bv = call(V, "GET", "/api/wall")
+        st_vw, _ = call(V, "POST", "/api/handover", {"open": "x"})
+        st_empty, _ = call(O1, "POST", "/api/handover", {"open": " ", "watching": "", "done": ""})
+        st_long, _ = call(O1, "POST", "/api/handover", {"open": "x" * 2001})
+        st_w, note = call(O1, "POST", "/api/handover", {"open": "Pump 2 liner swap not finished.", "watching": "SPP on the alarm well.", "done": "Reamed to bottom."})
+        _, bdg1 = call(V, "GET", "/api/badges")
+        _, wall1 = call(V, "GET", "/api/wall")
+        st_self, self_ = call(O1, "POST", "/api/handover/ack", {"id": note.get("id")})
+        st_ack, ack = call(O2, "POST", "/api/handover/ack", {"id": note.get("id"), "comment": "read at 06:05"})
+        st_ack2, _ = call(O2, "POST", "/api/handover/ack", {"id": note.get("id")})
+        st_nf, _ = call(O2, "POST", "/api/handover/ack", {"id": "H-none"})
+        _, bdg2 = call(V, "GET", "/api/badges")
+        _, hv = call(V, "GET", "/api/handover")
+        # wall passes
+        st_pv, _ = call(V, "POST", "/api/wall/pass", {"label": "doghouse"})
+        st_pnl, _ = call(A, "POST", "/api/wall/pass", {"label": " "})
+        st_pc, wp = call(A, "POST", "/api/wall/pass", {"label": "doghouse wall", "days": 30})
+        _, plist = call(A, "GET", "/api/wall/passes")
+        W = opener()
+        st_bad, _ = call(W, "POST", "/api/wall/enter", {"code": wp["id"] + ".wrong"})
+        st_in, _ = call(W, "POST", "/api/wall/enter", {"code": wp["code"]})
+        st_ww, wwall = call(W, "GET", "/api/wall")
+        _, wsess = call(W, "GET", "/api/session")
+        st_wo, _ = call(W, "GET", "/api/overview")
+        st_wh, _ = call(W, "GET", "/api/handover")
+        st_wpost, _ = call(W, "POST", "/api/handover", {"open": "x"})
+        st_rv, _ = call(A, "POST", "/api/wall/pass/revoke", {"id": wp["id"]})
+        st_after, _ = call(W, "GET", "/api/wall")
+        _, plist2 = call(A, "GET", "/api/wall/passes")
+        L = opener()
+        lock = [call(L, "POST", "/api/wall/enter", {"code": "W000000.nope"})[0] for _ in range(6)]
+    finally:
+        svc.stop()
+    audit = [e["action"] for e in ws.audit_log(5000)]
+    raw = (Path(wsp, "wall", "passes.json")).read_text(encoding="utf-8")
+    lines = Path(wsp, "handover", "notes.jsonl").read_text(encoding="utf-8").splitlines()
+    c2 = {
+        "roles": st_anon == 401 and st_v == 200 and len(bv["tiles"]) == 4 and st_vw == 403,
+        "note": st_empty == 400 and st_long == 400 and st_w == 200 and note["by"] == "op1" and note["ack"] is None
+                and note["board"]["SERIOUS"] == 2 and note["board"]["alarms_unacked"] == 1
+                and wall1["handover"]["id"] == note["id"] and wall1["handover"]["ack"] is None and bdg1["handover_unacked"] == 1,
+        "ack": st_self == 403 and "author cannot acknowledge" in self_["error"] and st_ack == 200 and ack["by"] == "op2"
+               and st_ack2 == 400 and st_nf == 404 and bdg2["handover_unacked"] == 0
+               and hv["notes"][0]["ack"]["by"] == "op2" and hv["notes"][0]["ack"]["comment"] == "read at 06:05",
+        "append_only": len(lines) == 2 and json.loads(lines[0])["kind"] == "note" and json.loads(lines[1])["kind"] == "ack",
+        "audit": all(a_ in audit for a_ in ("handover.write", "handover.ack", "wall.pass.create", "wall.pass.failed", "wall.pass.enter", "wall.pass.revoke")),
+    }
+    failed2 = [k for k, v in c2.items() if not v]
+    ok(not failed2,
+       "CD2 the shift handover: a viewer reads the wall and cannot write; an empty or over-long note refused; a signed note keeps its author, its "
+       "time and the board's counts (2 serious, 1 unacknowledged alarm), shows on the wall and badges the tab until acknowledged; the author "
+       "cannot acknowledge it, the next operator can (once, with a comment); notes and acknowledgements appended, never rewritten; every step in "
+       "the audit log"
+       + (f" [failed: {', '.join(failed2)}]" if failed2 else ""))
+
+    page_src = (pkg / "web" / "app.html").read_text(encoding="utf-8")
+    help_md = (pkg / "help" / "control-room.md").read_text(encoding="utf-8")
+    c3 = {
+        "pass": st_pv == 403 and st_pnl == 400 and st_pc == 200 and wp["code"].startswith(wp["id"] + ".") and wp["code"].split(".", 1)[1] not in raw
+                and "hash" not in json.dumps(plist) and plist["passes"][0]["status"] == "ACTIVE",
+        "wall_only": st_bad == 401 and st_in == 200 and st_ww == 200 and len(wwall["tiles"]) == 4 and wsess["user"] is None
+                     and wsess["wall_pass"]["label"] == "doghouse wall" and st_wo == 401 and st_wh == 401 and st_wpost == 401,
+        "revoke": st_rv == 200 and st_after == 401 and plist2["passes"][0]["status"] == "REVOKED" and lock[-1] == 429,
+        "page": all(x in page_src for x in ("VIEWS.wall", "VIEWS.control", "['#/control', 'Control room']", "/api/wall", "/api/handover/ack",
+                                             "NOT UPDATING", "wakeLock", "requestFullscreen", "SIGNED OUT", "This is a wall screen",
+                                             "crShape", "data-age0", "turns pages", "/api/wall/pass/revoke", "handover_unacked")),
+        "help": any(t["topic"] == "control-room" for t in H.topics()) and not H.check("control-room") and H.VIEW_TOPIC.get("control") == "control-room"
+                and "never called stale" in help_md.split("will not call a measurement")[1],
+    }
+    failed3 = [k for k, v in c3.items() if not v]
+    ok(not failed3,
+       "CD3 wall passes: issued by an administrator only and named for the screen, the code shown once and only its hash kept; a wrong pass "
+       "refused and locked out after repeats; a pass shows the wall and nothing else - no other page, no handover list, no writing - until it is "
+       "revoked; the page: the wall full screen with its shapes and words, the age of its last answer, NOT UPDATING across the screen when the "
+       "service stops, the screen kept awake, pages turned when the wells do not fit, the handover written and acknowledged; the help page"
+       + (f" [failed: {', '.join(failed3)}]" if failed3 else ""))
+
+
+
+def section_ce_well3d(tmp: str) -> None:
+    """Section CE - the well in 3D: the minimum-curvature method against a known arc and against the survey companies'
+    own positions, the shipped Volve excerpt checked against itself (frame, picks, surfaces, provenance), a workspace
+    well drawn from an attached survey among its formation-top contours with its recorder replay, the API and the tab."""
+    import base64 as _b64
+    import gzip as _gz
+    import hashlib as _hl
+    import http.cookiejar
+    import math as _m
+    import subprocess as _sp
+    import urllib.request
+    import urllib.error
+    from . import edr as E
+    from . import geodesy as G
+    from . import helplib as H
+    from . import shapefile as SF
+    from . import well3d as W
+    from .workspace import Workspace
+    from .service import Service, Users
+    pkg = Path(__file__).parent
+    d = Path(tmp, "ce"); d.mkdir()
+
+    R = 1000.0
+    q = W.min_curvature([0.0, _m.pi * R / 2], [0.0, 90.0], [0.0, 0.0])
+    p_arc = W.Path([0.0, _m.pi * R / 2], [0.0, 90.0], [0.0, 0.0], 0.0, 0.0, 0.0)
+    mid = p_arc.at(_m.pi * R / 4)
+    tr4 = W.read_witsml_trajectory(_gz.decompress((pkg / "reference" / "volve_traj_F-4.xml.gz").read_bytes()).decode())[0]
+    tr10 = W.read_witsml_trajectory(_gz.decompress((pkg / "reference" / "volve_traj_F-10.xml.gz").read_bytes()).decode())[0]
+    f5 = [t for n in ("36in", "17.5in", "12.25in", "8.5in") for t in W.read_witsml_trajectory(_gz.decompress((pkg / "reference" / f"volve_traj_F-5_{n}.xml.gz").read_bytes()).decode())]
+    rep5, _ = W.replay_from_stations(f5)
+    rep4, note4 = W.replay_from_stations([tr4])
+    csv_ok = "MD (ft),Inc (deg),Azi (deg)\n0,0,0\n1000,0,0\n2000,10,90\n3000,30,90\n4000,45,95\n5750,60,100\n"
+    refused = []
+    for bad in ("Depth,Angle\n0,0\n", "MD,Inc,Azi\n0,0,0\n0,1,2\n", "MD,Inc,Azi\n0,0,0\n100,x,2\n", "MD,Inc,Azi\n0,0,0\n100,200,2\n"):
+        try:
+            W.read_survey_csv(bad)
+        except ValueError:
+            refused.append(True)
+    st_csv = W.read_survey_csv(csv_ok, "ft")
+    c1 = {
+        "arc": abs(q["n"][1] - R) < 1e-9 and abs(q["v"][1] - R) < 1e-9 and abs(mid[1] - R * (1 - _m.cos(_m.pi / 4))) < 1e-6
+               and abs(mid[2] - R * _m.sin(_m.pi / 4)) < 1e-6 and abs(q["dls"][1] - 90.0 / (_m.pi * R / 2) * 30.0) < 1e-9,
+        "units": tr4["azi_ref"] == "true north" and abs(tr4["stations"][40]["incl"] - _m.degrees(0.446106147512686)) < 1e-9 and tr10["grid_cor_deg"] == -0.95
+                 and len(tr4["stations"]) == 87 and abs(st_csv[-1]["md"] - 5750 * 0.3048) < 1e-9,
+        "company": all(x is not None and x < 0.01 for x in (W.file_position_miss(tr4["stations"]), W.file_position_miss(tr10["stations"]),
+                                                            W.file_position_miss(f5[-1]["stations"]))),
+        "csv": len(refused) == 4,
+        "replay": rep5 is not None and len(rep5["t"]) == 100 and rep5["t"][0].startswith("2007-11-16") and rep5["t"][-1].startswith("2008-07-29")
+                  and rep4 is None and "1 distinct time" in note4 and all(m0 <= m1 or t0 < t1 for (t0, m0), (t1, m1) in zip(zip(rep5["t"], rep5["md"]), zip(rep5["t"][1:], rep5["md"][1:]))),
+    }
+    failed1 = [k for k, v in c1.items() if not v]
+    ok(not failed1,
+       "CE1 the minimum-curvature method: a quarter circle of 1,000 m radius ends 1,000 m out and 1,000 m down, and its middle lies on the arc "
+       "(292.9 m out, 707.1 m down), not on a chord; the WITSML reader takes radians and degrees alike; the paths computed from the Volve files "
+       "reproduce the TVD and displacements the survey companies wrote in them to under a centimetre; a survey table without its columns, with a "
+       "repeated depth or a bad row refused; the replay only from stations measured while drilling (F-5: 100 steps, Nov 2007 to Jul 2008; F-4's "
+       "one-time resurvey not replayed)"
+       + (f" [failed: {', '.join(failed1)}]" if failed1 else ""))
+
+    sc = W.volve_scene()
+    wl = {w["id"]: w for w in sc["wells"]}
+    picks = W.volve_picks()
+    dms = lambda dd, mm, ss: dd + mm / 60 + ss / 3600
+    frame = []
+    for wn, la, lo in (("NO 15/9-F-1 C", dms(58, 26, 29.907), dms(1, 53, 14.708)), ("NO 15/9-F-11 B", dms(58, 26, 29.960), dms(1, 53, 14.870))):
+        sb = next(p for p in picks if p["WELL"] == wn and p["PICKS"] == "Seabed")
+        a = G.unproject(float(sb["EASTING"]), float(sb["NORTHING"]), G.utm(31, True, "INTL1924"))
+        b = G.unproject(float(sb["EASTING"]), float(sb["NORTHING"]), G.utm(31, True, "WGS84"))
+        frame.append((G.geodesic_m(a[0], a[1], la, lo)["distance_m"], G.geodesic_m(b[0], b[1], la, lo)["distance_m"]))
+    prov = json.loads((pkg / "reference" / "volve_3d.provenance.json").read_text(encoding="utf-8"))
+    hashes_ok = all(_hl.sha256((pkg / "reference" / k).read_bytes()).hexdigest() == v["sha256"] for k, v in prov["outputs"].items())
+    pts = prov["outputs"]["volve_surfaces.json.gz"]["points"]
+    # the repository stores every text file with LF line ends (.gitattributes: text=auto eol=lf), so a package file holding
+    # CRLF would reach a checkout - and the published wheel - as different bytes from the ones hashed here
+    crlf = []
+    for dp, dn, fn in os.walk(pkg):
+        dn[:] = [x for x in dn if x not in ("__pycache__", "catalog_operator")]
+        for f_ in fn:
+            b_ = Path(dp, f_).read_bytes()
+            if b"\0" not in b_[:8000] and b"\r\n" in b_:
+                crlf.append(os.path.relpath(os.path.join(dp, f_), pkg))
+    c2 = {
+        "frame": all(a < 0.1 and 140 < b < 160 for a, b in frame) and "ED50" in sc["frame"]["crs"],
+        "picks": all(wl[k]["pick_miss_max_m"] < 0.3 and abs(wl[k]["ref_elev_m"] - 54.9) < 1e-6 for k in ("F-4", "F-5", "F-10"))
+                 and abs(wl["F-10"]["td_md"] - 5300.07) < 0.01 and len(sc["picks_only"]) == 31,
+        "azimuth": wl["F-4"]["azimuth"]["applied"] == "as grid" and "14 m" in wl["F-4"]["azimuth"]["note"] and wl["F-5"]["azimuth"]["file_reference"] == "grid north",
+        "surfaces": len(sc["surfaces"]) == 5 and all(s["affine_miss_m"] < 0.01 and s["tie_median_m"] < 7.0 and len(s["ties"]) >= 23 for s in sc["surfaces"])
+                    and all(pts[s["id"]] == [s["n_source_points"], s["n_points"]] for s in sc["surfaces"])
+                    and [s["id"] for s in sc["surfaces"]] == ["Ty", "SHETLAND", "BCU", "Hugin_Fm_Top", "Hugin_Base"],
+        "line_ends": not crlf and all(k.endswith(".gz") or b"\r\n" not in (pkg / "reference" / k).read_bytes() for k in prov["outputs"]),
+        "provenance": hashes_ok and "commercial purposes" in prov["license"] and "may not be sold" in prov["license"] and len(prov["sources"]) == 12
+                      and all(len(v["commit"]) == 40 for v in prov["sources"].values()) and "Equinor ASA" in sc["attribution"] and "not sold" in sc["attribution"],
+    }
+    failed2 = [k for k, v in c2.items() if not v]
+    ok(not failed2,
+       "CE2 the Volve excerpt checked against itself: its eastings and northings are ED50 zone 31N (the slots unproject within 0.1 m of the "
+       "operator's own log headers, about 150 m off if read as WGS84); the three paths land on Equinor's picks to under 0.3 m from a 54.9 m rig "
+       "floor; 15/9-F-4's azimuths read as grid, with the 14 m it would miss by otherwise stated; five surfaces drawn as delivered, every fourth "
+       "line, each within a few metres of its picks at the median; every shipped file hashed in its provenance with the pinned public copy it "
+       "came from (and no package file with CRLF line ends, which a checkout would change from the bytes hashed), the licence (commercial use, no sale) and the attribution on the page"
+       + (f" [failed: {', '.join(failed2)}]" if failed2 else ""))
+
+    ev = E.selftest_events(str(d))
+    wsp = str(d / "ws"); ws = Workspace.create(wsp, "Pad 3", actor="tester")
+    wb = ws.add_well_file(ev["path"], display="Well B", actor="tester"); wc = ws.add_well_file(ev["path"], display="Well C", actor="tester")
+    ws.set_well(wb["id"], surface={"lat": 31.42, "lon": -104.10, "datum": "NAD83", "elevation": {"ft": 3000.0, "datum": "KB", "source": "test plat"}}, actor="tester")
+    lat1, lon1 = 31.42, -104.10
+    kmlon = 1.0 / (111.195 * _m.cos(_m.radians(lat1)))
+    SF.write(str(d / "Top_A"), "PolyLine", [([[(lon1 + (i - 1.5) * kmlon + 0.0005 * k, lat1 - 0.05 + 0.001 * k) for k in range(100)]], {"ElevTop_ft": -1500.0 - 100.0 * i}) for i in range(4)],
+             [("ElevTop_ft", "F", 10, 1)], SF.NAD83_GEOGCS)
+    SF.write(str(d / "Faults"), "PolyLine", [([[(lon1 - 0.02 + 0.001 * j, lat1) for j in range(40)]], {"TYPE": "normal fault"})], [("TYPE", "C", 20, 0)], SF.NAD83_GEOGCS)
+    ws.add_layer(str(d / "Top_A.shp"), "Top A", actor="tester", formation="A", role="top", citation="test")
+    ws.add_layer(str(d / "Faults.shp"), "Faults", actor="tester", citation="test")
+    sc0 = W.workspace_scene(ws)                                   # no survey yet: the formation tops alone
+    att_refused = 0
+    for kw in ({"azimuth_reference": ""}, {"azimuth_reference": "grid"}, {"azimuth_reference": "magnetic"}, {"azimuth_reference": "true", "source": " "}):
+        k = {"azimuth_reference": "true", "source": "s", **kw}
+        try:
+            W.attach_trajectory(ws, wb["id"], "s.csv", csv_ok.encode(), "tester", k["azimuth_reference"], "ft", None, None, None, k["source"])
+        except ValueError:
+            att_refused += 1
+    ws.refresh_dashboard(actor="tester")
+    Users(str(Path(wsp, "users.json"))).add("v", "viewer-pass-1", "viewer"); Users(str(Path(wsp, "users.json"))).add("op", "op-pass-123", "operator")
+    svc = Service(wsp, port=0, scheduler=False).start()
+    base_url = svc.url.rstrip("/")
+
+    def opener(name, pw):
+        o = urllib.request.build_opener(urllib.request.HTTPCookieProcessor(http.cookiejar.CookieJar()))
+        rq = urllib.request.Request(base_url + "/api/login", method="POST", data=json.dumps({"name": name, "password": pw}).encode())
+        rq.add_header("Content-Type", "application/json"); rq.add_header("X-GEA-Action", "1"); o.open(rq, timeout=30)
+        return o
+
+    def call(o, method, path, body=None):
+        rq = urllib.request.Request(base_url + path, method=method, data=(json.dumps(body).encode() if body is not None else None))
+        if method == "POST":
+            rq.add_header("Content-Type", "application/json"); rq.add_header("X-GEA-Action", "1")
+        try:
+            with o.open(rq, timeout=120) as rr_:
+                return rr_.status, json.loads(rr_.read())
+        except urllib.error.HTTPError as e:
+            return e.code, json.loads(e.read())
+    try:
+        V, O = opener("v", "viewer-pass-1"), opener("op", "op-pass-123")
+        st_l, ls = call(V, "GET", "/api/3d/scenes")
+        st_v, vs = call(V, "GET", "/api/3d/scene/volve")
+        st_n, _ = call(V, "GET", "/api/3d/scene/nope")
+        body = {"filename": "survey.csv", "content_b64": _b64.b64encode(csv_ok.encode()).decode(), "azimuth_reference": "grid", "grid": "TX_CENTRAL",
+                "md_unit": "ft", "source": "test survey"}
+        st_pv, _ = call(V, "POST", f"/api/wells/{wb['id']}/trajectory", body)
+        st_b64, _ = call(O, "POST", f"/api/wells/{wb['id']}/trajectory", {**body, "content_b64": "%%%"})
+        st_at, at = call(O, "POST", f"/api/wells/{wb['id']}/trajectory", body)
+        st_w, wsc = call(V, "GET", "/api/3d/scene/workspace")
+    finally:
+        svc.stop()
+    kept = (Path(wsp, "wells", wb["id"], "trajectory", "survey.csv").read_bytes() == csv_ok.encode())
+    audit = [e["action"] for e in ws.audit_log(500)]
+    cli = _sp.run([sys.executable, "-m", "gea", "well3d", "--scene", "workspace", "--workspace", wsp], capture_output=True, text=True, cwd=str(pkg.parent))
+    page_src = (pkg / "web" / "app.html").read_text(encoding="utf-8")
+    help_md = (pkg / "help" / "well3d.md").read_text(encoding="utf-8")
+    w0 = (wsc.get("wells") or [{}])[0]
+    c3 = {
+        "attach": att_refused == 4 and st_pv == 403 and st_b64 == 400 and st_at == 200 and at["stations"] == 6 and at["grid"] == "TX_CENTRAL" and kept
+                  and at["sha256"] == _hl.sha256(csv_ok.encode()).hexdigest() and "well.trajectory" in audit,
+        "scene": st_l == 200 and [s["id"] for s in ls["scenes"]] == ["workspace", "volve"] and st_v == 200 and len(vs["wells"]) == 3 and st_n == 404
+                 and st_w == 200 and len(wsc["wells"]) == 1 and abs(w0["td_md"] - 1752.6) < 0.01 and "TX_CENTRAL grid north" in w0["azimuth"]["applied"]
+                 and w0["replay"]["kind"] == "drilling recorder" and "Gamma" in w0["channels"] and not w0["replay"]["t"][0].endswith("Z")
+                 and [n["id"] for n in wsc["not_placed"]] == [wc["id"]] and wsc["layers"][0]["lines"][0]["z"] == round(1500 * 0.3048, 1)
+                 and any(n["display"] == "Faults" and "no depth" in n["why"] for n in wsc["layers_not_drawn"]),
+        "tops_alone": sc0["wells"] == [] and [l["name"] for l in sc0["layers"]] == ["Top A"] and sc0["frame"]["area"] == "the placed wells"
+                      and {n["id"]: n["why"] for n in sc0["not_placed"]}[wb["id"]] == "no directional survey attached",
+        "cli": cli.returncode == 0 and "Well B: 6 stations" in cli.stdout and "contour lines at their elevations" in cli.stdout,
+        "page": all(x in page_src for x in ("VIEWS['3d']", "['#/3d', '3D']", "getContext('webgl'", "OES_element_index_uint", "w3Tube", "w3SurfMesh",
+                                             "/api/3d/scene/", "w3Play", "the recorder's clock", "ArrowLeft", "see-through", "/trajectory`", "w3Hit")),
+        "help": any(t["topic"] == "well3d" for t in H.topics()) and not H.check("well3d") and H.VIEW_TOPIC.get("3d") == "well3d"
+                and "a well with no survey" in help_md.split("will not call a measurement")[1],
+    }
+    failed3 = [k for k, v in c3.items() if not v]
+    ok(not failed3,
+       "CE3 a workspace well in 3D: a survey attached only with its azimuth reference stated (a grid named, a declination given) and its source, "
+       "kept byte for byte and audited; the path turned from TX_CENTRAL grid north by the convergence there, set below sea level from the "
+       "declared rig floor among the formation-top contours at their elevations (a fault trace listed as having no depth; the tops drawn alone "
+       "before any survey), the recorder's bit depth replayed on its own clock and its gamma along the path; a well without a position listed, not drawn; the API, the CLI and the 3D "
+       "tab (WebGL by hand, the replay, turn and zoom by pointer and keys); the help page"
+       + (f" [failed: {', '.join(failed3)}]" if failed3 else ""))
+
+
 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:
@@ -8532,6 +8884,8 @@ def main() -> int:
         section_ca_setup_checklist(tmp)
         section_cb_charts(tmp)
         section_cc_site_map(tmp)
+        section_cd_control_room(tmp)
+        section_ce_well3d(tmp)
     if _FAILS:
         print(f"[ACCEPTANCE] {len(_FAILS)} FAILURES ({_PASS} passed):")
         for f in _FAILS: