diff --git a/process/divertor.py b/process/divertor.py index 02bac08915..2be5806a49 100644 --- a/process/divertor.py +++ b/process/divertor.py @@ -32,8 +32,10 @@ def run(self, output: bool) -> None: :param output: indicate whether output should be written to the output file, or not :type output: boolean """ - - if pv.itart == 1 and (dv.i_hldiv == 0 or dv.i_hldiv == 1): + if dv.i_hldiv == 0 and output: + po.ovarre(self.outfile, "Divertor heat load (MW/m2)", "(hldiv)", dv.hldiv) + return + if dv.i_hldiv == 1: self.divtart( pv.rmajor, pv.rminor, @@ -61,692 +63,6 @@ def run(self, output: bool) -> None: ) return - # Scale geometric quantities - - # Perpendicular diffusivity in the plasma scrapeoff (m2/s) - # Assume ions transport 33% of electron power - - xperp = dv.xpertin * 1.33e0 - - # Reference null to strike distances - # Only set up for outer divertor for double-null - - # plsep = min(plsepo,pi*rminor) # Obscure reason to set a limit... - plsep = bv.plsepo - - # Scale plasma quantities - - delne = dv.prn1 * pv.dene * 1.0e-20 # scrapeoff density by main plasma - pwr = pv.pdivt # power flow to divertor (MW) - aionso = pv.m_fuel_amu # scrape-off layer ion mass - - # Divertor Zeff: either scaled or use input value - zeffso = (1.0e0 + 0.8e0 * (pv.zeff - 1.0e0)) if dv.divdum == 0 else dv.zeffdiv - - # Strike point field values - - bpstk = pv.bp * 0.45e0 - btstk = pv.bt * pv.rmajor / bv.rspo - rbpbt = bpstk / btstk - - # Parallel diffusivity in the plasma scrapeoff (m2/s) - - xpara = dv.xparain / zeffso - - # Null radius - - rnull = pv.rmajor - pv.rminor * pv.triang - - # Divertor area and radius ratio - - rsrd = (rnull + pv.rmajor + pv.rminor) / (rnull + bv.rspo) - diva = constants.pi * (rnull + bv.rspo) * plsep - adas = diva / pv.a_plasma_surface - - # Main plasma separatrix area to divertor (and power fraction) - # +PJK Is the 2 related to 2 divertors (i.e. double-null assumed)? - frgd = (pv.a_plasma_surface_outboard) / (2.0e0 * pv.a_plasma_surface) - # -PJK - # Power flow to divertor - - pdiv = pwr * dv.ksic / 2.0e0 - qdiv = pdiv / (pv.a_plasma_surface * frgd) - - # Connection length scalings - # (2.5 factor comes from normalization to ITER 1990) - - tconl = ( - 2.5e0 - * pv.rmajor - * pv.q95 - * (1.0e0 + 1.0e0 / (pv.q95 * pv.aspect) ** 2) ** 0.5 - ) - dtheta = plsep / pv.rminor - dconl = ( - 2.5e0 - * bv.rspo - * pv.q95 - * dtheta - * (1.0e0 + 1.0e0 / (pv.q95 * pv.aspect) ** 2) ** 0.5 - ) - rconl = dconl / tconl - - # Minimum strike angle - - minstang = 0.5e0 - - # Call divertor routine - - ( - delta, - delw, - dv.dendiv, - dv.densin, - gamdt, - dv.lamp, - omlarg, - dv.ppdiv, - dv.ppdivr, - dv.ptpdiv, - dv.tdiv, - dv.tsep, - ) = self.divert( - adas, - aionso, - dv.anginc, - delne, - dv.c1div, - dv.c2div, - dv.c3div, - dv.c4div, - dv.c5div, - dv.delld, - dv.fdfs, - dv.fififi, - frgd, - dv.frrp, - minstang, - dv.omegan, - qdiv, - pdiv, - rbpbt, - rconl, - pv.rmajor, - rsrd, - tconl, - xpara, - xperp, - ) - - # Heat load - - dv.hldiv = dv.ppdivr - - # Ratio of collision length to connection length - - dv.rlclolcn = 1.44e-3 * dv.tsep**2 / (delne * 15.0e0 * tconl) - - # output deleted as per discussion in #242. - # tnunn: 19/01/2022 - - def divert( - self, - adas, - m_ions_total_amu, - anginc, - delne, - c1div, - c2div, - c3div, - c4div, - c5div, - delld, - fdfs, - fififi, - frgd, - frrp, - minstang, - omegan, - qdiv, - pdiv, - rbpbt, - rconl, - rmaj, - rsrd, - tconl, - xpara, - xperp, - ): - """Harrison-Kukushkin analytic ITER divertor model - author: J Galambos, ORNL - author: P J Knight, CCFE, Culham Science Centre - - This subroutine performs the iteration described in M. Harrison's - and Kukushkin's analytic ITER divertor model. - Report ITER-IL-PH-13-9-e12 - - :param adas: divertor flux area / main plasma area (long separatrix) - :type adas: float - - :param m_ions_total_amu: ion mass (assumes fuel only) (AMU) - :type m_ions_total_amu: float - - - :param anginc: pol. angle of incidence of field line on plate (rad) - :type anginc: float - - :param c1div: fitting coefficient for plate temperature - :type c1div: float - - :param c2div: fitting coefficient for plate temperature - :type c2div: float - - :param c3div: fitting coefficient for heat load - :type c3div: float - - :param c4div: fitting coefficient for heat load - :type c4div: float - - :param c5div: fitting coefficient for 'omlarg' - :type c5div: float - - :param delld: coeff. for power distribution flow into scrapeoff - :type delld: float - - :param delne: scrapeoff density by main plasma (10**20 m-3) - :type delne: float - - :param fdfs: gradient ratio (private flux side/other side) in 'omlarg' - :type fdfs: float - - :param fififi: coeff. used in sheath energy transfer factor calc. - :type fififi: float - - :param frgd: separatrix area to divertor / total separatrix area - :type frgd: float - - :param frrp: fraction of radiated power to plate - :type frrp: float - - :param minstang: minimum strike angle (total) for heat flux calc. - :type minstang: float - - :param omegan: pressure ratio of (plate / main plasma) - :type omegan: float - - :param qdiv: heat flux across separatrix to divertor (MW/m2) - :type qdiv: float - - :param pdiv: power flow to plate (MW) - :type pdiv: float - - :param rbpbt: ratio of toroidal / poloidal field at strike point - :type rbpbt: float - - :param rconl: connection length ratio (divertor region/main plasma region) - :type rconl: float - - :param rmaj: major radius (m) - :type rmaj: float - - :param rsrd: ratio of separatrix radius / divertor radius - :type rsrd: float - - :param tconl: connection length along field line by main plasma (m) - :type tconl: float - - :param xpara: parallel diffusivity in the plasma scrapeoff (m2/s) - :type xpara: float - - :param xperp: perpend. diffusivity in the plasma scrapeoff (m2/s) - :type xperp: float - - :returns: - - delta (`float`) iteration relative error - - delw (`float`) energy flow thickness in scrape-off (m) - - dendiv (`float`) plasma density at divertor (10**20 m-3) - - densin (`float`) peak plasma density at divertor (on separatrix) (10**20 m-3) - - gamdt (`float`) plasma flow to plate (10**20/s) - - lamp (`float`) power flow width (m) - - omlarg (`float`) factor accounting for power flow to private flux region - - ppdiv (`float`) divertor heat load without radiation (MW/m2) - - ppdivr (`float`) divertor heat load with radiation (MW/m2) - - ptpdiv (`float`) peak plasma temperature at the divertor plate (eV) - - tdiv (`float`) temperature at the plate (eV) - - tsep (`float`) temperature at the separatrix (eV) - """ - - c27 = 0.2857143e0 - ei = 13.6e0 - epsilon = 0.001e0 - relerr = 1.0e-9 - - fprime = c5div * fdfs - facdenom = fprime * rsrd * (adas / frgd) ** 2 / rconl - facdenom = max(facdenom, 0.04e0) - omlarg = 1.0e0 / (rsrd * math.exp(-facdenom)) - omlarg = min(omlarg, 2.0e0) - coefl = 1.0e0 / delld + rconl / omlarg # little 'l' in Harrison model - - # Start iteration on 2 simultaneous equations (Newton's method) - - tdivges = 150.0e0 - tptsges = 0.9e0 - tdiv = tdivges - tpts = tptsges - - for _i in range(15): - # Find derivatives for Newton's method - - tptsp = tpts * (1.0e0 + epsilon) - deltx = tpts * epsilon - tdivp = tdiv * (1.0e0 + epsilon) - delty = tdiv * epsilon - - f1 = self.ftpts( - m_ions_total_amu, - coefl, - delne, - fififi, - omegan, - omlarg, - qdiv, - tconl, - xpara, - xperp, - tpts, - tdiv, - ) - f2 = self.ftdiv( - m_ions_total_amu, - coefl, - delne, - fififi, - omegan, - omlarg, - qdiv, - tconl, - xpara, - xperp, - tpts, - tdiv, - ) - - f1dx = ( - self.ftpts( - m_ions_total_amu, - coefl, - delne, - fififi, - omegan, - omlarg, - qdiv, - tconl, - xpara, - xperp, - tptsp, - tdiv, - ) - - f1 - ) / deltx - f1dy = ( - self.ftpts( - m_ions_total_amu, - coefl, - delne, - fififi, - omegan, - omlarg, - qdiv, - tconl, - xpara, - xperp, - tpts, - tdivp, - ) - - f1 - ) / delty - f2dx = ( - self.ftdiv( - m_ions_total_amu, - coefl, - delne, - fififi, - omegan, - omlarg, - qdiv, - tconl, - xpara, - xperp, - tptsp, - tdiv, - ) - - f2 - ) / deltx - f2dy = ( - self.ftdiv( - m_ions_total_amu, - coefl, - delne, - fififi, - omegan, - omlarg, - qdiv, - tconl, - xpara, - xperp, - tpts, - tdivp, - ) - - f2 - ) / delty - - denom = f1dx * f2dy - f1dy * f2dx - if denom == 0.0e0: - denom = 1.0e-10 - deltpts = (-f2dy * f1 + f1dy * f2) / denom - deltdiv = (f2dx * f1 - f1dx * f2) / denom - - # New guess - - tdiv = tdiv + deltdiv - tpts = tpts + deltpts - delta = abs(deltdiv / tdiv + deltpts / tpts) - - if delta < relerr: - break - - tdiv = max(tdiv, 0.1000e0) - tpts = max(tpts, 0.0010e0) - tpts = min(tpts, 0.9999e0) - - # Some other quantities - - ct = max(0.1e0, (c1div + c2div / (tdiv))) - ptpdiv = tdiv * ct - gamdiv = self.gammash(fififi, tdiv) # sheath coefficient - dendiv = delne / (omegan * tpts) - eier = self.erprcy( - tdiv, dendiv - ) # ionization + radiation energy / recycle event - - tsep = ( - 251.0e0 - * ( - (qdiv * tconl) ** 2 - / (c27 * xpara * (1.0e0 - tpts**3.5e0)) - * coefl - / (xperp * delne) - ) - ** 0.2222222e0 - ) - - cp = max(0.1e0, (c3div + c4div / (tdiv))) - angle = math.sin(anginc) * rbpbt - - if minstang != 0.0e0: - angle = max(angle, (minstang / 57.3e0)) - - ppdiv = ( - 2.48e2 - * (qdiv) ** 1.55556e0 - / (xperp * delne) ** 0.777778e0 - * (c27 * xpara) ** 0.2222222e0 - * tconl**0.555556e0 - * ((1.0e0 - tpts**3.5e0) / coefl) ** 0.222222e0 - / omlarg - * (1.0e0 + ei / (gamdiv * tdiv)) - / (1.0e0 + eier / (gamdiv * tdiv)) - * angle - * cp - ) - ppdivr = ppdiv * (1.0e0 + frrp * (eier - ei) / (gamdiv * tdiv)) - gamdt = 6.25e4 * ppdiv / (gamdiv * ptpdiv) - densin = omegan * tsep * delne / ptpdiv - delw = ( - 4.01e-3 - * (delne * xperp) ** 0.7777778e0 - * tconl**0.4444444e0 - * coefl**0.2222222e0 - / ( - (qdiv) ** 0.55555556e0 - * (c27 * xpara * (1.0e0 - tpts**3.5e0)) ** 0.22222e0 - ) - ) - lamp = pdiv * rsrd / (2.0e0 * constants.pi * rmaj * ppdiv) - - return ( - delta, - delw, - dendiv, - densin, - gamdt, - lamp, - omlarg, - ppdiv, - ppdivr, - ptpdiv, - tdiv, - tsep, - ) - - def erprcy(self, tdiv: float, ndiv: float) -> float: - """Function providing the (energy radiated + ionized) per neutral - recycle event from the Harrison / Kukushkin ITER model - author: J Galambos, ORNL - author: P J Knight, CCFE, Culham Science Centre - - - This function calculates the total energy (in eV) radiated and - ionized, per neutral recycle event, from the Harrison / Kukushkin - ITER model. - Report ITER-IL-PH-13-9-e12 - - :param tdiv: electron temperature at the plate (eV) - :type tdiv: float - - :param ndiv: electron density at the plate (10**20 m-3) - :type ndiv: float - - :returns: the total energy (in eV) radiated and ionized, per neutral recycle event, from the Harrison / Kukushkin ITER model. - :rtype: float - """ - - return max(17.5e0 + (5.0e0 + 37.5e0 / tdiv) * math.log10(10.0e0 / ndiv), 0.001) - - def ftdiv( - self, - m_ions_total_amu: float, - coefl: float, - delne: float, - fififi: float, - omegan: float, - omlarg: float, - qdiv: float, - tconl: float, - xpara: float, - xperp: float, - xx: float, - yy: float, - ) -> float: - """Function for divertor temperature solution - author: J Galambos, ORNL - author: P J Knight, CCFE, Culham Science Centre - - This function calculates an estimate for the divertor temperature (eV). - Report ITER-IL-PH-13-9-e12 - - :param m_ions_total_amu: ion mass (assumes fuel only) (AMU) - :type m_ions_total_amu: float - - :param coefl: little 'l' in Harrison model - :type coefl: float - - :param delne: scrapeoff density by main plasma (10**20 m-3) - :type delne: float - - :param fififi: coeff. used in sheath energy transfer factor calc. - :type fififi: float - - :param omegan: pressure ratio of (plate / main plasma) - :type omegan: float - - :param omlarg: factor accounting for power flow - :type omlarg: float - - :param qdiv: heat flux across separatrix to divertor (MW/m2) - :type qdiv: float - - :param tconl: connection length along field line by main plasma (m) - :type tconl: float - - :param xpara: parallel diffusivity in the plasma scrapeoff (m2/s) - :type xpara: float - - :param xperp: perpend. diffusivity in the plasma scrapeoff (m2/s) - :type xperp: float - - :param xx: T_plate / T_separatrix guess - :type xx: float - - :param yy: T_plate guess (eV) - :type yy: float - - :returns: an estimate for the divertor temperature (eV) - :rtype: float - """ - - c27 = 0.28571428e0 - - xxs = max(xx, 0.001e0) - xxs = min(xxs, 0.99999e0) - yys = max(yy, 0.1e0) - - dendiv = delne * omegan / xxs - gamdiv = self.gammash(fififi, yys) - eier = self.erprcy(yys, dendiv) - ff = ( - 20.16e0 - * m_ions_total_amu - * ( - (qdiv) ** 10 - * (c27 * xpara) ** 4 - / (xperp**5 * delne**14) - * tconl - * (1.0e0 - xxs**3.5e0) ** 4 - / coefl**4 - ) - ** 0.22222e0 - / (omegan * gamdiv * omlarg * (1.0e0 + eier / (gamdiv * yys))) ** 2 - ) - - return yys - ff - - def ftpts( - self, - m_ions_total_amu: float, - coefl: float, - delne: float, - fififi: float, - omegan: float, - omlarg: float, - qdiv: float, - tconl: float, - xpara: float, - xperp: float, - xx: float, - yy: float, - ) -> float: - """Function for divertor model temperature ratio solution - author: J Galambos, ORNL - author: P J Knight, CCFE, Culham Science Centre - - This function updates the divertor model temperature ratio solution. - Report ITER-IL-PH-13-9-e12 - - :param m_ions_total_amu: ion mass (assumes fuel only) (AMU) - :type m_ions_total_amu: float - - :param coefl: little 'l' in Harrison model - :type coefl: float - - :param delne: scrapeoff density by main plasma (10**20 m-3) - :type delne: float - - :param fififi: coeff. used in sheath energy transfer factor calc. - :type fififi: float - - :param omegan: pressure ratio of (plate / main plasma) - :type omegan: float - - :param omlarg: factor accounting for power flow - :type omlarg: float - - :param qdiv: heat flux across separatrix to divertor (MW/m2) - :type qdiv: float - - :param tconl: connection length along field line by main plasma (m) - :type tconl: float - - :param xpara: parallel diffusivity in the plasma scrapeoff (m2/s) - :type xpara: float - - :param xperp: perpend. diffusivity in the plasma scrapeoff (m2/s) - :type xperp: float - - :param xx: T_plate / T_separatrix guess - :type xx: float - - :param yy: T_plate guess (eV) - :type yy: float - - :returns: an estimate for the divertor temperature (eV) - :rtype: float - """ - - xxs = max(xx, 0.001e0) - xxs = min(xxs, 0.99999e0) - yys = max(yy, 0.1e0) - - dendiv = delne * omegan / xxs - gamdiv = self.gammash(fififi, yys) - eier = self.erprcy(yys, dendiv) - - ff = ( - xxs**3.5e0 - + 9.66e0 - * (xxs / m_ions_total_amu) ** 0.9e0 - * (xperp / (qdiv) ** 2) ** 0.8e0 - * coefl - / (2.0e0 * xpara / 7.0e0) - * tconl**0.2e0 - * (omlarg * omegan * gamdiv * (1.0e0 + eier / (gamdiv * yys))) ** 1.8e0 - * delne**2.6e0 - ) - - return 1.0e0 - ff - - def gammash(self, gcoef: float, tdiv: float) -> float: - """Function to provide the plasma sheath energy transfer coefficient - author: J Galambos, ORNL - author: P J Knight, CCFE, Culham Science Centre - - This function provides the plasma sheath energy transfer coefficient - from the Harrison / Kukushkin ITER model. - Report ITER-IL-PH-13-9-e12 - - :param gcoef: coefficient - :type gcoef: float - - :param tdiv: electron temperature at the plate (eV) - :type tdiv: float - - :returns: plasma sheath energy transfer coefficient - :rtype: float - """ - - return 8.3e0 - 6.0e0 * (0.07e0 - 0.18e0 * math.log10(3.0e0 * tdiv * gcoef)) - def divtart( self, rmajor: float, diff --git a/process/init.py b/process/init.py index a9acd9dc3a..c444238b44 100644 --- a/process/init.py +++ b/process/init.py @@ -283,7 +283,6 @@ def initialise_iterative_variables(): fortran.define_iteration_variables.init_itv_31() fortran.define_iteration_variables.init_itv_32() fortran.define_iteration_variables.init_itv_33() - fortran.define_iteration_variables.init_itv_34() fortran.define_iteration_variables.init_itv_35() fortran.define_iteration_variables.init_itv_36() fortran.define_iteration_variables.init_itv_37() diff --git a/process/input.py b/process/input.py index bfdcee41c7..ebf1df51e1 100644 --- a/process/input.py +++ b/process/input.py @@ -228,7 +228,6 @@ def __post_init__(self): "bootstrap_current_fraction_max": InputVariable( fortran.current_drive_variables, float, range=(-0.999, 0.999) ), - "bpsout": InputVariable(fortran.divertor_variables, float, range=(0.0, 10.0)), "breeder_f": InputVariable(fortran.fwbs_variables, float, range=(0.0, 1.0)), "breeder_multiplier": InputVariable( fortran.fwbs_variables, float, range=(0.0, 1.0) @@ -236,12 +235,6 @@ def __post_init__(self): "bz_channel_conduct_liq": InputVariable( fortran.fwbs_variables, float, range=(1e-06, 1000000.0) ), - "c1div": InputVariable(fortran.divertor_variables, float, range=(-100.0, 100.0)), - "c2div": InputVariable(fortran.divertor_variables, float, range=(-100.0, 100.0)), - "c3div": InputVariable(fortran.divertor_variables, float, range=(-100.0, 100.0)), - "c4div": InputVariable(fortran.divertor_variables, float, range=(-100.0, 100.0)), - "c5div": InputVariable(fortran.divertor_variables, float, range=(-100.0, 100.0)), - "c6div": InputVariable(fortran.divertor_variables, float, range=(-100.0, 100.0)), "casthi": InputVariable(fortran.tfcoil_variables, float, range=(0.0, 1.0)), "casthi_fraction": InputVariable(fortran.tfcoil_variables, float, range=(0.0, 1.0)), "casths": InputVariable(fortran.tfcoil_variables, float, range=(0.0, 1.0)), @@ -359,7 +352,6 @@ def __post_init__(self): "declfw": InputVariable(fortran.fwbs_variables, float, range=(0.01, 0.2)), "declshld": InputVariable(fortran.fwbs_variables, float, range=(0.01, 0.2)), "decomf": InputVariable(fortran.cost_variables, float, range=(0.0, 1.0)), - "delld": InputVariable(fortran.divertor_variables, float, range=(0.1, 2.0)), "denstl": InputVariable(fortran.fwbs_variables, float, range=(5000.0, 10000.0)), "dhecoil": InputVariable(fortran.tfcoil_variables, float, range=(0.0, 0.1)), "dintrt": InputVariable(fortran.cost_variables, float, range=(0.0, 0.1)), @@ -590,10 +582,8 @@ def __post_init__(self): "fcupfsu": InputVariable(fortran.pfcoil_variables, float, range=(0.0, 1.0)), "fcutfsu": InputVariable(fortran.tfcoil_variables, float, range=(0.0, 1.0)), "fdene": InputVariable(fortran.constraint_variables, float, range=(0.001, 10.0)), - "fdfs": InputVariable(fortran.divertor_variables, float, range=(0.0, 20.0)), "f_ster_div_single": InputVariable(fortran.fwbs_variables, float, range=(0.0, 1.0)), "fdiva": InputVariable(fortran.divertor_variables, float, range=(0.1, 2.0)), - "fdivcol": InputVariable(fortran.constraint_variables, float, range=(0.001, 10.0)), "fdivwet": InputVariable(fortran.stellarator_variables, float, range=(0.01, 1.0)), "fdtmp": InputVariable(fortran.constraint_variables, float, range=(0.001, 10.0)), "fecrh_ignition": InputVariable( @@ -607,7 +597,6 @@ def __post_init__(self): "fhldiv": InputVariable(fortran.constraint_variables, float, range=(0.001, 10.0)), "fhole": InputVariable(fortran.fwbs_variables, float, range=(0.0, 1.0)), "fhts": InputVariable(fortran.tfcoil_variables, float, range=(0.01, 1.0)), - "fififi": InputVariable(fortran.divertor_variables, float, range=(1e-06, 1.0)), "fiooic": InputVariable(fortran.constraint_variables, float, range=(0.001, 10.0)), "fipir": InputVariable(fortran.constraint_variables, float, range=(0.001, 10.0)), "fjohc": InputVariable(fortran.constraint_variables, float, range=(0.001, 10.0)), @@ -683,7 +672,6 @@ def __post_init__(self): "frholeg": InputVariable(fortran.tfcoil_variables, float, range=(0.01, 5.0)), "frminor": InputVariable(fortran.constraint_variables, float, range=(0.001, 10.0)), "frrmax": InputVariable(fortran.ife_variables, float, range=(1e-06, 1.0)), - "frrp": InputVariable(fortran.divertor_variables, float, range=(0.0, 1.0)), "fseppc": InputVariable( fortran.build_variables, float, range=(1000000.0, 1000000000.0) ), @@ -841,7 +829,6 @@ def __post_init__(self): ), "kappa": InputVariable(fortran.physics_variables, float, range=(0.99, 5.0)), "kappa95": InputVariable(fortran.physics_variables, float, range=(0.99, 5.0)), - "ksic": InputVariable(fortran.divertor_variables, float, range=(0.0, 2.0)), "layer_ins": InputVariable(fortran.tfcoil_variables, float, range=(0.0, 0.1)), "ld_ratio_cst": InputVariable(fortran.pfcoil_variables, float, range=(0.0, 5.0)), "len_fw_channel": InputVariable( @@ -936,7 +923,6 @@ def __post_init__(self): "f_z_cs_tf_internal": InputVariable( fortran.pfcoil_variables, float, range=(0.0, 2.0) ), - "omegan": InputVariable(fortran.divertor_variables, float, range=(0.1, 10.0)), "outgasfactor": InputVariable( fortran.vacuum_variables, float, range=(1e-06, 1000.0) ), @@ -1092,7 +1078,6 @@ def __post_init__(self): "rhopfbus": InputVariable(fortran.pfcoil_variables, float, range=(0.0, 1e-05)), "rinboard": InputVariable(fortran.build_variables, float, range=(0.1, 10.0)), "ripmax": InputVariable(fortran.tfcoil_variables, float, range=(0.1, 100.0)), - "rlenmax": InputVariable(fortran.divertor_variables, float, range=(0.0, 1.0)), "rmajor": InputVariable(fortran.physics_variables, float, range=(0.1, 50.0)), "robotics_h": InputVariable(fortran.buildings_variables, float, range=(1.0, 100.0)), "robotics_l": InputVariable( @@ -1429,9 +1414,7 @@ def __post_init__(self): ), "wsvfac": InputVariable(fortran.buildings_variables, float, range=(0.9, 3.0)), "xi_ebw": InputVariable(fortran.current_drive_variables, float, range=(0.0, 1.0)), - "xparain": InputVariable(fortran.divertor_variables, float, range=(0.01, 10000.0)), "xpertin": InputVariable(fortran.divertor_variables, float, range=(0.0, 10.0)), - "zeffdiv": InputVariable(fortran.divertor_variables, float, range=(0.01, 100.0)), "zeffmax": InputVariable(fortran.constraint_variables, float, range=(1.0, 10.0)), "blktmodel": InputVariable(fortran.fwbs_variables, int, choices=[0, 1]), "blkttype": InputVariable(fortran.fwbs_variables, int, choices=[1, 2, 3]), @@ -1444,7 +1427,6 @@ def __post_init__(self): fortran.rebco_variables, float, range=(1.0e6, 1.0e10) ), "cost_model": InputVariable(fortran.cost_variables, int, choices=[0, 1]), - "divdum": InputVariable(fortran.divertor_variables, int, choices=[0, 1]), "dwell_pump": InputVariable(fortran.vacuum_variables, int, choices=[0, 1, 2]), "i_fw_blkt_vv_shape": InputVariable(fortran.fwbs_variables, int, range=(1, 2)), "hcdportsize": InputVariable(fortran.fwbs_variables, int, range=(1, 2)), diff --git a/process/io/obsolete_vars.py b/process/io/obsolete_vars.py index 4fb5d0fbe6..f1926e71d7 100644 --- a/process/io/obsolete_vars.py +++ b/process/io/obsolete_vars.py @@ -273,6 +273,24 @@ "ipump": "i_fw_blkt_shared_coolant", "coolwh": "i_blkt_coolant_type", "iblanket_thickness": None, + "bpsout": None, + "c1div": None, + "c2div": None, + "c3div": None, + "c4div": None, + "c5div": None, + "c6div": None, + "delld": None, + "fdfs": None, + "fdivcol": None, + "fififi": None, + "frrp": None, + "ksic": None, + "omegan": None, + "rlenmax": None, + "xparain": None, + "zeffdiv": None, + "divdum": None, } OBS_VARS_HELP = { diff --git a/process/utilities/errorlist.json b/process/utilities/errorlist.json index a904309d60..790e476653 100644 --- a/process/utilities/errorlist.json +++ b/process/utilities/errorlist.json @@ -8,7 +8,7 @@ "comment2": [ "Increment n_errortypes if an error is added to this list" ], - "n_errortypes": 288, + "n_errortypes": 289, "errors": [ { "no": 1, @@ -1449,6 +1449,11 @@ "no": 288, "level": 3, "message": "PFCOIL: i_pf_location(i) should not be 1 if iohcl=0" + }, + { + "no": 289, + "level": 3, + "message": "CHECK: Constraint equation 22 has been deprecated." } ] } diff --git a/source/fortran/constraint_equations.f90 b/source/fortran/constraint_equations.f90 index 3858412ab6..e9032293e8 100755 --- a/source/fortran/constraint_equations.f90 +++ b/source/fortran/constraint_equations.f90 @@ -148,7 +148,7 @@ subroutine constraint_eqns(m,ieqn,cc,con,err,symbol,units) case (20); call constraint_eqn_020(tmp_cc, tmp_con, tmp_err, tmp_symbol, tmp_units) ! Equation for minor radius lower limit case (21); call constraint_eqn_021(tmp_cc, tmp_con, tmp_err, tmp_symbol, tmp_units) - ! Equation for divertor collision/connection length ratio upper limit + ! Obsolete case (22); call constraint_eqn_022(tmp_cc, tmp_con, tmp_err, tmp_symbol, tmp_units) ! Equation for conducting shell radius / rminor upper limit case (23); call constraint_eqn_023(tmp_cc, tmp_con, tmp_err, tmp_symbol, tmp_units) @@ -1195,19 +1195,22 @@ subroutine constraint_eqn_022(tmp_cc, tmp_con, tmp_err, tmp_symbol, tmp_units) !! fdivcol : input real : f-value for divertor collisionality !! rlenmax : input real : maximum value for length ratio (rlclolcn) !! rlclolcn : input real : ratio of collision length / connection length - use constraint_variables, only: fdivcol - use divertor_variables, only: rlenmax, rlclolcn + implicit none + real(dp), intent(out) :: tmp_cc real(dp), intent(out) :: tmp_con real(dp), intent(out) :: tmp_err character(len=1), intent(out) :: tmp_symbol character(len=10), intent(out) :: tmp_units + !! Constraints output - tmp_cc = rlclolcn / rlenmax - 1.0D0 * fdivcol - tmp_con = rlenmax * (1.0D0 - tmp_cc) - tmp_err = rlclolcn * tmp_cc - tmp_symbol = '<' + ! This constraint is depreciated + call report_error(289) + + tmp_con = 1.0D0 + tmp_err = 0.0D0 + tmp_symbol = '=' tmp_units = '' end subroutine constraint_eqn_022 diff --git a/source/fortran/constraint_variables.f90 b/source/fortran/constraint_variables.f90 index f96f696dea..ffb4df99ec 100644 --- a/source/fortran/constraint_variables.f90 +++ b/source/fortran/constraint_variables.f90 @@ -54,9 +54,6 @@ module constraint_variables !! f-value for density limit (`constraint equation 5`, `iteration variable 9`) !! (invalid if `ipedestal=3`) - real(dp) :: fdivcol - !! f-value for divertor collisionality (`constraint equation 22`, `iteration variable 34`) - real(dp) :: fdtmp !! f-value for first wall coolant temperature rise !! (`constraint equation 38`, `iteration variable 62`) @@ -332,7 +329,6 @@ subroutine init_constraint_variables fcpttf = 1.0D0 fr_conducting_wall = 1.0D0 fdene = 1.0D0 - fdivcol = 1.0D0 fdtmp = 1.0D0 fflutf = 1.0D0 ffuspow = 1.0D0 diff --git a/source/fortran/divertor_variables.f90 b/source/fortran/divertor_variables.f90 index c150eade29..654b5fd1c5 100644 --- a/source/fortran/divertor_variables.f90 +++ b/source/fortran/divertor_variables.f90 @@ -14,9 +14,6 @@ module divertor_variables public - real(dp) :: adas - !! area divertor / area main plasma (along separatrix) - real(dp) :: anginc !! angle of incidence of field line on plate (rad) @@ -29,49 +26,12 @@ module divertor_variables real(dp) :: betao !! poloidal plane angle between divertor plate and leg, outboard (rad) - real(dp) :: bpsout - !! reference B_p at outboard divertor strike point (T) - - real(dp) :: c1div - !! fitting coefficient to adjust ptpdiv, ppdiv - - real(dp) :: c2div - !! fitting coefficient to adjust ptpdiv, ppdiv - - real(dp) :: c3div - !! fitting coefficient to adjust ptpdiv, ppdiv - - real(dp) :: c4div - !! fitting coefficient to adjust ptpdiv, ppdiv - - real(dp) :: c5div - !! fitting coefficient to adjust ptpdiv, ppdiv - - real(dp) :: c6div - !! fitting coefficient to adjust ptpdiv, ppdiv - - real(dp) :: delld - !! coeff for power distribution along main plasma - - real(dp) :: dendiv - !! plasma density at divertor (10**20 /m3) - - real(dp) :: densin - !! density at plate (on separatrix) (10**20 /m3) - real(dp) :: divclfr !! divertor coolant fraction real(dp) :: divdens !! divertor structure density (kg/m3) - integer :: divdum - !! switch for divertor Zeff model: - !! - !! - =0 calc - !! - =1 input - !#TODO: switch name should be changed to i_ - real(dp) :: dz_divertor !! divertor structure vertical thickness (m) @@ -84,25 +44,12 @@ module divertor_variables real(dp) :: divsur !! divertor surface area (m2) - real(dp) :: fdfs - !! radial gradient ratio - real(dp) :: fdiva !! divertor area fudge factor (for ITER, Sept 1990) - real(dp) :: fhout - !! fraction of power to outboard divertor (for single null) - - real(dp) :: fififi - !! coefficient for gamdiv - !#TODO: what the hell is this variable name... - real(dp) :: flux_exp !! The plasma flux expansion in the divertor (default 2; Wade 2020) - real(dp) :: frrp - !! fraction of radiated power to plate - real(dp) :: hldiv !! divertor heat load (MW/m2) @@ -116,39 +63,13 @@ module divertor_variables real(dp) :: hldivlim !! heat load limit (MW/m2) - real(dp) :: ksic - !! power fraction for outboard double-null scrape-off plasma - - real(dp) :: lamp - !! power flow width (m) - - real(dp) :: minstang - !! minimum strike angle for heat flux calculation - - real(dp) :: omegan - !! pressure ratio (nT)_plasma / (nT)_scrape-off - - real(dp) :: omlarg - !! power spillage to private flux factor - - real(dp) :: ppdivr - !! peak heat load at plate (with radiation) (MW/m2) - real(dp) :: prn1 !! n-scrape-off / n-average plasma; (input for `ipedestal=0`, = nesep/dene if `ipedestal>=1`) - real(dp) :: ptpdiv - !! peak temperature at the plate (eV) real(dp) :: rconl !! connection length ratio, outboard side - real(dp) :: rlclolcn - !! ratio of collision length / connection length - - real(dp) :: rlenmax - !! maximum value for length ratio (rlclolcn) (`constraintg eqn 22`) - real(dp) :: rsrd !! effective separatrix/divertor radius ratio @@ -158,72 +79,35 @@ module divertor_variables real(dp) :: tdiv !! temperature at divertor (eV) (input for stellarator only, calculated for tokamaks) - real(dp) :: tsep - !! temperature at the separatrix (eV) - - real(dp) :: xparain - !! parallel heat transport coefficient (m2/s) - real(dp) :: xpertin !! perpendicular heat transport coefficient (m2/s) - real(dp) :: zeffdiv - !! Zeff in the divertor region (if `divdum/=0`) - contains subroutine init_divertor_variables !! Initialise divertor_variables implicit none - adas = 0.0D0 anginc = 0.262D0 beta_div = 1.0D0 betai = 1.0D0 betao = 1.0D0 - bpsout = 0.60D0 - c1div = 0.45D0 - c2div = -7.0D0 - c3div = 0.54D0 - c4div = -3.6D0 - c5div = 0.7D0 - c6div = 0.0D0 - delld = 1.0D0 - dendiv = 0.0D0 - densin = 0.0D0 divclfr = 0.3D0 divdens = 1.0D4 - divdum = 0 dz_divertor = 0.2D0 divmas = 0.0D0 divplt = 0.035D0 divsur = 0.0D0 - fdfs = 10.0D0 fdiva = 1.11D0 - fhout = 0.0D0 - fififi = 4.0D-3 flux_exp = 2.0D0 - frrp = 0.4D0 hldiv = 0.0D0 - i_hldiv = 0 + i_hldiv = 2 hldivlim = 5.0D0 - ksic = 0.8D0 - lamp = 0.0D0 - minstang = 0.0D0 - omegan = 1.0D0 - omlarg = 0.0D0 - ppdivr = 0.0D0 prn1 = 0.285D0 - ptpdiv = 0.0D0 rconl = 0.0D0 - rlclolcn = 0.0D0 - rlenmax = 0.5D0 rsrd = 0.0D0 tconl = 0.0D0 tdiv = 2.0D0 - tsep = 0.0D0 - xparain = 2.1D3 xpertin = 2.0D0 - zeffdiv = 1.0D0 end subroutine init_divertor_variables end module divertor_variables diff --git a/source/fortran/iteration_variables.f90 b/source/fortran/iteration_variables.f90 index 30e8d84628..34540907ed 100755 --- a/source/fortran/iteration_variables.f90 +++ b/source/fortran/iteration_variables.f90 @@ -788,28 +788,11 @@ end subroutine set_itv_33 !--------------------------------- - subroutine init_itv_34 - !!
  • (34) fdivcol (f-value for equation 22) - use numerics, only: lablxc, boundl, boundu - implicit none - lablxc(34) = 'fdivcol ' - boundl(34) = 0.001D0 - boundu(34) = 1.000D0 - end subroutine init_itv_34 - real(kind(1.d0)) function itv_34() - use constraint_variables, only: fdivcol implicit none - itv_34 = fdivcol + write(*,*) 'Iteration variable 34 is no longer in use.' end function itv_34 - subroutine set_itv_34(ratio) - use constraint_variables, only: fdivcol - implicit none - real(kind(1.d0)) :: ratio - fdivcol = ratio - end subroutine set_itv_34 - !--------------------------------- subroutine init_itv_35 @@ -4021,7 +4004,7 @@ subroutine loadxc case (31); xcm(i) = itv_31() case (32); xcm(i) = itv_32() case (33); xcm(i) = itv_33() - case (34); xcm(i) = itv_34() + case (34); case (35); xcm(i) = itv_35() case (36); xcm(i) = itv_36() case (37); xcm(i) = itv_37() @@ -4286,7 +4269,7 @@ subroutine convxc(xc,nn) case (31); call set_itv_31(ratio) case (32); call set_itv_32(ratio) case (33); call set_itv_33(ratio) - case (34); call set_itv_34(ratio) + case (34); case (35); call set_itv_35(ratio) case (36); call set_itv_36(ratio) case (37); call set_itv_37(ratio) diff --git a/source/fortran/numerics.f90 b/source/fortran/numerics.f90 index 55d872c6d9..707ed262e0 100755 --- a/source/fortran/numerics.f90 +++ b/source/fortran/numerics.f90 @@ -239,7 +239,7 @@ module numerics !!
  • (31) gapomin !!
  • (32) frminor (f-value for equation 21) !!
  • (33) fportsz (f-value for equation 20) - !!
  • (34) fdivcol (f-value for equation 22) + !!
  • (34) NOT USED !!
  • (35) fpeakb (f-value for equation 25) !!
  • (36) fbeta_max (f-value for equation 24) !!
  • (37) j_cs_flat_top_end diff --git a/tests/integration/ref_dicts.json b/tests/integration/ref_dicts.json index b4e902fdd7..ef037bbfab 100644 --- a/tests/integration/ref_dicts.json +++ b/tests/integration/ref_dicts.json @@ -209,7 +209,6 @@ ".false.", ".false." ], - "adas": 0.0, "adivflnc": 7.0, "adj_appr": [ "", @@ -1152,7 +1151,6 @@ 0.0, 0.0 ], - "bpsout": 0.6, "breeder": "'Orthosilicate'", "breeder_f": 0.5, "breeder_multiplier": 0.75, @@ -1178,7 +1176,6 @@ "len_blkt_inboard_channel_total": 0.0, "len_blkt_outboard_channel_total": 0.0, "c": "2.99792458e10_rkind", - "c1div": 0.45, "c21": 0.0, "c211": 0.0, "c212": 0.0, @@ -1271,11 +1268,6 @@ "c245": 0.0, "c25": 0.0, "c26": 0.0, - "c2div": -7.0, - "c3div": 0.54, - "c4div": -3.6, - "c5div": 0.7, - "c6div": 0.0, "c_": 299792500.0, "calpha": null, "capcost": 0.0, @@ -1615,18 +1607,15 @@ "deflect": 0.0, "degrad": 0.01745329251, "degree": "pi/180.0D0", - "delld": 1.0, "delta": null, "delta_eta": 0.0, "delta_var": null, - "dendiv": 0.0, "dene": 9.8e+19, "denh2o": 985.0, "denhe": 5.4, "nd_fuel_ions": 0.0, "denpbli": 9839.0, "densbreed": 0.0, - "densin": 0.0, "denstl": 7800.0, "denwc": 15630.0, "dewmkg": 0.0, @@ -1646,7 +1635,6 @@ "divclfr": 0.3, "divcst": 0.0, "divdens": 10000.0, - "divdum": 0.0, "dz_divertor": 0.2, "divlife": 0.0, "divmas": 0.0, @@ -1907,7 +1895,6 @@ "fcwr": 1.0, "fdene": 1.0, "f_deuterium": 0.5, - "fdfs": 10.0, "fdiags": [ "FLT_DEFAULT", "FLT_DEFAULT", @@ -1920,7 +1907,6 @@ ], "f_ster_div_single": 0.115, "fdiva": 1.11, - "fdivcol": 1.0, "fdivwet": 0.333333333333333, "fdtmp": 1.0, "fe2": null, @@ -1940,9 +1926,7 @@ "f_helium3": 0.0, "fhldiv": 1.0, "fhole": 0.0, - "fhout": 0.0, "fhts": 0.5, - "fififi": 0.004, "figmer": 0.0, "fileprefix": "", "fimp": [ @@ -2031,7 +2015,6 @@ "frholeg": 1.0, "frminor": 1.0, "frrmax": 1.0, - "frrp": 0.4, "fseppc": 350000000.0, "fstr_wp": 1.0, "fstrcase": 1.0, @@ -2889,7 +2872,6 @@ "kh2o": 0.651, "ki": null, "kron": "reshape((/unit, zero, zero, zero, unit, zero, zero, zero, unit/), (/3, 3/))", - "ksic": 0.8, "lablcc": [ "Betaconsistency", "Globalpowerbalanceconsistency", @@ -3198,7 +3180,6 @@ "lambda_q_VV": 6.92, "lambda_q_omp": 0.002, "lambdaio": 0.0, - "lamp": 0.0, "latentheat": 2257000.0, "layer_ins": 0.0, "lb": null, @@ -3306,7 +3287,6 @@ "mhd": null, "mi": 0.0, "minmax": 7.0, - "minstang": 0.0, "ml": null, "mlp1": null, "mode": null, @@ -3448,8 +3428,6 @@ "dr_cs": 0.811, "f_z_cs_tf_internal": 0.71, "om": null, - "omegan": 1.0, - "omlarg": 0.0, "one": 1.0, "opt_file": 15.0, "outfile": 0.0, @@ -3597,7 +3575,6 @@ "fusion_power": 0.0, "p_cs_resistive_flat_top": 0.0, "p_pf_coil_resistive_total_flat_top": 0.0, - "ppdivr": 0.0, "len_plasma_poloidal": 0.0, "ppump": 0.0, "ppumpmw": 0.0, @@ -3673,7 +3650,6 @@ 0.0, 0.0 ], - "ptpdiv": 0.0, "p_electron_transport_loss_mw": 0.0, "pden_electron_transport_loss_mw": 0.0, "p_ion_transport_loss_mw": 0.0, @@ -3941,7 +3917,6 @@ 0.0, 0.0 ], - "rconl": 0.0, "rcool": 0.005, "r_cryostat_inboard": 0.0, "reactor_clrnc": 4.0, @@ -4228,8 +4203,6 @@ "j_cs_critical_pulse_start": 0.0, "j_pf_wp_critical": 0.0, "rkind": "double", - "rlclolcn": 0.0, - "rlenmax": 0.5, "ind_plasma_internal_norm": 0.9, "ind_plasma": 0.0, "rmaj": null, @@ -4276,7 +4249,6 @@ "rsldi": 0.0, "rsldo": 0.0, "rspo": 0.0, - "rsrd": 0.0, "rtanbeam": 0.0, "rtanmax": 0.0, "run_tests": 0.0, @@ -7487,7 +7459,6 @@ "t_burn": 1000.0, "t_burn_0": 0.0, "tcomrepl": 0.5, - "tconl": 0.0, "tcoolin": 313.15, "temp_cp_average": 373.15, "tcpav2": 0.0, @@ -7638,7 +7609,6 @@ "trithtmw": 15.0, "tritprate": 0.0, "triv": 40000.0, - "tsep": 0.0, "ttarget": 5.0, "tturb": 0.0, "turbine_hall_h": 15.0, @@ -8555,7 +8525,6 @@ "xk0": null, "xlabel": "", "xlabel_2": "", - "xparain": 2100.0, "xpertin": 2.0, "yarc": [ 0.0, @@ -8579,7 +8548,6 @@ "zeff": 0.0, "zeff_div": 0.0, "zeffai": 0.0, - "zeffdiv": 1.0, "zeffmax": 3.6, "zero": 0.0, "zfxf": [ @@ -8783,7 +8751,6 @@ "acs": "Area of space inside conductor (m2)", "acstf": "Cable space area (per turn) [m2]\n Includes the area of voids and central helium channel", "active_constraints": "active_constraints(ipeqns) : Logical array showing which constraints are active", - "adas": "area divertor / area main plasma (along separatrix)", "adivflnc": "allowable divertor heat fluence (MW-yr/m2)", "adj_appr": "", "admv": "administration building volume (m3)", @@ -8936,7 +8903,6 @@ "bp_ratio_ob": "", "b_pf_coil_peak": "peak field at coil i (T)", "bpf2": "", - "bpsout": "reference B_p at outboard divertor strike point (T)", "breeder": "Breeder type (allowed values are Orthosilicate, Metatitanate or Zirconate)", "breeder_f": "Volume ratio: Li4SiO4/(Be12Ti+Li4SiO4) (`iteration variable 108`)", "breeder_multiplier": "combined breeder/multipler fraction of blanket by volume", @@ -8962,7 +8928,6 @@ "len_blkt_inboard_channel_total": "Inboard/outboard blanket flow lengths (m)", "len_blkt_outboard_channel_total": "Inboard/outboard blanket flow lengths (m)", "c": "", - "c1div": "fitting coefficient to adjust ptpdiv, ppdiv", "c21": "", "c211": "", "c212": "", @@ -9055,11 +9020,6 @@ "c245": "", "c25": "", "c26": "", - "c2div": "fitting coefficient to adjust ptpdiv, ppdiv", - "c3div": "fitting coefficient to adjust ptpdiv, ppdiv", - "c4div": "fitting coefficient to adjust ptpdiv, ppdiv", - "c5div": "fitting coefficient to adjust ptpdiv, ppdiv", - "c6div": "fitting coefficient to adjust ptpdiv, ppdiv", "c_": "", "calpha": "", "capcost": "total capital cost including interest (M$)", @@ -9264,18 +9224,15 @@ "deflect": "TF coil radial deflection (displacement) radial distribution [m]", "degrad": "degrees to radians, = pi/180", "degree": "", - "delld": "coeff for power distribution along main plasma", "delta": "", "delta_eta": "", "delta_var": "", - "dendiv": "plasma density at divertor (10**20 /m3)", "dene": "electron density (/m3) (`iteration variable 6`) (calculated if `ipedestal=3`)", "denh2o": "density of water (kg/m3)", "denhe": "", "nd_fuel_ions": "fuel ion density (/m3)", "denpbli": "", "densbreed": "density of breeder material (kg/m3) (`i_blanket_type=2` (KIT HCPB))", - "densin": "density at plate (on separatrix) (10**20 /m3)", "denstl": "density of steel (kg/m3)", "denwc": "density of tungsten carbide (kg/m3)", "dewmkg": "total mass of vacuum vessel + cryostat (kg) (calculated if blktmodel>0)", @@ -9295,7 +9252,6 @@ "divclfr": "divertor coolant fraction", "divcst": "divertor direct cost (M$)", "divdens": "divertor structure density (kg/m3)", - "divdum": "switch for divertor Zeff model:\n", "dz_divertor": "divertor structure vertical thickness (m)", "divlife": "Full power lifetime of divertor (y)", "divmas": "divertor plate mass (kg)", @@ -9520,11 +9476,9 @@ "fcwr": "f-value for conducting wall radius / rminor limit\n (`constraint equation 23`, `iteration variable 104`)", "fdene": "f-value for density limit (`constraint equation 5`, `iteration variable 9`)\n (invalid if `ipedestal=3`)", "f_deuterium": "deuterium fuel fraction", - "fdfs": "radial gradient ratio", "fdiags": "", "f_ster_div_single": "Solid angle fraction taken by one divertor", "fdiva": "divertor area fudge factor (for ITER, Sept 1990)", - "fdivcol": "f-value for divertor collisionality (`constraint equation 22`, `iteration variable 34`)", "fdivwet": "wetted fraction of the divertor area", "fdtmp": "f-value for first wall coolant temperature rise\n (`constraint equation 38`, `iteration variable 62`)", "fe2": "", @@ -9544,9 +9498,7 @@ "f_helium3": "helium-3 fuel fraction", "fhldiv": "f-value for divertor heat load (`constraint equation 18`, `iteration variable 27`)", "fhole": "area fraction taken up by other holes (IFE)", - "fhout": "fraction of power to outboard divertor (for single null)", "fhts": "technology adjustment factor for critical current density fit for isumat..=2\n Bi-2212 superconductor, to describe the level of technology assumed (i.e. to\n account for stress, fatigue, radiation, AC losses, joints or manufacturing\n variations; 1.0 would be very optimistic)", - "fififi": "coefficient for gamdiv", "figmer": "physics figure of merit (= plasma_currentaspect*sbar, where `sbar=1`)", "fileprefix": "input file prefix", "fimp": "", @@ -9620,7 +9572,6 @@ "frholeg": "Ouboard legs resistivity enhancement factor. Only used for `itart=1`.", "frminor": "f-value for minor radius limit (`constraint equation 21`, `iteration variable 32`)", "frrmax": "f-value for maximum IFE repetition rate (`constraint equation 50`, `iteration variable 86`)", - "frrp": "fraction of radiated power to plate", "fseppc": "Separation force in CS coil pre-compression structure", "fstr_wp": "f-value for maxiumum TF coil strain absolute value\n (`constraint equation 88`, `iteration variable 165`)", "fstrcase": "f-value for maximum TF coil case Tresca yield criterion\n (`constraint equation 31`, `iteration variable 48`)", @@ -9947,7 +9898,6 @@ "kh2o": "thermal conductivity of water (W/m/K)", "ki": "", "kron": "", - "ksic": "power fraction for outboard double-null scrape-off plasma", "lablcc": "lablcc(ipeqns) : labels describing constraint equations (corresponding itvs)\n\n\n\n", "lablmm": "lablmm(ipnfoms) : labels describing figures of merit:\n\n\n", "lablxc": "lablxc(ipnvars) : labels describing iteration variables