In [1]:

    

import yt


    

In [18]:

    

ds = yt.load("gen-mesh-one-material_out.e", step=-1)

ad = ds.all_data()

temp = ad[('connect1', 'temp')]
group1 = ad['group1']
group2 = ad['group2']

print(temp.max())
print(temp.min())
print(group1.max())
print(group1.min())
print(group2.max())
print(group2.min())

g1slice = yt.SlicePlot(ds, 'z', 'group1')
g2slice = yt.SlicePlot(ds, 'z', 'group2')


    

    




yt : [INFO     ] 2016-12-22 10:33:41,113 Loading coordinates
yt : [INFO     ] 2016-12-22 10:33:41,114 Loading connectivity
yt : [INFO     ] 2016-12-22 10:33:41,124 Parameters: current_time              = 2.0
yt : [INFO     ] 2016-12-22 10:33:41,125 Parameters: domain_dimensions         = [1 1 1]
yt : [INFO     ] 2016-12-22 10:33:41,126 Parameters: domain_left_edge          = [-0.6 -0.6  0. ]
yt : [INFO     ] 2016-12-22 10:33:41,128 Parameters: domain_right_edge         = [ 6.6  6.6  1. ]
yt : [INFO     ] 2016-12-22 10:33:41,129 Parameters: cosmological_simulation   = 0
yt : [INFO     ] 2016-12-22 10:33:41,137 Loading coordinates
yt : [INFO     ] 2016-12-22 10:33:41,138 Loading connectivity
yt : [INFO     ] 2016-12-22 10:33:41,269 xlim = -0.600000 6.600000
yt : [INFO     ] 2016-12-22 10:33:41,269 ylim = -0.600000 6.600000
yt : [INFO     ] 2016-12-22 10:33:41,271 xlim = -0.600000 6.600000
yt : [INFO     ] 2016-12-22 10:33:41,272 ylim = -0.600000 6.600000
yt : [INFO     ] 2016-12-22 10:33:41,275 Making a fixed resolution buffer of (('connect1', 'group1')) 800 by 800






    




959.3496444727219 dimensionless
900.0 dimensionless
8.138120750714538 dimensionless
8.137864325565221 dimensionless
0.8231674445177946 dimensionless
0.8230987236017343 dimensionless






    




yt : [INFO     ] 2016-12-22 10:33:41,526 xlim = -0.600000 6.600000
yt : [INFO     ] 2016-12-22 10:33:41,527 ylim = -0.600000 6.600000
yt : [INFO     ] 2016-12-22 10:33:41,528 xlim = -0.600000 6.600000
yt : [INFO     ] 2016-12-22 10:33:41,529 ylim = -0.600000 6.600000
yt : [INFO     ] 2016-12-22 10:33:41,531 Making a fixed resolution buffer of (('connect1', 'group2')) 800 by 800

In [16]:

    

g1slice.set_log('group1', False)
g1slice.set_zlim('group1', group1.min(), group1.max())


    

    
Out[16]:

In [17]:

    

g2slice.set_log('group2', False)
g2slice.set_zlim('group2', group2.min(), group2.max())


    

    
Out[17]:

In [19]:

    

ds = yt.load("../gen-mesh-one-material_out.e", step=-1)

ad = ds.all_data()

temp = ad[('connect1', 'temp')]
group1 = ad['group1']
group2 = ad['group2']

print(temp.max())
print(temp.min())
print(group1.max())
print(group1.min())
print(group2.max())
print(group2.min())

g1slice = yt.SlicePlot(ds, 'z', 'group1')
g2slice = yt.SlicePlot(ds, 'z', 'group2')


    

    




yt : [INFO     ] 2016-12-22 10:35:23,543 Loading coordinates
yt : [INFO     ] 2016-12-22 10:35:23,544 Loading connectivity
yt : [INFO     ] 2016-12-22 10:35:23,551 Parameters: current_time              = 2.0
yt : [INFO     ] 2016-12-22 10:35:23,552 Parameters: domain_dimensions         = [1 1 1]
yt : [INFO     ] 2016-12-22 10:35:23,554 Parameters: domain_left_edge          = [-0.6 -0.6  0. ]
yt : [INFO     ] 2016-12-22 10:35:23,555 Parameters: domain_right_edge         = [ 6.6  6.6  1. ]
yt : [INFO     ] 2016-12-22 10:35:23,557 Parameters: cosmological_simulation   = 0
yt : [INFO     ] 2016-12-22 10:35:23,558 Loading coordinates
yt : [INFO     ] 2016-12-22 10:35:23,560 Loading connectivity
yt : [INFO     ] 2016-12-22 10:35:23,697 xlim = -0.600000 6.600000
yt : [INFO     ] 2016-12-22 10:35:23,698 ylim = -0.600000 6.600000
yt : [INFO     ] 2016-12-22 10:35:23,699 xlim = -0.600000 6.600000
yt : [INFO     ] 2016-12-22 10:35:23,700 ylim = -0.600000 6.600000
yt : [INFO     ] 2016-12-22 10:35:23,702 Making a fixed resolution buffer of (('connect1', 'group1')) 800 by 800






    




900.4045415306368 dimensionless
900.0 dimensionless
8.097011032290823 dimensionless
8.097006592146112 dimensionless
0.8138626175116652 dimensionless
0.8138613746801414 dimensionless






    




yt : [INFO     ] 2016-12-22 10:35:23,956 xlim = -0.600000 6.600000
yt : [INFO     ] 2016-12-22 10:35:23,957 ylim = -0.600000 6.600000
yt : [INFO     ] 2016-12-22 10:35:23,958 xlim = -0.600000 6.600000
yt : [INFO     ] 2016-12-22 10:35:23,959 ylim = -0.600000 6.600000
yt : [INFO     ] 2016-12-22 10:35:23,961 Making a fixed resolution buffer of (('connect1', 'group2')) 800 by 800

In [20]:

    

g1slice.set_log('group1', False)
g1slice.set_zlim('group1', group1.min(), group1.max())


    

    
Out[20]:

In [21]:

    

g2slice.set_log('group2', False)
g2slice.set_zlim('group2', group2.min(), group2.max())


    

    
Out[21]:

In [ ]: