diff --git a/Makefile b/Makefile
index f55a1efec7688d807fde89681b38ddb90c7b51db..2012f489f280bf9e0ab358066303059fd20eb393 100644
--- a/Makefile
+++ b/Makefile
@@ -2,7 +2,7 @@
 
 include ${FSLCONFDIR}/default.mk
 
-OPTFLAGS = -O3
+OPTFLAGS = -O3 
 
 PROJNAME = melodic2
 
diff --git a/melreport.cc b/melreport.cc
index b49e489194eebba891b0f5f7cd5a6cadaa2f31da..ba4e32e5b8616f8f9a86b4951b87318c716f0613 100644
--- a/melreport.cc
+++ b/melreport.cc
@@ -133,17 +133,19 @@ namespace Melodic{
       
       {//plot frequency 
 	miscplot newplot;
+	RowVector empty(1);
+	empty = 0.0;
 	int fact = int(std::pow(10.0,int(std::log10(float(melodat.get_mix().Nrows())))));
       
 	if(opts.tr.value()>0.0)
-	  newplot.timeseries(melodat.get_fmix().Column(cnum).t(),
+	  newplot.timeseries(empty | melodat.get_fmix().Column(cnum).t(),
 			     report.appendDir(string("f")+
 					      num2str(cnum)+".png"),
 			     string("FFT of timecourse (in Hz / ") +num2str(fact)+")",
 			     fact/(opts.tr.value()*melodat.get_mix().Nrows()), 150,
 			     0,2);
 	else
-	  newplot.timeseries(melodat.get_fmix().Column(cnum).t(),
+	  newplot.timeseries(empty | melodat.get_fmix().Column(cnum).t(),
 			     report.appendDir(string("f")+
 					      num2str(cnum)+".png"),
 			     string(string("FFT of timecourse (in cycles); ")
@@ -154,7 +156,7 @@ namespace Melodic{
 				  +"* TR)/cycles"));
 	write_ascii_matrix(report.appendDir(string("f")
 					    +num2str(cnum)+".txt"),
-			   melodat.get_mix().Column(cnum));
+			   melodat.get_fmix().Column(cnum));
 	IChtml << "<A HREF=\"" << string("f")
 	  +num2str(cnum)+".txt" << "\"> ";
 	IChtml << "<img BORDER=0 SRC=\""