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")
# # 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"])
# # 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")
# # 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"])
def __init__(self,config=None): super(FourierModeEstimator, self).__init__() self._base_config = DottedConfiguration.create(module=__name__,defaultcfg="fmts.yml",cfg=config) self.initialze()