[yt-users] Colourbar tick labels

[Britton Smith]

Hi Libby,

Check out the powFormat in the plot_modifications repo.  It'll do that.

Britton

On Mon, Jun 6, 2011 at 11:36 AM, Elizabeth Harper-Clark <
h-clark at astro.utoronto.ca> wrote:

> Argh, sorry guys. First message didn't send as list didn't like the size of
> the attachment. here is the message:
>
>
> Hi guys,
>
> Quick question: how do I change how yt formats the number labels on the
> ticks on a colourbar? Rather than 0.0001, 0.001, 0.01, 0.1, and 1.0 x 1e28 I
> would rather the ticks were labelled 1e24, 1e25, 1e26, 1e27, and 1e28 (or
> even better 10 with power superscripted)!
>
> Thank you,
>
> Libby
>
>
>
>
> On 6 June 2011 13:32, Elizabeth Harper-Clark <h-clark at astro.utoronto.ca>wrote:
>
>> ps, here is my script:
>>
>>
>> import matplotlib
>> matplotlib.use("Agg")
>> from yt.mods import *
>> import pylab
>> import sys
>> import numpy as na
>>
>> def _MyRadius(field, data):
>>     center = data.get_field_parameter("center")
>>     dx = data["x"] - center[0]
>>     dy = data["y"] - center[1]
>>     dz = data["z"] - center[2]
>>     return na.sqrt(dx*dx + dy*dy + dz*dz)
>> add_field("Radius", function=_MyRadius, take_log=False, units='')
>>
>> def _ConvertAccel(data):
>>     return data.convert("Length") / (data.convert("Time"))**2.0
>> def _RadiationAccelerationMagnitude(field, data):
>>     return ( (data["RadAccel1"])**2.0 + \
>>              (data["RadAccel2"])**2.0 + \
>>              (data["RadAccel3"])**2.0 )**(1.0/2.0)
>>     #return ( (data["Density"])*0.0+1.e-28)
>> add_field("RadiationAccelerationMagnitude",
>> function=_RadiationAccelerationMagnitude,
>> take_log=False,convert_function=_ConvertAccel, units=r"\rm{cm}/\rm{s}^{2}")
>> def _RadiationForceMagnitude(field, data):
>>     return ( data["RadiationAccelerationMagnitude"] * data["Density"] *
>> data["CellVolume"])
>> add_field("RadiationForceMagnitude", function=_RadiationForceMagnitude,
>> take_log=False, units=r"\rm{dynes}")
>>
>> def _RadiationForce1(field, data):
>>     return ( (data["RadAccel1"] * data["Density"] * data["CellVolume"]))
>> add_field("RadiationForce1", function=_RadiationForce1, take_log=False,
>> units=r"\rm{dynes}")
>> def _RadiationForce2(field, data):
>>     return (( data["RadAccel2"] * data["Density"] * data["CellVolume"]))
>> add_field("RadiationForce2", function=_RadiationForce2, take_log=False,
>> units=r"\rm{dynes}")
>> def _RadiationForce3(field, data):
>>     return (( data["RadAccel3"] * data["Density"] * data["CellVolume"]))
>> add_field("RadiationForce3", function=_RadiationForce3, take_log=False,
>> units=r"\rm{dynes}")
>>
>>
>>
>> xc = 0.5
>> yc = 0.5
>> zc = 0.5
>> width = 179.2
>> frame_template = "aaSN/%s/RadPresProj_%s_%s"
>>
>> folder = sys.argv[1]
>> starti = int(sys.argv[2])
>> endi = int(sys.argv[3])
>>
>>
>> Q = 1.e47
>> RPSF = 2.0
>> nu = 21.62 * 1.6022e-12 # in ergs
>> print "Q:",Q
>> print "RPScaleFactor:", RPSF
>> print "nu:",nu
>>
>> rho_min = 1.e24
>> rho_max = 5.e27
>>
>> for i in range(starti, endi ,1):
>>     print "a, i=",i
>>     pf = load("SciNet/%s/DD%04i/data%04i" % (folder, i, i))
>>     if pf is None: continue
>>     sp = pf.h.sphere([xc,yc,zc], 1.0)
>>     v, c = pf.h.find_max("Density")
>>     print i, c
>>     print v,c
>>     pc = PlotCollection(pf, center=(xc,yc,zc))
>>     total = 0
>>     print "Initial time, code units =", pf["InitialTime"]
>>     for ax in range(0,1):
>>         p = pc.add_projection("RadiationForceMagnitude", ax,
>> weight_field="Density")
>>         p.set_log_field(True)
>>         p.set_cmap('Blues')
>>         p.set_zlim(1e24, 1e28)
>>         p.set_width(width,"pc")
>>         sp= pf.h.all_data()
>>         v1 = "RadiationForce2"
>>         v2 = "RadiationForce3"
>>         p.modify["quiver"](v1, v2, 12)
>>         if na.any(pf.h.grid_particle_count):
>>             colours = sp["ParticleMassMsun"]
>>             p.modify["particles"](1.0, p_size=10.0, col='r', minimum_mass=
>> 6e0)
>>             print "          XXXXX Particles present XXXXX"
>>         Timemyr = pf["InitialTime"]*pf["years"]/1.e6
>>         p.modify["text"]((0.006,0.899),"%6.1f Myr" % (Timemyr), text_args
>> = {'color':'w','fontsize':'30'})
>>         p.modify["text"]((0.005,0.9),"%6.1f Myr" % (Timemyr), text_args =
>> {'color':'b','fontsize':'30' })
>>         pc.save(frame_template % (sys.argv[1],ax,i), override=True,
>> format="eps")
>>
>>
>>
>>
>
>
> --
> Elizabeth Harper-Clark MA MSci PhD
>
> www.astro.utoronto.ca/~h-clark <http://www.astro.utoronto.ca/%7Eh-clark>
> h-clark at cita.utoronto.ca
> AIM: edphc1
> MSN: edphc1 at hotmail.com
> Skype: eharperclark
> Office phone: 416-978-5759
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20110606/f367bc33/attachment.htm>