From 6c1951c3fd25bf975b7a63edc59d140ca534010e Mon Sep 17 00:00:00 2001 From: caes Date: Thu, 10 Aug 2017 05:39:07 -0400 Subject: [PATCH] R plot scripts --- sed/plots.R | 76 +++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 76 insertions(+) create mode 100644 sed/plots.R diff --git a/sed/plots.R b/sed/plots.R new file mode 100644 index 0000000..8562276 --- /dev/null +++ b/sed/plots.R @@ -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)