From f48507c8e8b081482fc08cd5d0e5e7c678f9cdb5 Mon Sep 17 00:00:00 2001 From: caes Date: Mon, 28 Aug 2017 00:10:02 -0400 Subject: [PATCH] added note about SED cutoff for future program --- notes/SED_cutoff | 1 + scripts/cloudy/magdziarz.mpi_grid.in | 4 ++-- scripts/cloudy/mehdipour.mpi_grid.in | 4 ++-- sed/plots.R | 4 ++++ 4 files changed, 9 insertions(+), 4 deletions(-) create mode 100644 notes/SED_cutoff diff --git a/notes/SED_cutoff b/notes/SED_cutoff new file mode 100644 index 0000000..1841973 --- /dev/null +++ b/notes/SED_cutoff @@ -0,0 +1 @@ +fun = A nu^(hnu/kTcut) diff --git a/scripts/cloudy/magdziarz.mpi_grid.in b/scripts/cloudy/magdziarz.mpi_grid.in index 0ccae80..3334c8d 100644 --- a/scripts/cloudy/magdziarz.mpi_grid.in +++ b/scripts/cloudy/magdziarz.mpi_grid.in @@ -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 \ No newline at end of file +iterate to convergence max=15 \ No newline at end of file diff --git a/scripts/cloudy/mehdipour.mpi_grid.in b/scripts/cloudy/mehdipour.mpi_grid.in index 2254381..5df3023 100644 --- a/scripts/cloudy/mehdipour.mpi_grid.in +++ b/scripts/cloudy/mehdipour.mpi_grid.in @@ -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 \ No newline at end of file diff --git a/sed/plots.R b/sed/plots.R index 5e03996..b66405f 100644 --- a/sed/plots.R +++ b/sed/plots.R @@ -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)