<div dir="ltr">I've opened an issue, let's move discussion over there:<div><br></div><div><a href="https://github.com/yt-project/yt/issues/1552">https://github.com/yt-project/yt/issues/1552</a><br></div><div><br></div><div>Tazkera, you can subscribe to notifications on the issue by clicking the subscribe button.</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Sep 7, 2017 at 2:35 PM, John Zuhone <span dir="ltr"><<a href="mailto:jzuhone@gmail.com" target="_blank">jzuhone@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Hi Tazkera, <div><br></div><div>Thanks for uploading your file—it will make it far simpler to debug. I’ll have a look at this and get back to you. </div><div><br></div><div>Best,</div><div><br></div><div>John</div><div><br><div><blockquote type="cite"><div><div class="h5"><div>On Sep 7, 2017, at 3:24 PM, tazkera haque <<a href="mailto:h.tazkera@gmail.com" target="_blank">h.tazkera@gmail.com</a>> wrote:</div><br class="m_3537214411666842434Apple-interchange-newline"></div></div><div><div><div class="h5"><div dir="ltr"><div><div>HI People, <br><br></div>I am sorry to ask you a very novice question, but I am new to deal with FITS data and I wanted to load it into yt to do some processing. I uploaded my file here.<br><br></div><div>First of all, I am trying to load my first FITS data in yt and it returns me an error which I don't understand:<br><br><b>import yt</b><br><b>ds = yt.load("signal_wfilt_sn.fits"<wbr>)<br>ds.print_stats()<br><br></b></div><div><b>the error:<br><br></b></div><div><div class="m_3537214411666842434gmail-output_area" style="box-sizing:border-box;padding:0px;display:flex;font-family:'Helvetica Neue',Helvetica,Arial,sans-serif;font-size:14px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255)"><div class="m_3537214411666842434gmail-output_stderr m_3537214411666842434gmail-output_subarea m_3537214411666842434gmail-output_text m_3537214411666842434gmail-output_stream" style="box-sizing:border-box;overflow-x:auto;padding:0.4em;text-align:left;line-height:1.21429em;background-image:none;background-color:rgb(255,221,221);background-position:0% 0%;background-repeat:repeat repeat"><pre style="box-sizing:border-box;overflow:auto;font-family:monospace;font-size:inherit;display:block;padding:0px;margin:0px;line-height:inherit;word-break:break-all;background-color:transparent;border:0px none currentcolor;border-top-left-radius:0px;border-top-right-radius:0px;border-bottom-right-radius:0px;border-bottom-left-radius:0px;white-space:pre-wrap;vertical-align:baseline">yt : [WARNING ] 2017-09-07 15:12:42,963 Cannot find time
yt : [INFO ] 2017-09-07 15:12:42,965 Detected these axes: RA---TAN DEC--TAN
yt : [WARNING ] 2017-09-07 15:12:42,979 No length conversion provided. Assuming 1 = 1 cm.
yt : [INFO ] 2017-09-07 15:12:42,996 Parameters: current_time = 0.0
yt : [INFO ] 2017-09-07 15:12:42,997 Parameters: domain_dimensions = [1629 1629 1]
yt : [INFO ] 2017-09-07 15:12:43,001 Parameters: domain_left_edge = [ 0.5 0.5 0.5]
yt : [INFO ] 2017-09-07 15:12:43,004 Parameters: domain_right_edge = [ 1.62950000e+03 1.62950000e+03 1.50000000e+00]
yt : [INFO ] 2017-09-07 15:12:43,006 Parameters: cosmological_simulation = 0.0
</pre></div></div><div class="m_3537214411666842434gmail-output_area" style="box-sizing:border-box;padding:0px;display:flex;font-family:'Helvetica Neue',Helvetica,Arial,sans-serif;font-size:14px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:normal;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255)"><div class="m_3537214411666842434gmail-output_subarea m_3537214411666842434gmail-output_text m_3537214411666842434gmail-output_error" style="box-sizing:border-box;overflow-x:auto;padding:0.4em;text-align:left;line-height:1.21429em"><pre style="box-sizing:border-box;overflow:auto;font-family:monospace;font-size:inherit;display:block;padding:0px;margin:0px;line-height:inherit;word-break:break-all;background-color:transparent;border:0px none currentcolor;border-top-left-radius:0px;border-top-right-radius:0px;border-bottom-right-radius:0px;border-bottom-left-radius:0px;white-space:pre-wrap;vertical-align:baseline"><span class="m_3537214411666842434gmail-ansired" style="box-sizing:border-box;color:darkred">------------------------------<wbr>------------------------------<wbr>---------------</span>
<span class="m_3537214411666842434gmail-ansired" style="box-sizing:border-box;color:darkred">IndexError</span> Traceback (most recent call last)
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><ipython-input-38-<wbr>71d5999c6b09></span> in <span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><module><span class="m_3537214411666842434gmail-ansiblue" style="box-sizing:border-box;color:darkblue">()</span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 1</span> <span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">import</span> yt<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 2</span> ds2 <span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">=</span> yt<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">.</span>load<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">(</span><span class="m_3537214411666842434gmail-ansiblue" style="box-sizing:border-box;color:darkblue">"/home/trina/<wbr>Documents/ASTRO_732/A-01/<wbr>signal_wfilt_sn.fits"</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">)</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">----> 3<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"> </span>ds2<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">.</span>print_stats<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">(</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">)</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">/home/trina/anaconda2/lib/<wbr>python2.7/site-packages/yt/<wbr>utilities/parallel_tools/<wbr>parallel_analysis_interface.<wbr>pyc</span> in <span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue">root_only<span class="m_3537214411666842434gmail-ansiblue" style="box-sizing:border-box;color:darkblue">(*args, **kwargs)</span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 318</span> <span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">def</span> root_only<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">(</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">*</span>args<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">,</span> <span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">**</span>kwargs<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">)</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">:</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 319</span> <span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">if</span> <span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">not</span> parallel_capable<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">:</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">--> 320<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"> <span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">return</span> func<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">(</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">*</span>args<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">,</span> <span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">**</span>kwargs<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">)</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 321</span> comm <span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">=</span> _get_comm<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">(</span>args<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">)</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 322</span> rv <span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">=</span> None<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">/home/trina/anaconda2/lib/<wbr>python2.7/site-packages/yt/<wbr>data_objects/static_output.pyc</span> in <span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue">print_stats<span class="m_3537214411666842434gmail-ansiblue" style="box-sizing:border-box;color:darkblue">(self)</span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 455</span> <span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">@</span>parallel_root_only<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 456</span> <span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">def</span> print_stats<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">(</span>self<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">)</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">:</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">--> 457<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"> </span>self<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">.</span>index<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">.</span>print_stats<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">(</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">)</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 458</span> <span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 459</span> <span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">@</span>property<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">/home/trina/anaconda2/lib/<wbr>python2.7/site-packages/yt/<wbr>data_objects/static_output.pyc</span> in <span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue">index<span class="m_3537214411666842434gmail-ansiblue" style="box-sizing:border-box;color:darkblue">(self)</span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 417</span> <span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">raise</span> RuntimeError<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">(</span><span class="m_3537214411666842434gmail-ansiblue" style="box-sizing:border-box;color:darkblue">"You should not instantiate Dataset."</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">)</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 418</span> self._instantiated_index = self._index_class(
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">--> 419<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"> self, dataset_type=self.dataset_<wbr>type)
</span><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 420</span> <span class="m_3537214411666842434gmail-ansired" style="box-sizing:border-box;color:darkred"># Now we do things that we need an instantiated index for</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 421</span> <span class="m_3537214411666842434gmail-ansired" style="box-sizing:border-box;color:darkred"># ...first off, we create our field_info now.</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">/home/trina/anaconda2/lib/<wbr>python2.7/site-packages/yt/<wbr>frontends/fits/data_<wbr>structures.pyc</span> in <span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue">__init__<span class="m_3537214411666842434gmail-ansiblue" style="box-sizing:border-box;color:darkblue">(self, ds, dataset_type)</span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 90</span> self<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">.</span>_handle <span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">=</span> ds<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">.</span>_handle<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 91</span> self<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">.</span>float_type <span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">=</span> np<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">.</span>float64<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">---> 92<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"> </span>GridIndex<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">.</span>__init__<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">(</span>self<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">,</span>ds<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">,</span>dat<wbr>aset_type<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">)</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 93</span> <span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 94</span> <span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">def</span> _initialize_data_storage<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">(</span>self<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">)</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"><wbr>:</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">/home/trina/anaconda2/lib/<wbr>python2.7/site-packages/yt/<wbr>geometry/geometry_handler.pyc</span> in <span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue">__init__<span class="m_3537214411666842434gmail-ansiblue" style="box-sizing:border-box;color:darkblue">(self, ds, dataset_type)</span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 56</span> <span class="m_3537214411666842434gmail-ansired" style="box-sizing:border-box;color:darkred"># potentially quite expensive, and should be done with the indexing.</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 57</span> mylog<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">.</span>debug<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">(</span><span class="m_3537214411666842434gmail-ansiblue" style="box-sizing:border-box;color:darkblue">"Detecting fields."</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">)</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">---> 58<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"> </span>self<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">.</span>_detect_output_fields<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">(</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">)</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 59</span> <span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 60</span> <span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">def</span> _initialize_state_variables<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">(</span>se<wbr>lf<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">)</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">:</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">/home/trina/anaconda2/lib/<wbr>python2.7/site-packages/yt/<wbr>frontends/fits/data_<wbr>structures.pyc</span> in <span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue">_detect_output_fields<span class="m_3537214411666842434gmail-ansiblue" style="box-sizing:border-box;color:darkblue">(self)</span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 171</span> <span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">continue</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 172</span> <span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">if</span> self<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">.</span>_ensure_same_dims<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">(</span>hdu<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">)</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">:</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">--> 173<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"> </span>units <span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">=</span> self<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">.</span>_determine_image_units<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">(</span>hd<wbr>u<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">.</span>header<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">,</span> known_units<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">)</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 174</span> <span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">try</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">:</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 175</span> <span class="m_3537214411666842434gmail-ansired" style="box-sizing:border-box;color:darkred"># Grab field name from btype</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">/home/trina/anaconda2/lib/<wbr>python2.7/site-packages/yt/<wbr>frontends/fits/data_<wbr>structures.pyc</span> in <span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue">_determine_image_units<span class="m_3537214411666842434gmail-ansiblue" style="box-sizing:border-box;color:darkblue">(self, header, known_units)</span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 116</span> n <span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">+=</span> <span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue">1</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 117</span> <span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">if</span> n <span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">!=</span> len<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">(</span>units<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">)</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">:</span> field_units <span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">=</span> <span class="m_3537214411666842434gmail-ansiblue" style="box-sizing:border-box;color:darkblue">"dimensionless"</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">--> 118<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"> <span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">if</span> field_units<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">[</span><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue">0</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">]</span> <span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">==</span> <span class="m_3537214411666842434gmail-ansiblue" style="box-sizing:border-box;color:darkblue">"/"</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">:</span><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 119</span> field_units <span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">=</span> <span class="m_3537214411666842434gmail-ansiblue" style="box-sizing:border-box;color:darkblue">"1%s"</span> <span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)">%</span> field_units<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"> 120</span> <span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen">return</span> field_units<span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"></span>
<span class="m_3537214411666842434gmail-ansired" style="box-sizing:border-box;color:darkred">IndexError</span>: string index out of range<br><br></span></span></span></span></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:inherit;display:block;padding:0px;margin:0px;line-height:inherit;word-break:break-all;background-color:transparent;border:0px none currentcolor;border-top-left-radius:0px;border-top-right-radius:0px;border-bottom-right-radius:0px;border-bottom-left-radius:0px;white-space:pre-wrap;vertical-align:baseline"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"><span><span style="background-color:rgb(0,0,0)"><span></span></span>Can you Please suggest me how to load and work on this data?</span><br><br></span></span></span></span></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:inherit;display:block;padding:0px;margin:0px;line-height:inherit;word-break:break-all;background-color:transparent;border:0px none currentcolor;border-top-left-radius:0px;border-top-right-radius:0px;border-bottom-right-radius:0px;border-bottom-left-radius:0px;white-space:pre-wrap;vertical-align:baseline"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"><span>My goal is to create an image from the signal to noise data in the data file 'signal_wfilt_sn.fits' using only data<br>within 20 arcminutes of the center of the image. the pixel size is 3x3 arc second</span><br><br></span></span></span></span></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:inherit;display:block;padding:0px;margin:0px;line-height:inherit;word-break:break-all;background-color:transparent;border:0px none currentcolor;border-top-left-radius:0px;border-top-right-radius:0px;border-bottom-right-radius:0px;border-bottom-left-radius:0px;white-space:pre-wrap;vertical-align:baseline"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansicyan" style="box-sizing:border-box;color:steelblue"><span class="m_3537214411666842434gmail-ansigreen" style="box-sizing:border-box;color:darkgreen"><span class="m_3537214411666842434gmail-ansiyellow" style="box-sizing:border-box;color:rgb(196,160,0)"><span>I had an initial code in mind which will create a sphere with a radius of 20 arcminute:</span><br><br></span></span></span></span></span></span></span></span></span></span></span></span></span></span></span></span></span></span><br><b>import yt<br>ds = yt.load("signal_wfilt_sn.fits"<wbr>)<br>ds.print_stats()<br>sp = ds.sphere(ds.domain_center, (400, 'code_length'))<br></b><br></pre><pre style="box-sizing:border-box;overflow:auto;font-family:monospace;font-size:inherit;display:block;padding:0px;margin:0px;line-height:inherit;word-break:break-all;background-color:transparent;border:0px none currentcolor;border-top-left-radius:0px;border-top-right-radius:0px;border-bottom-right-radius:0px;border-bottom-left-radius:0px;white-space:pre-wrap;vertical-align:baseline">However, I first need to be able to load and explore the data. Your help and suggestion is greatly appreciated.<br><br></pre><pre style="box-sizing:border-box;overflow:auto;font-family:monospace;font-size:inherit;display:block;padding:0px;margin:0px;line-height:inherit;word-break:break-all;background-color:transparent;border:0px none currentcolor;border-top-left-radius:0px;border-top-right-radius:0px;border-bottom-right-radius:0px;border-bottom-left-radius:0px;white-space:pre-wrap;vertical-align:baseline">Thanks a lot<b><br></b></pre><pre style="box-sizing:border-box;overflow:auto;font-family:monospace;font-size:inherit;display:block;padding:0px;margin:0px;line-height:inherit;word-break:break-all;background-color:transparent;border:0px none currentcolor;border-top-left-radius:0px;border-top-right-radius:0px;border-bottom-right-radius:0px;border-bottom-left-radius:0px;white-space:pre-wrap;vertical-align:baseline"><br></pre><pre style="box-sizing:border-box;overflow:auto;font-family:monospace;font-size:inherit;display:block;padding:0px;margin:0px;line-height:inherit;word-break:break-all;background-color:transparent;border:0px none currentcolor;border-top-left-radius:0px;border-top-right-radius:0px;border-bottom-right-radius:0px;border-bottom-left-radius:0px;white-space:pre-wrap;vertical-align:baseline">Tazkera<br></pre><pre style="box-sizing:border-box;overflow:auto;font-family:monospace;font-size:inherit;display:block;padding:0px;margin:0px;line-height:inherit;word-break:break-all;background-color:transparent;border:0px none currentcolor;border-top-left-radius:0px;border-top-right-radius:0px;border-bottom-right-radius:0px;border-bottom-left-radius:0px;white-space:pre-wrap;vertical-align:baseline"><br></pre></div></div><br class="m_3537214411666842434gmail-Apple-interchange-newline"></div><br><br><div><b><br></b></div><div><br></div><div><br><br><br><br><br></div></div>
</div></div><span id="m_3537214411666842434cid:A32F5937-5A32-4303-881C-AEF5628D7EA6"><signal_wfilt_sn.fits></span>________<wbr>______________________________<wbr>_________<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" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br></div></blockquote></div><br></div></div><br>______________________________<wbr>_________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">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/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
<br></blockquote></div><br></div>