build

Build

Source code in process/models/build.py

class Build:
    def __init__(self):
        self.outfile = constants.NOUT
        self.mfile = constants.MFILE

    def run(self):
        self.calculate_radial_build(output=False)
        self.calculate_vertical_build(output=False)

        (
            current_drive_variables.radius_beam_tangency,
            current_drive_variables.radius_beam_tangency_max,
        ) = self.calculate_beam_port_size(
            f_radius_beam_tangency_rmajor=current_drive_variables.f_radius_beam_tangency_rmajor,
            rmajor=physics_variables.rmajor,
            n_tf_coils=tfcoil_variables.n_tf_coils,
            dx_tf_inboard_out_toroidal=tfcoil_variables.dx_tf_inboard_out_toroidal,
            dr_tf_outboard=build_variables.dr_tf_outboard,
            r_tf_outboard_mid=build_variables.r_tf_outboard_mid,
            dx_beam_duct=current_drive_variables.dx_beam_duct,
            dx_beam_shield=current_drive_variables.dx_beam_shield,
        )

    def calculate_beam_port_size(
        self,
        f_radius_beam_tangency_rmajor: float,
        rmajor: float,
        n_tf_coils: int,
        dx_tf_inboard_out_toroidal: float,
        dr_tf_outboard: float,
        r_tf_outboard_mid: float,
        dx_beam_duct: float,
        dx_beam_shield: float,
    ) -> tuple[float, float]:
        """Calculates the maximum possible tangency radius for adequate beam access.

        Parameters
        ----------
        f_radius_beam_tangency_rmajor : float
            Fraction of rmajor for beam tangency
        rmajor : float
            Major radius
        n_tf_coils : int
            Number of TF coils
        dx_tf_inboard_out_toroidal : float
            Toroidal width of outboard TF coil
        dr_tf_outboard : float
            Radial thickness of outboard TF coil leg
        r_tf_outboard_mid : float
            Major radius of centre of outboard TF coil
        dx_beam_duct : float
            Width of beam duct
        dx_beam_shield : float
            Shielding width on both sides of beam duct

        Returns
        -------
        tuple[float, float]
            Tuple containing (radius_beam_tangency, radius_beam_tangency_max)
        """

        # Have kept the single letter variable names to match the original code and documentation diagram.
        radius_beam_tangency = f_radius_beam_tangency_rmajor * rmajor

        omega = 2.0 * np.pi / n_tf_coils

        a = 0.5e0 * dx_tf_inboard_out_toroidal
        try:
            assert a < np.inf
        except AssertionError:
            logger.exception("a is inf. Kludging to 1e10.")
            a = 1e10

        b = dr_tf_outboard
        try:
            assert b < np.inf
        except AssertionError:
            logger.exception("b is inf. Kludging to 1e10.")
            b = 1e10

        c = dx_beam_duct + 2.0e0 * dx_beam_shield

        d = r_tf_outboard_mid - 0.5e0 * b
        try:
            assert d < np.inf
        except AssertionError:
            logger.exception("d is inf. Kludging to 1e10.")
            d = 1e10

        e = np.sqrt(a**2 + (d + b) ** 2)
        f = np.sqrt(a**2 + d**2)

        theta = omega - np.arctan(a / d)
        phi = theta - np.arcsin(a / e)

        g = np.sqrt(e**2 + f**2 - 2.0e0 * e * f * np.cos(phi))

        if g > c:
            h = np.sqrt(g**2 - c**2)
            alpha = np.arctan(h / c)
            eps = np.arcsin(e * np.sin(phi) / g) - alpha
            radius_beam_tangency_max = f * np.cos(eps) - 0.5e0 * c
        else:
            logger.error(
                f"Max beam tangency radius set =0 temporarily; change dx_beam_duct. {g=} {c=}"
            )
            radius_beam_tangency_max = 0.0e0

        return radius_beam_tangency, radius_beam_tangency_max

    def calculate_vertical_build(self, output: bool):
        """Determines the vertical build of the machine.

        This method calculates various parameters related to the vertical build of the machine,
        such as thicknesses, radii, and areas. Results can be outputted with the `output` flag.

        Parameters
        ----------
        output : bool
            Flag indicating whether to output results
        """
        # Set the X-point heights for the top and bottom of the plasma
        # Assumes top-down plasma symmetry
        build_variables.z_plasma_xpoint_upper = (
            physics_variables.rminor * physics_variables.kappa
        )
        build_variables.z_plasma_xpoint_lower = (
            physics_variables.rminor * physics_variables.kappa
        )

        if output:
            po.oheadr(self.outfile, "Vertical Build")

            po.ovarin(
                self.mfile,
                "Divertor null switch",
                "(i_single_null)",
                physics_variables.i_single_null,
            )

            if physics_variables.i_single_null == 0:
                po.ocmmnt(self.outfile, "Double null case")

                # Start at the top and work down.

                vertical_build_upper = (
                    buildings_variables.dz_tf_cryostat
                    + build_variables.dr_tf_inboard
                    + build_variables.dr_tf_shld_gap
                    + build_variables.dz_shld_thermal
                    + build_variables.dz_shld_vv_gap
                    + build_variables.dz_vv_upper
                    + build_variables.dz_shld_upper
                    + divertor_variables.dz_divertor
                    + build_variables.dz_xpoint_divertor
                    + build_variables.z_plasma_xpoint_upper
                )

                # To calculate vertical offset between TF coil centre and plasma centre
                vbuile1 = vertical_build_upper

                po.obuild(
                    self.outfile,
                    "Cryostat roof structure*",
                    buildings_variables.dz_tf_cryostat,
                    vertical_build_upper,
                    "(dz_tf_cryostat)",
                )
                po.ovarre(
                    self.mfile,
                    "Cryostat roof structure*",
                    "(dz_tf_cryostat)",
                    buildings_variables.dz_tf_cryostat,
                )
                vertical_build_upper = (
                    vertical_build_upper - buildings_variables.dz_tf_cryostat
                )

                # Top of TF coil
                tf_top = vertical_build_upper

                po.obuild(
                    self.outfile,
                    "TF coil",
                    build_variables.dr_tf_inboard,
                    vertical_build_upper,
                    "(dr_tf_inboard)",
                )
                vertical_build_upper = (
                    vertical_build_upper - build_variables.dr_tf_inboard
                )

                po.obuild(
                    self.outfile,
                    "Gap",
                    build_variables.dr_tf_shld_gap,
                    vertical_build_upper,
                    "(dr_tf_shld_gap)",
                )
                vertical_build_upper = (
                    vertical_build_upper - build_variables.dr_tf_shld_gap
                )

                po.obuild(
                    self.outfile,
                    "Thermal shield, vertical",
                    build_variables.dz_shld_thermal,
                    vertical_build_upper,
                    "(dz_shld_thermal)",
                )

                po.ovarre(
                    self.mfile,
                    "Thermal shield, vertical (m)",
                    "(dz_shld_thermal)",
                    build_variables.dz_shld_thermal,
                )
                vertical_build_upper = (
                    vertical_build_upper - build_variables.dz_shld_thermal
                )

                po.obuild(
                    self.outfile,
                    "Gap",
                    build_variables.dz_shld_vv_gap,
                    vertical_build_upper,
                    "(dz_shld_vv_gap)",
                )
                po.ovarre(
                    self.mfile,
                    "Vessel - TF coil vertical gap (m)",
                    "(dz_shld_vv_gap)",
                    build_variables.dz_shld_vv_gap,
                )
                vertical_build_upper = (
                    vertical_build_upper - build_variables.dz_shld_vv_gap
                )

                po.obuild(
                    self.outfile,
                    "Vacuum vessel (and shielding)",
                    build_variables.dz_vv_upper + build_variables.dz_shld_upper,
                    vertical_build_upper,
                    "(dz_vv_upper+dz_shld_upper)",
                )
                vertical_build_upper = (
                    vertical_build_upper
                    - build_variables.dz_vv_upper
                    - build_variables.dz_shld_upper
                )
                po.ovarre(
                    self.mfile,
                    "Topside vacuum vessel radial thickness (m)",
                    "(dz_vv_upper)",
                    build_variables.dz_vv_upper,
                )
                po.ovarre(
                    self.mfile,
                    "Top radiation shield thickness (m)",
                    "(dz_shld_upper)",
                    build_variables.dz_shld_upper,
                )

                po.obuild(
                    self.outfile,
                    "Divertor structure",
                    divertor_variables.dz_divertor,
                    vertical_build_upper,
                    "(dz_divertor)",
                )
                po.ovarre(
                    self.mfile,
                    "Divertor structure vertical thickness (m)",
                    "(dz_divertor)",
                    divertor_variables.dz_divertor,
                )
                vertical_build_upper = (
                    vertical_build_upper - divertor_variables.dz_divertor
                )

                po.obuild(
                    self.outfile,
                    "Top scrape-off",
                    build_variables.dz_xpoint_divertor,
                    vertical_build_upper,
                    "(dz_xpoint_divertor)",
                )
                po.ovarre(
                    self.mfile,
                    "Top scrape-off vertical thickness (m)",
                    "(dz_xpoint_divertor)",
                    build_variables.dz_xpoint_divertor,
                )
                vertical_build_upper = (
                    vertical_build_upper - build_variables.dz_xpoint_divertor
                )

                po.obuild(
                    self.outfile,
                    "Plasma upper X-point height (m)",
                    build_variables.z_plasma_xpoint_upper,
                    vertical_build_upper,
                    "(z_plasma_xpoint_upper)",
                )
                po.ovarre(
                    self.mfile,
                    "Plasma upper X-point height (m)",
                    "(z_plasma_xpoint_upper)",
                    build_variables.z_plasma_xpoint_upper,
                )
                vertical_build_upper = (
                    vertical_build_upper - build_variables.z_plasma_xpoint_upper
                )

                po.obuild(self.outfile, "Midplane", 0.0e0, vertical_build_upper)

                vertical_build_upper = (
                    vertical_build_upper - build_variables.z_plasma_xpoint_lower
                )
                po.obuild(
                    self.outfile,
                    "Plasma lower X-point height (m)",
                    build_variables.z_plasma_xpoint_lower,
                    vertical_build_upper,
                    "(z_plasma_xpoint_lower)",
                )
                po.ovarre(
                    self.mfile,
                    "Plasma lower X-point height (m)",
                    "(z_plasma_xpoint_lower)",
                    build_variables.z_plasma_xpoint_lower,
                )

                vertical_build_upper = (
                    vertical_build_upper - build_variables.dz_xpoint_divertor
                )
                po.obuild(
                    self.outfile,
                    "Lower scrape-off",
                    build_variables.dz_xpoint_divertor,
                    vertical_build_upper,
                    "(dz_xpoint_divertor)",
                )
                po.ovarre(
                    self.mfile,
                    "Bottom scrape-off vertical thickness (m)",
                    "(dz_xpoint_divertor)",
                    build_variables.dz_xpoint_divertor,
                )

                vertical_build_upper = (
                    vertical_build_upper - divertor_variables.dz_divertor
                )
                po.obuild(
                    self.outfile,
                    "Divertor structure",
                    divertor_variables.dz_divertor,
                    vertical_build_upper,
                    "(dz_divertor)",
                )
                po.ovarre(
                    self.mfile,
                    "Divertor structure vertical thickness (m)",
                    "(dz_divertor)",
                    divertor_variables.dz_divertor,
                )

                vertical_build_upper = (
                    vertical_build_upper - build_variables.dz_shld_lower
                )

                vertical_build_upper = vertical_build_upper - build_variables.dz_vv_lower
                po.obuild(
                    self.outfile,
                    "Vacuum vessel (and shielding)",
                    build_variables.dz_vv_lower + build_variables.dz_shld_lower,
                    vertical_build_upper,
                    "(dz_vv_lower+dz_shld_lower)",
                )
                po.ovarre(
                    self.mfile,
                    "Bottom radiation shield thickness (m)",
                    "(dz_shld_lower)",
                    build_variables.dz_shld_lower,
                )
                po.ovarre(
                    self.mfile,
                    "Underside vacuum vessel radial thickness (m)",
                    "(dz_vv_lower)",
                    build_variables.dz_vv_lower,
                )

                vertical_build_upper = (
                    vertical_build_upper - build_variables.dz_shld_vv_gap
                )
                po.obuild(
                    self.outfile,
                    "Gap",
                    build_variables.dz_shld_vv_gap,
                    vertical_build_upper,
                    "(dz_shld_vv_gap)",
                )

                vertical_build_upper = (
                    vertical_build_upper - build_variables.dz_shld_thermal
                )
                po.obuild(
                    self.outfile,
                    "Thermal shield, vertical",
                    build_variables.dz_shld_thermal,
                    vertical_build_upper,
                    "(dz_shld_thermal)",
                )

                vertical_build_upper = (
                    vertical_build_upper - build_variables.dr_tf_shld_gap
                )
                po.obuild(
                    self.outfile,
                    "Gap",
                    build_variables.dr_tf_shld_gap,
                    vertical_build_upper,
                    "(dr_tf_shld_gap)",
                )

                vertical_build_upper = (
                    vertical_build_upper - build_variables.dr_tf_inboard
                )
                po.obuild(
                    self.outfile,
                    "TF coil",
                    build_variables.dr_tf_inboard,
                    vertical_build_upper,
                    "(dr_tf_inboard)",
                )

                # Total height of TF coil
                tf_height = tf_top - vertical_build_upper
                # Inner vertical dimension of TF coil
                build_variables.dh_tf_inner_bore = (
                    tf_height - 2 * build_variables.dr_tf_inboard
                )

                vertical_build_upper = (
                    vertical_build_upper - buildings_variables.dz_tf_cryostat
                )
                po.obuild(
                    self.outfile,
                    "Cryostat floor structure**",
                    buildings_variables.dz_tf_cryostat,
                    vertical_build_upper,
                    "(dz_tf_cryostat)",
                )

                # To calculate vertical offset between TF coil centre and plasma centre
                build_variables.dz_tf_plasma_centre_offset = (
                    vbuile1 + vertical_build_upper
                ) / 2.0e0

                # End of Double null case
            else:
                po.ocmmnt(self.outfile, "Single null case")
                build_variables.dz_vv_upper = 0.5 * (
                    build_variables.dz_vv_upper + build_variables.dz_vv_lower
                )

                build_variables.dz_fw_upper = 0.5 * (
                    build_variables.dr_fw_inboard + build_variables.dr_fw_outboard
                )

                vbuild = (
                    buildings_variables.dz_tf_cryostat
                    + build_variables.dr_tf_inboard
                    + build_variables.dr_tf_shld_gap
                    + build_variables.dz_shld_thermal
                    + build_variables.dz_shld_vv_gap
                    + build_variables.dz_vv_upper
                    + build_variables.dr_shld_blkt_gap
                    + build_variables.dz_shld_upper
                    + build_variables.dz_blkt_upper
                    + build_variables.dz_fw_upper
                    + build_variables.dz_fw_plasma_gap
                    + build_variables.z_plasma_xpoint_upper
                )

                # To calculate vertical offset between TF coil centre and plasma centre
                vbuile1 = vbuild

                po.obuild(
                    self.outfile,
                    "Cryostat roof structure*",
                    buildings_variables.dz_tf_cryostat,
                    vbuild,
                    "(dz_tf_cryostat)",
                )
                po.ovarre(
                    self.mfile,
                    "Cryostat roof structure*",
                    "(dz_tf_cryostat)",
                    buildings_variables.dz_tf_cryostat,
                )
                vbuild = vbuild - buildings_variables.dz_tf_cryostat

                # Top of TF coil
                tf_top = vbuild

                po.obuild(
                    self.outfile,
                    "TF coil",
                    build_variables.dr_tf_inboard,
                    vbuild,
                    "(dr_tf_inboard)",
                )
                vbuild = vbuild - build_variables.dr_tf_inboard

                po.obuild(
                    self.outfile,
                    "Gap",
                    build_variables.dr_tf_shld_gap,
                    vbuild,
                    "(dr_tf_shld_gap)",
                )
                vbuild = vbuild - build_variables.dr_tf_shld_gap

                po.obuild(
                    self.outfile,
                    "Thermal shield, vertical",
                    build_variables.dz_shld_thermal,
                    vbuild,
                    "(dz_shld_thermal)",
                )
                po.ovarre(
                    self.mfile,
                    "Thermal shield, vertical (m)",
                    "(dz_shld_thermal)",
                    build_variables.dz_shld_thermal,
                )
                vbuild = vbuild - build_variables.dz_shld_thermal

                po.obuild(
                    self.outfile,
                    "Gap",
                    build_variables.dz_shld_vv_gap,
                    vbuild,
                    "(dz_shld_vv_gap)",
                )
                po.ovarre(
                    self.mfile,
                    "Vessel - TF coil vertical gap (m)",
                    "(dz_shld_vv_gap)",
                    build_variables.dz_shld_vv_gap,
                )
                vbuild = vbuild - build_variables.dz_shld_vv_gap

                po.obuild(
                    self.outfile,
                    "Vacuum vessel (and shielding)",
                    build_variables.dz_vv_upper + build_variables.dz_shld_upper,
                    vbuild,
                    "(dz_vv_upper+dz_shld_upper)",
                )
                vbuild = (
                    vbuild - build_variables.dz_vv_upper - build_variables.dz_shld_upper
                )
                po.ovarre(
                    self.mfile,
                    "Topside vacuum vessel radial thickness (m)",
                    "(dz_vv_upper)",
                    build_variables.dz_vv_upper,
                )
                po.ovarre(
                    self.mfile,
                    "Top radiation shield thickness (m)",
                    "(dz_shld_upper)",
                    build_variables.dz_shld_upper,
                )

                po.obuild(
                    self.outfile,
                    "Gap",
                    build_variables.dr_shld_blkt_gap,
                    vbuild,
                    "(dr_shld_blkt_gap)",
                )
                vbuild = vbuild - build_variables.dr_shld_blkt_gap

                po.obuild(
                    self.outfile,
                    "Top blanket",
                    build_variables.dz_blkt_upper,
                    vbuild,
                    "(dz_blkt_upper)",
                )
                po.ovarre(
                    self.mfile,
                    "Top blanket vertical thickness (m)",
                    "(dz_blkt_upper)",
                    build_variables.dz_blkt_upper,
                )
                vbuild = vbuild - build_variables.dz_blkt_upper

                dz_fw_upper = 0.5e0 * (
                    build_variables.dr_fw_inboard + build_variables.dr_fw_outboard
                )
                po.obuild(
                    self.outfile, "Top first wall", dz_fw_upper, vbuild, "(dz_fw_upper)"
                )
                po.ovarre(
                    self.mfile,
                    "Top first wall vertical thickness (m)",
                    "(dz_fw_upper)",
                    dz_fw_upper,
                )
                vbuild = vbuild - dz_fw_upper

                po.obuild(
                    self.outfile,
                    "Top scrape-off",
                    build_variables.dz_fw_plasma_gap,
                    vbuild,
                    "(dz_fw_plasma_gap)",
                )
                po.ovarre(
                    self.mfile,
                    "Top scrape-off vertical thickness (m)",
                    "(dz_fw_plasma_gap)",
                    build_variables.dz_fw_plasma_gap,
                )
                vbuild = vbuild - build_variables.dz_fw_plasma_gap

                po.obuild(
                    self.outfile,
                    "Plasma upper X-point height (m)",
                    build_variables.z_plasma_xpoint_upper,
                    vbuild,
                    "(z_plasma_xpoint_upper)",
                )
                po.ovarre(
                    self.mfile,
                    "Plasma upper X-point height (m)",
                    "(z_plasma_xpoint_upper)",
                    build_variables.z_plasma_xpoint_upper,
                )
                vbuild = vbuild - build_variables.z_plasma_xpoint_upper

                po.obuild(self.outfile, "Midplane", 0.0e0, vbuild)

                vbuild = vbuild - build_variables.z_plasma_xpoint_lower
                po.obuild(
                    self.outfile,
                    "Plasma lower X-point height (m)",
                    build_variables.z_plasma_xpoint_lower,
                    vbuild,
                    "(z_plasma_xpoint_lower)",
                )
                po.ovarre(
                    self.mfile,
                    "Plasma lower X-point height (m)",
                    "(z_plasma_xpoint_lower)",
                    build_variables.z_plasma_xpoint_lower,
                )

                vbuild = vbuild - build_variables.dz_xpoint_divertor
                po.obuild(
                    self.outfile,
                    "Lower scrape-off",
                    build_variables.dz_xpoint_divertor,
                    vbuild,
                    "(dz_xpoint_divertor)",
                )
                po.ovarre(
                    self.mfile,
                    "Bottom scrape-off vertical thickness (m)",
                    "(dz_xpoint_divertor)",
                    build_variables.dz_xpoint_divertor,
                )

                vbuild = vbuild - divertor_variables.dz_divertor
                po.obuild(
                    self.outfile,
                    "Divertor structure",
                    divertor_variables.dz_divertor,
                    vbuild,
                    "(dz_divertor)",
                )
                po.ovarre(
                    self.mfile,
                    "Divertor structure vertical thickness (m)",
                    "(dz_divertor)",
                    divertor_variables.dz_divertor,
                )

                vbuild = vbuild - build_variables.dz_shld_lower

                vbuild = vbuild - build_variables.dz_vv_lower
                po.obuild(
                    self.outfile,
                    "Vacuum vessel (and shielding)",
                    build_variables.dz_vv_lower + build_variables.dz_shld_lower,
                    vbuild,
                    "(dz_vv_lower+dz_shld_lower)",
                )
                po.ovarre(
                    self.mfile,
                    "Bottom radiation shield thickness (m)",
                    "(dz_shld_lower)",
                    build_variables.dz_shld_lower,
                )
                po.ovarre(
                    self.mfile,
                    "Underside vacuum vessel radial thickness (m)",
                    "(dz_vv_lower)",
                    build_variables.dz_vv_lower,
                )

                vbuild = vbuild - build_variables.dz_shld_vv_gap
                po.obuild(
                    self.outfile,
                    "Gap",
                    build_variables.dz_shld_vv_gap,
                    vbuild,
                    "(dz_shld_vv_gap)",
                )

                vbuild = vbuild - build_variables.dz_shld_thermal
                po.obuild(
                    self.outfile,
                    "Thermal shield, vertical",
                    build_variables.dz_shld_thermal,
                    vbuild,
                    "(dz_shld_thermal)",
                )

                vbuild = vbuild - build_variables.dr_tf_shld_gap
                po.obuild(
                    self.outfile,
                    "Gap",
                    build_variables.dr_tf_shld_gap,
                    vbuild,
                    "(dr_tf_shld_gap)",
                )

                vbuild = vbuild - build_variables.dr_tf_inboard
                po.obuild(
                    self.outfile,
                    "TF coil",
                    build_variables.dr_tf_inboard,
                    vbuild,
                    "(dr_tf_inboard)",
                )

                # Total height of TF coil
                tf_height = tf_top - vbuild
                # Inner vertical dimension of TF coil
                build_variables.dh_tf_inner_bore = (
                    tf_height - 2 * build_variables.dr_tf_inboard
                )

                vbuild = vbuild - buildings_variables.dz_tf_cryostat

                po.obuild(
                    self.outfile,
                    "Cryostat floor structure**",
                    buildings_variables.dz_tf_cryostat,
                    vbuild,
                    "(dz_tf_cryostat)",
                )

                # To calculate vertical offset between TF coil centre and plasma centre
                build_variables.dz_tf_plasma_centre_offset = (vbuile1 + vbuild) / 2.0e0

                # end of Single null case

            po.ovarre(
                self.mfile,
                "Ratio of Central Solenoid height to TF coil internal height",
                "(f_z_cs_tf_internal)",
                pfcoil_variables.f_z_cs_tf_internal,
            )
            po.ocmmnt(
                self.outfile,
                "\n*Cryostat roof allowance includes uppermost PF coil and outer thermal shield.\n*Cryostat floor allowance includes lowermost PF coil, outer thermal shield and gravity support.",
            )

        # Output the cdivertor geometry
        divht = self.divgeom(output)
        # Issue #481 Remove build_variables.vgaptf
        if build_variables.dz_xpoint_divertor < 0.00001e0:
            build_variables.dz_xpoint_divertor = divht

        # If build_variables.dz_xpoint_divertor /= 0 use the value set by the user.

        # Height to inside edge of TF coil. TF coils are assumed to be symmetrical.
        # Therefore this applies to single and double null cases.
        build_variables.z_tf_inside_half = (
            build_variables.z_plasma_xpoint_upper
            + build_variables.dz_xpoint_divertor
            + divertor_variables.dz_divertor
            + build_variables.dz_shld_lower
            + build_variables.dz_vv_lower
            + build_variables.dz_shld_vv_gap
            + build_variables.dz_shld_thermal
            + build_variables.dr_tf_shld_gap
        )

        #  Vertical locations of divertor coils
        if physics_variables.i_single_null == 0:
            build_variables.z_tf_top = (
                build_variables.z_tf_inside_half + build_variables.dr_tf_inboard
            )
            build_variables.dz_tf_upper_lower_midplane = 0.0e0
        else:
            build_variables.z_tf_top = (
                build_variables.dr_tf_inboard
                + build_variables.dr_tf_shld_gap
                + build_variables.dz_shld_thermal
                + build_variables.dz_shld_vv_gap
                + build_variables.dz_vv_upper
                + build_variables.dz_shld_upper
                + build_variables.dr_shld_blkt_gap
                + build_variables.dz_blkt_upper
                + 0.5e0
                * (build_variables.dr_fw_inboard + build_variables.dr_fw_outboard)
                + build_variables.dz_fw_plasma_gap
                + build_variables.z_plasma_xpoint_upper
            )
            build_variables.dz_tf_upper_lower_midplane = build_variables.z_tf_top - (
                build_variables.z_tf_inside_half + build_variables.dr_tf_inboard
            )

    def divgeom(self, output: bool):
        """Divertor geometry calculation

        This subroutine determines the divertor geometry.
        The inboard (i) and outboard (o) plasma surfaces
        are approximated by arcs, and followed past the X-point to
        determine the maximum height.
        TART option: Peng SOFT paper

        Parameters
        ----------
        output: bool

        Returns
        -------
        divht:
            divertor height (m)
        """
        if physics_variables.itart == 1:
            return 1.75e0 * physics_variables.rminor
        #  Conventional tokamak divertor model
        #  options for seperate upper and lower physics_variables.triangularity

        kap = physics_variables.kappa
        triu = physics_variables.triang
        tril = physics_variables.triang

        # Old method: assumes that divertor arms are continuations of arcs
        #
        # Outboard side
        # build_variables.plsepo = poloidal length along the separatrix from null to
        # strike point on outboard [default 1.5 m]
        # thetao = arc angle between the strike point and the null point
        #
        # xpointo = physics_variables.rmajor + 0.5e0*physics_variables.rminor*(kap**2 + tri**2 - 1.0e0) /     #     (1.0e0 - tri)
        # rprimeo = (xpointo - physics_variables.rmajor + physics_variables.rminor)
        # phio = asin(kap*physics_variables.rminor/rprimeo)
        # thetao = build_variables.plsepo/rprimeo
        #
        # Initial strike point
        #
        # yspointo = rprimeo * sin(thetao + phio)
        # xspointo = xpointo - rprimeo * cos(thetao + phio)
        #
        # Outboard strike point radius - normalized to ITER
        #
        # rstrko = xspointo + 0.14e0
        #
        # Uppermost divertor strike point (end of power decay)
        # anginc = angle of incidence of scrape-off field lines on the
        # divertor (rad)
        #
        # +**PJK 25/07/11 Changed sign of anginc contribution
        # yprimeb = soleno * cos(thetao + phio - anginc)
        #
        # divht = yprimeb + yspointo - kap*physics_variables.rminor

        # New method, assuming straight legs -- superceded by new method 26/5/2016
        # Assumed 90 degrees at X-pt -- wrong#
        #
        #  Find half-angle of outboard arc
        # denomo = (tril**2 + kap**2 - 1.0e0)/( 2.0e0*(1.0e0+tril) ) - tril
        # thetao = atan(kap/denomo)
        # Angle between horizontal and inner divertor leg
        # alphad = (pi/2.0e0) - thetao

        # Method 26/05/2016
        # Find radius of inner and outer plasma arcs

        rco = 0.5 * np.sqrt(
            (physics_variables.rminor**2 * ((tril + 1.0e0) ** 2 + kap**2) ** 2)
            / ((tril + 1.0e0) ** 2)
        )
        rci = 0.5 * np.sqrt(
            (physics_variables.rminor**2 * ((tril - 1.0e0) ** 2 + kap**2) ** 2)
            / ((tril - 1.0e0) ** 2)
        )

        # Find angles between vertical and legs
        # Inboard arc angle = outboard leg angle

        thetao = np.arcsin(1.0e0 - (physics_variables.rminor * (1.0e0 - tril)) / rci)

        # Outboard arc angle = inboard leg angle

        thetai = np.arcsin(1.0e0 - (physics_variables.rminor * (1.0e0 + tril)) / rco)

        #  Position of lower x-pt
        rxpt = physics_variables.rmajor - tril * physics_variables.rminor
        zxpt = -1.0e0 * kap * physics_variables.rminor

        # Position of inner strike point
        # rspi = rxpt - build_variables.plsepi*cos(alphad)
        # zspi = zxpt - build_variables.plsepi*sin(alphad)
        rspi = rxpt - build_variables.plsepi * np.cos(thetai)
        zspi = zxpt - build_variables.plsepi * np.sin(thetai)

        # Position of outer strike point
        # build_variables.rspo = rxpt + build_variables.plsepo*cos((pi/2.0e0)-alphad)
        # zspo = zxpt - build_variables.plsepo*sin((pi/2.0e0)-alphad)
        build_variables.rspo = rxpt + build_variables.plsepo * np.cos(thetao)
        zspo = zxpt - build_variables.plsepo * np.sin(thetao)

        # Position of inner plate ends
        # rplti = rspi - (build_variables.plleni/2.0e0)*sin(divertor_variables.betai + alphad - pi/2.0e0)
        # zplti = zspi + (build_variables.plleni/2.0e0)*cos(divertor_variables.betai + alphad - pi/2.0e0)
        # rplbi = rspi + (build_variables.plleni/2.0e0)*sin(divertor_variables.betai + alphad - pi/2.0e0)
        # zplbi = zspi - (build_variables.plleni/2.0e0)*cos(divertor_variables.betai + alphad - pi/2.0e0)
        rplti = rspi + (build_variables.plleni / 2.0e0) * np.cos(
            thetai + divertor_variables.betai
        )
        zplti = zspi + (build_variables.plleni / 2.0e0) * np.sin(
            thetai + divertor_variables.betai
        )
        rplbi = rspi - (build_variables.plleni / 2.0e0) * np.cos(
            thetai + divertor_variables.betai
        )
        zplbi = zspi - (build_variables.plleni / 2.0e0) * np.sin(
            thetai + divertor_variables.betai
        )

        # Position of outer plate ends
        # rplto = build_variables.rspo + (build_variables.plleno/2.0e0)*sin(divertor_variables.betao - alphad)
        # zplto = zspo + (build_variables.plleno/2.0e0)*cos(divertor_variables.betao - alphad)
        # rplbo = build_variables.rspo - (build_variables.plleno/2.0e0)*sin(divertor_variables.betao - alphad)
        # zplbo = zspo - (build_variables.plleno/2.0e0)*cos(divertor_variables.betao - alphad)
        rplto = build_variables.rspo - (build_variables.plleno / 2.0e0) * np.cos(
            thetao + divertor_variables.betao
        )
        zplto = zspo + (build_variables.plleno / 2.0e0) * np.sin(
            thetao + divertor_variables.betao
        )
        rplbo = build_variables.rspo + (build_variables.plleno / 2.0e0) * np.cos(
            thetao + divertor_variables.betao
        )
        zplbo = zspo - (build_variables.plleno / 2.0e0) * np.sin(
            thetao + divertor_variables.betao
        )

        divht = max(zplti, zplto) - min(zplbo, zplbi)

        if output:
            if divertor_variables.n_divertors == 1:
                po.oheadr(self.outfile, "Divertor build and plasma position")
                po.ocmmnt(self.outfile, "Divertor Configuration = Single Null Divertor")
                po.oblnkl(self.outfile)
                ptop_radial = physics_variables.rmajor - triu * physics_variables.rminor
                ptop_vertical = kap * physics_variables.rminor
                po.ovarrf(
                    self.outfile,
                    "Plasma top position, radial (m)",
                    "(ptop_radial)",
                    ptop_radial,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Plasma top position, vertical (m)",
                    "(ptop_vertical)",
                    ptop_vertical,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Plasma geometric centre, radial (m)",
                    "(rmajor.)",
                    physics_variables.rmajor,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Plasma geometric centre, vertical (m)",
                    "(0.0)",
                    0.0e0,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Plasma lower triangularity",
                    "(tril)",
                    tril,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Plasma elongation",
                    "(kappa.)",
                    kap,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "TF coil vertical offset (m)",
                    "(dz_tf_plasma_centre_offset)",
                    build_variables.dz_tf_plasma_centre_offset,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Plasma outer arc radius of curvature (m)",
                    "(rco)",
                    rco,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Plasma inner arc radius of curvature (m)",
                    "(rci)",
                    rci,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile, "Plasma lower X-pt, radial (m)", "(rxpt)", rxpt, "OP "
                )
                po.ovarrf(
                    self.outfile,
                    "Plasma lower X-pt, vertical (m)",
                    "(zxpt)",
                    zxpt,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Poloidal plane angle between vertical and inner leg (rad)",
                    "(thetai)",
                    thetai,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Poloidal plane angle between vertical and outer leg (rad)",
                    "(thetao)",
                    thetao,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Poloidal plane angle between inner leg and plate (rad)",
                    "(betai)",
                    divertor_variables.betai,
                )
                po.ovarrf(
                    self.outfile,
                    "Poloidal plane angle between outer leg and plate (rad)",
                    "(betao)",
                    divertor_variables.betao,
                )
                po.ovarrf(
                    self.outfile,
                    "Inner divertor leg poloidal length (m)",
                    "(plsepi)",
                    build_variables.plsepi,
                )
                po.ovarrf(
                    self.outfile,
                    "Outer divertor leg poloidal length (m)",
                    "(plsepo)",
                    build_variables.plsepo,
                )
                po.ovarrf(
                    self.outfile,
                    "Inner divertor plate length (m)",
                    "(plleni)",
                    build_variables.plleni,
                )
                po.ovarrf(
                    self.outfile,
                    "Outer divertor plate length (m)",
                    "(plleno)",
                    build_variables.plleno,
                )
                po.ovarrf(
                    self.outfile,
                    "Inner strike point, radial (m)",
                    "(rspi)",
                    rspi,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Inner strike point, vertical (m)",
                    "(zspi)",
                    zspi,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile, "Inner plate top, radial (m)", "(rplti)", rplti, "OP "
                )
                po.ovarrf(
                    self.outfile,
                    "Inner plate top, vertical (m)",
                    "(zplti)",
                    zplti,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Inner plate bottom, radial (m)",
                    "(rplbi)",
                    rplbi,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Inner plate bottom, vertical (m)",
                    "(zplbi)",
                    zplbi,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Outer strike point, radial (m)",
                    "(rspo)",
                    build_variables.rspo,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Outer strike point, vertical (m)",
                    "(zspo)",
                    zspo,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile, "Outer plate top, radial (m)", "(rplto)", rplto, "OP "
                )
                po.ovarrf(
                    self.outfile,
                    "Outer plate top, vertical (m)",
                    "(zplto)",
                    zplto,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Outer plate bottom, radial (m)",
                    "(rplbo)",
                    rplbo,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Outer plate bottom, vertical (m)",
                    "(zplbo)",
                    zplbo,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Calculated maximum divertor height (m)",
                    "(divht)",
                    divht,
                    "OP ",
                )

            elif divertor_variables.n_divertors == 2:
                po.oheadr(self.outfile, "Divertor build and plasma position")
                po.ocmmnt(self.outfile, "Divertor Configuration = Double Null Divertor")
                po.oblnkl(self.outfile)
                # Assume upper and lower divertors geometries are symmetric.
                ptop_radial = physics_variables.rmajor - triu * physics_variables.rminor
                ptop_vertical = kap * physics_variables.rminor
                po.ovarrf(
                    self.outfile,
                    "Plasma top position, radial (m)",
                    "(ptop_radial)",
                    ptop_radial,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Plasma top position, vertical (m)",
                    "(ptop_vertical)",
                    ptop_vertical,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Plasma geometric centre, radial (m)",
                    "(rmajor.)",
                    physics_variables.rmajor,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Plasma geometric centre, vertical (m)",
                    "(0.0)",
                    0.0e0,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Plasma physics_variables.triangularity",
                    "(tril)",
                    tril,
                    "OP ",
                )
                po.ovarrf(self.outfile, "Plasma elongation", "(kappa.)", kap, "OP ")
                po.ovarrf(
                    self.outfile,
                    "TF coil vertical offset (m)",
                    "(dz_tf_plasma_centre_offset)",
                    build_variables.dz_tf_plasma_centre_offset,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile, "Plasma upper X-pt, radial (m)", "(rxpt)", rxpt, "OP "
                )
                po.ovarrf(
                    self.outfile,
                    "Plasma upper X-pt, vertical (m)",
                    "(-zxpt)",
                    -zxpt,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Plasma outer arc radius of curvature (m)",
                    "(rco)",
                    rco,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Plasma inner arc radius of curvature (m)",
                    "(rci)",
                    rci,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile, "Plasma lower X-pt, radial (m)", "(rxpt)", rxpt, "OP "
                )
                po.ovarrf(
                    self.outfile,
                    "Plasma lower X-pt, vertical (m)",
                    "(zxpt)",
                    zxpt,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Poloidal plane angle between vertical and inner leg (rad)",
                    "(thetai)",
                    thetai,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Poloidal plane angle between vertical and outer leg (rad)",
                    "(thetao)",
                    thetao,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Poloidal plane angle between inner leg and plate (rad)",
                    "(betai)",
                    divertor_variables.betai,
                )
                po.ovarrf(
                    self.outfile,
                    "Poloidal plane angle between outer leg and plate (rad)",
                    "(betao)",
                    divertor_variables.betao,
                )
                po.ovarrf(
                    self.outfile,
                    "Inner divertor leg poloidal length (m)",
                    "(plsepi)",
                    build_variables.plsepi,
                )
                po.ovarrf(
                    self.outfile,
                    "Outer divertor leg poloidal length (m)",
                    "(plsepo)",
                    build_variables.plsepo,
                )
                po.ovarrf(
                    self.outfile,
                    "Inner divertor plate length (m)",
                    "(lleni)",
                    build_variables.plleni,
                )
                po.ovarrf(
                    self.outfile,
                    "Outer divertor plate length (m)",
                    "(plleno)",
                    build_variables.plleno,
                )
                po.ovarrf(
                    self.outfile,
                    "Upper inner strike point, radial (m)",
                    "(rspi)",
                    rspi,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Upper inner strike point, vertical (m)",
                    "(-zspi)",
                    -zspi,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Upper inner plate top, radial (m)",
                    "(rplti)",
                    rplti,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Upper inner plate top, vertical (m)",
                    "(-zplti)",
                    -zplti,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Upper inner plate bottom, radial (m)",
                    "(rplbi)",
                    rplbi,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Upper inner plate bottom, vertical (m)",
                    "(-zplbi)",
                    -zplbi,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Upper outer strike point, radial (m)",
                    "(rspo)",
                    build_variables.rspo,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Upper outer strike point, vertical (m)",
                    "(-zspo)",
                    -zspo,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Upper outer plate top, radial (m)",
                    "(rplto)",
                    rplto,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Upper outer plate top, vertical (m)",
                    "(-zplto)",
                    -zplto,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Upper outer plate bottom, radial (m)",
                    "(rplbo)",
                    rplbo,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Upper outer plate bottom, vertical (m)",
                    "(-zplbo)",
                    -zplbo,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Lower inner strike point, radial (m)",
                    "(rspi)",
                    rspi,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Lower inner strike point, vertical (m)",
                    "(zspi)",
                    zspi,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Lower inner plate top, radial (m)",
                    "(rplti)",
                    rplti,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Lower inner plate top, vertical (m)",
                    "(zplti)",
                    zplti,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Lower inner plate bottom, radial (m)",
                    "(rplbi)",
                    rplbi,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Lower inner plate bottom, vertical (m)",
                    "(zplbi)",
                    zplbi,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Lower outer strike point, radial (m)",
                    "(rspo)",
                    build_variables.rspo,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Lower outer strike point, vertical (m)",
                    "(zspo)",
                    zspo,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Lower outer plate top, radial (m)",
                    "(rplto)",
                    rplto,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Lower outer plate top, vertical (m)",
                    "(zplto)",
                    zplto,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Lower outer plate bottom, radial (m)",
                    "(rplbo)",
                    rplbo,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Lower outer plate bottom, vertical (m)",
                    "(zplbo)",
                    zplbo,
                    "OP ",
                )
                po.ovarrf(
                    self.outfile,
                    "Calculated maximum divertor height (m)",
                    "(divht)",
                    divht,
                    "OP ",
                )
            else:
                po.oheadr(self.outfile, "Divertor build and plasma position")
                po.ocmmnt(
                    self.outfile,
                    "ERROR: null value not supported, check i_single_null value.",
                )
        return divht

    def plasma_outboard_edge_toroidal_ripple(
        self,
        ripple_b_tf_plasma_edge_max: float,
        r_tf_outboard_mid: float,
        n_tf_coils: int,
        rmajor: float,
        rminor: float,
        r_tf_wp_inboard_inner,
        r_tf_wp_inboard_centre: float,
        r_tf_wp_inboard_outer: float,
        dx_tf_wp_primary_toroidal: float,
        i_tf_shape: int,
        i_tf_sup: int,
        dx_tf_wp_insulation: float,
        dx_tf_wp_insertion_gap: float,
    ) -> tuple[float, float, int]:
        """Plasma outboard toroidal field (TF) ripple calculation.

        This routine computes the TF ripple amplitude at the midplane outboard
        plasma edge and the minimum radius of the TF coil centre that would
        produce a specified maximum allowed ripple. The calculation uses
        fitted coefficients derived from numerical modelling (MAGINT) and
        includes a simplified analytical picture-frame coil model for
        i_tf_shape == 2.

        Parameters
        ----------
        ripple_b_tf_plasma_edge_max : float
            Maximum allowed ripple at plasma edge (percent)
        r_tf_outboard_mid : float
            Radius to the centre of the outboard TF coil leg (m)
        n_tf_coils : int
            Number of TF coils
        rmajor : float
            Plasma major radius (m)
        rminor : float
            Plasma minor radius (m)
        r_tf_wp_inboard_inner : float
            Inner winding-pack inboard radius (m)
        r_tf_wp_inboard_centre : float
            Centre winding-pack inboard radius (m)
        r_tf_wp_inboard_outer : float
            Outer winding-pack inboard radius (m)
        dx_tf_wp_primary_toroidal : float
            Primary toroidal winding-pack thickness (m)
        i_tf_shape : int
            TF coil shape switch (2 => picture-frame analytical model)
        i_tf_sup : int
            TF coil support flag (1 => superconducting)
        dx_tf_wp_insulation : float
            Winding-pack insulation thickness (m)
        dx_tf_wp_insertion_gap : float
            Winding-pack insertion gap (m)

        Returns
        -------
        tuple[float, float, int]
            Tuple containing:
            - ripple: Calculated ripple at plasma edge (percent)
            - r_tf_outboard_midmin: Minimum r_tf_outboard_mid that yields the specified maximum ripple (m)
            - flag: Applicability flag (0 = OK, non-zero = fitted-range concern)

        Notes
        -----
        - Fitted coefficients originate from parametric MAGINT runs (M. Kovari, 2014).
        - Picture-frame coil analytical model (Ken McClements, 2022) is used when
        `i_tf_shape == 2` and gives approximate results (within ~10% of numerical).
        - The routine sets an applicability flag when fitted-range assumptions are exceeded.
        """
        if i_tf_sup == 1:
            # Minimal inboard WP radius [m]
            r_wp_min = r_tf_wp_inboard_inner

            # Rectangular WP
            if tfcoil_variables.i_tf_wp_geom == 0:
                r_wp_max = r_wp_min

            # Double rectangle WP
            elif tfcoil_variables.i_tf_wp_geom == 1:
                r_wp_max = r_tf_wp_inboard_centre

            # Trapezoidal WP
            elif tfcoil_variables.i_tf_wp_geom == 2:
                r_wp_max = r_tf_wp_inboard_outer

            # Calculated maximum toroidal WP toroidal thickness [m]
            dx_tf_wp_conductor_max = dx_tf_wp_primary_toroidal - 2.0 * (
                dx_tf_wp_insulation + dx_tf_wp_insertion_gap
            )

        # Resistive magnet case
        else:
            # Radius used to define the dx_tf_wp_conductor_max [m]
            r_wp_max = r_tf_wp_inboard_outer
            # Calculated maximum toroidal WP toroidal thickness [m]
            dx_tf_wp_conductor_max = 2.0e0 * r_wp_max * np.tan(np.pi / n_tf_coils)

        flag = 0
        if i_tf_shape == 2:
            # Ken McClements ST picture frame coil analytical ripple calc
            # Calculated ripple for coil at r_tf_outboard_mid (%)
            ripple = 100.0e0 * ((rmajor + rminor) / r_tf_outboard_mid) ** (n_tf_coils)
            #  Calculated r_tf_outboard_mid to produce a ripple of amplitude ripple_b_tf_plasma_edge_max
            r_tf_outboard_midmin = (rmajor + rminor) / (
                (0.01e0 * ripple_b_tf_plasma_edge_max) ** (1.0e0 / n_tf_coils)
            )
        else:
            # Winding pack to iter-coil at plasma centre toroidal lenth ratio
            x = dx_tf_wp_conductor_max * n_tf_coils / rmajor

            # Fitting parameters
            c1 = 0.875e0 - 0.0557e0 * x
            c2 = 1.617e0 + 0.0832e0 * x

            #  Calculated ripple for coil at r_tf_outboard_mid (%)
            ripple = (
                100.0e0
                * c1
                * ((rmajor + rminor) / r_tf_outboard_mid) ** (n_tf_coils - c2)
            )

            #  Calculated r_tf_outboard_mid to produce a ripple of amplitude ripple_b_tf_plasma_edge_max
            base = 0.01 * ripple_b_tf_plasma_edge_max / c1
            # Avoid potential negative or complex result: kludge base to be
            # small and positive if required
            try:
                assert base > 1e-6
            except AssertionError:
                logger.exception("base is <= 1e-6. Kludging to 1e-6.")
                base = 1e-6

            r_tf_outboard_midmin = (rmajor + rminor) / (
                base ** (1.0 / (n_tf_coils - c2))
            )

            try:
                assert r_tf_outboard_midmin < np.inf
            except AssertionError:
                logger.exception(
                    "r_tf_outboard_midmin is inf. Kludging to a large value instead."
                )
                r_tf_outboard_midmin = (rmajor + rminor) * 3

            #  Notify via flag if a range of applicability is violated
            flag = 0
            if (x < 0.737e0) or (x > 2.95e0):
                flag = 1
            if (n_tf_coils < 16) or (n_tf_coils > 20):
                flag = 2
            if ((rmajor + rminor) / r_tf_outboard_mid < 0.7e0) or (
                (rmajor + rminor) / r_tf_outboard_mid > 0.8e0
            ):
                flag = 3

        return ripple, r_tf_outboard_midmin, flag

    def calculate_radial_build(self, output: bool):
        """This method determines the radial build of the machine.
        It calculates various parameters related to the build of the machine,
        such as thicknesses, radii, and areas.
        Results can be outputted with the `output` flag.

        Parameters
        ----------
        output : bool
            Flag indicating whether to output the results
        """

        if fwbs_variables.blktmodel > 0:
            build_variables.dr_blkt_inboard = (
                build_variables.blbuith
                + build_variables.blbmith
                + build_variables.blbpith
            )
            build_variables.dr_blkt_outboard = (
                build_variables.blbuoth
                + build_variables.blbmoth
                + build_variables.blbpoth
            )
            build_variables.dz_shld_upper = 0.5e0 * (
                build_variables.dr_shld_inboard + build_variables.dr_shld_outboard
            )

        #  Top/bottom blanket thickness
        build_variables.dz_blkt_upper = 0.5e0 * (
            build_variables.dr_blkt_inboard + build_variables.dr_blkt_outboard
        )

        if physics_variables.i_single_null == 1:
            #  Check if build_variables.dz_fw_plasma_gap has been set too small
            build_variables.dz_fw_plasma_gap = max(
                0.5e0
                * (
                    build_variables.dr_fw_plasma_gap_inboard
                    + build_variables.dr_fw_plasma_gap_outboard
                ),
                build_variables.dz_fw_plasma_gap,
            )

        # Calculate pre-compression structure thickness is build_variables.i_cs_precomp=1
        if build_variables.i_cs_precomp == 1 and build_variables.i_tf_inside_cs == 0:
            build_variables.dr_cs_precomp = build_variables.fseppc / (
                2.0e0
                * np.pi
                * build_variables.fcspc
                * build_variables.sigallpc
                * (
                    build_variables.dr_bore
                    + build_variables.dr_bore
                    + build_variables.dr_cs
                )
            )
        elif build_variables.i_cs_precomp == 1 and build_variables.i_tf_inside_cs == 1:
            build_variables.dr_cs_precomp = build_variables.fseppc / (
                2.0e0
                * np.pi
                * build_variables.fcspc
                * build_variables.sigallpc
                * (
                    2.0 * build_variables.dr_bore
                    + 2.0 * build_variables.dr_tf_inboard
                    + 2.0 * build_variables.dr_cs_tf_gap
                    + build_variables.dr_cs
                )
            )
        else:
            build_variables.dr_cs_precomp = 0.0e0

        # Issue #514 Radial dimensions of inboard leg
        # Calculate build_variables.dr_tf_inboard if tfcoil_variables.dr_tf_wp_with_insulation is an iteration variable (140)
        if 140 in numerics.ixc[0 : numerics.nvar]:
            build_variables.dr_tf_inboard = (
                tfcoil_variables.dr_tf_wp_with_insulation
                + tfcoil_variables.dr_tf_plasma_case
                + tfcoil_variables.dr_tf_nose_case
            )

        if build_variables.i_tf_inside_cs == 1:
            build_variables.r_tf_inboard_in = (
                build_variables.dr_bore
                # NOTE: dr_bore is just the hollow space, the
                # true dr_bore size used for flux calculations
                # is dr_bore + dr_tf_inboard + dr_cs_tf_gap
            )
        else:
            # Inboard side inner radius [m]
            build_variables.r_tf_inboard_in = (
                build_variables.dr_bore
                + build_variables.dr_cs
                + build_variables.dr_cs_precomp
                + build_variables.dr_cs_tf_gap
            )

        # Radial build to tfcoil middle [m]
        build_variables.r_tf_inboard_mid = (
            build_variables.r_tf_inboard_in + 0.5e0 * build_variables.dr_tf_inboard
        )

        # Radial build to tfcoil plasma facing side [m]
        build_variables.r_tf_inboard_out = (
            build_variables.r_tf_inboard_in + build_variables.dr_tf_inboard
        )

        # WP radial thickness [m]
        # Calculated only if not used as an iteration variable
        if 140 not in numerics.ixc[0 : numerics.nvar]:
            tfcoil_variables.dr_tf_wp_with_insulation = (
                build_variables.dr_tf_inboard
                - tfcoil_variables.dr_tf_plasma_case
                - tfcoil_variables.dr_tf_nose_case
            )

        # Radius of the centrepost at the top of the machine
        if physics_variables.itart == 1 and tfcoil_variables.i_tf_sup != 1:
            # build_variables.r_cp_top is set using the plasma shape
            if build_variables.i_r_cp_top == 0:
                build_variables.r_cp_top = (
                    physics_variables.rmajor
                    - physics_variables.rminor * physics_variables.triang
                    - (
                        build_variables.dr_tf_shld_gap
                        + build_variables.dr_shld_thermal_inboard
                        + build_variables.dr_shld_inboard
                        + build_variables.dr_shld_blkt_gap
                        + build_variables.dr_blkt_inboard
                        + build_variables.dr_fw_inboard
                        + 3.0e0 * build_variables.dr_fw_plasma_gap_inboard
                    )
                    + tfcoil_variables.drtop
                )

                # Notify user that build_variables.r_cp_top has been set to 1.01*build_variables.r_tf_inboard_out (lvl 2 error)
                if build_variables.r_cp_top < 1.01e0 * build_variables.r_tf_inboard_out:
                    logger.error(
                        "TF CP top radius (r_cp_top) replaced by 1.01*r_tf_inboard_out -> potential top rbuild issue"
                        f"{build_variables.r_cp_top=} {build_variables.r_tf_inboard_out=}"
                    )

                    # build_variables.r_cp_top correction
                    build_variables.r_cp_top = build_variables.r_tf_inboard_out * 1.01e0

                # Top and mid-plane TF coil CP radius ratio
                build_variables.f_r_cp = (
                    build_variables.r_cp_top / build_variables.r_tf_inboard_out
                )

            # User defined build_variables.r_cp_top
            elif build_variables.i_r_cp_top == 1:
                # Notify user that build_variables.r_cp_top has been set to 1.01*build_variables.r_tf_inboard_out (lvl 2 error)
                if build_variables.r_cp_top < 1.01e0 * build_variables.r_tf_inboard_out:
                    logger.error(
                        "TF CP top radius (r_cp_top) replaced by 1.01*r_tf_inboard_out -> potential top rbuild issue"
                        f"{build_variables.r_cp_top=} {build_variables.r_tf_inboard_out=}"
                    )

                    # build_variables.r_cp_top correction
                    build_variables.r_cp_top = build_variables.r_tf_inboard_out * 1.01e0

                # Top / mid-plane TF CP radius ratio
                build_variables.f_r_cp = (
                    build_variables.r_cp_top / build_variables.r_tf_inboard_out
                )

            # build_variables.r_cp_top set as a fraction of the outer TF midplane radius
            elif build_variables.i_r_cp_top == 2:
                build_variables.r_cp_top = (
                    build_variables.f_r_cp * build_variables.r_tf_inboard_out
                )

        else:  # End of physics_variables.itart == 1 .and. tfcoil_variables.i_tf_sup /= 1
            build_variables.r_cp_top = build_variables.r_tf_inboard_out

        if build_variables.i_r_cp_top != 0 and (
            build_variables.r_cp_top
            > physics_variables.rmajor
            - physics_variables.rminor * physics_variables.triang
            - (
                build_variables.dr_tf_shld_gap
                + build_variables.dr_shld_thermal_inboard
                + build_variables.dr_shld_inboard
                + build_variables.dr_shld_blkt_gap
                + build_variables.dr_blkt_inboard
                + build_variables.dr_fw_inboard
                + 3.0e0 * build_variables.dr_fw_plasma_gap_inboard
            )
            + tfcoil_variables.drtop
        ):
            logger.error(
                f"Top CP radius larger that its value determined with plasma shape {build_variables.r_cp_top=}"
            )
        if build_variables.i_tf_inside_cs == 1:
            #  Radial position of vacuum vessel [m]
            build_variables.r_vv_inboard_out = (
                build_variables.r_tf_inboard_out
                + build_variables.dr_cs
                + build_variables.dr_cs_tf_gap
                + build_variables.dr_cs_precomp
                + build_variables.dr_tf_shld_gap
                + build_variables.dr_shld_thermal_inboard
                + build_variables.dr_shld_vv_gap_inboard
                + build_variables.dr_vv_inboard
            )
        else:
            build_variables.r_vv_inboard_out = (
                build_variables.r_tf_inboard_out
                + build_variables.dr_tf_shld_gap
                + build_variables.dr_shld_thermal_inboard
                + build_variables.dr_shld_vv_gap_inboard
                + build_variables.dr_vv_inboard
            )
        # Radial position of the inner side of inboard neutronic shield [m]
        build_variables.r_sh_inboard_in = build_variables.r_vv_inboard_out

        # Radial position of the plasma facing side of inboard neutronic shield [m]
        build_variables.r_sh_inboard_out = (
            build_variables.r_sh_inboard_in + build_variables.dr_shld_inboard
        )

        #  Radial build to centre of plasma (should be equal to physics_variables.rmajor)
        build_variables.rbld = (
            build_variables.r_sh_inboard_out
            + build_variables.dr_shld_blkt_gap
            + build_variables.dr_blkt_inboard
            + build_variables.dr_fw_inboard
            + build_variables.dr_fw_plasma_gap_inboard
            + physics_variables.rminor
        )

        #  Radius to inner edge of inboard shield
        build_variables.r_shld_inboard_inner = (
            physics_variables.rmajor
            - physics_variables.rminor
            - build_variables.dr_fw_plasma_gap_inboard
            - build_variables.dr_fw_inboard
            - build_variables.dr_blkt_inboard
            - build_variables.dr_shld_inboard
        )

        #  Radius to outer edge of outboard shield
        build_variables.r_shld_outboard_outer = (
            physics_variables.rmajor
            + physics_variables.rminor
            + build_variables.dr_fw_plasma_gap_outboard
            + build_variables.dr_fw_outboard
            + build_variables.dr_blkt_outboard
            + build_variables.dr_shld_outboard
        )

        #  Thickness of outboard TF coil legs
        if tfcoil_variables.i_tf_sup != 1:
            build_variables.dr_tf_outboard = (
                build_variables.f_dr_tf_outboard_inboard * build_variables.dr_tf_inboard
            )
        else:
            build_variables.dr_tf_outboard = build_variables.dr_tf_inboard

        #  Radius to centre of outboard TF coil legs
        build_variables.r_tf_outboard_mid = (
            build_variables.r_shld_outboard_outer
            + build_variables.dr_shld_blkt_gap
            + build_variables.dr_vv_outboard
            + build_variables.gapomin
            + build_variables.dr_shld_thermal_outboard
            + build_variables.dr_tf_shld_gap
            + 0.5e0 * build_variables.dr_tf_outboard
        )

        # TF coil horizontal build_variables.dr_bore [m]
        build_variables.dr_tf_inner_bore = (
            build_variables.r_tf_outboard_mid - 0.5e0 * build_variables.dr_tf_outboard
        ) - (build_variables.r_tf_inboard_mid - 0.5e0 * build_variables.dr_tf_inboard)

        (
            tfcoil_variables.ripple_b_tf_plasma_edge,
            r_tf_outboard_midl,
            build_variables.ripflag,
        ) = self.plasma_outboard_edge_toroidal_ripple(
            ripple_b_tf_plasma_edge_max=tfcoil_variables.ripple_b_tf_plasma_edge_max,
            r_tf_outboard_mid=build_variables.r_tf_outboard_mid,
            n_tf_coils=tfcoil_variables.n_tf_coils,
            rmajor=physics_variables.rmajor,
            rminor=physics_variables.rminor,
            r_tf_wp_inboard_inner=superconducting_tf_coil_variables.r_tf_wp_inboard_inner,
            r_tf_wp_inboard_centre=superconducting_tf_coil_variables.r_tf_wp_inboard_centre,
            r_tf_wp_inboard_outer=superconducting_tf_coil_variables.r_tf_wp_inboard_outer,
            dx_tf_wp_primary_toroidal=tfcoil_variables.dx_tf_wp_primary_toroidal,
            i_tf_shape=tfcoil_variables.i_tf_shape,
            i_tf_sup=tfcoil_variables.i_tf_sup,
            dx_tf_wp_insulation=tfcoil_variables.dx_tf_wp_insulation,
            dx_tf_wp_insertion_gap=tfcoil_variables.dx_tf_wp_insertion_gap,
        )

        #  If the tfcoil_variables.ripple is too large then move the outboard TF coil leg
        if r_tf_outboard_midl > build_variables.r_tf_outboard_mid:
            build_variables.r_tf_outboard_mid = r_tf_outboard_midl
            build_variables.dr_shld_vv_gap_outboard = (
                build_variables.r_tf_outboard_mid
                - 0.5e0 * build_variables.dr_tf_outboard
                - build_variables.dr_vv_outboard
                - build_variables.r_shld_outboard_outer
                - build_variables.dr_shld_thermal_outboard
                - build_variables.dr_tf_shld_gap
                - build_variables.dr_shld_blkt_gap
            )
            build_variables.dr_tf_inner_bore = (
                build_variables.r_tf_outboard_mid
                - 0.5e0 * build_variables.dr_tf_outboard
            ) - (
                build_variables.r_tf_inboard_mid - 0.5e0 * build_variables.dr_tf_inboard
            )
        else:
            build_variables.dr_shld_vv_gap_outboard = build_variables.gapomin

        (
            tfcoil_variables.ripple_b_tf_plasma_edge,
            r_tf_outboard_midl,
            build_variables.ripflag,
        ) = self.plasma_outboard_edge_toroidal_ripple(
            ripple_b_tf_plasma_edge_max=tfcoil_variables.ripple_b_tf_plasma_edge_max,
            r_tf_outboard_mid=build_variables.r_tf_outboard_mid,
            n_tf_coils=tfcoil_variables.n_tf_coils,
            rmajor=physics_variables.rmajor,
            rminor=physics_variables.rminor,
            r_tf_wp_inboard_inner=superconducting_tf_coil_variables.r_tf_wp_inboard_inner,
            r_tf_wp_inboard_centre=superconducting_tf_coil_variables.r_tf_wp_inboard_centre,
            r_tf_wp_inboard_outer=superconducting_tf_coil_variables.r_tf_wp_inboard_outer,
            dx_tf_wp_primary_toroidal=tfcoil_variables.dx_tf_wp_primary_toroidal,
            i_tf_shape=tfcoil_variables.i_tf_shape,
            i_tf_sup=tfcoil_variables.i_tf_sup,
            dx_tf_wp_insulation=tfcoil_variables.dx_tf_wp_insulation,
            dx_tf_wp_insertion_gap=tfcoil_variables.dx_tf_wp_insertion_gap,
        )

        #

        if output:
            #  Print out device build

            po.oheadr(self.outfile, "Radial Build")

            if build_variables.ripflag != 0:
                po.ocmmnt(
                    self.outfile,
                    "(Ripple result may not be accurate, as the fit was outside",
                )
                po.ocmmnt(self.outfile, " its range of applicability.)")
                po.oblnkl(self.outfile)
                logger.warning(
                    "Ripple result may be inaccurate, as the fit has been extrapolated"
                )

                if build_variables.ripflag == 1:
                    diagnostic = (
                        tfcoil_variables.dx_tf_wp_primary_toroidal
                        * tfcoil_variables.n_tf_coils
                        / physics_variables.rmajor
                    )
                    logger.warning(
                        f"(TF coil ripple calculation) Dimensionless coil width X out of fitted range. {diagnostic=}"
                    )
                elif build_variables.ripflag == 2:
                    logger.warning(
                        f"(TF coil ripple calculation) No of TF coils not between 16 and 20 inclusive {tfcoil_variables.n_tf_coils=}"
                    )
                else:
                    diagnostic = (
                        physics_variables.rmajor + physics_variables.rminor
                    ) / build_variables.r_tf_outboard_mid

                    logger.warning(
                        f"(TF coil ripple calculation) (R+a)/rtot={diagnostic} out of fitted range."
                    )

            po.ovarin(
                self.outfile,
                "TF coil radial placement switch",
                "(i_tf_inside_cs)",
                build_variables.i_tf_inside_cs,
            )
            po.ovarrf(
                self.outfile,
                "Inboard build thickness (m)",
                "(dr_inboard_build)",
                physics_variables.rmajor - physics_variables.rminor,
                "OP ",
            )

            if build_variables.i_tf_inside_cs == 1:
                po.ocmmnt(
                    self.outfile,
                    (
                        "\n (The stated machine dr_bore size is just for the hollow space, "
                    ),
                )
                po.ocmmnt(
                    self.outfile,
                    (
                        "the true dr_bore size used for calculations is dr_bore + dr_tf_inboard + dr_cs_tf_gap)\n"
                    ),
                )
            if (
                build_variables.i_tf_inside_cs == 1
                and tfcoil_variables.i_tf_bucking >= 2
            ):
                po.ocmmnt(
                    self.outfile,
                    "(Bore hollow space has been filled with a solid metal cyclinder to act as wedge support)\n",
                )

            # an array that holds the following information
            # description, variable name, thickness, radius
            radial_build_data = []

            radius = 0.0e0
            radial_build_data.append(["Device centreline", None, 0.0, radius])
            if (
                build_variables.i_tf_inside_cs == 1
                and tfcoil_variables.i_tf_bucking >= 2
            ):
                radius = radius + build_variables.dr_bore

                radial_build_data.append([
                    "Machine dr_bore wedge support cylinder",
                    "dr_bore",
                    build_variables.dr_bore,
                    radius,
                ])
            elif (
                build_variables.i_tf_inside_cs == 1 and tfcoil_variables.i_tf_bucking < 2
            ):
                radius = radius + build_variables.dr_bore
                radial_build_data.append([
                    "Machine dr_bore hole",
                    "dr_bore",
                    build_variables.dr_bore,
                    radius,
                ])
            else:
                radius = radius + build_variables.dr_bore
                radial_build_data.append([
                    "Machine dr_bore",
                    "dr_bore",
                    build_variables.dr_bore,
                    radius,
                ])
            if build_variables.i_tf_inside_cs == 1:
                radius += build_variables.dr_tf_inboard
                radial_build_data.append([
                    "TF coil inboard leg (in dr_bore)",
                    "dr_tf_inboard",
                    build_variables.dr_tf_inboard,
                    radius,
                ])

                radius += build_variables.dr_cs_tf_gap
                radial_build_data.append([
                    "CS precompresion to TF coil radial gap",
                    "dr_cs_tf_gap",
                    build_variables.dr_cs_tf_gap,
                    radius,
                ])

            radius = radius + build_variables.dr_cs
            radial_build_data.append([
                "Central solenoid",
                "dr_cs",
                build_variables.dr_cs,
                radius,
            ])

            radius = radius + build_variables.dr_cs_precomp
            radial_build_data.append([
                "CS precompression",
                "dr_cs_precomp",
                build_variables.dr_cs_precomp,
                radius,
            ])
            if build_variables.i_tf_inside_cs == 0:
                radius = radius + build_variables.dr_cs_tf_gap
                radial_build_data.append([
                    "CS precompresion to TF coil radial gap",
                    "dr_cs_tf_gap",
                    build_variables.dr_cs_tf_gap,
                    radius,
                ])

                radius = radius + build_variables.dr_tf_inboard
                radial_build_data.append([
                    "TF coil inboard leg",
                    "dr_tf_inboard",
                    build_variables.dr_tf_inboard,
                    radius,
                ])

            radius = radius + build_variables.dr_tf_shld_gap
            radial_build_data.append([
                "TF coil inboard leg insulation gap",
                "dr_tf_shld_gap",
                build_variables.dr_tf_shld_gap,
                radius,
            ])

            radius = radius + build_variables.dr_shld_thermal_inboard
            radial_build_data.append([
                "Thermal shield, inboard",
                "dr_shld_thermal_inboard",
                build_variables.dr_shld_thermal_inboard,
                radius,
            ])

            radius = radius + build_variables.dr_shld_vv_gap_inboard
            radial_build_data.append([
                "Thermal shield to vessel radial gap",
                "dr_shld_vv_gap_inboard",
                build_variables.dr_shld_vv_gap_inboard,
                radius,
            ])

            radius += build_variables.dr_vv_inboard
            radial_build_data.append([
                "Inboard vacuum vessel",
                "dr_vv_inboard",
                build_variables.dr_vv_inboard,
                radius,
            ])

            radius += build_variables.dr_shld_inboard
            radial_build_data.append([
                "Inner radiation shield",
                "dr_shld_inboard",
                build_variables.dr_shld_inboard,
                radius,
            ])

            radius = radius + build_variables.dr_shld_blkt_gap
            radial_build_data.append([
                "Gap",
                "dr_shld_blkt_gap",
                build_variables.dr_shld_blkt_gap,
                radius,
            ])

            radius = radius + build_variables.dr_blkt_inboard
            radial_build_data.append([
                "Inboard blanket",
                "dr_blkt_inboard",
                build_variables.dr_blkt_inboard,
                radius,
            ])

            radius = radius + build_variables.dr_fw_inboard
            radial_build_data.append([
                "Inboard first wall",
                "dr_fw_inboard",
                build_variables.dr_fw_inboard,
                radius,
            ])

            radius = radius + build_variables.dr_fw_plasma_gap_inboard
            radial_build_data.append([
                "Inboard scrape-off",
                "dr_fw_plasma_gap_inboard",
                build_variables.dr_fw_plasma_gap_inboard,
                radius,
            ])

            radius = radius + physics_variables.rminor
            radial_build_data.append([
                "Plasma geometric centre",
                "rminor",
                physics_variables.rminor,
                radius,
            ])

            radius = radius + physics_variables.rminor
            radial_build_data.append([
                "Plasma outboard edge",
                "rminor",
                physics_variables.rminor,
                radius,
            ])

            radius = radius + build_variables.dr_fw_plasma_gap_outboard
            radial_build_data.append([
                "Outboard scrape-off",
                "dr_fw_plasma_gap_outboard",
                build_variables.dr_fw_plasma_gap_outboard,
                radius,
            ])

            radius = radius + build_variables.dr_fw_outboard
            radial_build_data.append([
                "Outboard first wall",
                "dr_fw_outboard",
                build_variables.dr_fw_outboard,
                radius,
            ])

            radius = radius + build_variables.dr_blkt_outboard
            radial_build_data.append([
                "Outboard blanket",
                "dr_blkt_outboard",
                build_variables.dr_blkt_outboard,
                radius,
            ])

            radius = radius + build_variables.dr_shld_blkt_gap
            radial_build_data.append([
                "Gap",
                "dr_shld_blkt_gap",
                build_variables.dr_shld_blkt_gap,
                radius,
            ])

            radius += build_variables.dr_shld_outboard
            radial_build_data.append([
                "Outer radiation shield",
                "dr_shld_outboard",
                build_variables.dr_shld_outboard,
                radius,
            ])

            radius += build_variables.dr_vv_outboard
            radial_build_data.append([
                "Outboard vacuum vessel",
                "dr_vv_outboard",
                build_variables.dr_vv_outboard,
                radius,
            ])

            radius = radius + build_variables.dr_shld_vv_gap_outboard
            radial_build_data.append([
                "Vessel to TF gap",
                "dr_shld_vv_gap_outboard",
                build_variables.dr_shld_vv_gap_outboard,
                radius,
            ])

            radius = radius + build_variables.dr_shld_thermal_outboard
            radial_build_data.append([
                "Ouboard thermal shield",
                "dr_shld_thermal_outboard",
                build_variables.dr_shld_thermal_outboard,
                radius,
            ])

            radius = radius + build_variables.dr_tf_shld_gap
            radial_build_data.append([
                "Gap",
                "dr_tf_shld_gap",
                build_variables.dr_tf_shld_gap,
                radius,
            ])

            radius = radius + build_variables.dr_tf_outboard
            radial_build_data.append([
                "TF coil outboard leg",
                "dr_tf_outboard",
                build_variables.dr_tf_outboard,
                radius,
            ])

            for description, variable, thickness, radius in radial_build_data:
                po.obuild(
                    self.outfile,
                    description,
                    thickness,
                    radius,
                    f"({variable})" if variable else "",
                )

            # use manual index to ensure count is contiguous in the event
            # of a `None` variable component
            index = 0
            for description, variable, thickness, radius in radial_build_data:
                if variable is None:
                    continue

                index += 1

                po.ovarre(
                    self.mfile,
                    f"{description} radial thickness (m)",
                    f"({variable})",
                    thickness,
                )

                po.ovarst(
                    self.mfile,
                    f"Radial build component {index}",
                    f"(radial_label({index}))",
                    f'"{variable}"',
                )
                po.ovarre(
                    self.mfile,
                    f"Radial build cumulative radius {index}",
                    f"(radial_cum({index}))",
                    radius,
                )

            if (current_drive_variables.i_hcd_primary in [5, 8]) or (
                current_drive_variables.i_hcd_secondary in [5, 8]
            ):
                po.ovarre(
                    self.mfile,
                    "Width of neutral beam duct where it passes between the TF coils (m)",
                    "(dx_beam_duct)",
                    current_drive_variables.dx_beam_duct,
                )

