Added on-axis temperature and density plotting

resources/default_plotsettings.yml

@@ -0,0 +1,305 @@
+#-----------------------------------------------------------------------
+# OpenPOPCON Plotting Settings
+#-----------------------------------------------------------------------
+
+yax: "nG" # y-axis variable; nG for greenwald fraction, n20 for density
+xax: "T_i_av" # x-axis variable; T_i_av for average ion temperature, 
+              #                  T_i_ax for maximum ion temperature,
+              #                  T_e_av for average electron temperature,
+              #                  T_e_ax for maximum electron temperature
+figsize: [8, 6] # figure size in inches
+
+#----------------------
+# Options per output
+#----------------------
+
+plotoptions: 
+    Paux: 
+        plot: True # do not plot this output
+        color: r # line color
+        linewidth: 3 # line width
+        label: $P_{aux}$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.2d' # format for the on-contour labels
+        spacing: manual # lin, log, or manual spacing for the contour levels
+        manuallevels: [0.1,1,10,50,100,500,2000] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 10 # number of contour levels
+    Pfusion: 
+        plot: True # do not plot this output
+        color: k # line color
+        linewidth: 2 # line width
+        label: $P_{fus}$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.2d' # format for the on-contour labels
+        spacing: manual # lin, log, or manual spacing for the contour levels
+        manuallevels: [100,200,500,1000,2000,4000] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 10 # number of contour levels
+    Pfusionheating: 
+        plot: False # do not plot this output
+        color: k # line color
+        linewidth: 2 # line width
+        label: $P_{fus}$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.2d' # format for the on-contour labels
+        spacing: manual # lin, log, or manual spacing for the contour levels
+        manuallevels: [100,200,500,1000,2000,4000] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 10 # number of contour levels
+    Psol:
+        plot: False # do not plot this output
+        color: k # line color
+        linewidth: 2 # line width
+        label: $P_{SOL}$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.2d' # format for the on-contour labels
+        spacing: manual # lin, log, or manual spacing for the contour levels
+        manuallevels: [10,20,50,80,100] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 10 # number of contour levels
+    Pohmic: 
+        plot: False # do not plot this output
+        color: k # line color
+        linewidth: 2 # line width
+        label: $P_{OH}$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.2f' # format for the on-contour labels
+        spacing: manual # lin, log, or manual spacing for the contour levels
+        manuallevels: [0.01,0.1,0.2,0.5,1.] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 10 # number of contour levels
+    Prad: 
+        plot: True # do not plot this output
+        color: purple # line color
+        linewidth: 2 # line width
+        label: $P_{rad}$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.2d' # format for the on-contour labels
+        spacing: manual # lin, log, or manual spacing for the contour levels
+        manuallevels: [10,50,100,200] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 3 # number of contour levels
+    Pheat: 
+        plot: False # do not plot this output
+        color: k # line color
+        linewidth: 2 # line width
+        label: $P_{heat}$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.2f' # format for the on-contour labels
+        spacing: lin # lin, log, or manual spacing for the contour levels
+        manuallevels: [] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 10 # number of contour levels
+    Ploss: 
+        plot: False # do not plot this output
+        color: k # line color
+        linewidth: 2 # line width
+        label: $P_{loss}$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.2f' # format for the on-contour labels
+        spacing: lin # lin, log, or manual spacing for the contour levels
+        manuallevels: [] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 10 # number of contour levels
+    Palpha: 
+        plot: False # do not plot this output
+        color: k # line color
+        linewidth: 2 # line width
+        label: $P_{\alpha}$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.2f' # format for the on-contour labels
+        spacing: lin # lin, log, or manual spacing for the contour levels
+        manuallevels: [] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 10 # number of contour levels
+    Pdd: 
+        plot: False # do not plot this output
+        color: k # line color
+        linewidth: 2 # line width
+        label: $P_{DD}$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.2f' # format for the on-contour labels
+        spacing: lin # lin, log, or manual spacing for the contour levels
+        manuallevels: [] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 10 # number of contour levels
+    Pdt: 
+        plot: False # do not plot this output
+        color: k # line color
+        linewidth: 2 # line width
+        label: $P_{DT}$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.2f' # format for the on-contour labels
+        spacing: lin # lin, log, or manual spacing for the contour levels
+        manuallevels: [] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 10 # number of contour levels
+    tauE: 
+        plot: False # do not plot this output
+        color: k # line color
+        linewidth: 2 # line width
+        label: $\tau_E$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.2f' # format for the on-contour labels
+        spacing: lin # lin, log, or manual spacing for the contour levels
+        manuallevels: [] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 10 # number of contour levels
+    Q: 
+        plot: True # do not plot this output
+        color: orange # line color
+        linewidth: 2 # line width
+        label: $Q$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.f' # format for the on-contour labels
+        spacing: manual # lin, log, or manual spacing for the contour levels
+        manuallevels: [0.5,1,2,5,10,25,50,100] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        min : -1 # only used if scale is specified
+        max : 2 # only used if scale is specified
+        levels: 4 # number of contour levels
+    H89: 
+        plot: False # do not plot this output
+        color: blue # line color
+        linewidth: 2 # line width
+        label: $H_{89}$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.2f' # format for the on-contour labels
+        spacing: lin # lin, log, or manual spacing for the contour levels
+        manuallevels: [] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 10 # number of contour levels
+    H98: 
+        plot: True # do not plot this output
+        color: blue # line color
+        linewidth: 2 # line width
+        label: $H_{98y2}$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.1f' # format for the on-contour labels
+        spacing: lin # lin, log, or manual spacing for the contour levels
+        manuallevels: [] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 5 # number of contour levels
+    vloop: 
+        plot: False # do not plot this output
+        color: k # line color
+        linewidth: 2 # line width
+        label: $V_{loop}$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.2f' # format for the on-contour labels
+        spacing: manual # lin, log, or manual spacing for the contour levels
+        manuallevels: [] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 10 # number of contour levels
+    betaN: 
+        plot: True # do not plot this output
+        color: forestgreen # line color
+        linewidth: 2 # line width
+        label: $\beta_N$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.d' # format for the on-contour labels
+        spacing: lin # lin, log, or manual spacing for the contour levels
+        manuallevels: [1,2,3] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 5 # number of contour levels
+    n_G_frac: 
+        plot: False # do not plot this output
+        color: k # line color
+        linewidth: 2 # line width
+        label: $n_G$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.2f' # format for the on-contour labels
+        spacing: lin # lin, log, or manual spacing for the contour levels
+        manuallevels: [] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 10 # number of contour levels
+    n20_max: 
+        plot: False # do not plot this output
+        color: k # line color
+        linewidth: 2 # line width
+        label: $n_{i20}(0)$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.2f' # format for the on-contour labels
+        spacing: lin # lin, log, or manual spacing for the contour levels
+        manuallevels: [] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 10 # number of contour levels
+    n20_avg: 
+        plot: False # do not plot this output
+        color: k # line color
+        linewidth: 2 # line width
+        label: $\langle n_{i20}\rangle$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.2f' # format for the on-contour labels
+        spacing: lin # lin, log, or manual spacing for the contour levels
+        manuallevels: [] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 10 # number of contour levels
+    T_i_max: 
+        plot: False # do not plot this output
+        color: k # line color
+        linewidth: 2 # line width
+        label: $T_i(0)$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.2f' # format for the on-contour labels
+        spacing: lin # lin, log, or manual spacing for the contour levels
+        manuallevels: [] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 10 # number of contour levels
+    T_i_avg: 
+        plot: False # do not plot this output
+        color: k # line color
+        linewidth: 2 # line width
+        label: $\langle T_i \rangle$ # label for the line
+        fontsize: 12 # font size for labels
+        fmt: '%.2f' # format for the on-contour labels
+        spacing: lin # lin, log, or manual spacing for the contour levels
+        manuallevels: [] # only used if manual spacing
+        scale: minmax # minmax or specified; below not used if manual spacing
+        # min : 0.0 # only used if scale is specified
+        # max : 1.0 # only used if scale is specified
+        levels: 10 # number of contour levels
+
+
+

resources/examples/MANTA/POPCON_input_example.yaml

@@ -10,22 +10,24 @@ delta: -0.5 # triangularity
 
 #----------------------
 # I_P and q* can be calculated from each other. Pick one to specify:
-# I_P: 10. # plasma current in MA
+I_P: 10. # plasma current in MA
 # or
-qstar: 2.5
+# qstar: 2.5
 #----------------------
 
 #----------------------
 # B_0 and B_coil can be calculated from each other. Pick one to specify:
-# B_0: 11.0 # on-axis magnetic field in Tesla
+B_0: 11.0 # on-axis magnetic field in Tesla
 # or
-B_coil: 24.0 # magnetic field at the coil in Tesla
-wall_thickness: 1.1 # wall thickness between magnet and separatrix in meters
+# B_coil: 24.0 # magnetic field at the coil in Tesla
+# wall_thickness: 1.1 # wall thickness between magnet and separatrix in meters
 #----------------------
 
+#----------------------
 M_i: 2.5 # average ion mass in amu
-tipeak_over_tepeak: 1.0 # ratio of ion to electron temperature at peak
+tipeak_over_tepeak: 1.1 # ratio of ion to electron temperature at peak
 fuel: 2 # 1 for D-D, 2 for D-T, 3 for D-He3
+#----------------------
 
 #----------------------
 # Impurity Fractions
