plot_proc.py fails if heating is not ecrh or nbi

[jonmaddock]

In GitLab by @asbrown on Oct 13, 2020, 15:00

In plot_proc.py the current drive function has options for ecrh and nbi, but if iefrf=4 (Ehst Lower Hybrid) is selected then plot_proc.py will fail as the array data is not created but later in the code data.append is used.

    if ecrh:
        data = [(pinjie, "Steady state auxiliary power", "MW"),
            ("pheat", "Power for heating only", "MW"),
            ("bootipf", "Bootstrap fraction", ""),
            ("faccd", "Auxiliary fraction", ""),
            ("facoh", "Inductive fraction", ""),
            ("powerht", "Plasma heating used for H factor", "MW"),
            (pdivr, r"$\frac{P_{\mathrm{div}}}{R_{0}}$", "MW m$^{-1}$"),
            (pdivnr, r"$\frac{P_{\mathrm{div}}}{<n> R_{0}}$",
             r"$\times 10^{-20}$ MW m$^{2}$"),
            (flh, r"$\frac{P_{\mathrm{div}}}{P_{\mathrm{LH}}}$", ""),
            (hstar, "H* (non-rad. corr.)", "")]

Another data array for Ehst Lower Hybrid to be made in the current drive function and an else statement to catch future current drive options.

[jonmaddock]

In GitLab by @skahn on Oct 14, 2020, 12:44

Out of curiosity, what current drive system have you been trying ?

[jonmaddock]

In GitLab by @asbrown on Oct 14, 2020, 15:14

iefrf=4 (Ehst Lower Hybrid), for STARFIRE

[jonmaddock]

In GitLab by @mkovari on Oct 15, 2020, 10:10

Adam, why??

[jonmaddock]

In GitLab by @asbrown on Oct 15, 2020, 10:44

Here are the pieces of the code that cause the problem:

def plot_current_drive_info(axis, mfile_data, scan):
...
    iefrf = mfile_data.data["iefrf"].get_scan(scan)
...
    if((iefrf == 5)or(iefrf==8)):
        nbi = True
        axis.text(-0.05, 1, 'Neutral Beam Current Drive:', ha='left', va='center')
    if((iefrf == 3)or(iefrf==7)or(iefrf==10)or(iefrf==11)):
        ecrh = True
        axis.text(-0.05, 1, 'Electron Cyclotron Current Drive:', ha='left', va='center')
...
    if ecrh:
        data = [(pinjie, "Steady state auxiliary power", "MW"),
...
    if nbi:
        data = [(pinjie, "Steady state auxiliary power", "MW"),
...
    if coe == 0.0 :   
        data.append(("", "", ""))

In the plot_current_drive_info function the current drive option (iefrf) is read in, which can be a value between 1 and 12.

If the iefrf is 3, 5, 7, 8, 10, 11 then either the variables nbi or ecrh are given a true flag. This means nbi and ecrh are both false wheniefrf is 1, 2, 4, 6, 9, or 12.

When either nbi or ecrh are set to true the array data is created. A subsequent append operation is done on data depending on the coe value. But, if iefrf is 1, 2, 4, 6, 9, 12 then data is not created and the append operation causes plot_proc.py to fail.

[grmtrkngtn]

This issue was resolved and merged, see pull request.