outfile = constants.NOUT instance-attribute

mfile = constants.MFILE instance-attribute

run()

Source code in process/models/build.py

def run(self):
    self.calculate_radial_build(output=False)
    self.calculate_vertical_build(output=False)

    (
        current_drive_variables.radius_beam_tangency,
        current_drive_variables.radius_beam_tangency_max,
    ) = self.calculate_beam_port_size(
        f_radius_beam_tangency_rmajor=current_drive_variables.f_radius_beam_tangency_rmajor,
        rmajor=physics_variables.rmajor,
        n_tf_coils=tfcoil_variables.n_tf_coils,
        dx_tf_inboard_out_toroidal=tfcoil_variables.dx_tf_inboard_out_toroidal,
        dr_tf_outboard=build_variables.dr_tf_outboard,
        r_tf_outboard_mid=build_variables.r_tf_outboard_mid,
        dx_beam_duct=current_drive_variables.dx_beam_duct,
        dx_beam_shield=current_drive_variables.dx_beam_shield,
    )

calculate_beam_port_size(f_radius_beam_tangency_rmajor, rmajor, n_tf_coils, dx_tf_inboard_out_toroidal, dr_tf_outboard, r_tf_outboard_mid, dx_beam_duct, dx_beam_shield)

Calculates the maximum possible tangency radius for adequate beam access.

