Flip order iK - it in output of S(q,t)

src/EXECS/Heis_Fourier_to_Sqt.cc

@@ -92,16 +92,18 @@ int main(int argc, char* argv[])
     if (FTim_outfile.fail()) ABACUSerror("Could not open FTim_outfile... ");
 
 
-    for (int iK = iKmin; iK <= iKmax; ++iK) {
+    for (int it = 0; it <= Npts_t; ++it) {
       FTre_outfile << setprecision(16);
       FTim_outfile << setprecision(16);
-      if (iK > iKmin) {
+
+      if (it > 0) {
 	FTre_outfile << endl;
 	FTim_outfile << endl;
       }
-      FTre_outfile << real(FT[0][iK]);
-      FTim_outfile << real(FT[0][iK]);
-      for (int it = 1; it <= Npts_t; ++it) {
+      FTre_outfile << real(FT[it][iKmin]);
+      FTim_outfile << real(FT[it][iKmin]);
+
+      for (int iK = iKmin + 1; iK <= iKmax; ++iK) {
 	FTre_outfile << "\t" << real(FT[it][iK]);
 	FTim_outfile << "\t" << imag(FT[it][iK]);
       }