Tutorial

We need some scripts from BHAc code


this python scripts are here(we need to put his in our ./bashrc)
export PYTHONPATH="$BHAC_DIR/tools/python:"

ipython --pylab

import read,amrplot
import matplotlib.pyplot as plt
import numpy as np

d1=read.load(0,file='output/data',type='vtu')

d1.rho

p1=amrplot.polyplot(d1.rho,d1,blocks=None,blockWidth = 10,blockHeight = 10,min=-6,max=0,cmap='jet',showcbar=True,clear=True)

p1=amrplot.polyplot(np.log10(d1.p/d1.rho),d1)

p1.ax.set_ybound([-50,50])

p1.ax.set_xbound([0,50])


p1.ax.set_title(r'$\log_{10}\rho$')
p1.save('logrho.png')
export LC_ALL=en_US.UTF-8
export LANG=en_US.UTF-8