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))