Entrained Material

In [1]:

    

%pylab ipympl
from ppmpy import ppm
import os

data_dir = '/data/ASDR/PPMstar/rpod2_PPMstar/YProfiles'
project = 'O-shell-M25'

mpl_logger = logging.getLogger('matplotlib')
mpl_logger.setLevel(logging.WARNING)


    

    




Populating the interactive namespace from numpy and matplotlib

In [2]:

    

ppm.set_YProf_path(os.path.join(data_dir,project))


    

In [3]:

    

D1 = ppm.yprofile('D1')
D2 = ppm.yprofile('D2')
D1.entrainment_rate(range(150, 332, 10), 7.5, 8.5, var = 'vxz',
                     criterion = 'min_grad', offset = -1.,
                     integrate_both_fluids = False, ifig0 = 18,
                     mdot_curve_label = r'D1 $\left( 768^3 \right)$')


    

    




Reading attributes from file  YProfile-01-0348.bobaaa
There are 349 YProfile files in the /data/ASDR/PPMstar/rpod2_PPMstar/YProfiles/O-shell-M25/D1 directory.
Ndump values range from 0 to 348
Time values range from 0.0 to 3480.0
Reading attributes from file  YProfile-01-0163.bobaaa
There are 164 YProfile files in the /data/ASDR/PPMstar/rpod2_PPMstar/YProfiles/O-shell-M25/D2 directory.
Ndump values range from 0 to 163
Time values range from 0.0 to 1640.0






    






 
 










    






 
 










    
Out[3]:





1.1567752336505737e-06

In [4]:

    

cases = ('D1', 'D8', 'D5', 'D6', 'D9', 'D10')
c0 = (241,154,142,155,123,78)  
ppm.plot_entr_v_ut(cases,c0,10,7.5,8.5,'tangential','max','$768^3$',ifig = 3,
                  integrate_both_fluids = False)
xlim(1.4,2.4);ylim(-6.5,-3.6)


    

    




Reading attributes from file  YProfile-01-0348.bobaaa
There are 349 YProfile files in the /data/ASDR/PPMstar/rpod2_PPMstar/YProfiles/O-shell-M25/D1 directory.
Ndump values range from 0 to 348
Time values range from 0.0 to 3480.0
Reading attributes from file  YProfile-01-0217.bobaaa
There are 218 YProfile files in the /data/ASDR/PPMstar/rpod2_PPMstar/YProfiles/O-shell-M25/D8 directory.
Ndump values range from 0 to 217
Time values range from 0.0 to 2170.0
Reading attributes from file  YProfile-01-0279.bobaaa
There are 280 YProfile files in the /data/ASDR/PPMstar/rpod2_PPMstar/YProfiles/O-shell-M25/D5 directory.
Ndump values range from 0 to 279
Time values range from 0.0 to 2790.0
Reading attributes from file  YProfile-01-0248.bobaaa
There are 249 YProfile files in the /data/ASDR/PPMstar/rpod2_PPMstar/YProfiles/O-shell-M25/D6 directory.
Ndump values range from 0 to 248
Time values range from 0.0 to 2480.0
Reading attributes from file  YProfile-01-0261.bobaaa
There are 262 YProfile files in the /data/ASDR/PPMstar/rpod2_PPMstar/YProfiles/O-shell-M25/D9 directory.
Ndump values range from 0 to 261
Time values range from 0.0 to 2610.0
Reading attributes from file  YProfile-01-0262.bobaaa
There are 263 YProfile files in the /data/ASDR/PPMstar/rpod2_PPMstar/YProfiles/O-shell-M25/D10 directory.
Ndump values range from 0 to 262
Time values range from 0.0 to 2620.0






    






 
 










    
Out[4]:





(-7, -3.5)

In [6]:

    

H1 = ppm.yprofile('/data/ppm_rpod2/YProfiles/agb-entrainment-convergence/H1')
ppm.compare_entrained_material([H1,D2],['O shell','AGB'], fname = 271)


    

    




Reading attributes from file  YProfile-01-0498.bobaaa
There are 499 YProfile files in the /data/ppm_rpod2/YProfiles/agb-entrainment-convergence/H1 directory.
Ndump values range from 0 to 498
Time values range from 0.0 to 29954.6






    






 
 

In [7]:

    

D1.prof_time([0,3,4,5,10,15,20],yaxis_thing = 'FV',
             logy=True,num_type='time',markevery = 50)


    

    






 
 










    




/usr/local/lib/python3.6/dist-packages/ppmpy/ppm.py:3113: RuntimeWarning: divide by zero encountered in log10
  pl.plot(Y,np.log10(y),utils.linestyle(i+ls_offset)[0],

In [ ]: