barium.py.svn-base
No OneTemporary
Actions

Subscribers

None

File Metadata

	#!/usr/bin/env python

	# data from Cescutti 2006, tab 3

	import Ptools as pt
	from numpy import *
	import scipy.interpolate
	from pNbody import myNumeric

	ms = array([10,12. ,15. ,30.,])
	Ba = array([0,9e-7 ,3e-8 ,1e-9])
	Eu = array([0,4.5e-8,3e-9 ,5e-10])


	pt.plot(ms,Ba,'o')
	#pt.plot(ms,Eu)

	# interpolation
	msi = array([10.,13.,15.,18.,20.,25.,40.,70.])

	tck = scipy.interpolate.fitpack.splrep(ms,Ba,k=1,s=0)
	Bai = scipy.interpolate.fitpack.splev(msi,tck)
	#Bai = myNumeric.lininterp1d(msi.astype(float32),ms.astype(float32),Ba.astype(float32))

	pt.plot(msi,Bai,':')

	pt.show()