Write data

fourier_transform.c

@@ -274,12 +274,41 @@ void write_fftw_data(){
   fits_write_subset(fptrimg, TDOUBLE, fpixel, lpixel, image_imag, &status);
   fits_close_file(fptrimg, &status);
 
+ #endif //FITSIO
+
+  for (int isector=0; isector<size; isector++)
+    {
+
+      MPI_Barrier(MPI_COMM_WORLD);
+      
+      if(isector == rank)
+	{
+
+	  printf("%d writing\n",isector);
+
+	 #ifdef FITSIO
+
+	  fpixel[0] = 1;
+	  fpixel[1] = isector*yaxis+1;
+	  lpixel[0] = xaxis;
+	  lpixel[1] = (isector+1)*yaxis;
+
+	  status = 0;
+	  fits_open_image(&fptreal, testfitsreal, READWRITE, &status);
+	  fits_write_subset(fptreal, TDOUBLE, fpixel, lpixel, image_real, &status);
+	  fits_close_file(fptreal, &status);
+
+	  status = 0;
+	  fits_open_image(&fptrimg, testfitsimag, READWRITE, &status);
+	  fits_write_subset(fptrimg, TDOUBLE, fpixel, lpixel, image_imag, &status);
+	  fits_close_file(fptrimg, &status);
+
 	 #endif //FITSIO
 
 	  file.pFilereal = fopen (out.fftfile2,"ab");
 	  file.pFileimg = fopen (out.fftfile3,"ab");
 
-  uint global_index = rank*(xaxis*yaxis)*sizeof(double);
+	  uint global_index = isector*(xaxis*yaxis)*sizeof(double);
 
 	  fseek(file.pFilereal, global_index, SEEK_SET);
 	  fwrite(image_real, xaxis*yaxis, sizeof(double), file.pFilereal);
@@ -288,6 +317,8 @@ void write_fftw_data(){
 
 	  fclose(file.pFilereal);
 	  fclose(file.pFileimg);
+	}
+    }
  
   MPI_Barrier(MPI_COMM_WORLD);