Homec4science
Diffusion pNbody (master) Browse

pNbody/pNbody/SSP/libRGBs.py
737d7b683637master

Last Change
Raw File

libRGBs.py
View Options

#!/usr/bin/env python
import scipy
import scipy.interpolate
from numpy import *
from pylab import *
import sys,os
from pNbody import *
#Should be in its own class...
Kroupa_meanstar = 0.4820
def Kroupa_IMF(m):
#Meanstar = 0.4820
A = 0.10629
if m <= 0.08:
return(A*m**-0.3)
elif m <= 0.5:
return(A*m**-1.8)
elif m <= 1:
return(A*m**-2.7)
else:
return(A*m**-2.3)
def Kroupa_IMFint(m):
A = 0.10629
f = zeros(size(m))
f[where(m<=0.08)] = A*m**0.7/0.7 - 0.01864976527083534
f[where((m > 0.08)*(m<=0.5))] = -A*m**-.8/.8 + 1.0021446918927897 + 0.0072656
f[where((m>0.5)*(m<=1))] = -A*m**-1.7/1.7 + 0.2031395463734218+0.778083245211
f[where(m>1)] = -A*m**-1.3/1.3 + A*1**-1.3/1.3+0.918699
return(f)
class NRGB:
def __init__(self,file):
self.fileMS = file+'MS.dat'
self.fileRGB = file+'rgb.dat'
self.ReadMS() #Read data file, create self.data
self.ReadRGB()
self.InterpolateMS()
self.InterpolateTRGB()
def ReadMS(self):
pass
def ReadRGB(self):
pass
def InterpolateMS(self):
pass
def InterpolateTRGB(self):
pass
Kroupa_meanstar = 0.4820
def Kroupa_IMF(m):
#Meanstar = 0.4820
A = 0.10629
if m <= 0.08:
return(A*m**-0.3)
elif m <= 0.5:
return(A*m**-1.8)
elif m <= 1:
return(A*m**-2.7)
else:
return(A*m**-2.3)
def Kroupa_IMFint(m):
A = 0.10629
f = zeros(size(m))
f[where(m<=0.08)] = A*m**0.7/0.7 - 0.01864976527083534
f[where((m > 0.08)*(m<=0.5))] = -A*m**-.8/.8 + 1.0021446918927897 + 0.0072656
f[where((m>0.5)*(m<=1))] = -A*m**-1.7/1.7 + 0.2031395463734218+0.778083245211
f[where(m>1)] = -A*m**-1.3/1.3 + A*1**-1.3/1.3+0.918699
return(f)
def RGBs(self,Zs,Ages):
Frac = Kroupa_IMFint(self.InterpolateTRGB.ev(Zs,Ages)) - Kroupa_IMFint(self.InterpolateMS.ev(Zs,Ages))
return(Frac/Kroupa_meanstar)
c4science · Help