Diffusion Analysis


In [1]:
%pylab ipympl
from ppmpy import ppm
import matplotlib.pyplot as pl
import os

project_dir = '/data/ASDR/PPMstar/rpod2_PPMstar/YProfiles/O-shell-M25/'
data_dir = '/data/ASDR/CSA/'
ppm.set_YProf_path(project_dir)


Populating the interactive namespace from numpy and matplotlib

In [4]:
%%capture  --no-stdout --no-display
D1 = ppm.yprofile('D1')
res = D1.Dinv(231, 331, plot_Dlt0=False, approx_D = True, linelims = [7.7,8.0],
        r0= 6,debug =False)
pl.xlim(7.6,8.3)
pl.ylim(7,15)


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
Out[4]:
(7, 15)

In [7]:
%%capture  --no-stdout --no-display
ppm.plot_diffusion_profiles('D2',data_dir + 'Stellar_models/O-shell-M25/M25Z0.02/',28900,
                        7.8489,(1,160),2.*460.,7.8282,10.**12.27,0.21,0.055,1.6)
pl.xlim(3,9)
pl.ylim(11,16)
pl.legend(loc=6)


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
1386 in profiles.index file ...
Found and load nearest profile for cycle 28900
reading profile/data/ASDR/CSA/Stellar_models/O-shell-M25/M25Z0.02//LOGS/log593.data ...
 reading ...100% 

Closing profile tool ...
Out[7]:
<matplotlib.legend.Legend at 0x7f9b24419780>

In [ ]: