Example #1
0
def main ():
    config = DottedConfiguration.fromfile ("../input/config.yaml")
    
    if len (sys.argv) < 2:
        print "Proper use is './merge_cdf.py number_of_cdfs_to_merge'"
        exit ()
    
    merge ("../output/" + config ["output.file"] + ".cdf", range (int (sys.argv [1])), range (int (config ["time.steps"] / config ["output.every"]) + 1), "../output/" + config ["output.merge_file"] + ".cdf")
Example #2
0
#
#  Created by Alexander Rudy on 2012-11-12.
#  Copyright 2012 Alexander Rudy. All rights reserved.
#

import numpy as np
import matplotlib.pyplot as plt

from pystellar.energy import epp, eCNO

from pystellar.opacity import OpacityTable
from pystellar.threading import ObjectThread

from pyshell.config import DottedConfiguration

Config = DottedConfiguration()
Config.dn = DottedConfiguration
Config.load("Star.yml")
X = 1.000
Y = 0.000
rho = 1
Opacity = ObjectThread(OpacityTable,
                       ikwargs=dict(fkey='GN93hz', X=X, Y=Y, snap=True),
                       locking=True,
                       timeout=None)
Opacity.start()

T = np.logspace(6.5, 7.5, 100)

epsilon_pp = epp(T=T, rho=rho, X=X, c=Config["Data.Energy"])
epsilon_CNO = eCNO(T=T, rho=rho, X=X, XCNO=0.01, c=Config["Data.Energy"])
Example #3
0
#  
#  Created by Alexander Rudy on 2012-11-12.
#  Copyright 2012 Alexander Rudy. All rights reserved.
# 

import numpy as np
import matplotlib.pyplot as plt

from pystellar.energy import epp, eCNO

from pystellar.opacity import OpacityTable
from pystellar.threading import ObjectThread

from pyshell.config import DottedConfiguration

Config = DottedConfiguration()
Config.dn = DottedConfiguration
Config.load("Star.yml")
X = 1.000
Y = 0.000
rho = 1
Opacity = ObjectThread(OpacityTable,ikwargs=dict(fkey='GN93hz',X=X,Y=Y,snap=True),locking=True,timeout=None)
Opacity.start()

T = np.logspace(6.5,7.5,100)

epsilon_pp = epp(T=T,rho=rho,X=X,c=Config["Data.Energy"])
epsilon_CNO = eCNO(T=T,rho=rho,X=X,XCNO=0.01,c=Config["Data.Energy"])

plt.loglog(T,epsilon_CNO+epsilon_pp,'-',label='Total')
plt.loglog(T,epsilon_pp,':',label="PP")
Example #4
0
#  
#  Created by Alexander Rudy on 2012-11-12.
#  Copyright 2012 Alexander Rudy. All rights reserved.
# 

import numpy as np
import matplotlib.pyplot as plt

from pystellar.energy import epp, eCNO

from pystellar.opacity import OpacityTable
from pystellar.threading import ObjectThread

from pyshell.config import DottedConfiguration

Config = DottedConfiguration()
Config.dn = DottedConfiguration
Config.load("pystellar/Star.yml")
Config.load("Star.yml")
X = 0.700
Y = 0.280
XCNO = 0.7 * (1 - X - Y)
rho = 80

T = np.logspace(6.8,7.8,100)

epsilon_pp = epp(T=T,rho=rho,X=X,c=Config["Data.Energy"])
epsilon_CNO = eCNO(T=T,rho=rho,X=X,XCNO=XCNO,c=Config["Data.Energy"])

T1 = 18 * np.power(10,6)
epsilon_pp1 = epp(T=T1,rho=rho,X=X,c=Config["Data.Energy"])
Example #5
0
 def __init__(self,config=None):
     super(FourierModeEstimator, self).__init__()
     self._base_config = DottedConfiguration.create(module=__name__,defaultcfg="fmts.yml",cfg=config)
     self.initialze()