

In [1]:

    
%pylab nbagg
import ppmpy.ppm as ppm
import os

data_dir = '/data/ppm_rpod2/YProfiles/O-shell-M25/'









    



Populating the interactive namespace from numpy and matplotlib



In [2]:

    
ppm.plot_mesa_time_evo('/data/ppm_rpod2/Stellar_models/O-shell-M25/M25Z0.02/LOGS',
                       range(550,560),7.5829245098141646E+006,ifig=21)









    














    











    



Using old star.logsa file ...
 reading ...100% 

reading /data/ppm_rpod2/Stellar_models/O-shell-M25/M25Z0.02/LOGS/log550.data ...
 reading ...100% 

reading /data/ppm_rpod2/Stellar_models/O-shell-M25/M25Z0.02/LOGS/log551.data ...
 reading ...100% 

Closing profile tool ...
reading /data/ppm_rpod2/Stellar_models/O-shell-M25/M25Z0.02/LOGS/log552.data ...
 reading ...100% 

Closing profile tool ...
reading /data/ppm_rpod2/Stellar_models/O-shell-M25/M25Z0.02/LOGS/log553.data ...
 reading ...100% 

Closing profile tool ...
reading /data/ppm_rpod2/Stellar_models/O-shell-M25/M25Z0.02/LOGS/log554.data ...
 reading ...100% 

Closing profile tool ...
reading /data/ppm_rpod2/Stellar_models/O-shell-M25/M25Z0.02/LOGS/log555.data ...
 reading ...100% 

Closing profile tool ...
reading /data/ppm_rpod2/Stellar_models/O-shell-M25/M25Z0.02/LOGS/log556.data ...
 reading ...100% 

Closing profile tool ...
reading /data/ppm_rpod2/Stellar_models/O-shell-M25/M25Z0.02/LOGS/log557.data ...
 reading ...100% 

Closing profile tool ...
reading /data/ppm_rpod2/Stellar_models/O-shell-M25/M25Z0.02/LOGS/log558.data ...
 reading ...100% 

Closing profile tool ...
reading /data/ppm_rpod2/Stellar_models/O-shell-M25/M25Z0.02/LOGS/log559.data ...
 reading ...100% 

Closing profile tool ...
Closing star.log  tool ...
Closing profile tool ...



In [3]:

    
ppm.set_YProf_path(data_dir)
ppm.v_evolution(['D1','D2','D5','D6','D8','D9','D10'],4.,8.,'radial','max', sparse = 10, ifig = 12)

Time Evolution