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)