Parameters:

Name	Type	Description	Default
f_radius_beam_tangency_rmajor	float	Fraction of rmajor for beam tangency	required
rmajor	float	Major radius	required
n_tf_coils	int	Number of TF coils	required
dx_tf_inboard_out_toroidal	float	Toroidal width of outboard TF coil	required
dr_tf_outboard	float	Radial thickness of outboard TF coil leg	required
r_tf_outboard_mid	float	Major radius of centre of outboard TF coil	required
dx_beam_duct	float	Width of beam duct	required
dx_beam_shield	float	Shielding width on both sides of beam duct	required

Returns:

Type	Description
tuple[float, float]	Tuple containing (radius_beam_tangency, radius_beam_tangency_max)

Source code in process/models/build.py

def calculate_beam_port_size(
    self,
    f_radius_beam_tangency_rmajor: float,
    rmajor: float,
    n_tf_coils: int,
    dx_tf_inboard_out_toroidal: float,
    dr_tf_outboard: float,
    r_tf_outboard_mid: float,
    dx_beam_duct: float,
    dx_beam_shield: float,
) -> tuple[float, float]:
    """Calculates the maximum possible tangency radius for adequate beam access.

    Parameters
    ----------
    f_radius_beam_tangency_rmajor : float
        Fraction of rmajor for beam tangency
    rmajor : float
        Major radius
    n_tf_coils : int
        Number of TF coils
    dx_tf_inboard_out_toroidal : float
        Toroidal width of outboard TF coil
    dr_tf_outboard : float
        Radial thickness of outboard TF coil leg
    r_tf_outboard_mid : float
        Major radius of centre of outboard TF coil
    dx_beam_duct : float
        Width of beam duct
    dx_beam_shield : float
        Shielding width on both sides of beam duct

    Returns
    -------
    tuple[float, float]
        Tuple containing (radius_beam_tangency, radius_beam_tangency_max)
    """

    # Have kept the single letter variable names to match the original code and documentation diagram.
    radius_beam_tangency = f_radius_beam_tangency_rmajor * rmajor

    omega = 2.0 * np.pi / n_tf_coils

    a = 0.5e0 * dx_tf_inboard_out_toroidal
    try:
        assert a < np.inf
    except AssertionError:
        logger.exception("a is inf. Kludging to 1e10.")
        a = 1e10

    b = dr_tf_outboard
    try:
        assert b < np.inf
    except AssertionError:
        logger.exception("b is inf. Kludging to 1e10.")
        b = 1e10

    c = dx_beam_duct + 2.0e0 * dx_beam_shield

    d = r_tf_outboard_mid - 0.5e0 * b
    try:
        assert d < np.inf
    except AssertionError:
        logger.exception("d is inf. Kludging to 1e10.")
        d = 1e10

    e = np.sqrt(a**2 + (d + b) ** 2)
    f = np.sqrt(a**2 + d**2)

    theta = omega - np.arctan(a / d)
    phi = theta - np.arcsin(a / e)

    g = np.sqrt(e**2 + f**2 - 2.0e0 * e * f * np.cos(phi))

    if g > c:
        h = np.sqrt(g**2 - c**2)
        alpha = np.arctan(h / c)
        eps = np.arcsin(e * np.sin(phi) / g) - alpha
        radius_beam_tangency_max = f * np.cos(eps) - 0.5e0 * c
    else:
        logger.error(
            f"Max beam tangency radius set =0 temporarily; change dx_beam_duct. {g=} {c=}"
        )
        radius_beam_tangency_max = 0.0e0

    return radius_beam_tangency, radius_beam_tangency_max

