In [1]:
import os
import io
import copy
import glob
import urllib
import numpy as np
import h5py
import matplotlib.pyplot as plt
import matplotlib.colors
import seaborn as sns
import astropy.units as u
import astropy.constants as const
from scipy.ndimage import gaussian_filter
from sunpy.map import Map,GenericMap
import synthesizAR
from synthesizAR.util import EMCube
from synthesizAR.instruments import InstrumentHinodeEIS
%matplotlib inline
/home/wtb2/anaconda3/envs/systematic-ar-study/lib/python3.5/site-packages/IPython/html.py:14: ShimWarning: The `IPython.html` package has been deprecated. You should import from `notebook` instead. `IPython.html.widgets` has moved to `ipywidgets`.
"`IPython.html.widgets` has moved to `ipywidgets`.", ShimWarning)
WARNING: AstropyDeprecationWarning: astropy.utils.compat.odict.OrderedDict is now deprecated - import OrderedDict from the collections module instead [astropy.utils.compat.odict]
/home/wtb2/anaconda3/envs/systematic-ar-study/lib/python3.5/site-packages/matplotlib/__init__.py:1401: UserWarning: This call to matplotlib.use() has no effect
because the backend has already been chosen;
matplotlib.use() must be called *before* pylab, matplotlib.pyplot,
or matplotlib.backends is imported for the first time.
warnings.warn(_use_error_msg)
In [2]:
base1 = '/data/datadrive1/ar_forward_modeling/systematic_ar_study/noaa1109_tn{}'
base2 = '/data/datadrive2/ar_viz/systematic_ar_study/noaa1109_tn{}/'
In [3]:
eis = InstrumentHinodeEIS([7.5e3,1.25e4]*u.s)
In [4]:
frequencies = [250,750,'750-ion',2500,5000]
In [5]:
temperature_bin_edges = 10.**(np.arange(5.6, 7.0, 0.05))*u.K
Iterate over all "true" emission measure distributions and time-average them over the given interval.
In [6]:
time_averaged_ems = {'{}'.format(freq):None for freq in frequencies}
In [7]:
for freq in frequencies:
print('tn = {} s'.format(freq))
if type(freq) == int:
base = base1
else:
base = base2
# setup field and observer objects
field = synthesizAR.Skeleton.restore(os.path.join(base.format(freq),'field_checkpoint'))
observer = synthesizAR.Observer(field,[eis],ds=field._convert_angle_to_length(0.4*u.arcsec))
observer.build_detector_files(base.format(freq))
# iterate over time
for time in eis.observing_time:
print('t = {}'.format(time))
emcube = observer.make_emission_measure_map(time,eis,temperature_bin_edges=temperature_bin_edges)
if time_averaged_ems['{}'.format(freq)] is None:
time_averaged_ems['{}'.format(freq)] = emcube
for m in time_averaged_ems['{}'.format(freq)]:
m.data /= eis.observing_time.shape[0]
else:
for m1,m2 in zip(time_averaged_ems['{}'.format(freq)],emcube):
m1.data += m2.data/eis.observing_time.shape[0]
tn = 250 s
No HMI fits file supplied. A new HMI map object will not be created.
yt : [INFO ] 2017-06-20 14:05:38,707 Parameters: current_time = 0.0
yt : [INFO ] 2017-06-20 14:05:38,709 Parameters: domain_dimensions = [96 96 96]
yt : [INFO ] 2017-06-20 14:05:38,710 Parameters: domain_left_edge = [ 8.88147851e+09 1.66017386e+09 7.23009049e+08]
yt : [INFO ] 2017-06-20 14:05:38,711 Parameters: domain_right_edge = [ 4.35656994e+10 3.67056887e+10 3.54274434e+10]
yt : [INFO ] 2017-06-20 14:05:38,713 Parameters: cosmological_simulation = 0.0
/home/wtb2/anaconda3/envs/systematic-ar-study/lib/python3.5/site-packages/scipy/interpolate/fitpack.py:300: RuntimeWarning: The maximal number of iterations (20) allowed for finding smoothing
spline with fp=s has been reached. Probable cause: s too small.
(abs(fp-s)/s>0.001)
warnings.warn(RuntimeWarning(_iermess[ier][0]))
/home/wtb2/anaconda3/envs/systematic-ar-study/lib/python3.5/site-packages/scipy/interpolate/fitpack.py:300: RuntimeWarning: A theoretically impossible result when finding a smoothing spline
with fp = s. Probable cause: s too small. (abs(fp-s)/s>0.001)
warnings.warn(RuntimeWarning(_iermess[ier][0]))
t = 7500.0 s
t = 7510.0 s
t = 7520.0 s
t = 7530.0 s
t = 7540.0 s
t = 7550.0 s
t = 7560.0 s
t = 7570.0 s
t = 7580.0 s
t = 7590.0 s
t = 7600.0 s
t = 7610.0 s
t = 7620.0 s
t = 7630.0 s
t = 7640.0 s
t = 7650.0 s
t = 7660.0 s
t = 7670.0 s
t = 7680.0 s
t = 7690.0 s
t = 7700.0 s
t = 7710.0 s
t = 7720.0 s
t = 7730.0 s
t = 7740.0 s
t = 7750.0 s
t = 7760.0 s
t = 7770.0 s
t = 7780.0 s
t = 7790.0 s
t = 7800.0 s
t = 7810.0 s
t = 7820.0 s
t = 7830.0 s
t = 7840.0 s
t = 7850.0 s
t = 7860.0 s
t = 7870.0 s
t = 7880.0 s
t = 7890.0 s
t = 7900.0 s
t = 7910.0 s
t = 7920.0 s
t = 7930.0 s
t = 7940.0 s
t = 7950.0 s
t = 7960.0 s
t = 7970.0 s
t = 7980.0 s
t = 7990.0 s
t = 8000.0 s
t = 8010.0 s
t = 8020.0 s
t = 8030.0 s
t = 8040.0 s
t = 8050.0 s
t = 8060.0 s
t = 8070.0 s
t = 8080.0 s
t = 8090.0 s
t = 8100.0 s
t = 8110.0 s
t = 8120.0 s
t = 8130.0 s
t = 8140.0 s
t = 8150.0 s
t = 8160.0 s
t = 8170.0 s
t = 8180.0 s
t = 8190.0 s
t = 8200.0 s
t = 8210.0 s
t = 8220.0 s
t = 8230.0 s
t = 8240.0 s
t = 8250.0 s
t = 8260.0 s
t = 8270.0 s
t = 8280.0 s
t = 8290.0 s
t = 8300.0 s
t = 8310.0 s
t = 8320.0 s
t = 8330.0 s
t = 8340.0 s
t = 8350.0 s
t = 8360.0 s
t = 8370.0 s
t = 8380.0 s
t = 8390.0 s
t = 8400.0 s
t = 8410.0 s
t = 8420.0 s
t = 8430.0 s
t = 8440.0 s
t = 8450.0 s
t = 8460.0 s
t = 8470.0 s
t = 8480.0 s
t = 8490.0 s
t = 8500.0 s
t = 8510.0 s
t = 8520.0 s
t = 8530.0 s
t = 8540.0 s
t = 8550.0 s
t = 8560.0 s
t = 8570.0 s
t = 8580.0 s
t = 8590.0 s
t = 8600.0 s
t = 8610.0 s
t = 8620.0 s
t = 8630.0 s
t = 8640.0 s
t = 8650.0 s
t = 8660.0 s
t = 8670.0 s
t = 8680.0 s
t = 8690.0 s
t = 8700.0 s
t = 8710.0 s
t = 8720.0 s
t = 8730.0 s
t = 8740.0 s
t = 8750.0 s
t = 8760.0 s
t = 8770.0 s
t = 8780.0 s
t = 8790.0 s
t = 8800.0 s
t = 8810.0 s
t = 8820.0 s
t = 8830.0 s
t = 8840.0 s
t = 8850.0 s
t = 8860.0 s
t = 8870.0 s
t = 8880.0 s
t = 8890.0 s
t = 8900.0 s
t = 8910.0 s
t = 8920.0 s
t = 8930.0 s
t = 8940.0 s
t = 8950.0 s
t = 8960.0 s
t = 8970.0 s
t = 8980.0 s
t = 8990.0 s
t = 9000.0 s
t = 9010.0 s
t = 9020.0 s
t = 9030.0 s
t = 9040.0 s
t = 9050.0 s
t = 9060.0 s
t = 9070.0 s
t = 9080.0 s
t = 9090.0 s
t = 9100.0 s
t = 9110.0 s
t = 9120.0 s
t = 9130.0 s
t = 9140.0 s
t = 9150.0 s
t = 9160.0 s
t = 9170.0 s
t = 9180.0 s
t = 9190.0 s
t = 9200.0 s
t = 9210.0 s
t = 9220.0 s
t = 9230.0 s
t = 9240.0 s
t = 9250.0 s
t = 9260.0 s
t = 9270.0 s
t = 9280.0 s
t = 9290.0 s
t = 9300.0 s
t = 9310.0 s
t = 9320.0 s
t = 9330.0 s
t = 9340.0 s
t = 9350.0 s
t = 9360.0 s
t = 9370.0 s
t = 9380.0 s
t = 9390.0 s
t = 9400.0 s
t = 9410.0 s
t = 9420.0 s
t = 9430.0 s
t = 9440.0 s
t = 9450.0 s
t = 9460.0 s
t = 9470.0 s
t = 9480.0 s
t = 9490.0 s
t = 9500.0 s
t = 9510.0 s
t = 9520.0 s
t = 9530.0 s
t = 9540.0 s
t = 9550.0 s
t = 9560.0 s
t = 9570.0 s
t = 9580.0 s
t = 9590.0 s
t = 9600.0 s
t = 9610.0 s
t = 9620.0 s
t = 9630.0 s
t = 9640.0 s
t = 9650.0 s
t = 9660.0 s
t = 9670.0 s
t = 9680.0 s
t = 9690.0 s
t = 9700.0 s
t = 9710.0 s
t = 9720.0 s
t = 9730.0 s
t = 9740.0 s
t = 9750.0 s
t = 9760.0 s
t = 9770.0 s
t = 9780.0 s
t = 9790.0 s
t = 9800.0 s
t = 9810.0 s
t = 9820.0 s
t = 9830.0 s
t = 9840.0 s
t = 9850.0 s
t = 9860.0 s
t = 9870.0 s
t = 9880.0 s
t = 9890.0 s
t = 9900.0 s
t = 9910.0 s
t = 9920.0 s
t = 9930.0 s
t = 9940.0 s
t = 9950.0 s
t = 9960.0 s
t = 9970.0 s
t = 9980.0 s
t = 9990.0 s
t = 10000.0 s
t = 10010.0 s
t = 10020.0 s
t = 10030.0 s
t = 10040.0 s
t = 10050.0 s
t = 10060.0 s
t = 10070.0 s
t = 10080.0 s
t = 10090.0 s
t = 10100.0 s
t = 10110.0 s
t = 10120.0 s
t = 10130.0 s
t = 10140.0 s
t = 10150.0 s
t = 10160.0 s
t = 10170.0 s
t = 10180.0 s
t = 10190.0 s
t = 10200.0 s
t = 10210.0 s
t = 10220.0 s
t = 10230.0 s
t = 10240.0 s
t = 10250.0 s
t = 10260.0 s
t = 10270.0 s
t = 10280.0 s
t = 10290.0 s
t = 10300.0 s
t = 10310.0 s
t = 10320.0 s
t = 10330.0 s
t = 10340.0 s
t = 10350.0 s
t = 10360.0 s
t = 10370.0 s
t = 10380.0 s
t = 10390.0 s
t = 10400.0 s
t = 10410.0 s
t = 10420.0 s
t = 10430.0 s
t = 10440.0 s
t = 10450.0 s
t = 10460.0 s
t = 10470.0 s
t = 10480.0 s
t = 10490.0 s
t = 10500.0 s
t = 10510.0 s
t = 10520.0 s
t = 10530.0 s
t = 10540.0 s
t = 10550.0 s
t = 10560.0 s
t = 10570.0 s
t = 10580.0 s
t = 10590.0 s
t = 10600.0 s
t = 10610.0 s
t = 10620.0 s
t = 10630.0 s
t = 10640.0 s
t = 10650.0 s
t = 10660.0 s
t = 10670.0 s
t = 10680.0 s
t = 10690.0 s
t = 10700.0 s
t = 10710.0 s
t = 10720.0 s
t = 10730.0 s
t = 10740.0 s
t = 10750.0 s
t = 10760.0 s
t = 10770.0 s
t = 10780.0 s
t = 10790.0 s
t = 10800.0 s
t = 10810.0 s
t = 10820.0 s
t = 10830.0 s
t = 10840.0 s
t = 10850.0 s
t = 10860.0 s
t = 10870.0 s
t = 10880.0 s
t = 10890.0 s
t = 10900.0 s
t = 10910.0 s
t = 10920.0 s
t = 10930.0 s
t = 10940.0 s
t = 10950.0 s
t = 10960.0 s
t = 10970.0 s
t = 10980.0 s
t = 10990.0 s
t = 11000.0 s
t = 11010.0 s
t = 11020.0 s
t = 11030.0 s
t = 11040.0 s
t = 11050.0 s
t = 11060.0 s
t = 11070.0 s
t = 11080.0 s
t = 11090.0 s
t = 11100.0 s
t = 11110.0 s
t = 11120.0 s
t = 11130.0 s
t = 11140.0 s
t = 11150.0 s
t = 11160.0 s
t = 11170.0 s
t = 11180.0 s
t = 11190.0 s
t = 11200.0 s
t = 11210.0 s
t = 11220.0 s
t = 11230.0 s
t = 11240.0 s
t = 11250.0 s
t = 11260.0 s
t = 11270.0 s
t = 11280.0 s
t = 11290.0 s
t = 11300.0 s
t = 11310.0 s
t = 11320.0 s
t = 11330.0 s
t = 11340.0 s
t = 11350.0 s
t = 11360.0 s
t = 11370.0 s
t = 11380.0 s
t = 11390.0 s
t = 11400.0 s
t = 11410.0 s
t = 11420.0 s
t = 11430.0 s
t = 11440.0 s
t = 11450.0 s
t = 11460.0 s
t = 11470.0 s
t = 11480.0 s
t = 11490.0 s
t = 11500.0 s
t = 11510.0 s
t = 11520.0 s
t = 11530.0 s
t = 11540.0 s
t = 11550.0 s
t = 11560.0 s
t = 11570.0 s
t = 11580.0 s
t = 11590.0 s
t = 11600.0 s
t = 11610.0 s
t = 11620.0 s
t = 11630.0 s
t = 11640.0 s
t = 11650.0 s
t = 11660.0 s
t = 11670.0 s
t = 11680.0 s
t = 11690.0 s
t = 11700.0 s
t = 11710.0 s
t = 11720.0 s
t = 11730.0 s
t = 11740.0 s
t = 11750.0 s
t = 11760.0 s
t = 11770.0 s
t = 11780.0 s
t = 11790.0 s
t = 11800.0 s
t = 11810.0 s
t = 11820.0 s
t = 11830.0 s
t = 11840.0 s
t = 11850.0 s
t = 11860.0 s
t = 11870.0 s
t = 11880.0 s
t = 11890.0 s
t = 11900.0 s
t = 11910.0 s
t = 11920.0 s
t = 11930.0 s
t = 11940.0 s
t = 11950.0 s
t = 11960.0 s
t = 11970.0 s
t = 11980.0 s
t = 11990.0 s
t = 12000.0 s
t = 12010.0 s
t = 12020.0 s
t = 12030.0 s
t = 12040.0 s
t = 12050.0 s
t = 12060.0 s
t = 12070.0 s
t = 12080.0 s
t = 12090.0 s
t = 12100.0 s
t = 12110.0 s
t = 12120.0 s
t = 12130.0 s
t = 12140.0 s
t = 12150.0 s
t = 12160.0 s
t = 12170.0 s
t = 12180.0 s
t = 12190.0 s
t = 12200.0 s
t = 12210.0 s
t = 12220.0 s
t = 12230.0 s
t = 12240.0 s
t = 12250.0 s
t = 12260.0 s
t = 12270.0 s
t = 12280.0 s
t = 12290.0 s
t = 12300.0 s
t = 12310.0 s
t = 12320.0 s
t = 12330.0 s
t = 12340.0 s
t = 12350.0 s
t = 12360.0 s
t = 12370.0 s
t = 12380.0 s
t = 12390.0 s
t = 12400.0 s
t = 12410.0 s
t = 12420.0 s
t = 12430.0 s
t = 12440.0 s
t = 12450.0 s
t = 12460.0 s
t = 12470.0 s
t = 12480.0 s
t = 12490.0 s
tn = 750 s
No HMI fits file supplied. A new HMI map object will not be created.
yt : [INFO ] 2017-06-20 14:08:19,303 Parameters: current_time = 0.0
yt : [INFO ] 2017-06-20 14:08:19,304 Parameters: domain_dimensions = [96 96 96]
yt : [INFO ] 2017-06-20 14:08:19,305 Parameters: domain_left_edge = [ 8.88147851e+09 1.66017386e+09 7.23009049e+08]
yt : [INFO ] 2017-06-20 14:08:19,306 Parameters: domain_right_edge = [ 4.35656994e+10 3.67056887e+10 3.54274434e+10]
yt : [INFO ] 2017-06-20 14:08:19,307 Parameters: cosmological_simulation = 0.0
t = 7500.0 s
t = 7510.0 s
t = 7520.0 s
t = 7530.0 s
t = 7540.0 s
t = 7550.0 s
t = 7560.0 s
t = 7570.0 s
t = 7580.0 s
t = 7590.0 s
t = 7600.0 s
t = 7610.0 s
t = 7620.0 s
t = 7630.0 s
t = 7640.0 s
t = 7650.0 s
t = 7660.0 s
t = 7670.0 s
t = 7680.0 s
t = 7690.0 s
t = 7700.0 s
t = 7710.0 s
t = 7720.0 s
t = 7730.0 s
t = 7740.0 s
t = 7750.0 s
t = 7760.0 s
t = 7770.0 s
t = 7780.0 s
t = 7790.0 s
t = 7800.0 s
t = 7810.0 s
t = 7820.0 s
t = 7830.0 s
t = 7840.0 s
t = 7850.0 s
t = 7860.0 s
t = 7870.0 s
t = 7880.0 s
t = 7890.0 s
t = 7900.0 s
t = 7910.0 s
t = 7920.0 s
t = 7930.0 s
t = 7940.0 s
t = 7950.0 s
t = 7960.0 s
t = 7970.0 s
t = 7980.0 s
t = 7990.0 s
t = 8000.0 s
t = 8010.0 s
t = 8020.0 s
t = 8030.0 s
t = 8040.0 s
t = 8050.0 s
t = 8060.0 s
t = 8070.0 s
t = 8080.0 s
t = 8090.0 s
t = 8100.0 s
t = 8110.0 s
t = 8120.0 s
t = 8130.0 s
t = 8140.0 s
t = 8150.0 s
t = 8160.0 s
t = 8170.0 s
t = 8180.0 s
t = 8190.0 s
t = 8200.0 s
t = 8210.0 s
t = 8220.0 s
t = 8230.0 s
t = 8240.0 s
t = 8250.0 s
t = 8260.0 s
t = 8270.0 s
t = 8280.0 s
t = 8290.0 s
t = 8300.0 s
t = 8310.0 s
t = 8320.0 s
t = 8330.0 s
t = 8340.0 s
t = 8350.0 s
t = 8360.0 s
t = 8370.0 s
t = 8380.0 s
t = 8390.0 s
t = 8400.0 s
t = 8410.0 s
t = 8420.0 s
t = 8430.0 s
t = 8440.0 s
t = 8450.0 s
t = 8460.0 s
t = 8470.0 s
t = 8480.0 s
t = 8490.0 s
t = 8500.0 s
t = 8510.0 s
t = 8520.0 s
t = 8530.0 s
t = 8540.0 s
t = 8550.0 s
t = 8560.0 s
t = 8570.0 s
t = 8580.0 s
t = 8590.0 s
t = 8600.0 s
t = 8610.0 s
t = 8620.0 s
t = 8630.0 s
t = 8640.0 s
t = 8650.0 s
t = 8660.0 s
t = 8670.0 s
t = 8680.0 s
t = 8690.0 s
t = 8700.0 s
t = 8710.0 s
t = 8720.0 s
t = 8730.0 s
t = 8740.0 s
t = 8750.0 s
t = 8760.0 s
t = 8770.0 s
t = 8780.0 s
t = 8790.0 s
t = 8800.0 s
t = 8810.0 s
t = 8820.0 s
t = 8830.0 s
t = 8840.0 s
t = 8850.0 s
t = 8860.0 s
t = 8870.0 s
t = 8880.0 s
t = 8890.0 s
t = 8900.0 s
t = 8910.0 s
t = 8920.0 s
t = 8930.0 s
t = 8940.0 s
t = 8950.0 s
t = 8960.0 s
t = 8970.0 s
t = 8980.0 s
t = 8990.0 s
t = 9000.0 s
t = 9010.0 s
t = 9020.0 s
t = 9030.0 s
t = 9040.0 s
t = 9050.0 s
t = 9060.0 s
t = 9070.0 s
t = 9080.0 s
t = 9090.0 s
t = 9100.0 s
t = 9110.0 s
t = 9120.0 s
t = 9130.0 s
t = 9140.0 s
t = 9150.0 s
t = 9160.0 s
t = 9170.0 s
t = 9180.0 s
t = 9190.0 s
t = 9200.0 s
t = 9210.0 s
t = 9220.0 s
t = 9230.0 s
t = 9240.0 s
t = 9250.0 s
t = 9260.0 s
t = 9270.0 s
t = 9280.0 s
t = 9290.0 s
t = 9300.0 s
t = 9310.0 s
t = 9320.0 s
t = 9330.0 s
t = 9340.0 s
t = 9350.0 s
t = 9360.0 s
t = 9370.0 s
t = 9380.0 s
t = 9390.0 s
t = 9400.0 s
t = 9410.0 s
t = 9420.0 s
t = 9430.0 s
t = 9440.0 s
t = 9450.0 s
t = 9460.0 s
t = 9470.0 s
t = 9480.0 s
t = 9490.0 s
t = 9500.0 s
t = 9510.0 s
t = 9520.0 s
t = 9530.0 s
t = 9540.0 s
t = 9550.0 s
t = 9560.0 s
t = 9570.0 s
t = 9580.0 s
t = 9590.0 s
t = 9600.0 s
t = 9610.0 s
t = 9620.0 s
t = 9630.0 s
t = 9640.0 s
t = 9650.0 s
t = 9660.0 s
t = 9670.0 s
t = 9680.0 s
t = 9690.0 s
t = 9700.0 s
t = 9710.0 s
t = 9720.0 s
t = 9730.0 s
t = 9740.0 s
t = 9750.0 s
t = 9760.0 s
t = 9770.0 s
t = 9780.0 s
t = 9790.0 s
t = 9800.0 s
t = 9810.0 s
t = 9820.0 s
t = 9830.0 s
t = 9840.0 s
t = 9850.0 s
t = 9860.0 s
t = 9870.0 s
t = 9880.0 s
t = 9890.0 s
t = 9900.0 s
t = 9910.0 s
t = 9920.0 s
t = 9930.0 s
t = 9940.0 s
t = 9950.0 s
t = 9960.0 s
t = 9970.0 s
t = 9980.0 s
t = 9990.0 s
t = 10000.0 s
t = 10010.0 s
t = 10020.0 s
t = 10030.0 s
t = 10040.0 s
t = 10050.0 s
t = 10060.0 s
t = 10070.0 s
t = 10080.0 s
t = 10090.0 s
t = 10100.0 s
t = 10110.0 s
t = 10120.0 s
t = 10130.0 s
t = 10140.0 s
t = 10150.0 s
t = 10160.0 s
t = 10170.0 s
t = 10180.0 s
t = 10190.0 s
t = 10200.0 s
t = 10210.0 s
t = 10220.0 s
t = 10230.0 s
t = 10240.0 s
t = 10250.0 s
t = 10260.0 s
t = 10270.0 s
t = 10280.0 s
t = 10290.0 s
t = 10300.0 s
t = 10310.0 s
t = 10320.0 s
t = 10330.0 s
t = 10340.0 s
t = 10350.0 s
t = 10360.0 s
t = 10370.0 s
t = 10380.0 s
t = 10390.0 s
t = 10400.0 s
t = 10410.0 s
t = 10420.0 s
t = 10430.0 s
t = 10440.0 s
t = 10450.0 s
t = 10460.0 s
t = 10470.0 s
t = 10480.0 s
t = 10490.0 s
t = 10500.0 s
t = 10510.0 s
t = 10520.0 s
t = 10530.0 s
t = 10540.0 s
t = 10550.0 s
t = 10560.0 s
t = 10570.0 s
t = 10580.0 s
t = 10590.0 s
t = 10600.0 s
t = 10610.0 s
t = 10620.0 s
t = 10630.0 s
t = 10640.0 s
t = 10650.0 s
t = 10660.0 s
t = 10670.0 s
t = 10680.0 s
t = 10690.0 s
t = 10700.0 s
t = 10710.0 s
t = 10720.0 s
t = 10730.0 s
t = 10740.0 s
t = 10750.0 s
t = 10760.0 s
t = 10770.0 s
t = 10780.0 s
t = 10790.0 s
t = 10800.0 s
t = 10810.0 s
t = 10820.0 s
t = 10830.0 s
t = 10840.0 s
t = 10850.0 s
t = 10860.0 s
t = 10870.0 s
t = 10880.0 s
t = 10890.0 s
t = 10900.0 s
t = 10910.0 s
t = 10920.0 s
t = 10930.0 s
t = 10940.0 s
t = 10950.0 s
t = 10960.0 s
t = 10970.0 s
t = 10980.0 s
t = 10990.0 s
t = 11000.0 s
t = 11010.0 s
t = 11020.0 s
t = 11030.0 s
t = 11040.0 s
t = 11050.0 s
t = 11060.0 s
t = 11070.0 s
t = 11080.0 s
t = 11090.0 s
t = 11100.0 s
t = 11110.0 s
t = 11120.0 s
t = 11130.0 s
t = 11140.0 s
t = 11150.0 s
t = 11160.0 s
t = 11170.0 s
t = 11180.0 s
t = 11190.0 s
t = 11200.0 s
t = 11210.0 s
t = 11220.0 s
t = 11230.0 s
t = 11240.0 s
t = 11250.0 s
t = 11260.0 s
t = 11270.0 s
t = 11280.0 s
t = 11290.0 s
t = 11300.0 s
t = 11310.0 s
t = 11320.0 s
t = 11330.0 s
t = 11340.0 s
t = 11350.0 s
t = 11360.0 s
t = 11370.0 s
t = 11380.0 s
t = 11390.0 s
t = 11400.0 s
t = 11410.0 s
t = 11420.0 s
t = 11430.0 s
t = 11440.0 s
t = 11450.0 s
t = 11460.0 s
t = 11470.0 s
t = 11480.0 s
t = 11490.0 s
t = 11500.0 s
t = 11510.0 s
t = 11520.0 s
t = 11530.0 s
t = 11540.0 s
t = 11550.0 s
t = 11560.0 s
t = 11570.0 s
t = 11580.0 s
t = 11590.0 s
t = 11600.0 s
t = 11610.0 s
t = 11620.0 s
t = 11630.0 s
t = 11640.0 s
t = 11650.0 s
t = 11660.0 s
t = 11670.0 s
t = 11680.0 s
t = 11690.0 s
t = 11700.0 s
t = 11710.0 s
t = 11720.0 s
t = 11730.0 s
t = 11740.0 s
t = 11750.0 s
t = 11760.0 s
t = 11770.0 s
t = 11780.0 s
t = 11790.0 s
t = 11800.0 s
t = 11810.0 s
t = 11820.0 s
t = 11830.0 s
t = 11840.0 s
t = 11850.0 s
t = 11860.0 s
t = 11870.0 s
t = 11880.0 s
t = 11890.0 s
t = 11900.0 s
t = 11910.0 s
t = 11920.0 s
t = 11930.0 s
t = 11940.0 s
t = 11950.0 s
t = 11960.0 s
t = 11970.0 s
t = 11980.0 s
t = 11990.0 s
t = 12000.0 s
t = 12010.0 s
t = 12020.0 s
t = 12030.0 s
t = 12040.0 s
t = 12050.0 s
t = 12060.0 s
t = 12070.0 s
t = 12080.0 s
t = 12090.0 s
t = 12100.0 s
t = 12110.0 s
t = 12120.0 s
t = 12130.0 s
t = 12140.0 s
t = 12150.0 s
t = 12160.0 s
t = 12170.0 s
t = 12180.0 s
t = 12190.0 s
t = 12200.0 s
t = 12210.0 s
t = 12220.0 s
t = 12230.0 s
t = 12240.0 s
t = 12250.0 s
t = 12260.0 s
t = 12270.0 s
t = 12280.0 s
t = 12290.0 s
t = 12300.0 s
t = 12310.0 s
t = 12320.0 s
t = 12330.0 s
t = 12340.0 s
t = 12350.0 s
t = 12360.0 s
t = 12370.0 s
t = 12380.0 s
t = 12390.0 s
t = 12400.0 s
t = 12410.0 s
t = 12420.0 s
t = 12430.0 s
t = 12440.0 s
t = 12450.0 s
t = 12460.0 s
t = 12470.0 s
t = 12480.0 s
t = 12490.0 s
tn = 750-ion s
No HMI fits file supplied. A new HMI map object will not be created.
yt : [INFO ] 2017-06-20 14:11:01,603 Parameters: current_time = 0.0
yt : [INFO ] 2017-06-20 14:11:01,603 Parameters: domain_dimensions = [96 96 96]
yt : [INFO ] 2017-06-20 14:11:01,604 Parameters: domain_left_edge = [ 8.88147851e+09 1.66017386e+09 7.23009049e+08]
yt : [INFO ] 2017-06-20 14:11:01,605 Parameters: domain_right_edge = [ 4.35656994e+10 3.67056887e+10 3.54274434e+10]
yt : [INFO ] 2017-06-20 14:11:01,606 Parameters: cosmological_simulation = 0.0
t = 7500.0 s
t = 7510.0 s
t = 7520.0 s
t = 7530.0 s
t = 7540.0 s
t = 7550.0 s
t = 7560.0 s
t = 7570.0 s
t = 7580.0 s
t = 7590.0 s
t = 7600.0 s
t = 7610.0 s
t = 7620.0 s
t = 7630.0 s
t = 7640.0 s
t = 7650.0 s
t = 7660.0 s
t = 7670.0 s
t = 7680.0 s
t = 7690.0 s
t = 7700.0 s
t = 7710.0 s
t = 7720.0 s
t = 7730.0 s
t = 7740.0 s
t = 7750.0 s
t = 7760.0 s
t = 7770.0 s
t = 7780.0 s
t = 7790.0 s
t = 7800.0 s
t = 7810.0 s
t = 7820.0 s
t = 7830.0 s
t = 7840.0 s
t = 7850.0 s
t = 7860.0 s
t = 7870.0 s
t = 7880.0 s
t = 7890.0 s
t = 7900.0 s
t = 7910.0 s
t = 7920.0 s
t = 7930.0 s
t = 7940.0 s
t = 7950.0 s
t = 7960.0 s
t = 7970.0 s
t = 7980.0 s
t = 7990.0 s
t = 8000.0 s
t = 8010.0 s
t = 8020.0 s
t = 8030.0 s
t = 8040.0 s
t = 8050.0 s
t = 8060.0 s
t = 8070.0 s
t = 8080.0 s
t = 8090.0 s
t = 8100.0 s
t = 8110.0 s
t = 8120.0 s
t = 8130.0 s
t = 8140.0 s
t = 8150.0 s
t = 8160.0 s
t = 8170.0 s
t = 8180.0 s
t = 8190.0 s
t = 8200.0 s
t = 8210.0 s
t = 8220.0 s
t = 8230.0 s
t = 8240.0 s
t = 8250.0 s
t = 8260.0 s
t = 8270.0 s
t = 8280.0 s
t = 8290.0 s
t = 8300.0 s
t = 8310.0 s
t = 8320.0 s
t = 8330.0 s
t = 8340.0 s
t = 8350.0 s
t = 8360.0 s
t = 8370.0 s
t = 8380.0 s
t = 8390.0 s
t = 8400.0 s
t = 8410.0 s
t = 8420.0 s
t = 8430.0 s
t = 8440.0 s
t = 8450.0 s
t = 8460.0 s
t = 8470.0 s
t = 8480.0 s
t = 8490.0 s
t = 8500.0 s
t = 8510.0 s
t = 8520.0 s
t = 8530.0 s
t = 8540.0 s
t = 8550.0 s
t = 8560.0 s
t = 8570.0 s
t = 8580.0 s
t = 8590.0 s
t = 8600.0 s
t = 8610.0 s
t = 8620.0 s
t = 8630.0 s
t = 8640.0 s
t = 8650.0 s
t = 8660.0 s
t = 8670.0 s
t = 8680.0 s
t = 8690.0 s
t = 8700.0 s
t = 8710.0 s
t = 8720.0 s
t = 8730.0 s
t = 8740.0 s
t = 8750.0 s
t = 8760.0 s
t = 8770.0 s
t = 8780.0 s
t = 8790.0 s
t = 8800.0 s
t = 8810.0 s
t = 8820.0 s
t = 8830.0 s
t = 8840.0 s
t = 8850.0 s
t = 8860.0 s
t = 8870.0 s
t = 8880.0 s
t = 8890.0 s
t = 8900.0 s
t = 8910.0 s
t = 8920.0 s
t = 8930.0 s
t = 8940.0 s
t = 8950.0 s
t = 8960.0 s
t = 8970.0 s
t = 8980.0 s
t = 8990.0 s
t = 9000.0 s
t = 9010.0 s
t = 9020.0 s
t = 9030.0 s
t = 9040.0 s
t = 9050.0 s
t = 9060.0 s
t = 9070.0 s
t = 9080.0 s
t = 9090.0 s
t = 9100.0 s
t = 9110.0 s
t = 9120.0 s
t = 9130.0 s
t = 9140.0 s
t = 9150.0 s
t = 9160.0 s
t = 9170.0 s
t = 9180.0 s
t = 9190.0 s
t = 9200.0 s
t = 9210.0 s
t = 9220.0 s
t = 9230.0 s
t = 9240.0 s
t = 9250.0 s
t = 9260.0 s
t = 9270.0 s
t = 9280.0 s
t = 9290.0 s
t = 9300.0 s
t = 9310.0 s
t = 9320.0 s
t = 9330.0 s
t = 9340.0 s
t = 9350.0 s
t = 9360.0 s
t = 9370.0 s
t = 9380.0 s
t = 9390.0 s
t = 9400.0 s
t = 9410.0 s
t = 9420.0 s
t = 9430.0 s
t = 9440.0 s
t = 9450.0 s
t = 9460.0 s
t = 9470.0 s
t = 9480.0 s
t = 9490.0 s
t = 9500.0 s
t = 9510.0 s
t = 9520.0 s
t = 9530.0 s
t = 9540.0 s
t = 9550.0 s
t = 9560.0 s
t = 9570.0 s
t = 9580.0 s
t = 9590.0 s
t = 9600.0 s
t = 9610.0 s
t = 9620.0 s
t = 9630.0 s
t = 9640.0 s
t = 9650.0 s
t = 9660.0 s
t = 9670.0 s
t = 9680.0 s
t = 9690.0 s
t = 9700.0 s
t = 9710.0 s
t = 9720.0 s
t = 9730.0 s
t = 9740.0 s
t = 9750.0 s
t = 9760.0 s
t = 9770.0 s
t = 9780.0 s
t = 9790.0 s
t = 9800.0 s
t = 9810.0 s
t = 9820.0 s
t = 9830.0 s
t = 9840.0 s
t = 9850.0 s
t = 9860.0 s
t = 9870.0 s
t = 9880.0 s
t = 9890.0 s
t = 9900.0 s
t = 9910.0 s
t = 9920.0 s
t = 9930.0 s
t = 9940.0 s
t = 9950.0 s
t = 9960.0 s
t = 9970.0 s
t = 9980.0 s
t = 9990.0 s
t = 10000.0 s
t = 10010.0 s
t = 10020.0 s
t = 10030.0 s
t = 10040.0 s
t = 10050.0 s
t = 10060.0 s
t = 10070.0 s
t = 10080.0 s
t = 10090.0 s
t = 10100.0 s
t = 10110.0 s
t = 10120.0 s
t = 10130.0 s
t = 10140.0 s
t = 10150.0 s
t = 10160.0 s
t = 10170.0 s
t = 10180.0 s
t = 10190.0 s
t = 10200.0 s
t = 10210.0 s
t = 10220.0 s
t = 10230.0 s
t = 10240.0 s
t = 10250.0 s
t = 10260.0 s
t = 10270.0 s
t = 10280.0 s
t = 10290.0 s
t = 10300.0 s
t = 10310.0 s
t = 10320.0 s
t = 10330.0 s
t = 10340.0 s
t = 10350.0 s
t = 10360.0 s
t = 10370.0 s
t = 10380.0 s
t = 10390.0 s
t = 10400.0 s
t = 10410.0 s
t = 10420.0 s
t = 10430.0 s
t = 10440.0 s
t = 10450.0 s
t = 10460.0 s
t = 10470.0 s
t = 10480.0 s
t = 10490.0 s
t = 10500.0 s
t = 10510.0 s
t = 10520.0 s
t = 10530.0 s
t = 10540.0 s
t = 10550.0 s
t = 10560.0 s
t = 10570.0 s
t = 10580.0 s
t = 10590.0 s
t = 10600.0 s
t = 10610.0 s
t = 10620.0 s
t = 10630.0 s
t = 10640.0 s
t = 10650.0 s
t = 10660.0 s
t = 10670.0 s
t = 10680.0 s
t = 10690.0 s
t = 10700.0 s
t = 10710.0 s
t = 10720.0 s
t = 10730.0 s
t = 10740.0 s
t = 10750.0 s
t = 10760.0 s
t = 10770.0 s
t = 10780.0 s
t = 10790.0 s
t = 10800.0 s
t = 10810.0 s
t = 10820.0 s
t = 10830.0 s
t = 10840.0 s
t = 10850.0 s
t = 10860.0 s
t = 10870.0 s
t = 10880.0 s
t = 10890.0 s
t = 10900.0 s
t = 10910.0 s
t = 10920.0 s
t = 10930.0 s
t = 10940.0 s
t = 10950.0 s
t = 10960.0 s
t = 10970.0 s
t = 10980.0 s
t = 10990.0 s
t = 11000.0 s
t = 11010.0 s
t = 11020.0 s
t = 11030.0 s
t = 11040.0 s
t = 11050.0 s
t = 11060.0 s
t = 11070.0 s
t = 11080.0 s
t = 11090.0 s
t = 11100.0 s
t = 11110.0 s
t = 11120.0 s
t = 11130.0 s
t = 11140.0 s
t = 11150.0 s
t = 11160.0 s
t = 11170.0 s
t = 11180.0 s
t = 11190.0 s
t = 11200.0 s
t = 11210.0 s
t = 11220.0 s
t = 11230.0 s
t = 11240.0 s
t = 11250.0 s
t = 11260.0 s
t = 11270.0 s
t = 11280.0 s
t = 11290.0 s
t = 11300.0 s
t = 11310.0 s
t = 11320.0 s
t = 11330.0 s
t = 11340.0 s
t = 11350.0 s
t = 11360.0 s
t = 11370.0 s
t = 11380.0 s
t = 11390.0 s
t = 11400.0 s
t = 11410.0 s
t = 11420.0 s
t = 11430.0 s
t = 11440.0 s
t = 11450.0 s
t = 11460.0 s
t = 11470.0 s
t = 11480.0 s
t = 11490.0 s
t = 11500.0 s
t = 11510.0 s
t = 11520.0 s
t = 11530.0 s
t = 11540.0 s
t = 11550.0 s
t = 11560.0 s
t = 11570.0 s
t = 11580.0 s
t = 11590.0 s
t = 11600.0 s
t = 11610.0 s
t = 11620.0 s
t = 11630.0 s
t = 11640.0 s
t = 11650.0 s
t = 11660.0 s
t = 11670.0 s
t = 11680.0 s
t = 11690.0 s
t = 11700.0 s
t = 11710.0 s
t = 11720.0 s
t = 11730.0 s
t = 11740.0 s
t = 11750.0 s
t = 11760.0 s
t = 11770.0 s
t = 11780.0 s
t = 11790.0 s
t = 11800.0 s
t = 11810.0 s
t = 11820.0 s
t = 11830.0 s
t = 11840.0 s
t = 11850.0 s
t = 11860.0 s
t = 11870.0 s
t = 11880.0 s
t = 11890.0 s
t = 11900.0 s
t = 11910.0 s
t = 11920.0 s
t = 11930.0 s
t = 11940.0 s
t = 11950.0 s
t = 11960.0 s
t = 11970.0 s
t = 11980.0 s
t = 11990.0 s
t = 12000.0 s
t = 12010.0 s
t = 12020.0 s
t = 12030.0 s
t = 12040.0 s
t = 12050.0 s
t = 12060.0 s
t = 12070.0 s
t = 12080.0 s
t = 12090.0 s
t = 12100.0 s
t = 12110.0 s
t = 12120.0 s
t = 12130.0 s
t = 12140.0 s
t = 12150.0 s
t = 12160.0 s
t = 12170.0 s
t = 12180.0 s
t = 12190.0 s
t = 12200.0 s
t = 12210.0 s
t = 12220.0 s
t = 12230.0 s
t = 12240.0 s
t = 12250.0 s
t = 12260.0 s
t = 12270.0 s
t = 12280.0 s
t = 12290.0 s
t = 12300.0 s
t = 12310.0 s
t = 12320.0 s
t = 12330.0 s
t = 12340.0 s
t = 12350.0 s
t = 12360.0 s
t = 12370.0 s
t = 12380.0 s
t = 12390.0 s
t = 12400.0 s
t = 12410.0 s
t = 12420.0 s
t = 12430.0 s
t = 12440.0 s
t = 12450.0 s
t = 12460.0 s
t = 12470.0 s
t = 12480.0 s
t = 12490.0 s
tn = 2500 s
No HMI fits file supplied. A new HMI map object will not be created.
yt : [INFO ] 2017-06-20 14:13:37,839 Parameters: current_time = 0.0
yt : [INFO ] 2017-06-20 14:13:37,840 Parameters: domain_dimensions = [96 96 96]
yt : [INFO ] 2017-06-20 14:13:37,842 Parameters: domain_left_edge = [ 8.88147851e+09 1.66017386e+09 7.23009049e+08]
yt : [INFO ] 2017-06-20 14:13:37,843 Parameters: domain_right_edge = [ 4.35656994e+10 3.67056887e+10 3.54274434e+10]
yt : [INFO ] 2017-06-20 14:13:37,844 Parameters: cosmological_simulation = 0.0
t = 7500.0 s
t = 7510.0 s
t = 7520.0 s
t = 7530.0 s
t = 7540.0 s
t = 7550.0 s
t = 7560.0 s
t = 7570.0 s
t = 7580.0 s
t = 7590.0 s
t = 7600.0 s
t = 7610.0 s
t = 7620.0 s
t = 7630.0 s
t = 7640.0 s
t = 7650.0 s
t = 7660.0 s
t = 7670.0 s
t = 7680.0 s
t = 7690.0 s
t = 7700.0 s
t = 7710.0 s
t = 7720.0 s
t = 7730.0 s
t = 7740.0 s
t = 7750.0 s
t = 7760.0 s
t = 7770.0 s
t = 7780.0 s
t = 7790.0 s
t = 7800.0 s
t = 7810.0 s
t = 7820.0 s
t = 7830.0 s
t = 7840.0 s
t = 7850.0 s
t = 7860.0 s
t = 7870.0 s
t = 7880.0 s
t = 7890.0 s
t = 7900.0 s
t = 7910.0 s
t = 7920.0 s
t = 7930.0 s
t = 7940.0 s
t = 7950.0 s
t = 7960.0 s
t = 7970.0 s
t = 7980.0 s
t = 7990.0 s
t = 8000.0 s
t = 8010.0 s
t = 8020.0 s
t = 8030.0 s
t = 8040.0 s
t = 8050.0 s
t = 8060.0 s
t = 8070.0 s
t = 8080.0 s
t = 8090.0 s
t = 8100.0 s
t = 8110.0 s
t = 8120.0 s
t = 8130.0 s
t = 8140.0 s
t = 8150.0 s
t = 8160.0 s
t = 8170.0 s
t = 8180.0 s
t = 8190.0 s
t = 8200.0 s
t = 8210.0 s
t = 8220.0 s
t = 8230.0 s
t = 8240.0 s
t = 8250.0 s
t = 8260.0 s
t = 8270.0 s
t = 8280.0 s
t = 8290.0 s
t = 8300.0 s
t = 8310.0 s
t = 8320.0 s
t = 8330.0 s
t = 8340.0 s
t = 8350.0 s
t = 8360.0 s
t = 8370.0 s
t = 8380.0 s
t = 8390.0 s
t = 8400.0 s
t = 8410.0 s
t = 8420.0 s
t = 8430.0 s
t = 8440.0 s
t = 8450.0 s
t = 8460.0 s
t = 8470.0 s
t = 8480.0 s
t = 8490.0 s
t = 8500.0 s
t = 8510.0 s
t = 8520.0 s
t = 8530.0 s
t = 8540.0 s
t = 8550.0 s
t = 8560.0 s
t = 8570.0 s
t = 8580.0 s
t = 8590.0 s
t = 8600.0 s
t = 8610.0 s
t = 8620.0 s
t = 8630.0 s
t = 8640.0 s
t = 8650.0 s
t = 8660.0 s
t = 8670.0 s
t = 8680.0 s
t = 8690.0 s
t = 8700.0 s
t = 8710.0 s
t = 8720.0 s
t = 8730.0 s
t = 8740.0 s
t = 8750.0 s
t = 8760.0 s
t = 8770.0 s
t = 8780.0 s
t = 8790.0 s
t = 8800.0 s
t = 8810.0 s
t = 8820.0 s
t = 8830.0 s
t = 8840.0 s
t = 8850.0 s
t = 8860.0 s
t = 8870.0 s
t = 8880.0 s
t = 8890.0 s
t = 8900.0 s
t = 8910.0 s
t = 8920.0 s
t = 8930.0 s
t = 8940.0 s
t = 8950.0 s
t = 8960.0 s
t = 8970.0 s
t = 8980.0 s
t = 8990.0 s
t = 9000.0 s
t = 9010.0 s
t = 9020.0 s
t = 9030.0 s
t = 9040.0 s
t = 9050.0 s
t = 9060.0 s
t = 9070.0 s
t = 9080.0 s
t = 9090.0 s
t = 9100.0 s
t = 9110.0 s
t = 9120.0 s
t = 9130.0 s
t = 9140.0 s
t = 9150.0 s
t = 9160.0 s
t = 9170.0 s
t = 9180.0 s
t = 9190.0 s
t = 9200.0 s
t = 9210.0 s
t = 9220.0 s
t = 9230.0 s
t = 9240.0 s
t = 9250.0 s
t = 9260.0 s
t = 9270.0 s
t = 9280.0 s
t = 9290.0 s
t = 9300.0 s
t = 9310.0 s
t = 9320.0 s
t = 9330.0 s
t = 9340.0 s
t = 9350.0 s
t = 9360.0 s
t = 9370.0 s
t = 9380.0 s
t = 9390.0 s
t = 9400.0 s
t = 9410.0 s
t = 9420.0 s
t = 9430.0 s
t = 9440.0 s
t = 9450.0 s
t = 9460.0 s
t = 9470.0 s
t = 9480.0 s
t = 9490.0 s
t = 9500.0 s
t = 9510.0 s
t = 9520.0 s
t = 9530.0 s
t = 9540.0 s
t = 9550.0 s
t = 9560.0 s
t = 9570.0 s
t = 9580.0 s
t = 9590.0 s
t = 9600.0 s
t = 9610.0 s
t = 9620.0 s
t = 9630.0 s
t = 9640.0 s
t = 9650.0 s
t = 9660.0 s
t = 9670.0 s
t = 9680.0 s
t = 9690.0 s
t = 9700.0 s
t = 9710.0 s
t = 9720.0 s
t = 9730.0 s
t = 9740.0 s
t = 9750.0 s
t = 9760.0 s
t = 9770.0 s
t = 9780.0 s
t = 9790.0 s
t = 9800.0 s
t = 9810.0 s
t = 9820.0 s
t = 9830.0 s
t = 9840.0 s
t = 9850.0 s
t = 9860.0 s
t = 9870.0 s
t = 9880.0 s
t = 9890.0 s
t = 9900.0 s
t = 9910.0 s
t = 9920.0 s
t = 9930.0 s
t = 9940.0 s
t = 9950.0 s
t = 9960.0 s
t = 9970.0 s
t = 9980.0 s
t = 9990.0 s
t = 10000.0 s
t = 10010.0 s
t = 10020.0 s
t = 10030.0 s
t = 10040.0 s
t = 10050.0 s
t = 10060.0 s
t = 10070.0 s
t = 10080.0 s
t = 10090.0 s
t = 10100.0 s
t = 10110.0 s
t = 10120.0 s
t = 10130.0 s
t = 10140.0 s
t = 10150.0 s
t = 10160.0 s
t = 10170.0 s
t = 10180.0 s
t = 10190.0 s
t = 10200.0 s
t = 10210.0 s
t = 10220.0 s
t = 10230.0 s
t = 10240.0 s
t = 10250.0 s
t = 10260.0 s
t = 10270.0 s
t = 10280.0 s
t = 10290.0 s
t = 10300.0 s
t = 10310.0 s
t = 10320.0 s
t = 10330.0 s
t = 10340.0 s
t = 10350.0 s
t = 10360.0 s
t = 10370.0 s
t = 10380.0 s
t = 10390.0 s
t = 10400.0 s
t = 10410.0 s
t = 10420.0 s
t = 10430.0 s
t = 10440.0 s
t = 10450.0 s
t = 10460.0 s
t = 10470.0 s
t = 10480.0 s
t = 10490.0 s
t = 10500.0 s
t = 10510.0 s
t = 10520.0 s
t = 10530.0 s
t = 10540.0 s
t = 10550.0 s
t = 10560.0 s
t = 10570.0 s
t = 10580.0 s
t = 10590.0 s
t = 10600.0 s
t = 10610.0 s
t = 10620.0 s
t = 10630.0 s
t = 10640.0 s
t = 10650.0 s
t = 10660.0 s
t = 10670.0 s
t = 10680.0 s
t = 10690.0 s
t = 10700.0 s
t = 10710.0 s
t = 10720.0 s
t = 10730.0 s
t = 10740.0 s
t = 10750.0 s
t = 10760.0 s
t = 10770.0 s
t = 10780.0 s
t = 10790.0 s
t = 10800.0 s
t = 10810.0 s
t = 10820.0 s
t = 10830.0 s
t = 10840.0 s
t = 10850.0 s
t = 10860.0 s
t = 10870.0 s
t = 10880.0 s
t = 10890.0 s
t = 10900.0 s
t = 10910.0 s
t = 10920.0 s
t = 10930.0 s
t = 10940.0 s
t = 10950.0 s
t = 10960.0 s
t = 10970.0 s
t = 10980.0 s
t = 10990.0 s
t = 11000.0 s
t = 11010.0 s
t = 11020.0 s
t = 11030.0 s
t = 11040.0 s
t = 11050.0 s
t = 11060.0 s
t = 11070.0 s
t = 11080.0 s
t = 11090.0 s
t = 11100.0 s
t = 11110.0 s
t = 11120.0 s
t = 11130.0 s
t = 11140.0 s
t = 11150.0 s
t = 11160.0 s
t = 11170.0 s
t = 11180.0 s
t = 11190.0 s
t = 11200.0 s
t = 11210.0 s
t = 11220.0 s
t = 11230.0 s
t = 11240.0 s
t = 11250.0 s
t = 11260.0 s
t = 11270.0 s
t = 11280.0 s
t = 11290.0 s
t = 11300.0 s
t = 11310.0 s
t = 11320.0 s
t = 11330.0 s
t = 11340.0 s
t = 11350.0 s
t = 11360.0 s
t = 11370.0 s
t = 11380.0 s
t = 11390.0 s
t = 11400.0 s
t = 11410.0 s
t = 11420.0 s
t = 11430.0 s
t = 11440.0 s
t = 11450.0 s
t = 11460.0 s
t = 11470.0 s
t = 11480.0 s
t = 11490.0 s
t = 11500.0 s
t = 11510.0 s
t = 11520.0 s
t = 11530.0 s
t = 11540.0 s
t = 11550.0 s
t = 11560.0 s
t = 11570.0 s
t = 11580.0 s
t = 11590.0 s
t = 11600.0 s
t = 11610.0 s
t = 11620.0 s
t = 11630.0 s
t = 11640.0 s
t = 11650.0 s
t = 11660.0 s
t = 11670.0 s
t = 11680.0 s
t = 11690.0 s
t = 11700.0 s
t = 11710.0 s
t = 11720.0 s
t = 11730.0 s
t = 11740.0 s
t = 11750.0 s
t = 11760.0 s
t = 11770.0 s
t = 11780.0 s
t = 11790.0 s
t = 11800.0 s
t = 11810.0 s
t = 11820.0 s
t = 11830.0 s
t = 11840.0 s
t = 11850.0 s
t = 11860.0 s
t = 11870.0 s
t = 11880.0 s
t = 11890.0 s
t = 11900.0 s
t = 11910.0 s
t = 11920.0 s
t = 11930.0 s
t = 11940.0 s
t = 11950.0 s
t = 11960.0 s
t = 11970.0 s
t = 11980.0 s
t = 11990.0 s
t = 12000.0 s
t = 12010.0 s
t = 12020.0 s
t = 12030.0 s
t = 12040.0 s
t = 12050.0 s
t = 12060.0 s
t = 12070.0 s
t = 12080.0 s
t = 12090.0 s
t = 12100.0 s
t = 12110.0 s
t = 12120.0 s
t = 12130.0 s
t = 12140.0 s
t = 12150.0 s
t = 12160.0 s
t = 12170.0 s
t = 12180.0 s
t = 12190.0 s
t = 12200.0 s
t = 12210.0 s
t = 12220.0 s
t = 12230.0 s
t = 12240.0 s
t = 12250.0 s
t = 12260.0 s
t = 12270.0 s
t = 12280.0 s
t = 12290.0 s
t = 12300.0 s
t = 12310.0 s
t = 12320.0 s
t = 12330.0 s
t = 12340.0 s
t = 12350.0 s
t = 12360.0 s
t = 12370.0 s
t = 12380.0 s
t = 12390.0 s
t = 12400.0 s
t = 12410.0 s
t = 12420.0 s
t = 12430.0 s
t = 12440.0 s
t = 12450.0 s
t = 12460.0 s
t = 12470.0 s
t = 12480.0 s
t = 12490.0 s
tn = 5000 s
No HMI fits file supplied. A new HMI map object will not be created.
yt : [INFO ] 2017-06-20 14:16:13,119 Parameters: current_time = 0.0
yt : [INFO ] 2017-06-20 14:16:13,119 Parameters: domain_dimensions = [96 96 96]
yt : [INFO ] 2017-06-20 14:16:13,120 Parameters: domain_left_edge = [ 8.88147851e+09 1.66017386e+09 7.23009049e+08]
yt : [INFO ] 2017-06-20 14:16:13,121 Parameters: domain_right_edge = [ 4.35656994e+10 3.67056887e+10 3.54274434e+10]
yt : [INFO ] 2017-06-20 14:16:13,123 Parameters: cosmological_simulation = 0.0
t = 7500.0 s
t = 7510.0 s
t = 7520.0 s
t = 7530.0 s
t = 7540.0 s
t = 7550.0 s
t = 7560.0 s
t = 7570.0 s
t = 7580.0 s
t = 7590.0 s
t = 7600.0 s
t = 7610.0 s
t = 7620.0 s
t = 7630.0 s
t = 7640.0 s
t = 7650.0 s
t = 7660.0 s
t = 7670.0 s
t = 7680.0 s
t = 7690.0 s
t = 7700.0 s
t = 7710.0 s
t = 7720.0 s
t = 7730.0 s
t = 7740.0 s
t = 7750.0 s
t = 7760.0 s
t = 7770.0 s
t = 7780.0 s
t = 7790.0 s
t = 7800.0 s
t = 7810.0 s
t = 7820.0 s
t = 7830.0 s
t = 7840.0 s
t = 7850.0 s
t = 7860.0 s
t = 7870.0 s
t = 7880.0 s
t = 7890.0 s
t = 7900.0 s
t = 7910.0 s
t = 7920.0 s
t = 7930.0 s
t = 7940.0 s
t = 7950.0 s
t = 7960.0 s
t = 7970.0 s
t = 7980.0 s
t = 7990.0 s
t = 8000.0 s
t = 8010.0 s
t = 8020.0 s
t = 8030.0 s
t = 8040.0 s
t = 8050.0 s
t = 8060.0 s
t = 8070.0 s
t = 8080.0 s
t = 8090.0 s
t = 8100.0 s
t = 8110.0 s
t = 8120.0 s
t = 8130.0 s
t = 8140.0 s
t = 8150.0 s
t = 8160.0 s
t = 8170.0 s
t = 8180.0 s
t = 8190.0 s
t = 8200.0 s
t = 8210.0 s
t = 8220.0 s
t = 8230.0 s
t = 8240.0 s
t = 8250.0 s
t = 8260.0 s
t = 8270.0 s
t = 8280.0 s
t = 8290.0 s
t = 8300.0 s
t = 8310.0 s
t = 8320.0 s
t = 8330.0 s
t = 8340.0 s
t = 8350.0 s
t = 8360.0 s
t = 8370.0 s
t = 8380.0 s
t = 8390.0 s
t = 8400.0 s
t = 8410.0 s
t = 8420.0 s
t = 8430.0 s
t = 8440.0 s
t = 8450.0 s
t = 8460.0 s
t = 8470.0 s
t = 8480.0 s
t = 8490.0 s
t = 8500.0 s
t = 8510.0 s
t = 8520.0 s
t = 8530.0 s
t = 8540.0 s
t = 8550.0 s
t = 8560.0 s
t = 8570.0 s
t = 8580.0 s
t = 8590.0 s
t = 8600.0 s
t = 8610.0 s
t = 8620.0 s
t = 8630.0 s
t = 8640.0 s
t = 8650.0 s
t = 8660.0 s
t = 8670.0 s
t = 8680.0 s
t = 8690.0 s
t = 8700.0 s
t = 8710.0 s
t = 8720.0 s
t = 8730.0 s
t = 8740.0 s
t = 8750.0 s
t = 8760.0 s
t = 8770.0 s
t = 8780.0 s
t = 8790.0 s
t = 8800.0 s
t = 8810.0 s
t = 8820.0 s
t = 8830.0 s
t = 8840.0 s
t = 8850.0 s
t = 8860.0 s
t = 8870.0 s
t = 8880.0 s
t = 8890.0 s
t = 8900.0 s
t = 8910.0 s
t = 8920.0 s
t = 8930.0 s
t = 8940.0 s
t = 8950.0 s
t = 8960.0 s
t = 8970.0 s
t = 8980.0 s
t = 8990.0 s
t = 9000.0 s
t = 9010.0 s
t = 9020.0 s
t = 9030.0 s
t = 9040.0 s
t = 9050.0 s
t = 9060.0 s
t = 9070.0 s
t = 9080.0 s
t = 9090.0 s
t = 9100.0 s
t = 9110.0 s
t = 9120.0 s
t = 9130.0 s
t = 9140.0 s
t = 9150.0 s
t = 9160.0 s
t = 9170.0 s
t = 9180.0 s
t = 9190.0 s
t = 9200.0 s
t = 9210.0 s
t = 9220.0 s
t = 9230.0 s
t = 9240.0 s
t = 9250.0 s
t = 9260.0 s
t = 9270.0 s
t = 9280.0 s
t = 9290.0 s
t = 9300.0 s
t = 9310.0 s
t = 9320.0 s
t = 9330.0 s
t = 9340.0 s
t = 9350.0 s
t = 9360.0 s
t = 9370.0 s
t = 9380.0 s
t = 9390.0 s
t = 9400.0 s
t = 9410.0 s
t = 9420.0 s
t = 9430.0 s
t = 9440.0 s
t = 9450.0 s
t = 9460.0 s
t = 9470.0 s
t = 9480.0 s
t = 9490.0 s
t = 9500.0 s
t = 9510.0 s
t = 9520.0 s
t = 9530.0 s
t = 9540.0 s
t = 9550.0 s
t = 9560.0 s
t = 9570.0 s
t = 9580.0 s
t = 9590.0 s
t = 9600.0 s
t = 9610.0 s
t = 9620.0 s
t = 9630.0 s
t = 9640.0 s
t = 9650.0 s
t = 9660.0 s
t = 9670.0 s
t = 9680.0 s
t = 9690.0 s
t = 9700.0 s
t = 9710.0 s
t = 9720.0 s
t = 9730.0 s
t = 9740.0 s
t = 9750.0 s
t = 9760.0 s
t = 9770.0 s
t = 9780.0 s
t = 9790.0 s
t = 9800.0 s
t = 9810.0 s
t = 9820.0 s
t = 9830.0 s
t = 9840.0 s
t = 9850.0 s
t = 9860.0 s
t = 9870.0 s
t = 9880.0 s
t = 9890.0 s
t = 9900.0 s
t = 9910.0 s
t = 9920.0 s
t = 9930.0 s
t = 9940.0 s
t = 9950.0 s
t = 9960.0 s
t = 9970.0 s
t = 9980.0 s
t = 9990.0 s
t = 10000.0 s
t = 10010.0 s
t = 10020.0 s
t = 10030.0 s
t = 10040.0 s
t = 10050.0 s
t = 10060.0 s
t = 10070.0 s
t = 10080.0 s
t = 10090.0 s
t = 10100.0 s
t = 10110.0 s
t = 10120.0 s
t = 10130.0 s
t = 10140.0 s
t = 10150.0 s
t = 10160.0 s
t = 10170.0 s
t = 10180.0 s
t = 10190.0 s
t = 10200.0 s
t = 10210.0 s
t = 10220.0 s
t = 10230.0 s
t = 10240.0 s
t = 10250.0 s
t = 10260.0 s
t = 10270.0 s
t = 10280.0 s
t = 10290.0 s
t = 10300.0 s
t = 10310.0 s
t = 10320.0 s
t = 10330.0 s
t = 10340.0 s
t = 10350.0 s
t = 10360.0 s
t = 10370.0 s
t = 10380.0 s
t = 10390.0 s
t = 10400.0 s
t = 10410.0 s
t = 10420.0 s
t = 10430.0 s
t = 10440.0 s
t = 10450.0 s
t = 10460.0 s
t = 10470.0 s
t = 10480.0 s
t = 10490.0 s
t = 10500.0 s
t = 10510.0 s
t = 10520.0 s
t = 10530.0 s
t = 10540.0 s
t = 10550.0 s
t = 10560.0 s
t = 10570.0 s
t = 10580.0 s
t = 10590.0 s
t = 10600.0 s
t = 10610.0 s
t = 10620.0 s
t = 10630.0 s
t = 10640.0 s
t = 10650.0 s
t = 10660.0 s
t = 10670.0 s
t = 10680.0 s
t = 10690.0 s
t = 10700.0 s
t = 10710.0 s
t = 10720.0 s
t = 10730.0 s
t = 10740.0 s
t = 10750.0 s
t = 10760.0 s
t = 10770.0 s
t = 10780.0 s
t = 10790.0 s
t = 10800.0 s
t = 10810.0 s
t = 10820.0 s
t = 10830.0 s
t = 10840.0 s
t = 10850.0 s
t = 10860.0 s
t = 10870.0 s
t = 10880.0 s
t = 10890.0 s
t = 10900.0 s
t = 10910.0 s
t = 10920.0 s
t = 10930.0 s
t = 10940.0 s
t = 10950.0 s
t = 10960.0 s
t = 10970.0 s
t = 10980.0 s
t = 10990.0 s
t = 11000.0 s
t = 11010.0 s
t = 11020.0 s
t = 11030.0 s
t = 11040.0 s
t = 11050.0 s
t = 11060.0 s
t = 11070.0 s
t = 11080.0 s
t = 11090.0 s
t = 11100.0 s
t = 11110.0 s
t = 11120.0 s
t = 11130.0 s
t = 11140.0 s
t = 11150.0 s
t = 11160.0 s
t = 11170.0 s
t = 11180.0 s
t = 11190.0 s
t = 11200.0 s
t = 11210.0 s
t = 11220.0 s
t = 11230.0 s
t = 11240.0 s
t = 11250.0 s
t = 11260.0 s
t = 11270.0 s
t = 11280.0 s
t = 11290.0 s
t = 11300.0 s
t = 11310.0 s
t = 11320.0 s
t = 11330.0 s
t = 11340.0 s
t = 11350.0 s
t = 11360.0 s
t = 11370.0 s
t = 11380.0 s
t = 11390.0 s
t = 11400.0 s
t = 11410.0 s
t = 11420.0 s
t = 11430.0 s
t = 11440.0 s
t = 11450.0 s
t = 11460.0 s
t = 11470.0 s
t = 11480.0 s
t = 11490.0 s
t = 11500.0 s
t = 11510.0 s
t = 11520.0 s
t = 11530.0 s
t = 11540.0 s
t = 11550.0 s
t = 11560.0 s
t = 11570.0 s
t = 11580.0 s
t = 11590.0 s
t = 11600.0 s
t = 11610.0 s
t = 11620.0 s
t = 11630.0 s
t = 11640.0 s
t = 11650.0 s
t = 11660.0 s
t = 11670.0 s
t = 11680.0 s
t = 11690.0 s
t = 11700.0 s
t = 11710.0 s
t = 11720.0 s
t = 11730.0 s
t = 11740.0 s
t = 11750.0 s
t = 11760.0 s
t = 11770.0 s
t = 11780.0 s
t = 11790.0 s
t = 11800.0 s
t = 11810.0 s
t = 11820.0 s
t = 11830.0 s
t = 11840.0 s
t = 11850.0 s
t = 11860.0 s
t = 11870.0 s
t = 11880.0 s
t = 11890.0 s
t = 11900.0 s
t = 11910.0 s
t = 11920.0 s
t = 11930.0 s
t = 11940.0 s
t = 11950.0 s
t = 11960.0 s
t = 11970.0 s
t = 11980.0 s
t = 11990.0 s
t = 12000.0 s
t = 12010.0 s
t = 12020.0 s
t = 12030.0 s
t = 12040.0 s
t = 12050.0 s
t = 12060.0 s
t = 12070.0 s
t = 12080.0 s
t = 12090.0 s
t = 12100.0 s
t = 12110.0 s
t = 12120.0 s
t = 12130.0 s
t = 12140.0 s
t = 12150.0 s
t = 12160.0 s
t = 12170.0 s
t = 12180.0 s
t = 12190.0 s
t = 12200.0 s
t = 12210.0 s
t = 12220.0 s
t = 12230.0 s
t = 12240.0 s
t = 12250.0 s
t = 12260.0 s
t = 12270.0 s
t = 12280.0 s
t = 12290.0 s
t = 12300.0 s
t = 12310.0 s
t = 12320.0 s
t = 12330.0 s
t = 12340.0 s
t = 12350.0 s
t = 12360.0 s
t = 12370.0 s
t = 12380.0 s
t = 12390.0 s
t = 12400.0 s
t = 12410.0 s
t = 12420.0 s
t = 12430.0 s
t = 12440.0 s
t = 12450.0 s
t = 12460.0 s
t = 12470.0 s
t = 12480.0 s
t = 12490.0 s
Visualize the results to make sure we've averaged correctly.
In [8]:
fig = plt.figure(figsize=(20,15))
plt.subplots_adjust(right=0.87)
cax = fig.add_axes([0.88, 0.12, 0.025, 0.75])
plt.subplots_adjust(hspace=0.1)
for i in range(time_averaged_ems['250'].temperature_bin_edges.shape[0]-1):
# apply a filter to the
tmp = time_averaged_ems['250'][i].submap(u.Quantity([250,500],u.arcsec),u.Quantity([150,400],u.arcsec))
tmp.data = gaussian_filter(tmp.data,
eis.channels[0]['gaussian_width']['x'].value
)
# set up axes properly and add plot
ax = fig.add_subplot(6,5,i+1,projection=tmp)
im = tmp.plot(axes=ax,
annotate=False,
cmap=matplotlib.cm.get_cmap('magma'),
norm=matplotlib.colors.SymLogNorm(1, vmin=1e25, vmax=1e29)
)
# set title and labels
ax.set_title(r'${t0:.2f}-{t1:.2f}$ {uni}'.format(t0=np.log10(tmp.meta['temp_a']),
t1=np.log10(tmp.meta['temp_b']),uni='K'))
if i<25:
ax.coords[0].set_ticklabel_visible(False)
else:
ax.set_xlabel(r'$x$ ({})'.format(u.Unit(tmp.meta['cunit1'])))
if i%5==0:
ax.set_ylabel(r'$y$ ({})'.format(u.Unit(tmp.meta['cunit2'])))
else:
ax.coords[1].set_ticklabel_visible(False)
cbar = fig.colorbar(im,cax=cax)
Now save the results to our local temporary data folder.
In [12]:
for key in time_averaged_ems:
time_averaged_ems[key].save('../data/em_cubes_true_tn{}_t7500-12500.h5'.format(key))
In [10]:
foo = EMCube.restore('../data/em_cubes_tn250_t7500-12500.h5')
In [11]:
fig = plt.figure(figsize=(20,15))
plt.subplots_adjust(right=0.87)
cax = fig.add_axes([0.88, 0.12, 0.025, 0.75])
plt.subplots_adjust(hspace=0.1)
for i in range(foo.temperature_bin_edges.shape[0]-1):
# apply a filter to the
tmp = foo[i].submap(u.Quantity([250,500],u.arcsec),u.Quantity([150,400],u.arcsec))
tmp.data = gaussian_filter(tmp.data,
eis.channels[0]['gaussian_width']['x'].value
)
# set up axes properly and add plot
ax = fig.add_subplot(6,5,i+1,projection=tmp)
im = tmp.plot(axes=ax,
annotate=False,
cmap=matplotlib.cm.get_cmap('magma'),
norm=matplotlib.colors.SymLogNorm(1, vmin=1e25, vmax=1e29)
)
# set title and labels
ax.set_title(r'${t0:.2f}-{t1:.2f}$ {uni}'.format(t0=np.log10(tmp.meta['temp_a']),
t1=np.log10(tmp.meta['temp_b']),uni='K'))
if i<25:
ax.coords[0].set_ticklabel_visible(False)
else:
ax.set_xlabel(r'$x$ ({})'.format(u.Unit(tmp.meta['cunit1'])))
if i%5==0:
ax.set_ylabel(r'$y$ ({})'.format(u.Unit(tmp.meta['cunit2'])))
else:
ax.coords[1].set_ticklabel_visible(False)
cbar = fig.colorbar(im,cax=cax)
In [ ]:
Content source: wtbarnes/loops-workshop-2017-talk
Similar notebooks: