7 plutodir = os.environ[
'PLUTO_DIR']
8 wdir = plutodir+
'/Test_Problems/HD/Viscosity/Flow_Past_Cylinder/'
11 D = pp.pload(10,datatype=
'hdf5',level=3, w_dir=wdir)
14 print 'data.0010.hdf5 not found .. Loading data.0000.hdf5 at level 0'
15 D = pp.pload(0,datatype=
'hdf5',w_dir=wdir)
18 I.pldisplay(D, D.rho, x1=D.x1, x2=D.x2,
19 polar=[
True,
True], fignum=1,
20 cbar=[
True,
'horizontal'],figsize=[8,8],
21 label1=
'X', label2=
'Y',title=
r'Density [$\rho$]')
23 ax1.axis([-10, 50,-20.0, 20.0])
26 I.oplotbox(D.AMRLevel,lrange=[0,3],geom=D.geometry)
28 plt.savefig(
'amr_flowcyc.png')