calculate_vertical_build(output)

Determines the vertical build of the machine.

This method calculates various parameters related to the vertical build of the machine, such as thicknesses, radii, and areas. Results can be outputted with the output flag.

Parameters:

Name	Type	Description	Default
output	bool	Flag indicating whether to output results	required

Source code in process/models/build.py

def calculate_vertical_build(self, output: bool):
    """Determines the vertical build of the machine.

    This method calculates various parameters related to the vertical build of the machine,
    such as thicknesses, radii, and areas. Results can be outputted with the `output` flag.

    Parameters
    ----------
    output : bool
        Flag indicating whether to output results
    """
    # Set the X-point heights for the top and bottom of the plasma
    # Assumes top-down plasma symmetry
    build_variables.z_plasma_xpoint_upper = (
        physics_variables.rminor * physics_variables.kappa
    )
    build_variables.z_plasma_xpoint_lower = (
        physics_variables.rminor * physics_variables.kappa
    )

    if output:
        po.oheadr(self.outfile, "Vertical Build")

        po.ovarin(
            self.mfile,
            "Divertor null switch",
            "(i_single_null)",
            physics_variables.i_single_null,
        )

        if physics_variables.i_single_null == 0:
            po.ocmmnt(self.outfile, "Double null case")

            # Start at the top and work down.

            vertical_build_upper = (
                buildings_variables.dz_tf_cryostat
                + build_variables.dr_tf_inboard
                + build_variables.dr_tf_shld_gap
                + build_variables.dz_shld_thermal
                + build_variables.dz_shld_vv_gap
                + build_variables.dz_vv_upper
                + build_variables.dz_shld_upper
                + divertor_variables.dz_divertor
                + build_variables.dz_xpoint_divertor
                + build_variables.z_plasma_xpoint_upper
            )

            # To calculate vertical offset between TF coil centre and plasma centre
            vbuile1 = vertical_build_upper

            po.obuild(
                self.outfile,
                "Cryostat roof structure*",
                buildings_variables.dz_tf_cryostat,
                vertical_build_upper,
                "(dz_tf_cryostat)",
            )
            po.ovarre(
                self.mfile,
                "Cryostat roof structure*",
                "(dz_tf_cryostat)",
                buildings_variables.dz_tf_cryostat,
            )
            vertical_build_upper = (
                vertical_build_upper - buildings_variables.dz_tf_cryostat
            )

            # Top of TF coil
            tf_top = vertical_build_upper

            po.obuild(
                self.outfile,
                "TF coil",
                build_variables.dr_tf_inboard,
                vertical_build_upper,
                "(dr_tf_inboard)",
            )
            vertical_build_upper = (
                vertical_build_upper - build_variables.dr_tf_inboard
            )

            po.obuild(
                self.outfile,
                "Gap",
                build_variables.dr_tf_shld_gap,
                vertical_build_upper,
                "(dr_tf_shld_gap)",
            )
            vertical_build_upper = (
                vertical_build_upper - build_variables.dr_tf_shld_gap
            )

            po.obuild(
                self.outfile,
                "Thermal shield, vertical",
                build_variables.dz_shld_thermal,
                vertical_build_upper,
                "(dz_shld_thermal)",
            )

            po.ovarre(
                self.mfile,
                "Thermal shield, vertical (m)",
                "(dz_shld_thermal)",
                build_variables.dz_shld_thermal,
            )
            vertical_build_upper = (
                vertical_build_upper - build_variables.dz_shld_thermal
            )

            po.obuild(
                self.outfile,
                "Gap",
                build_variables.dz_shld_vv_gap,
                vertical_build_upper,
                "(dz_shld_vv_gap)",
            )
            po.ovarre(
                self.mfile,
                "Vessel - TF coil vertical gap (m)",
                "(dz_shld_vv_gap)",
                build_variables.dz_shld_vv_gap,
            )
            vertical_build_upper = (
                vertical_build_upper - build_variables.dz_shld_vv_gap
            )

            po.obuild(
                self.outfile,
                "Vacuum vessel (and shielding)",
                build_variables.dz_vv_upper + build_variables.dz_shld_upper,
                vertical_build_upper,
                "(dz_vv_upper+dz_shld_upper)",
            )
            vertical_build_upper = (
                vertical_build_upper
                - build_variables.dz_vv_upper
                - build_variables.dz_shld_upper
            )
            po.ovarre(
                self.mfile,
                "Topside vacuum vessel radial thickness (m)",
                "(dz_vv_upper)",
                build_variables.dz_vv_upper,
            )
            po.ovarre(
                self.mfile,
                "Top radiation shield thickness (m)",
                "(dz_shld_upper)",
                build_variables.dz_shld_upper,
            )

            po.obuild(
                self.outfile,
                "Divertor structure",
                divertor_variables.dz_divertor,
                vertical_build_upper,
                "(dz_divertor)",
            )
            po.ovarre(
                self.mfile,
                "Divertor structure vertical thickness (m)",
                "(dz_divertor)",
                divertor_variables.dz_divertor,
            )
            vertical_build_upper = (
                vertical_build_upper - divertor_variables.dz_divertor
            )

            po.obuild(
                self.outfile,
                "Top scrape-off",
                build_variables.dz_xpoint_divertor,
                vertical_build_upper,
                "(dz_xpoint_divertor)",
            )
            po.ovarre(
                self.mfile,
                "Top scrape-off vertical thickness (m)",
                "(dz_xpoint_divertor)",
                build_variables.dz_xpoint_divertor,
            )
            vertical_build_upper = (
                vertical_build_upper - build_variables.dz_xpoint_divertor
            )

            po.obuild(
                self.outfile,
                "Plasma upper X-point height (m)",
                build_variables.z_plasma_xpoint_upper,
                vertical_build_upper,
                "(z_plasma_xpoint_upper)",
            )
            po.ovarre(
                self.mfile,
                "Plasma upper X-point height (m)",
                "(z_plasma_xpoint_upper)",
                build_variables.z_plasma_xpoint_upper,
            )
            vertical_build_upper = (
                vertical_build_upper - build_variables.z_plasma_xpoint_upper
            )

            po.obuild(self.outfile, "Midplane", 0.0e0, vertical_build_upper)

            vertical_build_upper = (
                vertical_build_upper - build_variables.z_plasma_xpoint_lower
            )
            po.obuild(
                self.outfile,
                "Plasma lower X-point height (m)",
                build_variables.z_plasma_xpoint_lower,
                vertical_build_upper,
                "(z_plasma_xpoint_lower)",
            )
            po.ovarre(
                self.mfile,
                "Plasma lower X-point height (m)",
                "(z_plasma_xpoint_lower)",
                build_variables.z_plasma_xpoint_lower,
            )

            vertical_build_upper = (
                vertical_build_upper - build_variables.dz_xpoint_divertor
            )
            po.obuild(
                self.outfile,
                "Lower scrape-off",
                build_variables.dz_xpoint_divertor,
                vertical_build_upper,
                "(dz_xpoint_divertor)",
            )
            po.ovarre(
                self.mfile,
                "Bottom scrape-off vertical thickness (m)",
                "(dz_xpoint_divertor)",
                build_variables.dz_xpoint_divertor,
            )

            vertical_build_upper = (
                vertical_build_upper - divertor_variables.dz_divertor
            )
            po.obuild(
                self.outfile,
                "Divertor structure",
                divertor_variables.dz_divertor,
                vertical_build_upper,
                "(dz_divertor)",
            )
            po.ovarre(
                self.mfile,
                "Divertor structure vertical thickness (m)",
                "(dz_divertor)",
                divertor_variables.dz_divertor,
            )

            vertical_build_upper = (
                vertical_build_upper - build_variables.dz_shld_lower
            )

            vertical_build_upper = vertical_build_upper - build_variables.dz_vv_lower
            po.obuild(
                self.outfile,
                "Vacuum vessel (and shielding)",
                build_variables.dz_vv_lower + build_variables.dz_shld_lower,
                vertical_build_upper,
                "(dz_vv_lower+dz_shld_lower)",
            )
            po.ovarre(
                self.mfile,
                "Bottom radiation shield thickness (m)",
                "(dz_shld_lower)",
                build_variables.dz_shld_lower,
            )
            po.ovarre(
                self.mfile,
                "Underside vacuum vessel radial thickness (m)",
                "(dz_vv_lower)",
                build_variables.dz_vv_lower,
            )

            vertical_build_upper = (
                vertical_build_upper - build_variables.dz_shld_vv_gap
            )
            po.obuild(
                self.outfile,
                "Gap",
                build_variables.dz_shld_vv_gap,
                vertical_build_upper,
                "(dz_shld_vv_gap)",
            )

            vertical_build_upper = (
                vertical_build_upper - build_variables.dz_shld_thermal
            )
            po.obuild(
                self.outfile,
                "Thermal shield, vertical",
                build_variables.dz_shld_thermal,
                vertical_build_upper,
                "(dz_shld_thermal)",
            )

            vertical_build_upper = (
                vertical_build_upper - build_variables.dr_tf_shld_gap
            )
            po.obuild(
                self.outfile,
                "Gap",
                build_variables.dr_tf_shld_gap,
                vertical_build_upper,
                "(dr_tf_shld_gap)",
            )

            vertical_build_upper = (
                vertical_build_upper - build_variables.dr_tf_inboard
            )
            po.obuild(
                self.outfile,
                "TF coil",
                build_variables.dr_tf_inboard,
                vertical_build_upper,
                "(dr_tf_inboard)",
            )

            # Total height of TF coil
            tf_height = tf_top - vertical_build_upper
            # Inner vertical dimension of TF coil
            build_variables.dh_tf_inner_bore = (
                tf_height - 2 * build_variables.dr_tf_inboard
            )

            vertical_build_upper = (
                vertical_build_upper - buildings_variables.dz_tf_cryostat
            )
            po.obuild(
                self.outfile,
                "Cryostat floor structure**",
                buildings_variables.dz_tf_cryostat,
                vertical_build_upper,
                "(dz_tf_cryostat)",
            )

            # To calculate vertical offset between TF coil centre and plasma centre
            build_variables.dz_tf_plasma_centre_offset = (
                vbuile1 + vertical_build_upper
            ) / 2.0e0

            # End of Double null case
        else:
            po.ocmmnt(self.outfile, "Single null case")
            build_variables.dz_vv_upper = 0.5 * (
                build_variables.dz_vv_upper + build_variables.dz_vv_lower
            )

            build_variables.dz_fw_upper = 0.5 * (
                build_variables.dr_fw_inboard + build_variables.dr_fw_outboard
            )

            vbuild = (
                buildings_variables.dz_tf_cryostat
                + build_variables.dr_tf_inboard
                + build_variables.dr_tf_shld_gap
                + build_variables.dz_shld_thermal
                + build_variables.dz_shld_vv_gap
                + build_variables.dz_vv_upper
                + build_variables.dr_shld_blkt_gap
                + build_variables.dz_shld_upper
                + build_variables.dz_blkt_upper
                + build_variables.dz_fw_upper
                + build_variables.dz_fw_plasma_gap
                + build_variables.z_plasma_xpoint_upper
            )

            # To calculate vertical offset between TF coil centre and plasma centre
            vbuile1 = vbuild

            po.obuild(
                self.outfile,
                "Cryostat roof structure*",
                buildings_variables.dz_tf_cryostat,
                vbuild,
                "(dz_tf_cryostat)",
            )
            po.ovarre(
                self.mfile,
                "Cryostat roof structure*",
                "(dz_tf_cryostat)",
                buildings_variables.dz_tf_cryostat,
            )
            vbuild = vbuild - buildings_variables.dz_tf_cryostat

            # Top of TF coil
            tf_top = vbuild

            po.obuild(
                self.outfile,
                "TF coil",
                build_variables.dr_tf_inboard,
                vbuild,
                "(dr_tf_inboard)",
            )
            vbuild = vbuild - build_variables.dr_tf_inboard

            po.obuild(
                self.outfile,
                "Gap",
                build_variables.dr_tf_shld_gap,
                vbuild,
                "(dr_tf_shld_gap)",
            )
            vbuild = vbuild - build_variables.dr_tf_shld_gap

            po.obuild(
                self.outfile,
                "Thermal shield, vertical",
                build_variables.dz_shld_thermal,
                vbuild,
                "(dz_shld_thermal)",
            )
            po.ovarre(
                self.mfile,
                "Thermal shield, vertical (m)",
                "(dz_shld_thermal)",
                build_variables.dz_shld_thermal,
            )
            vbuild = vbuild - build_variables.dz_shld_thermal

            po.obuild(
                self.outfile,
                "Gap",
                build_variables.dz_shld_vv_gap,
                vbuild,
                "(dz_shld_vv_gap)",
            )
            po.ovarre(
                self.mfile,
                "Vessel - TF coil vertical gap (m)",
                "(dz_shld_vv_gap)",
                build_variables.dz_shld_vv_gap,
            )
            vbuild = vbuild - build_variables.dz_shld_vv_gap

            po.obuild(
                self.outfile,
                "Vacuum vessel (and shielding)",
                build_variables.dz_vv_upper + build_variables.dz_shld_upper,
                vbuild,
                "(dz_vv_upper+dz_shld_upper)",
            )
            vbuild = (
                vbuild - build_variables.dz_vv_upper - build_variables.dz_shld_upper
            )
            po.ovarre(
                self.mfile,
                "Topside vacuum vessel radial thickness (m)",
                "(dz_vv_upper)",
                build_variables.dz_vv_upper,
            )
            po.ovarre(
                self.mfile,
                "Top radiation shield thickness (m)",
                "(dz_shld_upper)",
                build_variables.dz_shld_upper,
            )

            po.obuild(
                self.outfile,
                "Gap",
                build_variables.dr_shld_blkt_gap,
                vbuild,
                "(dr_shld_blkt_gap)",
            )
            vbuild = vbuild - build_variables.dr_shld_blkt_gap

            po.obuild(
                self.outfile,
                "Top blanket",
                build_variables.dz_blkt_upper,
                vbuild,
                "(dz_blkt_upper)",
            )
            po.ovarre(
                self.mfile,
                "Top blanket vertical thickness (m)",
                "(dz_blkt_upper)",
                build_variables.dz_blkt_upper,
            )
            vbuild = vbuild - build_variables.dz_blkt_upper

            dz_fw_upper = 0.5e0 * (
                build_variables.dr_fw_inboard + build_variables.dr_fw_outboard
            )
            po.obuild(
                self.outfile, "Top first wall", dz_fw_upper, vbuild, "(dz_fw_upper)"
            )
            po.ovarre(
                self.mfile,
                "Top first wall vertical thickness (m)",
                "(dz_fw_upper)",
                dz_fw_upper,
            )
            vbuild = vbuild - dz_fw_upper

            po.obuild(
                self.outfile,
                "Top scrape-off",
                build_variables.dz_fw_plasma_gap,
                vbuild,
                "(dz_fw_plasma_gap)",
            )
            po.ovarre(
                self.mfile,
                "Top scrape-off vertical thickness (m)",
                "(dz_fw_plasma_gap)",
                build_variables.dz_fw_plasma_gap,
            )
            vbuild = vbuild - build_variables.dz_fw_plasma_gap

            po.obuild(
                self.outfile,
                "Plasma upper X-point height (m)",
                build_variables.z_plasma_xpoint_upper,
                vbuild,
                "(z_plasma_xpoint_upper)",
            )
            po.ovarre(
                self.mfile,
                "Plasma upper X-point height (m)",
                "(z_plasma_xpoint_upper)",
                build_variables.z_plasma_xpoint_upper,
            )
            vbuild = vbuild - build_variables.z_plasma_xpoint_upper

            po.obuild(self.outfile, "Midplane", 0.0e0, vbuild)

            vbuild = vbuild - build_variables.z_plasma_xpoint_lower
            po.obuild(
                self.outfile,
                "Plasma lower X-point height (m)",
                build_variables.z_plasma_xpoint_lower,
                vbuild,
                "(z_plasma_xpoint_lower)",
            )
            po.ovarre(
                self.mfile,
                "Plasma lower X-point height (m)",
                "(z_plasma_xpoint_lower)",
                build_variables.z_plasma_xpoint_lower,
            )

            vbuild = vbuild - build_variables.dz_xpoint_divertor
            po.obuild(
                self.outfile,
                "Lower scrape-off",
                build_variables.dz_xpoint_divertor,
                vbuild,
                "(dz_xpoint_divertor)",
            )
            po.ovarre(
                self.mfile,
                "Bottom scrape-off vertical thickness (m)",
                "(dz_xpoint_divertor)",
                build_variables.dz_xpoint_divertor,
            )

            vbuild = vbuild - divertor_variables.dz_divertor
            po.obuild(
                self.outfile,
                "Divertor structure",
                divertor_variables.dz_divertor,
                vbuild,
                "(dz_divertor)",
            )
            po.ovarre(
                self.mfile,
                "Divertor structure vertical thickness (m)",
                "(dz_divertor)",
                divertor_variables.dz_divertor,
            )

            vbuild = vbuild - build_variables.dz_shld_lower

            vbuild = vbuild - build_variables.dz_vv_lower
            po.obuild(
                self.outfile,
                "Vacuum vessel (and shielding)",
                build_variables.dz_vv_lower + build_variables.dz_shld_lower,
                vbuild,
                "(dz_vv_lower+dz_shld_lower)",
            )
            po.ovarre(
                self.mfile,
                "Bottom radiation shield thickness (m)",
                "(dz_shld_lower)",
                build_variables.dz_shld_lower,
            )
            po.ovarre(
                self.mfile,
                "Underside vacuum vessel radial thickness (m)",
                "(dz_vv_lower)",
                build_variables.dz_vv_lower,
            )

            vbuild = vbuild - build_variables.dz_shld_vv_gap
            po.obuild(
                self.outfile,
                "Gap",
                build_variables.dz_shld_vv_gap,
                vbuild,
                "(dz_shld_vv_gap)",
            )

            vbuild = vbuild - build_variables.dz_shld_thermal
            po.obuild(
                self.outfile,
                "Thermal shield, vertical",
                build_variables.dz_shld_thermal,
                vbuild,
                "(dz_shld_thermal)",
            )

            vbuild = vbuild - build_variables.dr_tf_shld_gap
            po.obuild(
                self.outfile,
                "Gap",
                build_variables.dr_tf_shld_gap,
                vbuild,
                "(dr_tf_shld_gap)",
            )

            vbuild = vbuild - build_variables.dr_tf_inboard
            po.obuild(
                self.outfile,
                "TF coil",
                build_variables.dr_tf_inboard,
                vbuild,
                "(dr_tf_inboard)",
            )

            # Total height of TF coil
            tf_height = tf_top - vbuild
            # Inner vertical dimension of TF coil
            build_variables.dh_tf_inner_bore = (
                tf_height - 2 * build_variables.dr_tf_inboard
            )

            vbuild = vbuild - buildings_variables.dz_tf_cryostat

            po.obuild(
                self.outfile,
                "Cryostat floor structure**",
                buildings_variables.dz_tf_cryostat,
                vbuild,
                "(dz_tf_cryostat)",
            )

            # To calculate vertical offset between TF coil centre and plasma centre
            build_variables.dz_tf_plasma_centre_offset = (vbuile1 + vbuild) / 2.0e0

            # end of Single null case

        po.ovarre(
            self.mfile,
            "Ratio of Central Solenoid height to TF coil internal height",
            "(f_z_cs_tf_internal)",
            pfcoil_variables.f_z_cs_tf_internal,
        )
        po.ocmmnt(
            self.outfile,
            "\n*Cryostat roof allowance includes uppermost PF coil and outer thermal shield.\n*Cryostat floor allowance includes lowermost PF coil, outer thermal shield and gravity support.",
        )

    # Output the cdivertor geometry
    divht = self.divgeom(output)
    # Issue #481 Remove build_variables.vgaptf
    if build_variables.dz_xpoint_divertor < 0.00001e0:
        build_variables.dz_xpoint_divertor = divht

    # If build_variables.dz_xpoint_divertor /= 0 use the value set by the user.

    # Height to inside edge of TF coil. TF coils are assumed to be symmetrical.
    # Therefore this applies to single and double null cases.
    build_variables.z_tf_inside_half = (
        build_variables.z_plasma_xpoint_upper
        + build_variables.dz_xpoint_divertor
        + divertor_variables.dz_divertor
        + build_variables.dz_shld_lower
        + build_variables.dz_vv_lower
        + build_variables.dz_shld_vv_gap
        + build_variables.dz_shld_thermal
        + build_variables.dr_tf_shld_gap
    )

    #  Vertical locations of divertor coils
    if physics_variables.i_single_null == 0:
        build_variables.z_tf_top = (
            build_variables.z_tf_inside_half + build_variables.dr_tf_inboard
        )
        build_variables.dz_tf_upper_lower_midplane = 0.0e0
    else:
        build_variables.z_tf_top = (
            build_variables.dr_tf_inboard
            + build_variables.dr_tf_shld_gap
            + build_variables.dz_shld_thermal
            + build_variables.dz_shld_vv_gap
            + build_variables.dz_vv_upper
            + build_variables.dz_shld_upper
            + build_variables.dr_shld_blkt_gap
            + build_variables.dz_blkt_upper
            + 0.5e0
            * (build_variables.dr_fw_inboard + build_variables.dr_fw_outboard)
            + build_variables.dz_fw_plasma_gap
            + build_variables.z_plasma_xpoint_upper
        )
        build_variables.dz_tf_upper_lower_midplane = build_variables.z_tf_top - (
            build_variables.z_tf_inside_half + build_variables.dr_tf_inboard
        )

