Plotting Radial Buckets



In [1]:

    
%pylab nbagg
import ppmpy.ppm as ppm
import ppmpy.rprofile as rprofile
import os









    



Populating the interactive namespace from numpy and matplotlib



In [2]:

    
data_dir = '/data/ppm_rpod2/RProfiles/'
!ls /data/ppm_rpod2/RProfiles/









    



AGBTP_M2.0Z1.e-5  C-ingestion  O-shell-M25  O-shell-mixing  sakurai



In [3]:

    
project = 'AGBTP_M2.0Z1.e-5'
!ls /data/ppm_rpod2/RProfiles/AGBTP_M2.0Z1.e-5/









    



F3  F4	F4e  F4f  F5  F9



In [4]:

    
data_path = os.path.join(data_dir,project,'F4')
print(data_path)
rp_set = rprofile.rprofile_set(data_path)
rp = rp_set.get_dump(100)
rp.get('radtop')









    



/data/ppm_rpod2/RProfiles/AGBTP_M2.0Z1.e-5/F4






    Out[4]:





29.009326934814453



In [5]:

    
rp.get('radbase')









    Out[5]:





11.18260383605957



In [6]:

    
ppm.plot_boundary_evolution(data_path, 27., 31., 700,sparse=20)



In [7]:

    
r_lim = (27, 30.5)
ut_lim = (1e-3, 19.)
dutdr_lim = (-0.099, 0.02)
dump = 560
hist_dmin = dump - 5
hist_dmax = dump + 5
ppm.upper_bound_ut(data_path,dump, hist_dmin, hist_dmax,r1 = r_lim[0],r2 = 31, derivative = False,
              ylims = [1e-3,19.])



In [8]:

    
ppm.upper_bound_ut(data_path,dump, hist_dmin, hist_dmax,r1 = r_lim[0],r2 = 31, derivative = True,
              ylims = (-0.099, 0.02))



In [9]:

    
ppm.plot_boundary_evolution(data_path,None,None,r_int=True,r_ref=13.,gamma = 5./3.,insert = True,
                       lims = [0., 2e3,27.6, 28.4],sparse=20)









    



/home/user/PyPPM/ppmpy/ppm.py:5983: RuntimeWarning: invalid value encountered in true_divide
  A = p/rho**gamma
/home/user/PyPPM/ppmpy/ppm.py:5992: RuntimeWarning: invalid value encountered in greater
  idx_top = np.argmax(rel_diff > threshold)






    














    











    



/usr/local/lib/python3.5/dist-packages/matplotlib/figure.py:1743: UserWarning: This figure includes Axes that are not compatible with tight_layout, so its results might be incorrect.
  warnings.warn("This figure includes Axes that are not "



In [ ]: