#!/home/hlrb2/pr95he/lu79qey3/yt-1.6/bin/python2.6

# This script shows the issue with the colorbar
# in phase plots.
# 
#
# Luigi Iapichino, February 2010.

# importing modules
import yt.mods, code
import math,os
from yt.mods import *
#from yt.lagos import hop
#import yt.lagos
import numpy as na
#import yt.raven as raven
from yt.raven import *


data_dir="/ptmp2/pr95he/lu79qey3/lss-sgs/64cubed/sgs1/files/outdir/"
# paths for ENZO data

fn = "RedshiftOutput0014"   # parameter file to load

##############################################################################

red = []

pf = load(data_dir+fn)

red = pf["CosmologyCurrentRedshift"]

print " \n "
print "Analyzing data:"
print "Redshift = "+str(red)
print " \n "

dd = pf.h.all_data() # This is an object that describes the entire box
pc = raven.PlotCollection(pf, center=[0.5,0.5,0.5])

pc.add_phase_object(dd, ["Density", "Temperature", "CellMassMsun"],
    cmap="Paired",weight=None,lazy_reader=True)

pc.save("test-colorbar",override=True)