divgeom(output)

Divertor geometry calculation

This subroutine determines the divertor geometry. The inboard (i) and outboard (o) plasma surfaces are approximated by arcs, and followed past the X-point to determine the maximum height. TART option: Peng SOFT paper

Parameters:

Name	Type	Description	Default
output	bool		required

Returns:

Name	Type	Description
divht		divertor height (m)

Source code in process/models/build.py

def divgeom(self, output: bool):
    """Divertor geometry calculation

    This subroutine determines the divertor geometry.
    The inboard (i) and outboard (o) plasma surfaces
    are approximated by arcs, and followed past the X-point to
    determine the maximum height.
    TART option: Peng SOFT paper

    Parameters
    ----------
    output: bool

    Returns
    -------
    divht:
        divertor height (m)
    """
    if physics_variables.itart == 1:
        return 1.75e0 * physics_variables.rminor
    #  Conventional tokamak divertor model
    #  options for seperate upper and lower physics_variables.triangularity

    kap = physics_variables.kappa
    triu = physics_variables.triang
    tril = physics_variables.triang

    # Old method: assumes that divertor arms are continuations of arcs
    #
    # Outboard side
    # build_variables.plsepo = poloidal length along the separatrix from null to
    # strike point on outboard [default 1.5 m]
    # thetao = arc angle between the strike point and the null point
    #
    # xpointo = physics_variables.rmajor + 0.5e0*physics_variables.rminor*(kap**2 + tri**2 - 1.0e0) /     #     (1.0e0 - tri)
    # rprimeo = (xpointo - physics_variables.rmajor + physics_variables.rminor)
    # phio = asin(kap*physics_variables.rminor/rprimeo)
    # thetao = build_variables.plsepo/rprimeo
    #
    # Initial strike point
    #
    # yspointo = rprimeo * sin(thetao + phio)
    # xspointo = xpointo - rprimeo * cos(thetao + phio)
    #
    # Outboard strike point radius - normalized to ITER
    #
    # rstrko = xspointo + 0.14e0
    #
    # Uppermost divertor strike point (end of power decay)
    # anginc = angle of incidence of scrape-off field lines on the
    # divertor (rad)
    #
    # +**PJK 25/07/11 Changed sign of anginc contribution
    # yprimeb = soleno * cos(thetao + phio - anginc)
    #
    # divht = yprimeb + yspointo - kap*physics_variables.rminor

    # New method, assuming straight legs -- superceded by new method 26/5/2016
    # Assumed 90 degrees at X-pt -- wrong#
    #
    #  Find half-angle of outboard arc
    # denomo = (tril**2 + kap**2 - 1.0e0)/( 2.0e0*(1.0e0+tril) ) - tril
    # thetao = atan(kap/denomo)
    # Angle between horizontal and inner divertor leg
    # alphad = (pi/2.0e0) - thetao

    # Method 26/05/2016
    # Find radius of inner and outer plasma arcs

    rco = 0.5 * np.sqrt(
        (physics_variables.rminor**2 * ((tril + 1.0e0) ** 2 + kap**2) ** 2)
        / ((tril + 1.0e0) ** 2)
    )
    rci = 0.5 * np.sqrt(
        (physics_variables.rminor**2 * ((tril - 1.0e0) ** 2 + kap**2) ** 2)
        / ((tril - 1.0e0) ** 2)
    )

    # Find angles between vertical and legs
    # Inboard arc angle = outboard leg angle

    thetao = np.arcsin(1.0e0 - (physics_variables.rminor * (1.0e0 - tril)) / rci)

    # Outboard arc angle = inboard leg angle

    thetai = np.arcsin(1.0e0 - (physics_variables.rminor * (1.0e0 + tril)) / rco)

    #  Position of lower x-pt
    rxpt = physics_variables.rmajor - tril * physics_variables.rminor
    zxpt = -1.0e0 * kap * physics_variables.rminor

    # Position of inner strike point
    # rspi = rxpt - build_variables.plsepi*cos(alphad)
    # zspi = zxpt - build_variables.plsepi*sin(alphad)
    rspi = rxpt - build_variables.plsepi * np.cos(thetai)
    zspi = zxpt - build_variables.plsepi * np.sin(thetai)

    # Position of outer strike point
    # build_variables.rspo = rxpt + build_variables.plsepo*cos((pi/2.0e0)-alphad)
    # zspo = zxpt - build_variables.plsepo*sin((pi/2.0e0)-alphad)
    build_variables.rspo = rxpt + build_variables.plsepo * np.cos(thetao)
    zspo = zxpt - build_variables.plsepo * np.sin(thetao)

    # Position of inner plate ends
    # rplti = rspi - (build_variables.plleni/2.0e0)*sin(divertor_variables.betai + alphad - pi/2.0e0)
    # zplti = zspi + (build_variables.plleni/2.0e0)*cos(divertor_variables.betai + alphad - pi/2.0e0)
    # rplbi = rspi + (build_variables.plleni/2.0e0)*sin(divertor_variables.betai + alphad - pi/2.0e0)
    # zplbi = zspi - (build_variables.plleni/2.0e0)*cos(divertor_variables.betai + alphad - pi/2.0e0)
    rplti = rspi + (build_variables.plleni / 2.0e0) * np.cos(
        thetai + divertor_variables.betai
    )
    zplti = zspi + (build_variables.plleni / 2.0e0) * np.sin(
        thetai + divertor_variables.betai
    )
    rplbi = rspi - (build_variables.plleni / 2.0e0) * np.cos(
        thetai + divertor_variables.betai
    )
    zplbi = zspi - (build_variables.plleni / 2.0e0) * np.sin(
        thetai + divertor_variables.betai
    )

    # Position of outer plate ends
    # rplto = build_variables.rspo + (build_variables.plleno/2.0e0)*sin(divertor_variables.betao - alphad)
    # zplto = zspo + (build_variables.plleno/2.0e0)*cos(divertor_variables.betao - alphad)
    # rplbo = build_variables.rspo - (build_variables.plleno/2.0e0)*sin(divertor_variables.betao - alphad)
    # zplbo = zspo - (build_variables.plleno/2.0e0)*cos(divertor_variables.betao - alphad)
    rplto = build_variables.rspo - (build_variables.plleno / 2.0e0) * np.cos(
        thetao + divertor_variables.betao
    )
    zplto = zspo + (build_variables.plleno / 2.0e0) * np.sin(
        thetao + divertor_variables.betao
    )
    rplbo = build_variables.rspo + (build_variables.plleno / 2.0e0) * np.cos(
        thetao + divertor_variables.betao
    )
    zplbo = zspo - (build_variables.plleno / 2.0e0) * np.sin(
        thetao + divertor_variables.betao
    )

    divht = max(zplti, zplto) - min(zplbo, zplbi)

    if output:
        if divertor_variables.n_divertors == 1:
            po.oheadr(self.outfile, "Divertor build and plasma position")
            po.ocmmnt(self.outfile, "Divertor Configuration = Single Null Divertor")
            po.oblnkl(self.outfile)
            ptop_radial = physics_variables.rmajor - triu * physics_variables.rminor
            ptop_vertical = kap * physics_variables.rminor
            po.ovarrf(
                self.outfile,
                "Plasma top position, radial (m)",
                "(ptop_radial)",
                ptop_radial,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Plasma top position, vertical (m)",
                "(ptop_vertical)",
                ptop_vertical,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Plasma geometric centre, radial (m)",
                "(rmajor.)",
                physics_variables.rmajor,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Plasma geometric centre, vertical (m)",
                "(0.0)",
                0.0e0,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Plasma lower triangularity",
                "(tril)",
                tril,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Plasma elongation",
                "(kappa.)",
                kap,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "TF coil vertical offset (m)",
                "(dz_tf_plasma_centre_offset)",
                build_variables.dz_tf_plasma_centre_offset,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Plasma outer arc radius of curvature (m)",
                "(rco)",
                rco,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Plasma inner arc radius of curvature (m)",
                "(rci)",
                rci,
                "OP ",
            )
            po.ovarrf(
                self.outfile, "Plasma lower X-pt, radial (m)", "(rxpt)", rxpt, "OP "
            )
            po.ovarrf(
                self.outfile,
                "Plasma lower X-pt, vertical (m)",
                "(zxpt)",
                zxpt,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Poloidal plane angle between vertical and inner leg (rad)",
                "(thetai)",
                thetai,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Poloidal plane angle between vertical and outer leg (rad)",
                "(thetao)",
                thetao,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Poloidal plane angle between inner leg and plate (rad)",
                "(betai)",
                divertor_variables.betai,
            )
            po.ovarrf(
                self.outfile,
                "Poloidal plane angle between outer leg and plate (rad)",
                "(betao)",
                divertor_variables.betao,
            )
            po.ovarrf(
                self.outfile,
                "Inner divertor leg poloidal length (m)",
                "(plsepi)",
                build_variables.plsepi,
            )
            po.ovarrf(
                self.outfile,
                "Outer divertor leg poloidal length (m)",
                "(plsepo)",
                build_variables.plsepo,
            )
            po.ovarrf(
                self.outfile,
                "Inner divertor plate length (m)",
                "(plleni)",
                build_variables.plleni,
            )
            po.ovarrf(
                self.outfile,
                "Outer divertor plate length (m)",
                "(plleno)",
                build_variables.plleno,
            )
            po.ovarrf(
                self.outfile,
                "Inner strike point, radial (m)",
                "(rspi)",
                rspi,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Inner strike point, vertical (m)",
                "(zspi)",
                zspi,
                "OP ",
            )
            po.ovarrf(
                self.outfile, "Inner plate top, radial (m)", "(rplti)", rplti, "OP "
            )
            po.ovarrf(
                self.outfile,
                "Inner plate top, vertical (m)",
                "(zplti)",
                zplti,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Inner plate bottom, radial (m)",
                "(rplbi)",
                rplbi,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Inner plate bottom, vertical (m)",
                "(zplbi)",
                zplbi,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Outer strike point, radial (m)",
                "(rspo)",
                build_variables.rspo,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Outer strike point, vertical (m)",
                "(zspo)",
                zspo,
                "OP ",
            )
            po.ovarrf(
                self.outfile, "Outer plate top, radial (m)", "(rplto)", rplto, "OP "
            )
            po.ovarrf(
                self.outfile,
                "Outer plate top, vertical (m)",
                "(zplto)",
                zplto,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Outer plate bottom, radial (m)",
                "(rplbo)",
                rplbo,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Outer plate bottom, vertical (m)",
                "(zplbo)",
                zplbo,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Calculated maximum divertor height (m)",
                "(divht)",
                divht,
                "OP ",
            )

        elif divertor_variables.n_divertors == 2:
            po.oheadr(self.outfile, "Divertor build and plasma position")
            po.ocmmnt(self.outfile, "Divertor Configuration = Double Null Divertor")
            po.oblnkl(self.outfile)
            # Assume upper and lower divertors geometries are symmetric.
            ptop_radial = physics_variables.rmajor - triu * physics_variables.rminor
            ptop_vertical = kap * physics_variables.rminor
            po.ovarrf(
                self.outfile,
                "Plasma top position, radial (m)",
                "(ptop_radial)",
                ptop_radial,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Plasma top position, vertical (m)",
                "(ptop_vertical)",
                ptop_vertical,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Plasma geometric centre, radial (m)",
                "(rmajor.)",
                physics_variables.rmajor,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Plasma geometric centre, vertical (m)",
                "(0.0)",
                0.0e0,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Plasma physics_variables.triangularity",
                "(tril)",
                tril,
                "OP ",
            )
            po.ovarrf(self.outfile, "Plasma elongation", "(kappa.)", kap, "OP ")
            po.ovarrf(
                self.outfile,
                "TF coil vertical offset (m)",
                "(dz_tf_plasma_centre_offset)",
                build_variables.dz_tf_plasma_centre_offset,
                "OP ",
            )
            po.ovarrf(
                self.outfile, "Plasma upper X-pt, radial (m)", "(rxpt)", rxpt, "OP "
            )
            po.ovarrf(
                self.outfile,
                "Plasma upper X-pt, vertical (m)",
                "(-zxpt)",
                -zxpt,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Plasma outer arc radius of curvature (m)",
                "(rco)",
                rco,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Plasma inner arc radius of curvature (m)",
                "(rci)",
                rci,
                "OP ",
            )
            po.ovarrf(
                self.outfile, "Plasma lower X-pt, radial (m)", "(rxpt)", rxpt, "OP "
            )
            po.ovarrf(
                self.outfile,
                "Plasma lower X-pt, vertical (m)",
                "(zxpt)",
                zxpt,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Poloidal plane angle between vertical and inner leg (rad)",
                "(thetai)",
                thetai,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Poloidal plane angle between vertical and outer leg (rad)",
                "(thetao)",
                thetao,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Poloidal plane angle between inner leg and plate (rad)",
                "(betai)",
                divertor_variables.betai,
            )
            po.ovarrf(
                self.outfile,
                "Poloidal plane angle between outer leg and plate (rad)",
                "(betao)",
                divertor_variables.betao,
            )
            po.ovarrf(
                self.outfile,
                "Inner divertor leg poloidal length (m)",
                "(plsepi)",
                build_variables.plsepi,
            )
            po.ovarrf(
                self.outfile,
                "Outer divertor leg poloidal length (m)",
                "(plsepo)",
                build_variables.plsepo,
            )
            po.ovarrf(
                self.outfile,
                "Inner divertor plate length (m)",
                "(lleni)",
                build_variables.plleni,
            )
            po.ovarrf(
                self.outfile,
                "Outer divertor plate length (m)",
                "(plleno)",
                build_variables.plleno,
            )
            po.ovarrf(
                self.outfile,
                "Upper inner strike point, radial (m)",
                "(rspi)",
                rspi,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Upper inner strike point, vertical (m)",
                "(-zspi)",
                -zspi,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Upper inner plate top, radial (m)",
                "(rplti)",
                rplti,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Upper inner plate top, vertical (m)",
                "(-zplti)",
                -zplti,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Upper inner plate bottom, radial (m)",
                "(rplbi)",
                rplbi,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Upper inner plate bottom, vertical (m)",
                "(-zplbi)",
                -zplbi,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Upper outer strike point, radial (m)",
                "(rspo)",
                build_variables.rspo,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Upper outer strike point, vertical (m)",
                "(-zspo)",
                -zspo,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Upper outer plate top, radial (m)",
                "(rplto)",
                rplto,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Upper outer plate top, vertical (m)",
                "(-zplto)",
                -zplto,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Upper outer plate bottom, radial (m)",
                "(rplbo)",
                rplbo,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Upper outer plate bottom, vertical (m)",
                "(-zplbo)",
                -zplbo,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Lower inner strike point, radial (m)",
                "(rspi)",
                rspi,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Lower inner strike point, vertical (m)",
                "(zspi)",
                zspi,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Lower inner plate top, radial (m)",
                "(rplti)",
                rplti,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Lower inner plate top, vertical (m)",
                "(zplti)",
                zplti,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Lower inner plate bottom, radial (m)",
                "(rplbi)",
                rplbi,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Lower inner plate bottom, vertical (m)",
                "(zplbi)",
                zplbi,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Lower outer strike point, radial (m)",
                "(rspo)",
                build_variables.rspo,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Lower outer strike point, vertical (m)",
                "(zspo)",
                zspo,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Lower outer plate top, radial (m)",
                "(rplto)",
                rplto,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Lower outer plate top, vertical (m)",
                "(zplto)",
                zplto,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Lower outer plate bottom, radial (m)",
                "(rplbo)",
                rplbo,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Lower outer plate bottom, vertical (m)",
                "(zplbo)",
                zplbo,
                "OP ",
            )
            po.ovarrf(
                self.outfile,
                "Calculated maximum divertor height (m)",
                "(divht)",
                divht,
                "OP ",
            )
        else:
            po.oheadr(self.outfile, "Divertor build and plasma position")
            po.ocmmnt(
                self.outfile,
                "ERROR: null value not supported, check i_single_null value.",
            )
    return divht

plasma_outboard_edge_toroidal_ripple(ripple_b_tf_plasma_edge_max, r_tf_outboard_mid, n_tf_coils, rmajor, rminor, r_tf_wp_inboard_inner, r_tf_wp_inboard_centre, r_tf_wp_inboard_outer, dx_tf_wp_primary_toroidal, i_tf_shape, i_tf_sup, dx_tf_wp_insulation, dx_tf_wp_insertion_gap)

Plasma outboard toroidal field (TF) ripple calculation.

This routine computes the TF ripple amplitude at the midplane outboard plasma edge and the minimum radius of the TF coil centre that would produce a specified maximum allowed ripple. The calculation uses fitted coefficients derived from numerical modelling (MAGINT) and includes a simplified analytical picture-frame coil model for i_tf_shape == 2.

Parameters:

Name	Type	Description	Default
ripple_b_tf_plasma_edge_max	float	Maximum allowed ripple at plasma edge (percent)	required
r_tf_outboard_mid	float	Radius to the centre of the outboard TF coil leg (m)	required
n_tf_coils	int	Number of TF coils	required
rmajor	float	Plasma major radius (m)	required
rminor	float	Plasma minor radius (m)	required
r_tf_wp_inboard_inner	float	Inner winding-pack inboard radius (m)	required
r_tf_wp_inboard_centre	float	Centre winding-pack inboard radius (m)	required
r_tf_wp_inboard_outer	float	Outer winding-pack inboard radius (m)	required
dx_tf_wp_primary_toroidal	float	Primary toroidal winding-pack thickness (m)	required
i_tf_shape	int	TF coil shape switch (2 => picture-frame analytical model)	required
i_tf_sup	int	TF coil support flag (1 => superconducting)	required
dx_tf_wp_insulation	float	Winding-pack insulation thickness (m)	required
dx_tf_wp_insertion_gap	float	Winding-pack insertion gap (m)	required

