def runid0000(self):
diff = self.data.diff
prettyID = [r'$D^{pinch}_j$', r'$D^{std}_{j}$']
title = r"$\bf{D}_j$ Considering Non-diffusive Effects"
adjustments = {}
caption = [
r"Main ion particle diffusion coefficient in the edge for",
r"DIII-D shot " + str(self.shotid) + "." + str(self.tid) +
" considering IOL and pinch",
r"velocity, $D^{pinch}_{j}$, compared with standard",
r"definition, $D^{std}_j.$"
]
prettyCompare(
('rhor', self.data.rhor),
[(prettyID[0], diff[0]), (prettyID[1], diff[1])],
# yrange=[-.5,.5],
datalabels=[prettyID[0], prettyID[1]],
title=title,
ylabel=[r"$D_j$", r"$\left[\frac{m^2}{s}\right]$"],
caption=caption,
adjust=adjustments,
# textSpace=.025,
# marginBottom=.2,
# yLabelAdj=-.6,
# marginLeft=.15,
size=(6, 9))
print self.data.rhor
adjustments = {}
caption = [
"Main ion particle diffusion coefficient in the core for",
"DIII-D shot 000000.0000 considering IOL and pinch velocity"
]
prettyCompare(
('rhor', self.data.rhor[:-30]),
[(prettyID[0], diff[0][:-30])],
corePlot=True,
# yrange=[-.5,.5],
datalabels=[prettyID[0]],
title=title,
ylabel=[r"$D_j$", r"$\left[\frac{m^2}{s}\right]$"],
caption=caption,
adjust=adjustments,
# textSpace=.02,
size=(6, 9))
def runid1570(self):
data=self.data
prettyID=[r'$\nu_{d,j}$']
title=r"Inferred $\nu_{d,j}$"
adjustments={}
caption=[r"Inferred $\nu_{d,j}$ for DIII-D shot "+str(self.shotid)+"."+str(self.tid),
"using pertrubation theory"]
prettyCompare(('rhor',self.data.rhor ),[(prettyID[0],data.nudrag)],
yrange=[0.,2.E4],
datalabels=[prettyID[0]],
title=title,
ylabel=[r"$\nu_{d,j}$",r'$\left[s^{-1}\right]$'],
caption=caption,
adjust=adjustments,
textSpace=.025,
marginBottom=.18,
yLabelAdj=-1.2,
marginLeft=.2,
size=(6,9),
toSN=4)
def runid1950(self):
chis, qs = self.chis, self.qs
prettyID = [
r'$Q^{cond}_j=Q^{total}_{j}$', # 0
r'$Q^{cond}_j=Q^{total}_j$ w/ IOL', # 1
r'$Q^{cond}_j=Q^{total}_j-Q^{conv}_j$', # 2
r'$Q^{cond}_j=Q^{total}_j-Q^{conv}_j-Q^{heatin}_j$', #3
r'$Q^{cond}_j=Q^{total}_j-Q^{conv}_j-Q^{heatin}_j-Q^{visc}_j$'
] #4
title = r"Comparison of Ion Heat Diffusivity Considering Various Non-diffusive Mechanisms"
adjustments = {0: .6, 2: -.60, 3: -1.9, 4: .25}
caption = [
"Main ion heat diffusivities in the edge for DIII-D shot 166606.1950 correcting for IOL, convective heat flux,",
"work done on the pressure tensor, and viscous heating (assuming 10% asymmetry)"
]
# caption='''
# Main ion heat diffusivities in the edge
# for DIII-D shot # 166606.1950 correcting
# for IOL,convective heat flux, work done
# on the pressure tensor, and viscous
# heating (assuming 10% asymmetry).'''
prettyCompare(('rhor', self.data.rhor), [(prettyID[0], chis[0]),
(prettyID[1], chis[1]),
(prettyID[2], chis[2]),
(prettyID[3], chis[3]),
(prettyID[4], chis[4])],
yrange=[-.4, .1],
datalabels=[
prettyID[0], prettyID[1], prettyID[2], prettyID[3],
prettyID[4]
],
title=title,
ylabel=[r"$\chi_j$", r'$\left[\frac{m^2}{s}\right]$'],
caption=caption,
adjust=adjustments,
xTickremove=None,
textSpace=.095,
marginBottom=.15,
marginLeft=0.075,
size=(16, 12))
adjustments = {}
caption = [
"Main ion heat diffusivities in the core for DIII-D shot 166606.1950 correcting for IOL, convective heat flux,",
"work done on the pressure tensor, and viscous heating (assuming 10% asymmetry)"
]
prettyCompare(
('rhor', self.data.rhor[:-30]), [(prettyID[0], chis[0][:-30]),
(prettyID[1], chis[1][:-30]),
(prettyID[2], chis[2][:-30]),
(prettyID[3], chis[3][:-30]),
(prettyID[4], chis[4][:-30])],
corePlot=True,
yrange=[-1., 5.],
datalabels=[
prettyID[0], prettyID[1], prettyID[2], prettyID[3], prettyID[4]
],
title=title,
ylabel=[r"$\chi_j$", r'$\left[\frac{m^2}{s}\right]$'],
caption=caption,
adjust=adjustments,
textSpace=.22)
prettyID = [
r'$Q^{total}_j$', r'$Q^{total}_j$ w/ IOL', r'$Q^{conv15}_j$',
r'$Q^{heatin}_j$', r'$Q^{visc}_j$'
]
caption = [
"Heat flux terms for DIII-D shot 166606.1950 for the edge plasma"
]
title = r"$\bf{Q}$ term comparison"
adjustments = {1: .95, 3: .75}
prettyCompare(('rhor', self.data.rhor), [(prettyID[0], qs[0]),
(prettyID[1], qs[1]),
(prettyID[2], qs[2]),
(prettyID[3], qs[3]),
(prettyID[4], qs[4])],
yrange=[-2.5E4, 2.5E4],
datalabels=[
prettyID[0], prettyID[1], prettyID[2], prettyID[3],
prettyID[4]
],
title=title,
ylabel=[r"$q$", r"$\left[\frac{W}{m^2}\right]$"],
caption=caption,
adjust=adjustments,
textSpace=.02)
prettyID = [
r'$Q^{total}_j$', r'$Q^{total}_j$ w/ IOL', r'$Q^{conv15}_j$',
r'$Q^{heatin}_j$', r'$Q^{visc}_j$'
]
caption = [
"Heat flux terms for DIII-D shot 166606.1950 for the core plasma"
]
title = r"$\bf{Q}$ term comparison"
adjustments = {}
prettyCompare(
('rhor', self.data.rhor[:-30]), [(prettyID[0], qs[0][:-30]),
(prettyID[1], qs[1][:-30]),
(prettyID[2], qs[2][:-30]),
(prettyID[3], qs[3][:-30]),
(prettyID[4], qs[4][:-30])],
corePlot=True,
yrange=[-1.E5, 1.E4],
datalabels=[
prettyID[0], prettyID[1], prettyID[2], prettyID[3], prettyID[4]
],
title=title,
ylabel=[r"$q$", r"$\left[\frac{W}{m^2}\right]$"],
caption=caption,
adjust=adjustments,
textSpace=.1)
def runid1560(self):
chis, qs = self.chis, self.qs
prettyID = [
r'$Q^{cond}_j=Q^{total}_{j}$', r'$Q^{cond}_j=Q^{total}_j$ w/ IOL',
r'$Q^{cond}_j=Q^{total}_j-Q^{conv}_j$'
] #,
# r'$Q^{cond}_j=Q^{total}_j-Q^{conv}_j-Q^{heatin}_j$',
# r'$Q^{cond}_j=Q^{total}_j-Q^{conv}_j-Q^{heatin}_j-Q^{visc}_j$']
title = r"Comparison of Ion Heat Diffusivity Considering Various Non-diffusive Mechanisms"
adjustments = {1: 1.4} #,2:2.8,3:-1.4}
caption = [
"Main ion heat diffusivities in the edge for DIII-D shot 118890.1560 correcting for IOL, convective heat flux",
"work done on the pressure tensor, and viscous heating (assuming 10% asymmetry)"
]
prettyCompare(
('rhor', self.data.rhor),
[(prettyID[0], chis[0]), (prettyID[1], chis[1]),
(prettyID[2], chis[2])
], #(prettyID[3],chis[3]),(prettyID[4],chis[4])],
yrange=[0., 5.],
datalabels=[prettyID[0], prettyID[1],
prettyID[2]], #,prettyID[3],prettyID[4]],
title=title,
ylabel=[r"$\chi_j$", r'', r'$\left[\frac{m^2}{s}\right]$'],
caption=caption,
adjust=adjustments,
xTickremove=None,
textSpace=.045,
marginBottom=.15,
marginLeft=0.1,
size=(16, 12))
adjustments = {}
caption = [
"Main ion heat diffusivities in the core for DIII-D shot 118890.1560 correcting for IOL, convective heat flux",
"work done on the pressure tensor, and viscous heating (assuming 10% asymmetry)"
]
prettyCompare(
('rhor', self.data.rhor[:-30]), [(prettyID[0], chis[0][:-30]),
(prettyID[1], chis[1][:-30]),
(prettyID[2], chis[2][:-30]),
(prettyID[3], chis[3][:-30]),
(prettyID[4], chis[4][:-30])],
corePlot=True,
yrange=[0., 2.5],
datalabels=[
prettyID[0], prettyID[1], prettyID[2], prettyID[3], prettyID[4]
],
title=title,
ylabel=[r"$\chi_j$", r'$\left[\frac{m^2}{s}\right]$'],
caption=caption,
adjust=adjustments,
textSpace=.22)
prettyID = [
r'$Q^{total}_j$', r'$Q^{total}_j$ w/ IOL', r'$Q^{conv15}_j$',
r'$Q^{heatin}_j$', r'$Q^{visc}_j$'
]
caption = [
"Heat flux terms for DIII-D shot 118890.1560 for the edge plasma"
]
title = r"$\bf{Q}$ term comparison"
adjustments = {}
prettyCompare(('rhor', self.data.rhor), [(prettyID[0], qs[0]),
(prettyID[1], qs[1]),
(prettyID[2], qs[2]),
(prettyID[3], qs[3]),
(prettyID[4], qs[4])],
yrange=[-0.5E4, 2.5E4],
datalabels=[
prettyID[0], prettyID[1], prettyID[2], prettyID[3],
prettyID[4]
],
title=title,
ylabel=[r"$q$", r"$\left[\frac{W}{m^2}\right]$"],
caption=caption,
adjust=adjustments,
textSpace=.02)
prettyID = [
r'$Q^{total}_j$', r'$Q^{total}_j$ w/ IOL', r'$Q^{conv15}_j$',
r'$Q^{heatin}_j$', r'$Q^{visc}_j$'
]
caption = [
"Heat flux terms for DIII-D shot 118890.1560 for the core plasma"
]
title = r"$\bf{Q}$ term comparison"
adjustments = {}
prettyCompare(
('rhor', self.data.rhor[:-30]), [(prettyID[0], qs[0][:-30]),
(prettyID[1], qs[1][:-30]),
(prettyID[2], qs[2][:-30]),
(prettyID[3], qs[3][:-30]),
(prettyID[4], qs[4][:-30])],
corePlot=True,
yrange=[-1.2E4, 2.5E4],
datalabels=[
prettyID[0], prettyID[1], prettyID[2], prettyID[3], prettyID[4]
],
title=title,
ylabel=[r"$q$", r"$\left[\frac{W}{m^2}\right]$"],
caption=caption,
adjust=adjustments,
textSpace=.1)