Science/cloudy-agn
1
0
Fork 0
mirror of https://asciireactor.com/otho/cloudy-agn.git synced 2024-12-04 18:45:07 +00:00
Code Issues Packages Projects Releases Wiki Activity

R plot scripts

Browse Source
This commit is contained in:
caes 2017-08-10 05:39:07 -04:00
parent 82fb8eb922
commit 6c1951c3fd
1 changed files with 76 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

76
sed/plots.R Normal file
View File

@ -0,0 +1,76 @@
data = read.table("mehdipour2013.tab")
plot (data,log="xy",type="l")
data=read.table("magdziarz1997.tab")
plot (data,log="xy",type="l")
data = read.table("mehdipour2013.tab")
plot (data,log="xy",type="l",xlim=c(.001,300),ylim=c(10e11,10e13))
# mehdipour axes
abline(v=0.001)
abline(h=40e11)
data = read.table("mehdipour2013.tab")
plot (data,log="xy",type="l")
# mehdipour axes
abline(v=0.001)
abline(h=40e11)
data=read.table("magdziarz1997.tab")
plot (data,log="xy",type="l",xlim=c(.001,300),ylim=c(1e-3,1e-1))
#magdziarz boundaries
abline(v=.0004)
abline(v=.03)
abline(v=.4)
abline(v=3)
abline(v=50)
abline(v=500)
#low uv peak
abline(h=0.04)
abline(v=0.008)
#high uv peak
abline(h=0.024)
abline(v=0.012)
#xray peak
abline(h=0.03)
abline(v=60)
coords=read.table("magdziarz1997_samples.tab")
points(coords)
hc=9.113e2
inci = read.table("cont_combined_inci",header=TRUE)
diffuse = read.table("cont_combined_diffuse",header=TRUE)
# To plot inci with error bars for bin size
plot(Enr,nFn,log="x",xlim=c(20,25),type="l")
arrows(Enr-d.anu.,nFn,Enr+d.anu.,nFn,length=0.05,angle=90,code=3)
arrows(Enr-d.anu.,nFn+5,Enr+d.anu.,nFn-5,length=0.05,angle=90,code=3)
arrows(Enr-d.anu.,nFn+5,Enr+d.anu.,nFn+5,length=0.05,angle=90,code=3)
arrows(Enr-d.anu.,nFn+50,Enr+d.anu.,nFn+50,length=0.05,angle=90,code=3)
arrows(Enr-d.anu.,nFn+1e9,Enr+d.anu.,nFn+1e9,length=0.05,angle=90,code=3)
plot(energy.Ryd,Total,log="x",type="l",xlim=c(1e-2,5))
arrows(energy.Ryd-d.anu.,Total,energy.Ryd+d.anu.,Total,length=0.05,angle=90,code=3)
plot(energy.Ryd,ConEmitLocal,log="x",type="l",xlim=c(1e-2,5))
lambda=hc/energy.Ryd
width_lambda=lambda*d.anu./energy.Ryd
plot(lambda,Total,log="x",type="l",ylim=c(0,0.2),xlim=c(2000,9400))
arrows(lambda-width_lambda/2,Total,lambda+width_lambda/2,Total,length=0.05,angle=90,code=3)