Returns:

Type	Description
tuple[float, float, int]	Tuple containing: - ripple: Calculated ripple at plasma edge (percent) - r_tf_outboard_midmin: Minimum r_tf_outboard_mid that yields the specified maximum ripple (m) - flag: Applicability flag (0 = OK, non-zero = fitted-range concern)

Notes

• Fitted coefficients originate from parametric MAGINT runs (M. Kovari, 2014).
• Picture-frame coil analytical model (Ken McClements, 2022) is used when i_tf_shape == 2 and gives approximate results (within ~10% of numerical).
• The routine sets an applicability flag when fitted-range assumptions are exceeded.

Source code in process/models/build.py

def plasma_outboard_edge_toroidal_ripple(
    self,
    ripple_b_tf_plasma_edge_max: float,
    r_tf_outboard_mid: float,
    n_tf_coils: int,
    rmajor: float,
    rminor: float,
    r_tf_wp_inboard_inner,
    r_tf_wp_inboard_centre: float,
    r_tf_wp_inboard_outer: float,
    dx_tf_wp_primary_toroidal: float,
    i_tf_shape: int,
    i_tf_sup: int,
    dx_tf_wp_insulation: float,
    dx_tf_wp_insertion_gap: float,
) -> tuple[float, float, int]:
    """Plasma outboard toroidal field (TF) ripple calculation.

    This routine computes the TF ripple amplitude at the midplane outboard
    plasma edge and the minimum radius of the TF coil centre that would
    produce a specified maximum allowed ripple. The calculation uses
    fitted coefficients derived from numerical modelling (MAGINT) and
    includes a simplified analytical picture-frame coil model for
    i_tf_shape == 2.

    Parameters
    ----------
    ripple_b_tf_plasma_edge_max : float
        Maximum allowed ripple at plasma edge (percent)
    r_tf_outboard_mid : float
        Radius to the centre of the outboard TF coil leg (m)
    n_tf_coils : int
        Number of TF coils
    rmajor : float
        Plasma major radius (m)
    rminor : float
        Plasma minor radius (m)
    r_tf_wp_inboard_inner : float
        Inner winding-pack inboard radius (m)
    r_tf_wp_inboard_centre : float
        Centre winding-pack inboard radius (m)
    r_tf_wp_inboard_outer : float
        Outer winding-pack inboard radius (m)
    dx_tf_wp_primary_toroidal : float
        Primary toroidal winding-pack thickness (m)
    i_tf_shape : int
        TF coil shape switch (2 => picture-frame analytical model)
    i_tf_sup : int
        TF coil support flag (1 => superconducting)
    dx_tf_wp_insulation : float
        Winding-pack insulation thickness (m)
    dx_tf_wp_insertion_gap : float
        Winding-pack insertion gap (m)

    Returns
    -------
    tuple[float, float, int]
        Tuple containing:
        - ripple: Calculated ripple at plasma edge (percent)
        - r_tf_outboard_midmin: Minimum r_tf_outboard_mid that yields the specified maximum ripple (m)
        - flag: Applicability flag (0 = OK, non-zero = fitted-range concern)

    Notes
    -----
    - Fitted coefficients originate from parametric MAGINT runs (M. Kovari, 2014).
    - Picture-frame coil analytical model (Ken McClements, 2022) is used when
    `i_tf_shape == 2` and gives approximate results (within ~10% of numerical).
    - The routine sets an applicability flag when fitted-range assumptions are exceeded.
    """
    if i_tf_sup == 1:
        # Minimal inboard WP radius [m]
        r_wp_min = r_tf_wp_inboard_inner

        # Rectangular WP
        if tfcoil_variables.i_tf_wp_geom == 0:
            r_wp_max = r_wp_min

        # Double rectangle WP
        elif tfcoil_variables.i_tf_wp_geom == 1:
            r_wp_max = r_tf_wp_inboard_centre

        # Trapezoidal WP
        elif tfcoil_variables.i_tf_wp_geom == 2:
            r_wp_max = r_tf_wp_inboard_outer

        # Calculated maximum toroidal WP toroidal thickness [m]
        dx_tf_wp_conductor_max = dx_tf_wp_primary_toroidal - 2.0 * (
            dx_tf_wp_insulation + dx_tf_wp_insertion_gap
        )

    # Resistive magnet case
    else:
        # Radius used to define the dx_tf_wp_conductor_max [m]
        r_wp_max = r_tf_wp_inboard_outer
        # Calculated maximum toroidal WP toroidal thickness [m]
        dx_tf_wp_conductor_max = 2.0e0 * r_wp_max * np.tan(np.pi / n_tf_coils)

    flag = 0
    if i_tf_shape == 2:
        # Ken McClements ST picture frame coil analytical ripple calc
        # Calculated ripple for coil at r_tf_outboard_mid (%)
        ripple = 100.0e0 * ((rmajor + rminor) / r_tf_outboard_mid) ** (n_tf_coils)
        #  Calculated r_tf_outboard_mid to produce a ripple of amplitude ripple_b_tf_plasma_edge_max
        r_tf_outboard_midmin = (rmajor + rminor) / (
            (0.01e0 * ripple_b_tf_plasma_edge_max) ** (1.0e0 / n_tf_coils)
        )
    else:
        # Winding pack to iter-coil at plasma centre toroidal lenth ratio
        x = dx_tf_wp_conductor_max * n_tf_coils / rmajor

        # Fitting parameters
        c1 = 0.875e0 - 0.0557e0 * x
        c2 = 1.617e0 + 0.0832e0 * x

        #  Calculated ripple for coil at r_tf_outboard_mid (%)
        ripple = (
            100.0e0
            * c1
            * ((rmajor + rminor) / r_tf_outboard_mid) ** (n_tf_coils - c2)
        )

        #  Calculated r_tf_outboard_mid to produce a ripple of amplitude ripple_b_tf_plasma_edge_max
        base = 0.01 * ripple_b_tf_plasma_edge_max / c1
        # Avoid potential negative or complex result: kludge base to be
        # small and positive if required
        try:
            assert base > 1e-6
        except AssertionError:
            logger.exception("base is <= 1e-6. Kludging to 1e-6.")
            base = 1e-6

        r_tf_outboard_midmin = (rmajor + rminor) / (
            base ** (1.0 / (n_tf_coils - c2))
        )

        try:
            assert r_tf_outboard_midmin < np.inf
        except AssertionError:
            logger.exception(
                "r_tf_outboard_midmin is inf. Kludging to a large value instead."
            )
            r_tf_outboard_midmin = (rmajor + rminor) * 3

        #  Notify via flag if a range of applicability is violated
        flag = 0
        if (x < 0.737e0) or (x > 2.95e0):
            flag = 1
        if (n_tf_coils < 16) or (n_tf_coils > 20):
            flag = 2
        if ((rmajor + rminor) / r_tf_outboard_mid < 0.7e0) or (
            (rmajor + rminor) / r_tf_outboard_mid > 0.8e0
        ):
            flag = 3

    return ripple, r_tf_outboard_midmin, flag

calculate_radial_build(output)

This method determines the radial build of the machine. It calculates various parameters related to the build of the machine, such as thicknesses, radii, and areas. Results can be outputted with the output flag.

Parameters:

Name	Type	Description	Default
output	bool	Flag indicating whether to output the results	required

Source code in process/models/build.py