@@ -43,23 +45,54 @@ Zeff_target: 2.1 # average effective charge of plasma
 impurity: 3 # 0 for He, 1 for Ne, 2 for Ar, 3 for Kr, 4 for Xe, 5 for W
 #----------------------
 
+#----------------------
 scalinglaw: "H_NT23" # Name of tau_E scaling law as defined in scaling_laws.yaml
 H_fac: 1.0 # Factor by which to multiply the scaling law
 nr: 100 # number of radial grid points
+#----------------------
+
+#----------------------
+# Profiles
+# If you specify a filename, the profiles will be read from that file and normalized.
+# Otherwise, the normalized profiles will be parabolic of the form
+# f(r) = (1 - offset)*(1 - r**alpha1)**alpha2 + offset
+# j profile comes from the equilibrium, the others come from the profs file.
+gfilename: "./resources/examples/MANTA/gMANTA"
+profsfilename: "./resources/examples/MANTA/profsMANTA.csv"
+
+# j_alpha1: 2.
+# j_alpha2: 3.
+# j_offset: 0.
+
+# ne_alpha1: 2.
+# ne_alpha2: 1.5
+# ne_offset: 0.
+
+# ni_alpha1: 2.
+# ni_alpha2: 1.5
+# ni_offset: 0.
+
+# Ti_alpha1: 2.
+# Ti_alpha2: 1.5
+# Ti_offset: 0.
+
+# Te_alpha1: 2.
+# Te_alpha2: 1.5
+# Te_offset: 0.
+#----------------------
+
 
 #----------------------
 # Algorithm settings
 #----------------------
 Nn: 20 # resolution of densities
 NTi: 20 # resolution of ion temperatures
-gfilename: "./resources/examples/MANTA/gMANTA"
-profsfilename: "./resources/examples/MANTA/profsMANTA.csv"
 nmax_frac: 1.2 # maximum density as a fraction of the Greenwald density
 nmin_frac: 0.3 # minimum density as a fraction of the Greenwald density
 Tmax_keV: 11. # maximum volume-averaged temperature in keV
 Tmin_keV: 3 # minimum volume-averaged temperature in keV
-maxit: 1000 # maximum number of iterations
-accel: 1.1 # acceleration factor
+maxit: 100000 # maximum number of iterations
+accel: 1.5 # acceleration factor
 err: 1e-6 # error tolerance
 verbosity: 1 # 0 for no output, 1 for minimal output, 2 for verbose output
 parallel: False