[yt-users] Code units have not been defined in Flash dataset
Alex Hill
ashill at haverford.edu
Sat Jan 28 18:38:47 PST 2017
Thanks, Nathan.
I can confirm that the fix in your pull request fixes the issue.
I can also confirm that this workaround works on an unpatched installation of yt.
Thanks!
Cheers,
Alex
---------
Alex Hill
Visiting Assistant Professor
Office: KINSC Link L106
Department of Astronomy, Haverford College
370 Lancaster Ave, Haverford, PA 19041 USA
phone: +1 484 297 2136
email/iMessage: ashill at haverford.edu
> On Jan 28, 2017, at 19:32, Nathan Goldbaum <nathan12343 at gmail.com> wrote:
>
> I should say that as a workaround you can do something like:
>
> left_edge = ds.arr([-500, -500,- 20], 'pc')
> right_edge = ds.arr([500, 500, 20], 'pc')
> b = ds.box(left_edge, right_edge)
>
> I'm already planning to release yt 3.3.4 soon, I'll try to make sure that this fix is part of that release.
>
> Thanks again for the bug report!
>
> -Nathan
>
> On Fri, Jan 27, 2017 at 3:50 PM, Nathan Goldbaum <nathan12343 at gmail.com> wrote:
> Hi Alex,
>
> Thanks for the report! I was able to reproduce the issue and I think I have a fix:
>
> https://bitbucket.org/yt_analysis/yt/pull-requests/2508
>
> If you could give that a test that would be great. No worries if you're not sure how to do that, just wanted an extra bit of confirmation :)
>
> -Nathan
>
>
>
> On Fri, Jan 27, 2017 at 3:32 PM, Alex Hill <ashill at haverford.edu> wrote:
> Hi,
>
> I’ve encountered a new error when specifying a box using units other than code units. This is a Flash 4 plot file in cgs units.
>
> Example:
>
> import yt
> yt.__version__
>
> ‘3.3.3’
>
> d = yt.load('sn34_smd75_sub1_hdf5_plt_cnt_0030')
> b = d.box([-500,-500,-20]*yt.units.pc.in_cgs(),[500,500,20]*yt.units.pc.in_cgs())
> fp = yt.PhasePlot(b,'n','pk',['cell_volume'],weight_field=None, fractional=True)
>
> I get the following error:
>
> UnitParseError: Code units have not been defined.
> Try creating the array or quantity using ds.arr or ds.quan instead.
>
> If I instead do
>
> b = d.box([-1.54e+21, -1.54e+21, -6.17e+19], [1.54e+21, 1.54e+21, 6.17e+19])
>
> or
>
> left = d.domain_left_edge
> right = d.domain_right_edge
> left[2] = -6.17e19
> right[2] = 6.17e19
> b = d.box(left, right)
>
> I do not get any error. If I specify yt.units.cm, I do get the error.
>
> This was not an issue on the same data file a few months ago. I think I have updated yt in the meantime, but I’m not sure what the version was.
>
> Is this a bug, or am I doing something wrong?
>
> Here is the unit registry for my data set:
>
> reg = d.unit_registry
>
> for un in reg.keys():
> if un.startswith('code_'):
> fmt_tup = (un, reg.lut[un][0], str(reg.lut[un][1]))
> print ("Unit name: {:<15}\nCGS conversion: {:<15}\nDimensions: {:<15}\n".format(*fmt_tup))
>
> Unit name: code_magnetic
> CGS conversion: 1.0
> Dimensions: sqrt((mass))/(sqrt((length))*(time))
>
> Unit name: code_time
> CGS conversion: 1.0
> Dimensions: (time)
>
> Unit name: code_temperature
> CGS conversion: 1.0
> Dimensions: (temperature)
>
> Unit name: code_pressure
> CGS conversion: 1.0
> Dimensions: (mass)/((length)*(time)**2)
>
> Unit name: code_mass
> CGS conversion: 1.0
> Dimensions: (mass)
>
> Unit name: code_velocity
> CGS conversion: 1.0
> Dimensions: (length)/(time)
>
> Unit name: code_length
> CGS conversion: 1.0
> Dimensions: (length)
>
> Unit name: code_density
> CGS conversion: 1.0
> Dimensions: (mass)/(length)**3
>
> Unit name: code_metallicity
> CGS conversion: 1.0
> Dimensions: 1
>
> Cheers,
> Alex
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
-------------- next part --------------
A non-text attachment was scrubbed...
Name: smime.p7s
Type: application/pkcs7-signature
Size: 3575 bytes
Desc: not available
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20170128/57240396/attachment.bin>
More information about the yt-users
mailing list