9 plutodir = os.environ[
'PLUTO_DIR']
10 wdir = plutodir+
'/Test_Problems/HD/Stellar_Wind/'
11 nlinf = pp.nlast_info(w_dir=wdir,datatype=
'float')
13 D = pp.pload(nlinf[
'nlast'],w_dir=wdir,datatype=
'float')
16 I.pldisplay(D, log10(D.rho),x1=D.x1,x2=D.x2,label1=
'x',label2=
'y',
17 title=
r'Log Density $\rho$ [Stellar Wind]',cbar=(
True,
'vertical'),figsize=[8,12])
22 Xmesh, Ymesh = meshgrid(D.x1.T,D.x2.T)
23 xcong = T.congrid(Xmesh,newdims,method=
'linear')
24 ycong = T.congrid(Ymesh,newdims,method=
'linear')
25 velxcong = T.congrid(D.vx1.T,newdims,method=
'linear')
26 velycong = T.congrid(D.vx2.T,newdims,method=
'linear')
27 gca().quiver(xcong, ycong, velxcong, velycong,color=
'w')
29 savefig(
'stellar_wind_1.png')