Science/cloudy-agn
1
0
Fork 0
mirror of https://asciireactor.com/otho/cloudy-agn.git synced 2025-10-31 13:48:05 +00:00
Code Issues Packages Projects Releases Wiki Activity

added note about SED cutoff for future program

Browse Source
This commit is contained in:
caes 2017-08-28 00:10:02 -04:00
parent 7e7879e969
commit f48507c8e8
4 changed files with 9 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
notes/SED_cutoff Normal file
View File

@ -0,0 +1 @@
fun = A nu^(hnu/kTcut)

4
scripts/cloudy/magdziarz.mpi_grid.in
View File

@ -6,7 +6,7 @@ title magdziarz 1997 hden_vs_phi .25 dex
//
//save last incident continuum "incident_continuum"
//save line labels "line_labels"
save continuum bins "continuum_bins"
//save continuum bins "continuum_bins"
//save diffuse continuum "diffuse_continuum"
save emitted continuum "emitted_continuum"
//
@ -371,4 +371,4 @@ grid range from 17 to 24 with .25 dex
////// Stopping condition and run
//
stop total column density 23
iterate to convergence max=40
iterate to convergence max=15

4
scripts/cloudy/mehdipour.mpi_grid.in
View File

@ -6,7 +6,7 @@ title mehdipour ngc5548 4thdex grid
//
//save last incident continuum "incident_continuum"
//save line labels "line_labels"
save continuum bins "continuum_bins"
//save continuum bins "continuum_bins"
//save diffuse continuum "diffuse_continuum"
save emitted continuum "emitted_continuum"
//
@ -371,4 +371,4 @@ grid range from 17 to 24 with .25 dex
////// Stopping condition and run
//
stop total column density 23
iterate to convergence max=40
iterate to convergence max=15

4
sed/plots.R
View File

@ -189,6 +189,10 @@ cut -f1,3 magdziarz_continuum_bins_angstroms > bins
#R code
cont = read.table("cont",header=TRUE)
bins = read.table("bins",header=TRUE)
plot(cont,log="xy",type="b",xlim=c(1000,1200),ylim=c(1e7,1e9));abline(v=1420)
# old stuff
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)