{
  "$schema": "https://ui.shadcn.com/schema/registry-item.json",
  "name": "subsea-manifold",
  "type": "registry:component",
  "title": "Subsea production manifold",
  "description": "Cluster-type subsea production manifold: two headers (production and test) with weld-capped ends rising to vertical outlet hubs, well slots on both long sides, each with a vertical jumper hub in a soft-landing funnel, a branch across both headers and an ROV-operated trunnion ball valve on each drop (gearbox and class-4 bucket). Header transmitters, saddles, overtrawlable tubular protection structure with sloped members, roof grating over the valves and open strips over the hubs, lift points, grab bars, anodes and an ROV-readable nameplate. Sits on a skirted mudmat; y = 0 is the seabed. Hub ports carry no flange rating (clamp hubs).",
  "dependencies": [
    "three",
    "@react-three/fiber",
    "@react-three/drei"
  ],
  "registryDependencies": [
    "https://deck-twin-library.demoareeb.workers.dev/r/core.json"
  ],
  "files": [
    {
      "path": "drilling/subsea-manifold.tsx",
      "type": "registry:file",
      "target": "components/deck-twin/drilling/subsea-manifold.tsx",
      "content": "// Subsea production manifold (API 17P / 17D): dual headers, valved well-slot branches with ROV buckets, vertical jumper hubs, overtrawlable protection structure, mudmat.\nimport * as THREE from 'three'\nimport { Kit, type MatKey, type V3 } from '../lib/kit'\nimport { paint } from '../lib/materials'\nimport { KitMesh, useKit, type GroupProps } from '../lib/mesh'\nimport { grating, member, pipeRun, tagPlate, weldCap } from '../lib/parts'\nimport { port, type Port } from '../lib/ports'\nimport { faceToFace, pipeID, pipeOD, pipeWall, type Rating } from '../lib/standards'\nimport { anode, grabBar, hub, mudmat, subseaBall } from './_drill'\n\nexport interface SubseaManifoldProps {\n  /** well slots (split over both long sides) */\n  slots?: number\n  /** slot branch size, inches */\n  slotNps?: number\n  /** header size, inches */\n  headerNps?: number\n  rating?: Rating\n  tag?: string\n}\n\nexport const defaults = { slots: 4, slotNps: 6, headerNps: 10, rating: 1500 as Rating, tag: 'PM-01' }\n\ntype P = typeof defaults\n\nconst HUB_LEN = 0.45\n\nfunction layout(p: Omit<P, 'tag'>) {\n  const n = Math.max(2, Math.min(8, Math.round(p.slots)))\n  const perSide = Math.ceil(n / 2)\n  const pitch = 2.2\n  const Lm = Math.max(6.5, perSide * pitch + 3.2)\n  const zH = 2.0\n  const Wm = 2 * (zH + 1.5)\n  const yHd = 1.0\n  const zHd = 0.55\n  const Lv = faceToFace('ball-full', p.slotNps, p.rating)\n  const yB = yHd + pipeOD(p.headerNps) / 2 + Lv + 0.45\n  const yHub = yB + 0.75\n  const xOut = Lm / 2 - 1.1\n  const x0 = -Lm / 2 + 0.9\n  const slotsAt = Array.from({ length: n }, (_, i) => {\n    const side = i % 2 === 0 ? -1 : 1\n    const j = Math.floor(i / 2)\n    const x = x0 + 0.7 + j * pitch + (side > 0 ? 0.3 : 0)\n    return { x, side, face: [x, yHub + HUB_LEN, side * zH] as V3 }\n  })\n  const outlets = [\n    { id: 'header-a', z: -zHd, face: [xOut, yHd + 0.9 + HUB_LEN, -zHd] as V3 },\n    { id: 'header-b', z: zHd, face: [xOut, yHd + 0.9 + HUB_LEN, zHd] as V3 },\n  ]\n  return { n, Lm, Wm, yHd, zHd, Lv, yB, zH, yHub, xOut, x0, slotsAt, outlets, yRoof: yHub + HUB_LEN + 0.55 }\n}\n\nexport function ports(props: SubseaManifoldProps = {}): Port[] {\n  const p = { ...defaults, ...props }\n  const l = layout(p)\n  return [\n    ...l.slotsAt.map((s, i) => port(`slot-${i + 1}`, `Well slot ${i + 1} hub`, s.face, [0, 1, 0], p.slotNps)),\n    ...l.outlets.map((o) => port(o.id, o.id === 'header-a' ? 'Production header hub' : 'Test / second header hub', o.face, [0, 1, 0], p.headerNps)),\n  ]\n}\n\nexport function SubseaManifold({ slots = defaults.slots, slotNps = defaults.slotNps, headerNps = defaults.headerNps, rating = defaults.rating, tag = defaults.tag, ...group }: SubseaManifoldProps & GroupProps) {\n  const data = useKit((k) => build(k, { slots, slotNps, headerNps, rating, tag }), [slots, slotNps, headerNps, rating, tag])\n  return (\n    <group {...group}>\n      <KitMesh data={data} />\n    </group>\n  )\n}\n\nexport default SubseaManifold\n\nfunction build(k: Kit, p: P) {\n  const l = layout(p)\n  const { Lm, Wm, yHd, zHd, Lv, yB, zH, yHub, xOut, x0, slotsAt, outlets, yRoof } = l\n  const Y: MatKey = paint('#e3b512', { weather: 1 })\n  const PIPE: MatKey = paint('#d9dcd8', { weather: 0.8 })\n  const t = mudmat(k, Lm, Wm, { mat: Y, skirt: 0.4 })\n  const hOD = pipeOD(p.headerNps)\n  const sOD = pipeOD(p.slotNps)\n\n  // headers on saddles, blind (weld-capped) at the -X end, rising to vertical outlet hubs at +X\n  for (const o of outlets) {\n    pipeRun(k, [[x0, yHd, o.z], [xOut, yHd, o.z], [xOut, yHd + 0.9, o.z]], p.headerNps, { mat: PIPE, schedule: '80' })\n    weldCap(k, p.headerNps, [x0, yHd, o.z], '-x', PIPE)\n    hub(k, [xOut, yHd + 0.9, o.z], 'y', pipeID(p.headerNps, '80') / 2, HUB_LEN, PIPE)\n    for (let x = x0 + 0.5; x < xOut - 0.4; x += 1.6) {\n      k.box(Y, [0.2, yHd - hOD / 2 - t, 0.3], [x, (yHd - hOD / 2 + t) / 2, o.z])\n      k.lathe(Y, [[hOD / 2 + 0.01, -0.1], [hOD / 2 + 0.03, -0.1], [hOD / 2 + 0.03, 0.1], [hOD / 2 + 0.01, 0.1]], [x, yHd, o.z], 'x', 32, Math.PI / 2, Math.PI)\n    }\n    // header pressure / temperature transmitter\n    k.cyl(PIPE, 0.05, 0.16, [xOut - 0.8, yHd + hOD / 2 + 0.08, o.z], 'y', 16)\n    k.rbox('stainless', [0.14, 0.2, 0.14], 0.02, [xOut - 0.8, yHd + hOD / 2 + 0.26, o.z])\n  }\n\n  // slot branches: hub → drop → run across to both header drops; each drop has an ROV-operated ball valve\n  const R = Math.max(0.12, (pipeID(p.slotNps, '80') / 2) * 2.4)\n  for (const s of slotsAt) {\n    const zs = s.side * zH\n    pipeRun(k, [[s.x, yHub, zs], [s.x, yB, zs], [s.x, yB, -s.side * zHd], [s.x, yB - 0.2, -s.side * zHd]], p.slotNps, { mat: PIPE, schedule: '80' })\n    // tee for the near-header drop\n    k.cyl(PIPE, sOD * 0.62, sOD * 1.5, [s.x, yB, s.side * zHd], 'z', 24)\n    k.cyl(PIPE, sOD * 0.6, 0.2, [s.x, yB - 0.1, s.side * zHd], 'y', 24)\n    for (const z of [zHd, -zHd]) {\n      const yv = yB - 0.2\n      k.with(new THREE.Matrix4().makeBasis(new THREE.Vector3(0, 0, 1), new THREE.Vector3(0, -1, 0), new THREE.Vector3(1, 0, 0)).setPosition(s.x, yv, z), () => subseaBall(k, p.slotNps, p.rating, Lv, { mat: PIPE }))\n      k.pipe(PIPE, [s.x, yv - Lv + 0.002, z], [s.x, yHd + hOD * 0.3, z], sOD, pipeWall(p.slotNps, '80'))\n    }\n    hub(k, [s.x, yHub, zs], 'y', pipeID(p.slotNps, '80') / 2, HUB_LEN, PIPE)\n    // hub support, guide posts and soft-landing funnel for the jumper connector\n    k.box(Y, [0.3, yB - sOD / 2 - t, 0.3], [s.x, (yB - sOD / 2 + t) / 2, zs])\n    const Rf = R + 0.2\n    const yf = yHub + HUB_LEN - 0.12\n    for (let i = 0; i < 4; i++) {\n      const a = (i / 4) * Math.PI * 2 + Math.PI / 4\n      member(k, [s.x + Math.cos(a) * Rf, t, zs + Math.sin(a) * Rf], [s.x + Math.cos(a) * Rf, yf, zs + Math.sin(a) * Rf], 'CHS 60.3x5', Y)\n    }\n    k.lathe(Y, [[Rf - 0.02, yf], [Rf, yf], [Rf + 0.3, yf + 0.25], [Rf + 0.28, yf + 0.25], [Rf - 0.02, yf]], [s.x, 0, zs], 'y', 36)\n  }\n\n  // protection structure: corner and intermediate columns, roof beams, sloped overtrawlable members, roof grating over the valves\n  const fx = Lm / 2 - 0.15\n  const fz = Wm / 2 - 0.15\n  const rx = fx - 0.6\n  const rz = zH + 0.75\n  const cols: number[] = []\n  for (let x = -rx; x <= rx + 1e-6; x += (2 * rx) / Math.max(2, Math.round((2 * rx) / 2.4))) cols.push(x)\n  for (const x of cols)\n    for (const sz of [-1, 1]) {\n      member(k, [x, t, sz * rz], [x, yRoof, sz * rz], 'CHS 219.1x10', Y)\n      member(k, [x * (fx / rx), t + 0.05, sz * fz], [x, yRoof, sz * rz], 'CHS 168.3x8', Y)\n    }\n  for (const sx of [-1, 1]) {\n    member(k, [sx * fx, t + 0.05, 0], [sx * rx, yRoof, 0], 'CHS 168.3x8', Y)\n    member(k, [sx * rx, yRoof, -rz], [sx * rx, yRoof, rz], 'CHS 219.1x10', Y)\n  }\n  for (const sz of [-1, 1]) member(k, [-rx, yRoof, sz * rz], [rx, yRoof, sz * rz], 'CHS 219.1x10', Y)\n  for (const x of cols.slice(1, -1)) member(k, [x, yRoof, -rz], [x, yRoof, rz], 'HEA 200', Y)\n  // roof grating over the header / valve zone, open strips over the slot hubs for jumper landing\n  const zIn = zH - 0.75\n  for (let i = 0; i < cols.length - 2; i++) {\n    const xa = cols[i]\n    const xb = cols[i + 1]\n    grating(k, [(xa + xb) / 2, yRoof + 0.13, 0], xb - xa - 0.12, 2 * zIn)\n    for (const sz of [-1, 1]) member(k, [xa, yRoof + 0.06, sz * zIn], [xb, yRoof + 0.06, sz * zIn], 'HEA 120', Y)\n  }\n  for (const sz of [-1, 1]) grabBar(k, [-rx * 0.6, yRoof - 0.45, sz * rz], [rx * 0.6, yRoof - 0.45, sz * rz], [0, 0, sz])\n  // lift points on the roof corners\n  for (const sx of [-1, 1])\n    for (const sz of [-1, 1]) {\n      k.box(Y, [0.4, 0.28, 0.05], [sx * rx, yRoof + 0.24, sz * rz])\n      k.torus('nut', 0.07, 0.016, [sx * rx, yRoof + 0.3, sz * rz], 'z', 6, 16)\n    }\n  for (const x of cols) for (const sz of [-1, 1]) anode(k, [x, t + 0.7, sz * (rz + 0.17)], 'y', 0.6)\n  k.box(Y, [1.4, 0.45, 0.03], [0, yRoof - 0.9, rz + 0.13])\n  tagPlate(k, p.tag, [0, yRoof - 0.9, rz + 0.147], 'z', 1.2)\n}\n"
    },
    {
      "path": "drilling/_drill.ts",
      "type": "registry:file",
      "target": "components/deck-twin/drilling/_drill.ts",
      "content": "// Shared drilling & subsea parts: API 6A flanges and gate valves, hammer-lug covers, convex hulls, ROV buckets, hubs and clamps, mudmats.\nimport * as THREE from 'three'\nimport { mat4, type Kit, type MatKey, type V3 } from '../lib/kit'\nimport { flange, liftingLug } from '../lib/parts'\nimport { flange as asmeFlange, inch, nearest, pipeOD, type Rating } from '../lib/standards'\n\n// ------------------------------------------------------------------ API 6A (typical 6B / 6BX dims, inches: OD, thickness, bolt circle, bolts, bolt dia)\n\nexport type Psi = 5000 | 10000 | 15000\n\ntype Row = [od: number, t: number, bc: number, n: number, bd: number]\n\nconst API: Record<Psi, Record<number, Row>> = {\n  5000: {\n    2.0625: [8.5, 1.81, 6.5, 8, 0.875], 2.5625: [9.62, 1.94, 7.5, 8, 1.0], 3.125: [10.5, 2.19, 8.0, 8, 1.125],\n    4.0625: [12.25, 2.44, 9.5, 8, 1.25], 5.125: [14.75, 3.19, 11.5, 8, 1.5], 7.0625: [15.5, 3.62, 12.5, 12, 1.375],\n  },\n  10000: {\n    2.0625: [7.88, 1.73, 6.25, 8, 0.75], 2.5625: [9.12, 2.02, 7.25, 8, 0.875], 3.0625: [10.62, 2.3, 8.5, 8, 1.0],\n    4.0625: [12.44, 2.77, 10.19, 8, 1.125], 5.125: [14.06, 3.12, 11.81, 12, 1.125], 7.0625: [18.88, 4.06, 15.88, 12, 1.5],\n  },\n  15000: {\n    2.0625: [8.75, 2.0, 6.88, 8, 0.875], 2.5625: [10.0, 2.25, 7.88, 8, 1.0], 3.0625: [11.31, 2.53, 9.06, 8, 1.125],\n    4.0625: [14.19, 3.09, 11.44, 8, 1.375], 5.125: [16.5, 3.88, 13.63, 12, 1.5], 7.0625: [19.88, 4.69, 16.88, 16, 1.5],\n  },\n}\n\n// typical API 6A slab-gate face-to-face, inches\nconst GATE: Record<Psi, Record<number, number>> = {\n  5000: { 2.0625: 14.62, 2.5625: 16.62, 3.125: 18.62, 4.0625: 21.62, 5.125: 28.62, 7.0625: 29.0 },\n  10000: { 2.0625: 20.5, 2.5625: 22.25, 3.0625: 24.38, 4.0625: 26.38, 5.125: 29.0, 7.0625: 35.0 },\n  15000: { 2.0625: 19.0, 2.5625: 21.0, 3.0625: 23.56, 4.0625: 29.0, 5.125: 35.0, 7.0625: 40.0 },\n}\n\nexport interface ApiFlange {\n  bore: number\n  od: number\n  t: number\n  bolts: number\n  boltDia: number\n  boltCircle: number\n  ring: number\n}\n\nexport function apiFlange(bore: number, psi: Psi = 10000): ApiFlange {\n  const tbl = API[psi] ?? API[10000]\n  const b = nearest(Object.keys(tbl).map(Number), bore)\n  const [od, t, bc, n, bd] = tbl[b]\n  return { bore: b, od: inch(od), t: inch(t), bolts: n, boltDia: inch(bd), boltCircle: inch(bc), ring: inch(bc - bd * 1.9) }\n}\n\nexport function apiGateLength(bore: number, psi: Psi = 10000): number {\n  const tbl = GATE[psi] ?? GATE[10000]\n  return inch(tbl[nearest(Object.keys(tbl).map(Number), bore)])\n}\n\n/** ring-gasket half gap at an RTJ joint */\nexport const G = 0.003\nconst rBore = (d: ApiFlange) => d.bore * 0.0127\n\n/** Hollow tube along local Y. */\nexport function tubeY(k: Kit, mat: MatKey, ro: number, ri: number, y0: number, y1: number, seg = 28) {\n  k.lathe(mat, [[ri, y0], [ro, y0], [ro, y1], [ri, y1]], [0, 0, 0], 'y', seg)\n  k.lathe('pipeBore', [[ri, Math.max(y0, y1)], [ri, Math.min(y0, y1)]], [0, 0, 0], 'y', seg)\n}\n\n/** Integral API flange, ring face at local y = face, body on side s (+1 above, -1 below). */\nexport function apiFlangeAt(k: Kit, d: ApiFlange, face: number, s: 1 | -1, mat: MatKey) {\n  const f0 = face + s * G\n  const f1 = face + s * (G + d.t)\n  const lo = Math.min(f0, f1)\n  const hi = Math.max(f0, f1)\n  const ch = Math.min(0.006, d.t * 0.12)\n  const rb = rBore(d)\n  const R = d.od / 2\n  const prof: [number, number][] = s > 0 ? [[rb, lo], [R, lo], [R, hi - ch], [R - ch, hi], [rb, hi]] : [[rb, lo], [R - ch, lo], [R, lo + ch], [R, hi], [rb, hi]]\n  k.lathe(mat, prof, [0, 0, 0], 'y', d.bolts >= 12 ? 40 : 28)\n}\n\n/** Studs, nuts and ring gasket for two API flanges meeting at local y = yj. */\nexport function apiJoint(k: Kit, d: ApiFlange, yj: number) {\n  const nh = d.boltDia * 0.95\n  const so = d.boltDia * 0.4\n  const len = 2 * (G + d.t) + 2 * nh + 2 * so\n  k.boltCircle('boltXylan', 'nut', d.bolts, d.boltCircle / 2, [0, yj - len / 2, 0], 'y', d.boltDia, len, [so + nh / 2, len - so - nh / 2])\n  k.torus('stainless', d.ring / 2, Math.max(0.004, d.boltDia * 0.18), [0, yj, 0], 'y', 4, 28)\n}\n\n/** Flanged spool between API faces y = 0 and y = L (local Y). */\nexport function apiSpool(k: Kit, d: ApiFlange, L: number, mat: MatKey) {\n  apiFlangeAt(k, d, 0, 1, mat)\n  apiFlangeAt(k, d, L, -1, mat)\n  tubeY(k, mat, d.od * 0.3, rBore(d), G + d.t - 0.002, L - G - d.t + 0.002, 24)\n}\n\n/** Open API face at local y = 0 (flange body below): machined face and ring groove; optional studs for a studded outlet. */\nexport function apiStuddedFace(k: Kit, d: ApiFlange, studs = true) {\n  const nh = d.boltDia * 0.95\n  k.lathe('machined', [[d.ring / 2 + d.boltDia * 0.6, -G + 0.0008], [rBore(d), -G + 0.0008]], [0, 0, 0], 'y', 32)\n  k.torus('pipeBore', d.ring / 2, Math.max(0.003, d.boltDia * 0.16), [0, -G + 0.001, 0], 'y', 6, 32)\n  if (studs) k.boltCircle('boltXylan', 'nut', d.bolts, d.boltCircle / 2, [0, -d.boltDia * 1.2, 0], 'y', d.boltDia, d.boltDia * 1.2 + G * 2 + d.t + nh * 1.4, [d.boltDia * 1.2 + G * 2 + d.t + nh * 0.5])\n}\n\nexport type GateOp = { kind: 'handwheel'; D: number; mat?: MatKey } | { kind: 'hydraulic'; r: number; len: number }\n\n/** API 6A slab gate valve between flange faces y = 0 and y = L, bonnet along local +Z. */\nexport function apiGate(k: Kit, d: ApiFlange, L: number, op: GateOp, mat: MatKey = 'valveBody') {\n  apiFlangeAt(k, d, 0, 1, mat)\n  apiFlangeAt(k, d, L, -1, mat)\n  const rb = rBore(d)\n  tubeY(k, mat, d.od * 0.3, rb, G + d.t - 0.002, L - G - d.t + 0.002, 24)\n  const w = d.od * 0.8\n  const dz = d.od * 0.92\n  const hb = Math.max(0.05, Math.min(d.od * 1.0, L - 2 * (G + d.t) - 0.02))\n  const yc = L / 2\n  k.rbox(mat, [w, hb, dz], Math.min(0.025, w * 0.1), [0, yc, 0])\n  k.cyl('brass', 0.008, 0.03, [w / 2 + 0.012, yc + hb * 0.25, 0], 'x', 8)\n  k.hex('brass', 0.017, 0.012, [w / 2 + 0.029, yc + hb * 0.25, 0], 'x')\n  const rbf = Math.min(w * 0.5, hb * 0.5)\n  let z = dz / 2\n  const tbf = Math.max(0.03, d.t * 0.6)\n  k.cyl(mat, rbf, tbf, [0, yc, z + tbf / 2], 'z', 28)\n  k.around(8, rbf * 0.8, [0, yc, z + tbf + 0.008], 'z', (p) => k.hex('nut', Math.max(0.022, d.boltDia * 1.3), 0.016, p, 'z'), Math.PI / 8)\n  z += tbf\n  const bl = 0.07 + d.od * 0.18\n  const rbo = rbf * 0.66\n  k.cyl(mat, rbo, bl, [0, yc, z + bl / 2], 'z', 24)\n  z += bl\n  if (op.kind === 'handwheel') {\n    const sr = Math.max(0.011, d.od * 0.045)\n    k.cyl(mat, rbo * 0.62, 0.04, [0, yc, z + 0.02], 'z', 18)\n    k.cyl('stainless', sr, 0.12, [0, yc, z + 0.06], 'z', 12)\n    k.handwheel(op.mat ?? 'actuatorRed', op.D, [0, yc, z + 0.1], 'z', 4)\n    k.hex('nut', sr * 2.4, 0.02, [0, yc, z + 0.12], 'z')\n    return z + 0.13\n  }\n  k.cyl(mat, rbo * 0.85, 0.05, [0, yc, z + 0.025], 'z', 24)\n  z += 0.05\n  const { r, len } = op\n  k.cyl('actuatorYellow', r * 1.08, 0.03, [0, yc, z + 0.015], 'z', 32)\n  k.cyl('actuatorYellow', r, len, [0, yc, z + 0.03 + len / 2], 'z', 32)\n  k.cyl('actuatorYellow', r * 1.08, 0.035, [0, yc, z + 0.03 + len + 0.0175], 'z', 32)\n  k.around(4, r * 1.04, [0, yc, 0], 'z', (p) => {\n    k.rod('stainless', [p[0], p[1], z], [p[0], p[1], z + 0.065 + len + 0.012], 0.007, 8)\n    k.hex('nut', 0.02, 0.012, [p[0], p[1], z + 0.065 + len + 0.006], 'z')\n  }, Math.PI / 4)\n  const zt = z + 0.065 + len\n  k.cyl('stainless', 0.02, 0.1, [0, yc, zt + 0.05], 'z', 14)\n  k.cyl('safetyGreen', 0.023, 0.025, [0, yc, zt + 0.03], 'z', 14)\n  k.cyl('nameplate', r * 1.002, 0.04, [0, yc, z + 0.03 + len * 0.62], 'z', 32)\n  return zt + 0.1\n}\n\n/** 1/2\" needle valve on a boss along +axis from base, T-handle to the side. Returns the outlet end. */\nexport function needleValve(k: Kit, base: V3, axis: 'x' | 'y' | 'z' | '-x' | '-z', handle: MatKey = 'actuatorRed'): V3 {\n  const v: V3 = axis === 'y' ? [0, 1, 0] : axis === 'z' ? [0, 0, 1] : axis === '-z' ? [0, 0, -1] : axis === 'x' ? [1, 0, 0] : [-1, 0, 0]\n  const p = (t: number): V3 => [base[0] + v[0] * t, base[1] + v[1] * t, base[2] + v[2] * t]\n  k.hex('stainless', 0.032, 0.02, p(0.01), axis)\n  k.box('stainless', [0.036, 0.036, 0.036], p(0.042))\n  const side: V3 = axis === 'y' ? [1, 0, 0] : [0, 1, 0]\n  const q = p(0.042)\n  k.rod('stainless', q, [q[0] + side[0] * 0.05, q[1] + side[1] * 0.05, q[2] + side[2] * 0.05], 0.006, 8)\n  k.box(handle, axis === 'y' ? [0.012, 0.012, 0.06] : [0.06, 0.012, 0.012], [q[0] + side[0] * 0.055, q[1] + side[1] * 0.055, q[2] + side[2] * 0.055])\n  k.rod('stainless', p(0.06), p(0.105), 0.008, 8)\n  return p(0.105)\n}\n\n/** Threaded hammer-lug cover (mud pump valve cover, union nut) on local +Y from y = 0: body, three lugs, centre boss. */\nexport function hammerCap(k: Kit, mat: MatKey, R: number, h: number, lugs = 3) {\n  k.cyl(mat, R, h, [0, h / 2, 0], 'y', 32)\n  k.cyl(mat, R * 0.96, 0.006, [0, h + 0.003, 0], 'y', 32)\n  for (let i = 0; i < lugs; i++) {\n    const a = (i / lugs) * Math.PI * 2\n    k.box(mat, [R * 0.42, h * 0.8, R * 0.32], [Math.cos(a) * (R + R * 0.14), h * 0.45, Math.sin(a) * (R + R * 0.14)], new THREE.Quaternion().setFromAxisAngle(new THREE.Vector3(0, 1, 0), -a))\n  }\n  k.cyl(mat, R * 0.35, h * 0.3, [0, h + h * 0.15, 0], 'y', 20)\n}\n\n// ------------------------------------------------------------------ 2D helpers\n\n/** Convex hull (CCW) of circles sampled at n points each, as a THREE.Shape in XY. */\nexport function hullOfCircles(circles: [number, number, number][], n = 32): THREE.Shape {\n  const hull = hullPoints(circles, n)\n  const sh = new THREE.Shape()\n  hull.forEach(([x, y], i) => (i ? sh.lineTo(x, y) : sh.moveTo(x, y)))\n  sh.closePath()\n  return sh\n}\n\n/** Plate band of `thick` × `depth` (along +Z from z0) following the hull outline of the circles: guard rims, belts. */\nexport function hullBand(k: Kit, mat: MatKey, circles: [number, number, number][], z0: number, depth: number, thick: number, n = 48) {\n  const h = hullPoints(circles, n)\n  for (let i = 0; i < h.length; i++) {\n    const [x0, y0] = h[i]\n    const [x1, y1] = h[(i + 1) % h.length]\n    const len = Math.hypot(x1 - x0, y1 - y0)\n    if (len < 1e-5) continue\n    const nx = (y1 - y0) / len\n    const ny = -(x1 - x0) / len\n    const q = new THREE.Quaternion().setFromAxisAngle(new THREE.Vector3(0, 0, 1), Math.atan2(y1 - y0, x1 - x0))\n    k.box(mat, [len + thick * 0.6, thick, depth], [(x0 + x1) / 2 - (nx * thick) / 2, (y0 + y1) / 2 - (ny * thick) / 2, z0 + depth / 2], q)\n  }\n}\n\n/** Flat convex face (triangle fan) over the hull of the circles at z = z0, facing +Z (or -Z when flip). UVs scaled by uvs per metre. */\nexport function hullFace(circles: [number, number, number][], z0: number, flip = false, uvs = 1, n = 48): THREE.BufferGeometry {\n  const h = hullPoints(circles, n)\n  const cx = h.reduce((a, p) => a + p[0], 0) / h.length\n  const cy = h.reduce((a, p) => a + p[1], 0) / h.length\n  const pos: number[] = [cx, cy, z0]\n  const uv: number[] = [cx * uvs, cy * uvs]\n  for (const [x, y] of h) pos.push(x, y, z0), uv.push(x * uvs, y * uvs)\n  const idx: number[] = []\n  for (let i = 0; i < h.length; i++) {\n    const a = 1 + i\n    const b = 1 + ((i + 1) % h.length)\n    if (flip) idx.push(0, b, a)\n    else idx.push(0, a, b)\n  }\n  const g = new THREE.BufferGeometry()\n  g.setAttribute('position', new THREE.Float32BufferAttribute(pos, 3))\n  g.setAttribute('uv', new THREE.Float32BufferAttribute(uv, 2))\n  g.setIndex(idx)\n  g.computeVertexNormals()\n  return g\n}\n\nexport function hullPoints(circles: [number, number, number][], n = 32): [number, number][] {\n  const pts: [number, number][] = []\n  for (const [x, y, r] of circles) for (let i = 0; i < n; i++) pts.push([x + Math.cos((i / n) * Math.PI * 2) * r, y + Math.sin((i / n) * Math.PI * 2) * r])\n  pts.sort((a, b) => a[0] - b[0] || a[1] - b[1])\n  const cross = (o: number[], a: number[], b: number[]) => (a[0] - o[0]) * (b[1] - o[1]) - (a[1] - o[1]) * (b[0] - o[0])\n  const lo: [number, number][] = []\n  for (const p of pts) {\n    while (lo.length >= 2 && cross(lo[lo.length - 2], lo[lo.length - 1], p) <= 0) lo.pop()\n    lo.push(p)\n  }\n  const up: [number, number][] = []\n  for (let i = pts.length - 1; i >= 0; i--) {\n    const p = pts[i]\n    while (up.length >= 2 && cross(up[up.length - 2], up[up.length - 1], p) <= 0) up.pop()\n    up.push(p)\n  }\n  return lo.slice(0, -1).concat(up.slice(0, -1))\n}\n\n/** Matrix placing a shape drawn in XY onto the plane z = z0 (extrusion runs +Z). */\nexport const atZ = (z0: number) => mat4([0, 0, z0])\n\n// ------------------------------------------------------------------ subsea\n\n/** Mudmat at seabed y = 0: plate with edge skirts and stiffeners, top at y = t. */\nexport function mudmat(k: Kit, sx: number, sz: number, opts: { skirt?: number; mat?: MatKey; t?: number; padeyes?: boolean } = {}) {\n  const mat = opts.mat ?? 'equipmentPaint'\n  const t = opts.t ?? 0.05\n  const sk = opts.skirt ?? 0.25\n  // skirts penetrate the seabed; plate rests on y = 0\n  k.box(mat, [sx, t, sz], [0, t / 2, 0])\n  for (const s of [-1, 1]) {\n    k.box(mat, [sx, sk, 0.02], [0, -sk / 2 + 0.001, s * (sz / 2 - 0.01)])\n    k.box(mat, [0.02, sk, sz], [s * (sx / 2 - 0.01), -sk / 2 + 0.001, 0])\n  }\n  // top stiffeners (flat bar grid)\n  const nx = Math.max(2, Math.round(sx / 1.2))\n  const nz = Math.max(2, Math.round(sz / 1.2))\n  for (let i = 1; i < nx; i++) k.box(mat, [0.02, 0.12, sz - 0.04], [-sx / 2 + (sx * i) / nx, t + 0.06, 0])\n  for (let i = 1; i < nz; i++) k.box(mat, [sx - 0.04, 0.12, 0.02], [0, t + 0.06, -sz / 2 + (sz * i) / nz])\n  // water evacuation holes with short vent pipes\n  for (const sx2 of [-1, 1]) for (const sz2 of [-1, 1]) k.cyl(mat, 0.05, 0.1, [sx2 * sx * 0.25, t + 0.05, sz2 * sz * 0.25], 'y', 14, undefined, true)\n  if (opts.padeyes !== false) for (const sx2 of [-1, 1]) for (const sz2 of [-1, 1]) liftingLug(k, [sx2 * (sx / 2 - 0.2), t, sz2 * (sz / 2 - 0.2)], 'x', 1.6, mat)\n  return t\n}\n\n/** Sacrificial Al-Zn-In anode on two stand-off cores, long axis along `along`, centre at c. */\nexport function anode(k: Kit, c: V3, along: 'x' | 'z' | 'y' = 'x', len = 0.6) {\n  const size: V3 = along === 'x' ? [len, 0.1, 0.12] : along === 'z' ? [0.12, 0.1, len] : [0.12, len, 0.1]\n  k.rbox('aluminium', size, 0.02, c)\n  const off = len * 0.55\n  for (const s of [-1, 1]) {\n    const e: V3 = along === 'x' ? [c[0] + s * off, c[1], c[2]] : along === 'z' ? [c[0], c[1], c[2] + s * off] : [c[0], c[1] + s * off, c[2]]\n    k.box('stainless', [0.03, 0.03, 0.03], e)\n  }\n}\n\n/** API 17D class-4 style ROV torque bucket (open funnel, docking slots) on a back plate, opening facing +axis from base. */\nexport function rovBucket(k: Kit, base: V3, axis: 'x' | 'z' | '-x' | '-z' | 'y', mat: MatKey = 'actuatorYellow') {\n  k.at(base, axis, () => {\n    k.box(mat, [0.32, 0.02, 0.32], [0, 0.01, 0])\n    k.lathe(mat, [[0.07, 0.02], [0.085, 0.02], [0.085, 0.14], [0.11, 0.2], [0.1, 0.2], [0.075, 0.145], [0.07, 0.145]], [0, 0, 0], 'y', 28)\n    k.lathe('pipeBore', [[0.07, 0.145], [0.07, 0.025]], [0, 0, 0], 'y', 24)\n    k.cyl('stainless', 0.025, 0.04, [0, 0.04, 0], 'y', 6)\n    // docking slots and position indicator\n    for (const s of [-1, 1]) k.box('nut', [0.016, 0.05, 0.02], [s * 0.086, 0.17, 0])\n    k.box('nut', [0.05, 0.004, 0.012], [0, 0.022, 0.13])\n  })\n}\n\n/** Subsea hub (API 17R style) on local +Y from base along `axis`: neck, tapered hub profile, seal face, optional protection cap. Face at base + len. */\nexport function hub(k: Kit, base: V3, axis: 'x' | 'z' | '-x' | '-z' | 'y', boreR: number, len = 0.35, mat: MatKey = 'equipmentPaint', cap = false) {\n  const R = Math.max(0.12, boreR * 2.4)\n  k.at(base, axis, () => {\n    k.lathe(mat, [[boreR, 0], [R * 0.62, 0], [R * 0.62, len * 0.5], [R * 0.75, len * 0.55], [R, len * 0.78], [R, len - 0.01], [R * 0.95, len], [boreR, len]], [0, 0, 0], 'y', 40)\n    k.lathe('pipeBore', [[boreR, len], [boreR, -0.05]], [0, 0, 0], 'y', 24)\n    k.lathe('machined', [[R * 0.9, len + 0.001], [boreR, len + 0.001]], [0, 0, 0], 'y', 32)\n    if (cap) {\n      k.lathe('actuatorYellow', [[R * 1.04, len * 0.7], [R * 1.12, len * 0.7], [R * 1.12, len + 0.12], [0, len + 0.12]], [0, 0, 0], 'y', 32)\n      k.cyl('nut', 0.03, 0.08, [0, len + 0.16, 0], 'y', 6)\n    }\n  })\n  return R\n}\n\n// ------------------------------------------------------------------ chokes and frames\n\n/** Matrix whose local +Y runs along `flow` and local +Z along `stem`, origin at o. */\nexport function flowFrame(o: V3, flow: V3, stem: V3) {\n  const y = new THREE.Vector3(...flow).normalize()\n  const z = new THREE.Vector3(...stem).normalize()\n  const x = new THREE.Vector3().crossVectors(y, z).normalize()\n  return new THREE.Matrix4().makeBasis(x, y, z).setPosition(...o)\n}\n\n/** Angle choke body centre to inlet face (along -X) and outlet face (along -Y). */\nexport function chokeDims(dIn: ApiFlange, dOut: ApiFlange) {\n  const B = Math.max(dIn.od, dOut.od) * 0.78\n  return { B, inlet: B / 2 + 0.05 + dIn.t + G, outlet: B / 2 + 0.09 + dOut.t + G }\n}\n\n/** API 6A / 16C angle choke at the origin: inlet face -X, outlet face -Y, bonnet and operator along +X. Returns the operator's far x and its hose port. */\nexport function angleChoke(k: Kit, dIn: ApiFlange, dOut: ApiFlange, op: 'handwheel' | 'hydraulic', mat: MatKey = 'valveBody'): { end: number; hose: V3 | null } {\n  const { B, inlet, outlet } = chokeDims(dIn, dOut)\n  k.rbox(mat, [B, B, B * 0.92], Math.min(0.03, B * 0.12), [0, 0, 0])\n  const ri = rBore(dIn)\n  const ro = rBore(dOut)\n  k.at([0, 0, 0], '-x', () => tubeY(k, mat, dIn.od * 0.3, ri, B / 2 - 0.01, inlet - dIn.t - G + 0.002, 24))\n  k.at([0, 0, 0], '-y', () => tubeY(k, mat, dOut.od * 0.32, ro, B / 2 - 0.01, outlet - dOut.t - G + 0.002, 24))\n  k.at([-inlet, 0, 0], '-x', () => apiFlangeAt(k, dIn, 0, -1, mat))\n  k.at([0, -outlet, 0], '-y', () => apiFlangeAt(k, dOut, 0, -1, mat))\n  k.disc('internalsSteel', ri * 1.02, [-B / 2 + 0.012, 0, 0], '-x', 24)\n  k.disc('internalsSteel', ro * 1.02, [0, -B / 2 + 0.012, 0], '-y', 24)\n  // bonnet nut with hammer lugs, wear-sleeve bleed plug\n  k.hex('brass', 0.022, 0.016, [0, B * 0.25, B * 0.46 + 0.008], 'z')\n  let x = B / 2\n  const rn = B * 0.42\n  k.cyl(mat, rn * 0.9, 0.02, [x + 0.01, 0, 0], 'x', 24)\n  x += 0.02\n  k.cyl(mat, rn, 0.07, [x + 0.035, 0, 0], 'x', 24)\n  for (let i = 0; i < 3; i++) {\n    const a = (i / 3) * Math.PI * 2 + Math.PI / 6\n    k.box(mat, [0.06, B * 0.2, B * 0.13], [x + 0.035, Math.cos(a) * (rn + B * 0.08), Math.sin(a) * (rn + B * 0.08)], new THREE.Quaternion().setFromEuler(new THREE.Euler(-a, 0, 0)))\n  }\n  x += 0.07\n  k.cyl(mat, rn * 0.55, 0.1, [x + 0.05, 0, 0], 'x', 20)\n  x += 0.1\n  // position indicator sleeve with scale ring\n  k.cyl('stainless', rn * 0.4, 0.07, [x + 0.035, 0, 0], 'x', 18)\n  k.cyl('brass', rn * 0.42, 0.012, [x + 0.02, 0, 0], 'x', 18)\n  k.box('actuatorRed', [0.05, 0.004, 0.012], [x + 0.035, rn * 0.4 + 0.002, 0])\n  x += 0.07\n  if (op === 'handwheel') {\n    k.cyl('stainless', 0.014, 0.1, [x + 0.05, 0, 0], 'x', 10)\n    k.handwheel('actuatorRed', Math.max(0.3, B * 1.3), [x + 0.07, 0, 0], 'x', 4)\n    return { end: x + 0.1, hose: null }\n  }\n  // hydraulic rotary-gear actuator: gear housing, double-acting cylinder, stroke-position transmitter, manual override\n  k.cyl(mat, rn * 0.35, 0.05, [x + 0.025, 0, 0], 'x', 14)\n  x += 0.05\n  k.rbox('actuatorYellow', [0.2, 0.2, 0.22], 0.02, [x + 0.1, 0, 0])\n  k.cyl('actuatorYellow', 0.07, 0.42, [x + 0.1, 0.0, 0.0], 'z', 24)\n  for (const s of [-1, 1]) {\n    k.cyl('actuatorYellow', 0.078, 0.03, [x + 0.1, 0, s * 0.21], 'z', 24)\n    k.cyl('brass', 0.012, 0.035, [x + 0.1, 0.07, s * 0.19], 'y', 10)\n  }\n  k.rbox('stainless', [0.08, 0.1, 0.08], 0.01, [x + 0.25, 0.04, 0])\n  k.cyl('brass', 0.012, 0.03, [x + 0.25, -0.025, 0], 'y', 8)\n  k.hex('nut', 0.04, 0.03, [x + 0.1, 0.115, 0], 'y')\n  k.box('nameplate', [0.08, 0.05, 0.002], [x + 0.1, 0, 0.112])\n  return { end: x + 0.3, hose: [x + 0.1, 0.105, 0.19] }\n}\n\n/** ROV grab bar: round bar on two stand-offs between a and b, stand-off along `out`. */\nexport function grabBar(k: Kit, a: V3, b: V3, out: V3, mat: MatKey = 'actuatorYellow') {\n  const o = (p: V3, s: number): V3 => [p[0] + out[0] * s, p[1] + out[1] * s, p[2] + out[2] * s]\n  k.rod(mat, o(a, 0.12), o(b, 0.12), 0.019, 10)\n  for (const p of [a, b]) k.rod(mat, p, o(p, 0.12), 0.015, 8)\n}\n\n/**\n * Subsea trunnion ball valve between flange faces at local y = 0 and y = L (flow along local +Y), stem and ROV\n * gearbox along local +Z, class-4 bucket opening along local +X. Flanges are ASME B16.5 with studs.\n */\nexport function subseaBall(k: Kit, nps: number, rating: Rating, L: number, opts: { mat?: MatKey; bucket?: boolean } = {}) {\n  const mat = opts.mat ?? 'equipmentPaint'\n  const f = asmeFlange(nps, rating)\n  const ro = pipeOD(nps) / 2\n  const R = Math.max(ro * 1.55, f.od * 0.5)\n  const fl0 = flange(k, nps, rating, [0, 0, 0], '-y', true)\n  const fl1 = flange(k, nps, rating, [0, L, 0], 'y', true)\n  k.lathe(mat, [[ro * 0.9, fl0], [ro * 1.15, fl0], [R * 0.9, L / 2 - R * 0.75], [R, L / 2 - R * 0.35], [R, L / 2 + R * 0.35], [R * 0.9, L / 2 + R * 0.75], [ro * 1.15, L - fl1], [ro * 0.9, L - fl1]], [0, 0, 0], 'y', 40)\n  // body joint flanges with studs\n  for (const y of [L / 2 - R * 0.55, L / 2 + R * 0.55]) {\n    k.cyl(mat, R * 1.08, 0.05, [0, y, 0], 'y', 40)\n    k.around(12, R * 1.0, [0, y + 0.035, 0], 'y', (p) => k.hex('nut', 0.05, 0.03, p, 'y'), Math.PI / 12)\n  }\n  // trunnion boss, stem housing, gearbox, bucket\n  k.cyl(mat, R * 0.45, R * 0.4, [0, L / 2, -R - R * 0.15], 'z', 28)\n  k.cyl(mat, R * 0.45, R * 0.5, [0, L / 2, R + R * 0.2], 'z', 28)\n  const zg = R + R * 0.45\n  k.rbox('actuatorYellow', [R * 1.5, R * 1.2, R * 0.9], 0.04, [0, L / 2, zg + R * 0.45])\n  k.cyl('actuatorYellow', R * 0.35, R * 0.5, [R * 0.75, L / 2, zg + R * 0.45], 'x', 24)\n  k.box('nut', [0.12, 0.02, 0.012], [0, L / 2 + R * 0.6, zg + R * 0.9 + 0.006])\n  if (opts.bucket !== false) rovBucket(k, [R * 1.05 + 0.25, L / 2, zg + R * 0.45], 'x')\n  k.rod('actuatorYellow', [R * 1.0, L / 2, zg + R * 0.45], [R * 1.05 + 0.25, L / 2, zg + R * 0.45], 0.05, 16)\n  return { R, top: zg + R * 0.9 }\n}\n"
    }
  ],
  "categories": [
    "drilling"
  ],
  "meta": {
    "export": "SubseaManifold",
    "anchor": "base",
    "standards": [
      "API 17P",
      "API 17D",
      "API 17R",
      "ISO 13628-15",
      "DNV-RP-B401"
    ]
  }
}