MBH1s = 10**array([float(i) for i in WM[:,10]])
MBH2s = 10**array([float(i) for i in WM[:,12]])
from rhomodels import NukerModelGenRho
for galaxy in range(15,16):
print 'galaxy ',galaxy+1, ' of ',len(WM)
name = names[galaxy]
alpha = alphas[galaxy]
beta = betas[galaxy]
gamma = gammas[galaxy]
M2L = M2Ls[galaxy]
MBH_Msun = MBH1s[galaxy]
rb = rbs[galaxy]
mub = mubs[galaxy]
rho0 = findrho0(rb,M2L,mub)
model1 = NukerModelGenRho('{0}'.format(name),alpha,beta,gamma,rb,rho0,MBH_Msun,GENERATE,memo = False)
model1.getrho()
Mencgood,psigood,Jc2good,ggood,Ggood,fgood = getrate(model1)
utest = arange(-7,1,0.05)
utest = 10**utest
d = dgdlnrp(utest,model1)
pklwrite('{0}/rrate.pkl'.format(model1.directory),utest)
pklwrite('{0}/rate.pkl'.format(model1.directory),d)
#fromfileplot(model1.name,model1.directory)
#fromfileplot(model2.name,model2.directory)
alpha = 1.0
beta = 4.0
gamma = 1.5
r0pc = 10**2.38#1.0
M2L = 6.27
mub = 19.98
rho0 = findrho0(r0pc,M2L,mub)#1e5
MBH_Msun = 10**6.04#1e3
name = 'NGC4467_test'
Menc,psi,Jc2,g,G,f,rate = 0,1,2,3,4,5,6
from rhomodels import NukerModelRho
model = NukerModelRho('{0}1'.format(name),alpha,beta,gamma,r0pc,rho0,MBH_Msun,GENERATE,memo = False)
name1 = str(model).split(' ')[0][11:]
Mencgood,psigood,Jc2good,ggood,Ggood,fgood,rategood = getrate(model)#,partial = {Menc:"OFF",psi:"OFF",Jc2:"OFF",g:"OFF",G:"OFF",f:"OFF",rate:"ON"})
from rhomodels import NukerModelGenRho
model1 = NukerModelGenRho('{0}2'.format(name),alpha,beta,gamma,r0pc,rho0,MBH_Msun,GENERATE,memo = False)
name2 = str(model1).split(' ')[0][11:]
model1.getrho()
Mencgood1,psigood1,Jc2good1,ggood1,Ggood1,fgood1,rategood1 = getrate(model1)#,partial = {Menc:"OFF",psi:"OFF",Jc2:"OFF",g:"OFF",G:"OFF",f:"OFF",rate:"ON"})
rtest2 = arange(-20,20,0.01)
rtest2 = 10**rtest2
plt.figure()
plt.loglog(rtest2,model.rho(rtest2))
plt.loglog(rtest2,model1.rho(rtest2))
plt.ylabel(r'$\rho$')
plt.xlabel(r'$r$')
plt.title('{0} vs {1}'.format(name1,name2))
plt.savefig('NukerRhoGals/compimages/newpsi_{0}_{1}_rho.png'.format(name1,name2))