def calculate_radial_build(self, output: bool):
    """This method determines the radial build of the machine.
    It calculates various parameters related to the build of the machine,
    such as thicknesses, radii, and areas.
    Results can be outputted with the `output` flag.

    Parameters
    ----------
    output : bool
        Flag indicating whether to output the results
    """

    if fwbs_variables.blktmodel > 0:
        build_variables.dr_blkt_inboard = (
            build_variables.blbuith
            + build_variables.blbmith
            + build_variables.blbpith
        )
        build_variables.dr_blkt_outboard = (
            build_variables.blbuoth
            + build_variables.blbmoth
            + build_variables.blbpoth
        )
        build_variables.dz_shld_upper = 0.5e0 * (
            build_variables.dr_shld_inboard + build_variables.dr_shld_outboard
        )

    #  Top/bottom blanket thickness
    build_variables.dz_blkt_upper = 0.5e0 * (
        build_variables.dr_blkt_inboard + build_variables.dr_blkt_outboard
    )

    if physics_variables.i_single_null == 1:
        #  Check if build_variables.dz_fw_plasma_gap has been set too small
        build_variables.dz_fw_plasma_gap = max(
            0.5e0
            * (
                build_variables.dr_fw_plasma_gap_inboard
                + build_variables.dr_fw_plasma_gap_outboard
            ),
            build_variables.dz_fw_plasma_gap,
        )

    # Calculate pre-compression structure thickness is build_variables.i_cs_precomp=1
    if build_variables.i_cs_precomp == 1 and build_variables.i_tf_inside_cs == 0:
        build_variables.dr_cs_precomp = build_variables.fseppc / (
            2.0e0
            * np.pi
            * build_variables.fcspc
            * build_variables.sigallpc
            * (
                build_variables.dr_bore
                + build_variables.dr_bore
                + build_variables.dr_cs
            )
        )
    elif build_variables.i_cs_precomp == 1 and build_variables.i_tf_inside_cs == 1:
        build_variables.dr_cs_precomp = build_variables.fseppc / (
            2.0e0
            * np.pi
            * build_variables.fcspc
            * build_variables.sigallpc
            * (
                2.0 * build_variables.dr_bore
                + 2.0 * build_variables.dr_tf_inboard
                + 2.0 * build_variables.dr_cs_tf_gap
                + build_variables.dr_cs
            )
        )
    else:
        build_variables.dr_cs_precomp = 0.0e0

    # Issue #514 Radial dimensions of inboard leg
    # Calculate build_variables.dr_tf_inboard if tfcoil_variables.dr_tf_wp_with_insulation is an iteration variable (140)
    if 140 in numerics.ixc[0 : numerics.nvar]:
        build_variables.dr_tf_inboard = (
            tfcoil_variables.dr_tf_wp_with_insulation
            + tfcoil_variables.dr_tf_plasma_case
            + tfcoil_variables.dr_tf_nose_case
        )

    if build_variables.i_tf_inside_cs == 1:
        build_variables.r_tf_inboard_in = (
            build_variables.dr_bore
            # NOTE: dr_bore is just the hollow space, the
            # true dr_bore size used for flux calculations
            # is dr_bore + dr_tf_inboard + dr_cs_tf_gap
        )
    else:
        # Inboard side inner radius [m]
        build_variables.r_tf_inboard_in = (
            build_variables.dr_bore
            + build_variables.dr_cs
            + build_variables.dr_cs_precomp
            + build_variables.dr_cs_tf_gap
        )

    # Radial build to tfcoil middle [m]
    build_variables.r_tf_inboard_mid = (
        build_variables.r_tf_inboard_in + 0.5e0 * build_variables.dr_tf_inboard
    )

    # Radial build to tfcoil plasma facing side [m]
    build_variables.r_tf_inboard_out = (
        build_variables.r_tf_inboard_in + build_variables.dr_tf_inboard
    )

    # WP radial thickness [m]
    # Calculated only if not used as an iteration variable
    if 140 not in numerics.ixc[0 : numerics.nvar]:
        tfcoil_variables.dr_tf_wp_with_insulation = (
            build_variables.dr_tf_inboard
            - tfcoil_variables.dr_tf_plasma_case
            - tfcoil_variables.dr_tf_nose_case
        )

    # Radius of the centrepost at the top of the machine
    if physics_variables.itart == 1 and tfcoil_variables.i_tf_sup != 1:
        # build_variables.r_cp_top is set using the plasma shape
        if build_variables.i_r_cp_top == 0:
            build_variables.r_cp_top = (
                physics_variables.rmajor
                - physics_variables.rminor * physics_variables.triang
                - (
                    build_variables.dr_tf_shld_gap
                    + build_variables.dr_shld_thermal_inboard
                    + build_variables.dr_shld_inboard
                    + build_variables.dr_shld_blkt_gap
                    + build_variables.dr_blkt_inboard
                    + build_variables.dr_fw_inboard
                    + 3.0e0 * build_variables.dr_fw_plasma_gap_inboard
                )
                + tfcoil_variables.drtop
            )

            # Notify user that build_variables.r_cp_top has been set to 1.01*build_variables.r_tf_inboard_out (lvl 2 error)
            if build_variables.r_cp_top < 1.01e0 * build_variables.r_tf_inboard_out:
                logger.error(
                    "TF CP top radius (r_cp_top) replaced by 1.01*r_tf_inboard_out -> potential top rbuild issue"
                    f"{build_variables.r_cp_top=} {build_variables.r_tf_inboard_out=}"
                )

                # build_variables.r_cp_top correction
                build_variables.r_cp_top = build_variables.r_tf_inboard_out * 1.01e0

            # Top and mid-plane TF coil CP radius ratio
            build_variables.f_r_cp = (
                build_variables.r_cp_top / build_variables.r_tf_inboard_out
            )

        # User defined build_variables.r_cp_top
        elif build_variables.i_r_cp_top == 1:
            # Notify user that build_variables.r_cp_top has been set to 1.01*build_variables.r_tf_inboard_out (lvl 2 error)
            if build_variables.r_cp_top < 1.01e0 * build_variables.r_tf_inboard_out:
                logger.error(
                    "TF CP top radius (r_cp_top) replaced by 1.01*r_tf_inboard_out -> potential top rbuild issue"
                    f"{build_variables.r_cp_top=} {build_variables.r_tf_inboard_out=}"
                )

                # build_variables.r_cp_top correction
                build_variables.r_cp_top = build_variables.r_tf_inboard_out * 1.01e0

            # Top / mid-plane TF CP radius ratio
            build_variables.f_r_cp = (
                build_variables.r_cp_top / build_variables.r_tf_inboard_out
            )

        # build_variables.r_cp_top set as a fraction of the outer TF midplane radius
        elif build_variables.i_r_cp_top == 2:
            build_variables.r_cp_top = (
                build_variables.f_r_cp * build_variables.r_tf_inboard_out
            )

    else:  # End of physics_variables.itart == 1 .and. tfcoil_variables.i_tf_sup /= 1
        build_variables.r_cp_top = build_variables.r_tf_inboard_out

    if build_variables.i_r_cp_top != 0 and (
        build_variables.r_cp_top
        > physics_variables.rmajor
        - physics_variables.rminor * physics_variables.triang
        - (
            build_variables.dr_tf_shld_gap
            + build_variables.dr_shld_thermal_inboard
            + build_variables.dr_shld_inboard
            + build_variables.dr_shld_blkt_gap
            + build_variables.dr_blkt_inboard
            + build_variables.dr_fw_inboard
            + 3.0e0 * build_variables.dr_fw_plasma_gap_inboard
        )
        + tfcoil_variables.drtop
    ):
        logger.error(
            f"Top CP radius larger that its value determined with plasma shape {build_variables.r_cp_top=}"
        )
    if build_variables.i_tf_inside_cs == 1:
        #  Radial position of vacuum vessel [m]
        build_variables.r_vv_inboard_out = (
            build_variables.r_tf_inboard_out
            + build_variables.dr_cs
            + build_variables.dr_cs_tf_gap
            + build_variables.dr_cs_precomp
            + build_variables.dr_tf_shld_gap
            + build_variables.dr_shld_thermal_inboard
            + build_variables.dr_shld_vv_gap_inboard
            + build_variables.dr_vv_inboard
        )
    else:
        build_variables.r_vv_inboard_out = (
            build_variables.r_tf_inboard_out
            + build_variables.dr_tf_shld_gap
            + build_variables.dr_shld_thermal_inboard
            + build_variables.dr_shld_vv_gap_inboard
            + build_variables.dr_vv_inboard
        )
    # Radial position of the inner side of inboard neutronic shield [m]
    build_variables.r_sh_inboard_in = build_variables.r_vv_inboard_out

    # Radial position of the plasma facing side of inboard neutronic shield [m]
    build_variables.r_sh_inboard_out = (
        build_variables.r_sh_inboard_in + build_variables.dr_shld_inboard
    )

    #  Radial build to centre of plasma (should be equal to physics_variables.rmajor)
    build_variables.rbld = (
        build_variables.r_sh_inboard_out
        + build_variables.dr_shld_blkt_gap
        + build_variables.dr_blkt_inboard
        + build_variables.dr_fw_inboard
        + build_variables.dr_fw_plasma_gap_inboard
        + physics_variables.rminor
    )

    #  Radius to inner edge of inboard shield
    build_variables.r_shld_inboard_inner = (
        physics_variables.rmajor
        - physics_variables.rminor
        - build_variables.dr_fw_plasma_gap_inboard
        - build_variables.dr_fw_inboard
        - build_variables.dr_blkt_inboard
        - build_variables.dr_shld_inboard
    )

    #  Radius to outer edge of outboard shield
    build_variables.r_shld_outboard_outer = (
        physics_variables.rmajor
        + physics_variables.rminor
        + build_variables.dr_fw_plasma_gap_outboard
        + build_variables.dr_fw_outboard
        + build_variables.dr_blkt_outboard
        + build_variables.dr_shld_outboard
    )

    #  Thickness of outboard TF coil legs
    if tfcoil_variables.i_tf_sup != 1:
        build_variables.dr_tf_outboard = (
            build_variables.f_dr_tf_outboard_inboard * build_variables.dr_tf_inboard
        )
    else:
        build_variables.dr_tf_outboard = build_variables.dr_tf_inboard

    #  Radius to centre of outboard TF coil legs
    build_variables.r_tf_outboard_mid = (
        build_variables.r_shld_outboard_outer
        + build_variables.dr_shld_blkt_gap
        + build_variables.dr_vv_outboard
        + build_variables.gapomin
        + build_variables.dr_shld_thermal_outboard
        + build_variables.dr_tf_shld_gap
        + 0.5e0 * build_variables.dr_tf_outboard
    )

    # TF coil horizontal build_variables.dr_bore [m]
    build_variables.dr_tf_inner_bore = (
        build_variables.r_tf_outboard_mid - 0.5e0 * build_variables.dr_tf_outboard
    ) - (build_variables.r_tf_inboard_mid - 0.5e0 * build_variables.dr_tf_inboard)

    (
        tfcoil_variables.ripple_b_tf_plasma_edge,
        r_tf_outboard_midl,
        build_variables.ripflag,
    ) = self.plasma_outboard_edge_toroidal_ripple(
        ripple_b_tf_plasma_edge_max=tfcoil_variables.ripple_b_tf_plasma_edge_max,
        r_tf_outboard_mid=build_variables.r_tf_outboard_mid,
        n_tf_coils=tfcoil_variables.n_tf_coils,
        rmajor=physics_variables.rmajor,
        rminor=physics_variables.rminor,
        r_tf_wp_inboard_inner=superconducting_tf_coil_variables.r_tf_wp_inboard_inner,
        r_tf_wp_inboard_centre=superconducting_tf_coil_variables.r_tf_wp_inboard_centre,
        r_tf_wp_inboard_outer=superconducting_tf_coil_variables.r_tf_wp_inboard_outer,
        dx_tf_wp_primary_toroidal=tfcoil_variables.dx_tf_wp_primary_toroidal,
        i_tf_shape=tfcoil_variables.i_tf_shape,
        i_tf_sup=tfcoil_variables.i_tf_sup,
        dx_tf_wp_insulation=tfcoil_variables.dx_tf_wp_insulation,
        dx_tf_wp_insertion_gap=tfcoil_variables.dx_tf_wp_insertion_gap,
    )

    #  If the tfcoil_variables.ripple is too large then move the outboard TF coil leg
    if r_tf_outboard_midl > build_variables.r_tf_outboard_mid:
        build_variables.r_tf_outboard_mid = r_tf_outboard_midl
        build_variables.dr_shld_vv_gap_outboard = (
            build_variables.r_tf_outboard_mid
            - 0.5e0 * build_variables.dr_tf_outboard
            - build_variables.dr_vv_outboard
            - build_variables.r_shld_outboard_outer
            - build_variables.dr_shld_thermal_outboard
            - build_variables.dr_tf_shld_gap
            - build_variables.dr_shld_blkt_gap
        )
        build_variables.dr_tf_inner_bore = (
            build_variables.r_tf_outboard_mid
            - 0.5e0 * build_variables.dr_tf_outboard
        ) - (
            build_variables.r_tf_inboard_mid - 0.5e0 * build_variables.dr_tf_inboard
        )
    else:
        build_variables.dr_shld_vv_gap_outboard = build_variables.gapomin

    (
        tfcoil_variables.ripple_b_tf_plasma_edge,
        r_tf_outboard_midl,
        build_variables.ripflag,
    ) = self.plasma_outboard_edge_toroidal_ripple(
        ripple_b_tf_plasma_edge_max=tfcoil_variables.ripple_b_tf_plasma_edge_max,
        r_tf_outboard_mid=build_variables.r_tf_outboard_mid,
        n_tf_coils=tfcoil_variables.n_tf_coils,
        rmajor=physics_variables.rmajor,
        rminor=physics_variables.rminor,
        r_tf_wp_inboard_inner=superconducting_tf_coil_variables.r_tf_wp_inboard_inner,
        r_tf_wp_inboard_centre=superconducting_tf_coil_variables.r_tf_wp_inboard_centre,
        r_tf_wp_inboard_outer=superconducting_tf_coil_variables.r_tf_wp_inboard_outer,
        dx_tf_wp_primary_toroidal=tfcoil_variables.dx_tf_wp_primary_toroidal,
        i_tf_shape=tfcoil_variables.i_tf_shape,
        i_tf_sup=tfcoil_variables.i_tf_sup,
        dx_tf_wp_insulation=tfcoil_variables.dx_tf_wp_insulation,
        dx_tf_wp_insertion_gap=tfcoil_variables.dx_tf_wp_insertion_gap,
    )

    #

    if output:
        #  Print out device build

        po.oheadr(self.outfile, "Radial Build")

        if build_variables.ripflag != 0:
            po.ocmmnt(
                self.outfile,
                "(Ripple result may not be accurate, as the fit was outside",
            )
            po.ocmmnt(self.outfile, " its range of applicability.)")
            po.oblnkl(self.outfile)
            logger.warning(
                "Ripple result may be inaccurate, as the fit has been extrapolated"
            )

            if build_variables.ripflag == 1:
                diagnostic = (
                    tfcoil_variables.dx_tf_wp_primary_toroidal
                    * tfcoil_variables.n_tf_coils
                    / physics_variables.rmajor
                )
                logger.warning(
                    f"(TF coil ripple calculation) Dimensionless coil width X out of fitted range. {diagnostic=}"
                )
            elif build_variables.ripflag == 2:
                logger.warning(
                    f"(TF coil ripple calculation) No of TF coils not between 16 and 20 inclusive {tfcoil_variables.n_tf_coils=}"
                )
            else:
                diagnostic = (
                    physics_variables.rmajor + physics_variables.rminor
                ) / build_variables.r_tf_outboard_mid

                logger.warning(
                    f"(TF coil ripple calculation) (R+a)/rtot={diagnostic} out of fitted range."
                )

        po.ovarin(
            self.outfile,
            "TF coil radial placement switch",
            "(i_tf_inside_cs)",
            build_variables.i_tf_inside_cs,
        )
        po.ovarrf(
            self.outfile,
            "Inboard build thickness (m)",
            "(dr_inboard_build)",
            physics_variables.rmajor - physics_variables.rminor,
            "OP ",
        )

        if build_variables.i_tf_inside_cs == 1:
            po.ocmmnt(
                self.outfile,
                (
                    "\n (The stated machine dr_bore size is just for the hollow space, "
                ),
            )
            po.ocmmnt(
                self.outfile,
                (
                    "the true dr_bore size used for calculations is dr_bore + dr_tf_inboard + dr_cs_tf_gap)\n"
                ),
            )
        if (
            build_variables.i_tf_inside_cs == 1
            and tfcoil_variables.i_tf_bucking >= 2
        ):
            po.ocmmnt(
                self.outfile,
                "(Bore hollow space has been filled with a solid metal cyclinder to act as wedge support)\n",
            )

        # an array that holds the following information
        # description, variable name, thickness, radius
        radial_build_data = []

        radius = 0.0e0
        radial_build_data.append(["Device centreline", None, 0.0, radius])
        if (
            build_variables.i_tf_inside_cs == 1
            and tfcoil_variables.i_tf_bucking >= 2
        ):
            radius = radius + build_variables.dr_bore

            radial_build_data.append([
                "Machine dr_bore wedge support cylinder",
                "dr_bore",
                build_variables.dr_bore,
                radius,
            ])
        elif (
            build_variables.i_tf_inside_cs == 1 and tfcoil_variables.i_tf_bucking < 2
        ):
            radius = radius + build_variables.dr_bore
            radial_build_data.append([
                "Machine dr_bore hole",
                "dr_bore",
                build_variables.dr_bore,
                radius,
            ])
        else:
            radius = radius + build_variables.dr_bore
            radial_build_data.append([
                "Machine dr_bore",
                "dr_bore",
                build_variables.dr_bore,
                radius,
            ])
        if build_variables.i_tf_inside_cs == 1:
            radius += build_variables.dr_tf_inboard
            radial_build_data.append([
                "TF coil inboard leg (in dr_bore)",
                "dr_tf_inboard",
                build_variables.dr_tf_inboard,
                radius,
            ])

            radius += build_variables.dr_cs_tf_gap
            radial_build_data.append([
                "CS precompresion to TF coil radial gap",
                "dr_cs_tf_gap",
                build_variables.dr_cs_tf_gap,
                radius,
            ])

        radius = radius + build_variables.dr_cs
        radial_build_data.append([
            "Central solenoid",
            "dr_cs",
            build_variables.dr_cs,
            radius,
        ])

        radius = radius + build_variables.dr_cs_precomp
        radial_build_data.append([
            "CS precompression",
            "dr_cs_precomp",
            build_variables.dr_cs_precomp,
            radius,
        ])
        if build_variables.i_tf_inside_cs == 0:
            radius = radius + build_variables.dr_cs_tf_gap
            radial_build_data.append([
                "CS precompresion to TF coil radial gap",
                "dr_cs_tf_gap",
                build_variables.dr_cs_tf_gap,
                radius,
            ])

            radius = radius + build_variables.dr_tf_inboard
            radial_build_data.append([
                "TF coil inboard leg",
                "dr_tf_inboard",
                build_variables.dr_tf_inboard,
                radius,
            ])

        radius = radius + build_variables.dr_tf_shld_gap
        radial_build_data.append([
            "TF coil inboard leg insulation gap",
            "dr_tf_shld_gap",
            build_variables.dr_tf_shld_gap,
            radius,
        ])

        radius = radius + build_variables.dr_shld_thermal_inboard
        radial_build_data.append([
            "Thermal shield, inboard",
            "dr_shld_thermal_inboard",
            build_variables.dr_shld_thermal_inboard,
            radius,
        ])

        radius = radius + build_variables.dr_shld_vv_gap_inboard
        radial_build_data.append([
            "Thermal shield to vessel radial gap",
            "dr_shld_vv_gap_inboard",
            build_variables.dr_shld_vv_gap_inboard,
            radius,
        ])

        radius += build_variables.dr_vv_inboard
        radial_build_data.append([
            "Inboard vacuum vessel",
            "dr_vv_inboard",
            build_variables.dr_vv_inboard,
            radius,
        ])

        radius += build_variables.dr_shld_inboard
        radial_build_data.append([
            "Inner radiation shield",
            "dr_shld_inboard",
            build_variables.dr_shld_inboard,
            radius,
        ])

        radius = radius + build_variables.dr_shld_blkt_gap
        radial_build_data.append([
            "Gap",
            "dr_shld_blkt_gap",
            build_variables.dr_shld_blkt_gap,
            radius,
        ])

        radius = radius + build_variables.dr_blkt_inboard
        radial_build_data.append([
            "Inboard blanket",
            "dr_blkt_inboard",
            build_variables.dr_blkt_inboard,
            radius,
        ])

        radius = radius + build_variables.dr_fw_inboard
        radial_build_data.append([
            "Inboard first wall",
            "dr_fw_inboard",
            build_variables.dr_fw_inboard,
            radius,
        ])

        radius = radius + build_variables.dr_fw_plasma_gap_inboard
        radial_build_data.append([
            "Inboard scrape-off",
            "dr_fw_plasma_gap_inboard",
            build_variables.dr_fw_plasma_gap_inboard,
            radius,
        ])

        radius = radius + physics_variables.rminor
        radial_build_data.append([
            "Plasma geometric centre",
            "rminor",
            physics_variables.rminor,
            radius,
        ])

        radius = radius + physics_variables.rminor
        radial_build_data.append([
            "Plasma outboard edge",
            "rminor",
            physics_variables.rminor,
            radius,
        ])

        radius = radius + build_variables.dr_fw_plasma_gap_outboard
        radial_build_data.append([
            "Outboard scrape-off",
            "dr_fw_plasma_gap_outboard",
            build_variables.dr_fw_plasma_gap_outboard,
            radius,
        ])

        radius = radius + build_variables.dr_fw_outboard
        radial_build_data.append([
            "Outboard first wall",
            "dr_fw_outboard",
            build_variables.dr_fw_outboard,
            radius,
        ])

        radius = radius + build_variables.dr_blkt_outboard
        radial_build_data.append([
            "Outboard blanket",
            "dr_blkt_outboard",
            build_variables.dr_blkt_outboard,
            radius,
        ])

        radius = radius + build_variables.dr_shld_blkt_gap
        radial_build_data.append([
            "Gap",
            "dr_shld_blkt_gap",
            build_variables.dr_shld_blkt_gap,
            radius,
        ])

        radius += build_variables.dr_shld_outboard
        radial_build_data.append([
            "Outer radiation shield",
            "dr_shld_outboard",
            build_variables.dr_shld_outboard,
            radius,
        ])

        radius += build_variables.dr_vv_outboard
        radial_build_data.append([
            "Outboard vacuum vessel",
            "dr_vv_outboard",
            build_variables.dr_vv_outboard,
            radius,
        ])

        radius = radius + build_variables.dr_shld_vv_gap_outboard
        radial_build_data.append([
            "Vessel to TF gap",
            "dr_shld_vv_gap_outboard",
            build_variables.dr_shld_vv_gap_outboard,
            radius,
        ])

        radius = radius + build_variables.dr_shld_thermal_outboard
        radial_build_data.append([
            "Ouboard thermal shield",
            "dr_shld_thermal_outboard",
            build_variables.dr_shld_thermal_outboard,
            radius,
        ])

        radius = radius + build_variables.dr_tf_shld_gap
        radial_build_data.append([
            "Gap",
            "dr_tf_shld_gap",
            build_variables.dr_tf_shld_gap,
            radius,
        ])

        radius = radius + build_variables.dr_tf_outboard
        radial_build_data.append([
            "TF coil outboard leg",
            "dr_tf_outboard",
            build_variables.dr_tf_outboard,
            radius,
        ])

        for description, variable, thickness, radius in radial_build_data:
            po.obuild(
                self.outfile,
                description,
                thickness,
                radius,
                f"({variable})" if variable else "",
            )

        # use manual index to ensure count is contiguous in the event
        # of a `None` variable component
        index = 0
        for description, variable, thickness, radius in radial_build_data:
            if variable is None:
                continue

            index += 1

            po.ovarre(
                self.mfile,
                f"{description} radial thickness (m)",
                f"({variable})",
                thickness,
            )

            po.ovarst(
                self.mfile,
                f"Radial build component {index}",
                f"(radial_label({index}))",
                f'"{variable}"',
            )
            po.ovarre(
                self.mfile,
                f"Radial build cumulative radius {index}",
                f"(radial_cum({index}))",
                radius,
            )

        if (current_drive_variables.i_hcd_primary in [5, 8]) or (
            current_drive_variables.i_hcd_secondary in [5, 8]
        ):
            po.ovarre(
                self.mfile,
                "Width of neutral beam duct where it passes between the TF coils (m)",
                "(dx_beam_duct)",
                current_drive_variables.dx_beam_duct,
            )