<div><div dir="auto">Hi Tazkera,</div><div dir="auto"><br></div><div dir="auto">Maybe you want to use units_override, which will let you specify units for your data if the ones yt infers from the output file are wrong for some reason:</div><div dir="auto"><br></div><div dir="auto"><a href="http://yt-project.org/doc/analyzing/units/comoving_units_and_code_units.html#Overriding-Code-Unit-Definitions">http://yt-project.org/doc/analyzing/units/comoving_units_and_code_units.html#Overriding-Code-Unit-Definitions</a><br></div><div dir="auto"><br></div><div dir="auto">That will let you avoid all of these ad-hoc unit conversions and work inside the unit system.</div><div dir="auto"><br></div><div dir="auto">Nathan</div><br><div class="gmail_quote"><div>On Thu, Oct 5, 2017 at 4:59 PM tazkera haque <<a href="mailto:h.tazkera@gmail.com">h.tazkera@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div><div><div>Dear yt people, <br><br></div>I am working on a
tipsy output file from a n-body simulation. My file contains some scaled
units for fields like density, coordinates, mass etc. To make a
physical sense of the data, I had to convert those fields units to some
physical units (c.g.s) and I derived new fields having cgs unit. for
example:<br><br><b>dd = ds.all_data()<br><br># derive a new domain width with conversion factor (6Mpc*current_time/hubble_constant) <br>def _NewDomainWidth(field, data):<br> return ds.domain_width*((3.086e18*6e6*0.909090893796)/(0.7))<br><br>ds.add_field(('gas', 'PhysicalDomainWidth'), function=_NewDomainWidth, units="code_length")</b><br></div><div><br></div><div><b># derive a new domain center with conversion factor (6Mpc*current_time/hubble_constant) <br>def _NewDomainCenter(field, data):<br> return ds.domain_center*((3.086e18*6e6*0.909090893796)/(0.7))<br><br>ds.add_field(('gas', 'PhysicalDomainCenter'), function=_NewDomainCenter, units="code_length")</b><br></div><br></div><div>note that here "code_length" is actually in "cm". but I could not just do dd["Gas",'Coordinates'].in_units('cm'),
because this conversion does not take into account some constant
corrections and returns the wrong length magnitude.<br><br></div>My goal was to make sliceplot and projectionplot of a newly derived field :<br><br><b># derive a new electron pressure field<br>def _ElectronPressure(field, data):<br> return data['Gas','Density']*(5.49e-7)*data['Gas','Temperature']<br><br>ds.add_field(('gas', 'ElectronPressure'), function=_ElectronPressure, units="(code_mass*K)/(code_length**3)")<br><br><br></b></div><b># integration of pressure field along line of sight<br></b><div><b>yt.ProjectionPlot(ds, 'z', ('gas', 'ElectronPressure'), width = dd['gas', 'PhysicalDomainWidth'], <br> center = dd['gas', 'PhysicalDomainCenter'] )</b></div><div><br></div><div><b>and it is returning me <br></b></div><pre style="box-sizing:border-box;overflow:auto;font-family:monospace;font-size:14px;display:block;padding:0px;margin:0px;line-height:inherit;word-break:break-all;color:rgb(0,0,0);background-color:rgb(255,255,255);border-width:0px;border-style:none;border-color:currentcolor;border-radius:0px;white-space:pre-wrap;vertical-align:baseline;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:left;text-indent:0px;text-transform:none;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansired" style="box-sizing:border-box;color:darkred">ValueError</span> Traceback (most recent call last)
<span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><ipython-input-8-fa351b6c09d1></span> in <span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><module><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansiblue" style="box-sizing:border-box;color:darkblue">()</span>
<span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 1</span> yt.ProjectionPlot(ds, 'z', ('gas', 'ElectronPressure'), width = dd['gas', 'PhysicalDomainWidth'],
<span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">----> 2<span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"> center = dd['gas', 'PhysicalDomainCenter'] )</span></span></span></pre><pre style="box-sizing:border-box;overflow:auto;font-family:monospace;font-size:14px;display:block;padding:0px;margin:0px;line-height:inherit;word-break:break-all;color:rgb(0,0,0);background-color:rgb(255,255,255);border-width:0px;border-style:none;border-color:currentcolor;border-radius:0px;white-space:pre-wrap;vertical-align:baseline;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:left;text-indent:0px;text-transform:none;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"><span style="color:rgb(0,0,0)"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansired" style="box-sizing:border-box">ValueError</span>: operands could not be broadcast together with shapes (39678,) (54608,) </span><br><br></span></span></span></span></span></span></span></span></span></span></span></span></span></pre><pre style="box-sizing:border-box;overflow:auto;font-family:monospace;font-size:14px;display:block;padding:0px;margin:0px;line-height:inherit;word-break:break-all;background-color:rgb(255,255,255);border-width:0px;border-style:none;border-color:currentcolor;border-radius:0px;white-space:pre-wrap;vertical-align:baseline;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:left;text-indent:0px;text-transform:none;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><span style="color:rgb(0,0,0)"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansiyellow" style="box-sizing:border-box">Now I cannot use width =ds.domain_width and center = ds.domain_center for the same unit error. I was wondering how I can make the projectionplot for my newly derived field and newly derived domain width and domain center.<br></span></span></span></span></span></span></span></span></span></span></span></span></span></span></pre><pre style="box-sizing:border-box;overflow:auto;font-family:monospace;font-size:14px;display:block;padding:0px;margin:0px;line-height:inherit;word-break:break-all;color:rgb(0,0,0);background-color:rgb(255,255,255);border-width:0px;border-style:none;border-color:currentcolor;border-radius:0px;white-space:pre-wrap;vertical-align:baseline;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:left;text-indent:0px;text-transform:none;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"><br><br></span></span></span></span></span></span></span></span></span></span></span></span></span></pre><pre style="box-sizing:border-box;overflow:auto;font-family:monospace;font-size:14px;display:block;padding:0px;margin:0px;line-height:inherit;word-break:break-all;color:rgb(0,0,0);background-color:rgb(255,255,255);border-width:0px;border-style:none;border-color:currentcolor;border-radius:0px;white-space:pre-wrap;vertical-align:baseline;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:left;text-indent:0px;text-transform:none;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"><span style="color:rgb(0,0,0)">Your suggestion is most welcome</span><br></span></span></span></span></span></span></span></span></span></span></span></span></span></pre><pre style="box-sizing:border-box;overflow:auto;font-family:monospace;font-size:14px;display:block;padding:0px;margin:0px;line-height:inherit;word-break:break-all;color:rgb(0,0,0);background-color:rgb(255,255,255);border-width:0px;border-style:none;border-color:currentcolor;border-radius:0px;white-space:pre-wrap;vertical-align:baseline;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:left;text-indent:0px;text-transform:none;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"><br></span></span></span></span></span></span></span></span></span></span></span></span></span></pre><pre style="box-sizing:border-box;overflow:auto;font-family:monospace;font-size:14px;display:block;padding:0px;margin:0px;line-height:inherit;word-break:break-all;background-color:rgb(255,255,255);border-width:0px;border-style:none;border-color:currentcolor;border-radius:0px;white-space:pre-wrap;vertical-align:baseline;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:left;text-indent:0px;text-transform:none;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><span style="color:rgb(0,0,0)"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansiyellow" style="box-sizing:border-box">Best<br></span></span></span></span></span></span></span></span></span></span></span></span></span></span></pre><pre style="box-sizing:border-box;overflow:auto;font-family:monospace;font-size:14px;display:block;padding:0px;margin:0px;line-height:inherit;word-break:break-all;background-color:rgb(255,255,255);border-width:0px;border-style:none;border-color:currentcolor;border-radius:0px;white-space:pre-wrap;vertical-align:baseline;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:left;text-indent:0px;text-transform:none;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><span style="color:rgb(0,0,0)"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansiyellow" style="box-sizing:border-box">Tazkera<br></span></span></span></span></span></span></span></span></span></span></span></span></span></span></pre><pre style="box-sizing:border-box;overflow:auto;font-family:monospace;font-size:14px;display:block;padding:0px;margin:0px;line-height:inherit;word-break:break-all;background-color:rgb(255,255,255);border-width:0px;border-style:none;border-color:currentcolor;border-radius:0px;white-space:pre-wrap;vertical-align:baseline;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:left;text-indent:0px;text-transform:none;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><span style="color:rgb(0,0,0)"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansicyan" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansigreen" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-m_6044548144393904699gmail-ansiyellow" style="box-sizing:border-box"><span class="m_4703861014847096559gmail-HOEnZb"><font color="#888888"></font></span></span></span></span></span></span></span></span></span></span></span></span></span></span></span></pre></div>
_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
</blockquote></div></div>