Paste #405

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
from yt.mods import *
import yt.extensions.volume_rendering as vr

pf = load("DD1701/DD1701")

L = [0.1, 0.3, 0.5]
mi, ma = -32.0, -25.0

W = 1.0
c = [0.5,0.5,0.5]

tf = vr.ColorTransferFunction((mi, ma), nbins=1024)
tf.add_layers(10, colormap='bds_highcontrast')

vp = pf.h.volume_rendering(L, W, c, 1024, tf)
vp.ray_cast()

image = (vp.image - vp.image.min()) / (vp.image.max() - vp.image.min())
vr.plot_rgb(image, 'PROC%03i' % MPI.COMM_WORLD.rank)