{
  "$schema": "https://ui.shadcn.com/schema/registry-item.json",
  "name": "subsea-tree",
  "type": "registry:component",
  "title": "Subsea Xmas tree (vertical)",
  "description": "Vertical subsea tree on a permanent guide base: skirted mudmat, four guide posts on a 6 ft radius, conductor and 18-3/4\" wellhead housings, hydraulic collet wellhead connector with lock indicators, valve block with fail-safe hydraulic actuators (PMV, PSV, AMV, PWV, AWV) and ROV overrides, production wing to a retrievable choke module (stepping actuator, ROV handle and bucket) and a horizontal flowline hub, annulus wing to an annulus hub, tree cap with T-handle. Tree frame with guide funnels, ROV panel (buckets, hot stabs, gauges), subsea control module with hydraulic tubing, grab bars, anodes and nameplate. y = 0 is the seabed. Hub ports carry no flange rating.",
  "dependencies": [
    "three",
    "@react-three/fiber",
    "@react-three/drei"
  ],
  "registryDependencies": [
    "https://deck-twin-library.demoareeb.workers.dev/r/core.json"
  ],
  "files": [
    {
      "path": "drilling/subsea-tree.tsx",
      "type": "registry:file",
      "target": "components/deck-twin/drilling/subsea-tree.tsx",
      "content": "// Vertical subsea Xmas tree (API 17D / ISO 13628-4) on a permanent guide base: wellhead connector, valve block with hydraulic actuators, choke module, flowline and annulus hubs, tree cap, ROV panel, SCM.\nimport { Kit, type MatKey, type V3 } from '../lib/kit'\nimport { paint } from '../lib/materials'\nimport { KitMesh, useKit, type GroupProps } from '../lib/mesh'\nimport { gauge, member, pipeRun, tagPlate } from '../lib/parts'\nimport { port, type Port } from '../lib/ports'\nimport { pipeID, pipeOD } from '../lib/standards'\nimport { anode, grabBar, hub, mudmat, rovBucket } from './_drill'\n\nexport interface SubseaTreeProps {\n  /** production bore / flowline hub, inches */\n  productionNps?: number\n  /** annulus hub, inches */\n  annulusNps?: number\n  /** retrievable choke module on the production wing */\n  choke?: boolean\n  tag?: string\n}\n\nexport const defaults = { productionNps: 5, annulusNps: 2, choke: true, tag: 'XT-S01' }\n\ntype P = typeof defaults\n\n// typical 5\" x 2\" vertical tree on an 18-3/4\" wellhead: guide posts on a 6 ft radius, ~5 m tall\nconst GP = 1.29\nconst yWH = 1.6\nconst yC0 = 1.75\nconst yC1 = 2.75\nconst yB1 = 4.05\n\nfunction layout(p: Omit<P, 'tag'>) {\n  const yW = 3.45\n  const xHub = 2.0\n  const yHubP = 2.1\n  return {\n    yW,\n    prod: [xHub + 0.4, yHubP, 0] as V3,\n    ann: [-xHub - 0.25, 3.2, 0] as V3,\n    xHub,\n    yHubP,\n    boreP: pipeID(p.productionNps, '160') / 2,\n    boreA: pipeID(p.annulusNps, '160') / 2,\n  }\n}\n\nexport function ports(props: SubseaTreeProps = {}): Port[] {\n  const p = { ...defaults, ...props }\n  const l = layout(p)\n  return [port('production', 'Production flowline hub', l.prod, [1, 0, 0], p.productionNps), port('annulus', 'Annulus hub', l.ann, [-1, 0, 0], p.annulusNps, undefined, 'utility')]\n}\n\nexport function SubseaTree({ productionNps = defaults.productionNps, annulusNps = defaults.annulusNps, choke = defaults.choke, tag = defaults.tag, ...group }: SubseaTreeProps & GroupProps) {\n  const data = useKit((k) => build(k, { productionNps, annulusNps, choke, tag }), [productionNps, annulusNps, choke, tag])\n  return (\n    <group {...group}>\n      <KitMesh data={data} />\n    </group>\n  )\n}\n\nexport default SubseaTree\n\n/** Fail-safe-close hydraulic gate valve actuator from base along `axis`: bonnet, spring cylinder, cap, ROV override bucket. */\nfunction actuator(k: Kit, base: V3, axis: 'x' | '-x' | 'z' | '-z', r: number, len: number, mat: MatKey, override = true) {\n  k.at(base, axis, () => {\n    k.cyl('valveBody', r * 0.7, 0.12, [0, 0.06, 0], 'y', 24)\n    k.cyl('valveBody', r * 0.95, 0.04, [0, 0.14, 0], 'y', 28)\n    k.cyl(mat, r, len, [0, 0.16 + len / 2, 0], 'y', 32)\n    k.cyl(mat, r * 1.06, 0.04, [0, 0.16 + len, 0], 'y', 32)\n    k.cyl('nameplate', r * 1.003, 0.05, [0, 0.16 + len * 0.6, 0], 'y', 32)\n    k.cyl('stainless', 0.025, 0.08, [0, 0.2 + len + 0.04, 0], 'y', 10)\n    k.cyl('safetyGreen', 0.03, 0.03, [0, 0.2 + len + 0.03, 0], 'y', 12)\n  })\n  if (override) {\n    const v = axis === 'x' ? [1, 0, 0] : axis === '-x' ? [-1, 0, 0] : axis === 'z' ? [0, 0, 1] : [0, 0, -1]\n    rovBucket(k, [base[0] + v[0] * (0.3 + len), base[1], base[2] + v[2] * (0.3 + len)], axis, mat)\n  }\n}\n\nfunction build(k: Kit, p: P) {\n  const l = layout(p)\n  const { yW, xHub, yHubP, boreP, boreA } = l\n  const Y: MatKey = paint('#e3b512', { weather: 1 })\n  const B: MatKey = paint('#c7cbcd', { weather: 0.8 })\n  const A: MatKey = paint('#e3b512', { weather: 0.6 })\n\n  // permanent guide base on a mudmat, guide posts, conductor and wellhead housings\n  const t = mudmat(k, 4.8, 4.8, { mat: Y, skirt: 0.35 })\n  for (const sx of [-1, 1])\n    for (const sz of [-1, 1]) {\n      member(k, [sx * GP, t, sz * GP], [sx * GP, 3.6, sz * GP], 'CHS 219.1x10', Y)\n      k.dome(Y, 0.11, [sx * GP, 3.6, sz * GP], 'y', Math.PI / 2, 1.6, 16)\n      member(k, [sx * GP, t + 0.4, sz * GP], [0, t + 0.4, sz * 0.5], 'HEA 200', Y)\n    }\n  k.lathe(Y, [[0.46, t], [1.0, t], [1.0, t + 0.6], [0.5, t + 0.6]], [0, 0, 0], 'y', 48)\n  k.lathe('flangeSteel', [[0.4, t], [0.457, t], [0.457, yWH], [0.4, yWH]], [0, 0, 0], 'y', 48)\n  k.lathe('flangeSteel', [[0.3, yWH - 0.05], [0.42, yWH - 0.05], [0.42, yC0 + 0.1], [0.3, yC0 + 0.1]], [0, 0, 0], 'y', 48)\n\n  // wellhead connector: hydraulic collet body, lock-indicator rods, unlock / lock lines\n  k.lathe(B, [[0.35, yC0], [0.62, yC0], [0.66, yC0 + 0.1], [0.66, yC1 - 0.15], [0.55, yC1 - 0.05], [0.45, yC1]], [0, 0, 0], 'y', 56)\n  for (let i = 0; i < 6; i++) {\n    const a = (i / 6) * Math.PI * 2\n    k.rod('stainless', [Math.cos(a) * 0.67, yC0 + 0.4, Math.sin(a) * 0.67], [Math.cos(a) * 0.67, yC1 + 0.05, Math.sin(a) * 0.67], 0.018, 8)\n  }\n\n  // tree body (master valve block) with production bore up the middle\n  const bw = 0.95\n  k.rbox(B, [bw, yB1 - yC1, bw], 0.04, [0, (yC1 + yB1) / 2, 0])\n  // production master (PMV) and production swab (PSV) actuators along +Z / -Z, annulus master (AMV) along -Z\n  actuator(k, [0, 3.05, bw / 2], 'z', 0.2, 0.5, A, false)\n  actuator(k, [0, 3.75, -bw / 2], '-z', 0.17, 0.45, A)\n  actuator(k, [-0.22, 3.05, -bw / 2], '-z', 0.15, 0.4, A, false)\n\n  // production wing: wing block, PWV actuator, crossover valve, then the choke module and the flowline hub\n  const xw = bw / 2\n  k.rbox(B, [0.55, 0.5, 0.5], 0.03, [xw + 0.27, yW, 0])\n  actuator(k, [xw + 0.27, yW, 0.25], 'z', 0.17, 0.45, A)\n  const prodOD = pipeOD(p.productionNps) * 1.3\n  if (p.choke) {\n    const xc = xw + 0.95\n    k.pipe(B, [xw + 0.54, yW, 0], [xc - 0.3, yW, 0], prodOD, prodOD * 0.2)\n    // choke module: body, insert clamp, stepping actuator, ROV handle and bucket\n    k.rbox(B, [0.6, 0.6, 0.6], 0.04, [xc, yW, 0])\n    k.cyl(B, 0.32, 0.12, [xc, yW + 0.36, 0], 'y', 32)\n    k.lathe('actuatorYellow', [[0.34, yW + 0.36], [0.4, yW + 0.36], [0.4, yW + 0.5], [0.34, yW + 0.5]], [xc, 0, 0], 'y', 32)\n    k.cyl(A, 0.24, 0.5, [xc, yW + 0.75, 0], 'y', 32)\n    k.cyl(A, 0.26, 0.05, [xc, yW + 1.02, 0], 'y', 32)\n    k.rod(Y, [xc - 0.25, yW + 1.12, 0], [xc + 0.25, yW + 1.12, 0], 0.025, 10)\n    for (const s of [-1, 1]) k.rod(Y, [xc + s * 0.25, yW + 1.12, 0], [xc + s * 0.25, yW + 1.04, 0], 0.02, 8)\n    rovBucket(k, [xc, yW + 0.75, 0.25], 'z', A)\n    pipeRun(k, [[xc, yW - 0.3, 0], [xc, yHubP, 0], [xHub, yHubP, 0]], p.productionNps, { mat: B, schedule: '160' })\n  } else {\n    pipeRun(k, [[xw + 0.54, yW, 0], [xw + 1.1, yW, 0], [xw + 1.1, yHubP, 0], [xHub, yHubP, 0]], p.productionNps, { mat: B, schedule: '160' })\n  }\n  hub(k, [xHub, yHubP, 0], 'x', boreP, 0.4, B)\n  k.box(Y, [0.3, yHubP - 0.2 - t, 0.5], [xHub - 0.15, (yHubP - 0.2 + t) / 2, 0])\n\n  // annulus wing to the annulus hub (-X)\n  pipeRun(k, [[-bw / 2, 3.2, 0], [-l.xHub, 3.2, 0]], p.annulusNps, { mat: B, schedule: '160' })\n  k.rbox(B, [0.4, 0.4, 0.4], 0.03, [-bw / 2 - 0.3, 3.2, 0])\n  actuator(k, [-bw / 2 - 0.3, 3.2, -0.2], '-z', 0.13, 0.35, A, false)\n  hub(k, [-l.xHub, 3.2, 0], '-x', boreA, 0.25, B)\n  member(k, [-l.xHub + 0.12, t, 0], [-l.xHub + 0.12, 3.2 - 0.13, 0], 'CHS 114.3x6.3', Y)\n  k.box(Y, [0.2, 0.03, 0.26], [-l.xHub + 0.12, 3.2 - 0.12, 0])\n  member(k, [-l.xHub + 0.12, 2.4, 0], [-GP, 2.4, -GP], 'CHS 88.9x6.3', Y)\n\n  // tree cap with ROV T-handle and lift eye\n  k.lathe(Y, [[0, yB1], [0.42, yB1], [0.42, yB1 + 0.35], [0.36, yB1 + 0.45], [0, yB1 + 0.45]], [0, 0, 0], 'y', 40)\n  k.rod(Y, [0, yB1 + 0.45, 0], [0, yB1 + 0.75, 0], 0.04, 10)\n  k.rod(Y, [-0.22, yB1 + 0.75, 0], [0.22, yB1 + 0.75, 0], 0.03, 10)\n  k.torus('nut', 0.06, 0.015, [0, yB1 + 0.85, 0], 'z', 6, 16)\n\n  // tree frame: guide funnels on the posts, columns, top ring beams, struts to the connector\n  const yT = 4.75\n  for (const sx of [-1, 1])\n    for (const sz of [-1, 1]) {\n      const x = sx * GP\n      const z = sz * GP\n      k.lathe(Y, [[0.13, 1.9], [0.15, 1.9], [0.38, 1.4], [0.36, 1.4], [0.13, 1.9]], [x, 0, z], 'y', 28)\n      k.lathe(Y, [[0.125, 1.9], [0.16, 1.9], [0.16, 3.9], [0.125, 3.9]], [x, 0, z], 'y', 28)\n      member(k, [x, 3.9, z], [x, yT, z], 'CHS 219.1x10', Y)\n      member(k, [x, 2.4, z], [sx * 0.47, 2.4, sz * 0.47], 'CHS 114.3x6.3', Y)\n    }\n  for (const s of [-1, 1]) {\n    member(k, [-GP, yT, s * GP], [GP, yT, s * GP], 'CHS 168.3x8', Y)\n    member(k, [s * GP, yT, -GP], [s * GP, yT, GP], 'CHS 168.3x8', Y)\n  }\n  for (const sx of [-1, 1]) for (const sz of [-1, 1]) member(k, [sx * GP, yT, sz * GP], [sx * 0.45, yB1 - 0.2, sz * 0.45], 'CHS 114.3x6.3', Y)\n\n  // ROV panel on the +Z face: buckets for the manual overrides, hot-stab receptacles, gauges\n  const zp = GP + 0.02\n  k.box(Y, [2 * GP, 1.3, 0.03], [0, 3.2, zp])\n  for (const x of [-0.8, -0.4, 0.2]) {\n    rovBucket(k, [x, 3.5, zp + 0.015], 'z', A)\n    k.box('nut', [0.18, 0.05, 0.003], [x, 3.25, zp + 0.017])\n  }\n  for (const x of [0.6, 0.75]) {\n    k.cyl('stainless', 0.05, 0.12, [x, 3.45, zp + 0.06], 'z', 16)\n    k.cyl('actuatorRed', 0.055, 0.03, [x, 3.45, zp + 0.13], 'z', 16)\n  }\n  gauge(k, [0.68, 2.85, zp + 0.06], 'z', 0.14, 'bar', 700)\n  gauge(k, [-0.4, 2.85, zp + 0.06], 'z', 0.14, 'bar', 700)\n\n  // subsea control module on its mounting base (-Z side), hydraulic jumpers to the actuators\n  const zs = -GP + 0.35\n  const xS = 0.65\n  member(k, [-GP, 2.2, -GP], [GP, 2.2, -GP], 'HEA 200', Y)\n  k.box(Y, [1.0, 0.1, 0.7], [xS, 2.35, zs])\n  k.rbox(A, [0.8, 1.0, 0.55], 0.05, [xS, 2.9, zs])\n  k.rod(Y, [xS - 0.25, 3.4, zs], [xS + 0.25, 3.4, zs], 0.03, 10)\n  for (const s of [-1, 1]) k.rod(Y, [xS + s * 0.25, 3.4, zs], [xS + s * 0.25, 3.5, zs], 0.025, 8)\n  k.box('nameplate', [0.3, 0.15, 0.003], [xS, 3.0, zs - 0.277])\n  const tubes: V3[] = [[0, 3.05, bw / 2 + 0.08], [0, 3.75, -bw / 2 - 0.08], [-0.22, 3.05, -bw / 2 - 0.08], [xw + 0.27, yW, 0.33]]\n  tubes.forEach((e, i) => k.tube('tubing', [[xS + (i - 1.5) * 0.08, 3.4, zs + 0.2], [xS + (i - 1.5) * 0.08, 3.6 + i * 0.03, zs + 0.2], [e[0], 3.6 + i * 0.03, zs + 0.2], [e[0], 3.6 + i * 0.03, e[2]], e], 0.008, 6, 0.06))\n\n  for (const sx of [-1, 1]) grabBar(k, [sx * (GP + 0.1), 4.2, -0.7], [sx * (GP + 0.1), 4.2, 0.7], [sx, 0, 0])\n  for (const sx of [-1, 1]) for (const sz of [-1, 1]) anode(k, [sx * (GP + 0.2), 2.8, sz * GP], 'y', 0.5)\n  k.box(Y, [1.0, 0.36, 0.02], [0, 4.35, zp])\n  tagPlate(k, p.tag, [0, 4.35, zp + 0.012], 'z', 0.9)\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": "SubseaTree",
    "anchor": "base",
    "standards": [
      "API 17D",
      "ISO 13628-4",
      "API 17H (ROV interfaces)",
      "API 6A"
    ]
  }
}