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)