Changes in parameters

fourier_transform.c

@@ -104,7 +104,7 @@ void write_fftw_data(){
                 for (int iv=0; iv<yaxis; iv++)
                 for (int iu=0; iu<xaxis; iu++)
                 {
-                          long global_index = (iu + (iv+isector*yaxis)*xaxis + iw*grid_size_x*grid_size_y)*sizeof(double);
+                          long global_index = (iu + (iv+isector*yaxis)*xaxis + iw*param.grid_size_x*param.grid_size_y)*sizeof(double);
                           long index = iu + iv*xaxis + iw*xaxis*yaxis;
                           fseek(file.pFilereal, global_index, SEEK_SET);
                           fwrite(&gridss_real[index], 1, sizeof(double), file.pFilereal);
@@ -113,7 +113,7 @@ void write_fftw_data(){
                 for (int iv=0; iv<yaxis; iv++)
                 for (int iu=0; iu<xaxis; iu++)
                 {
-                          long global_index = (iu + (iv+isector*yaxis)*xaxis + iw*grid_size_x*grid_size_y)*sizeof(double);
+                          long global_index = (iu + (iv+isector*yaxis)*xaxis + iw*param.grid_size_x*param.grid_size_y)*sizeof(double);
                           long index = iu + iv*xaxis + iw*xaxis*yaxis;
                           fseek(file.pFileimg, global_index, SEEK_SET);
                           fwrite(&gridss_img[index], 1, sizeof(double), file.pFileimg);

gridding.c

@@ -333,22 +333,22 @@ void write_grided_data()
                       gridss_real[i] = gridss[2*i];
                       gridss_img[i] = gridss[2*i+1];
               }
-              if (num_w_planes > 1)
+              if (param.num_w_planes > 1)
               {
-                      for (int iw=0; iw<num_w_planes; iw++)
+                      for (int iw=0; iw<param.num_w_planes; iw++)
                         for (int iv=0; iv<yaxis; iv++)
                           for (int iu=0; iu<xaxis; iu++)
                           {
-                               long global_index = (iu + (iv+isector*yaxis)*xaxis + iw*grid_size_x*grid_size_y)*sizeof(double);
+                               long global_index = (iu + (iv+isector*yaxis)*xaxis + iw*param.grid_size_x*param.grid_size_y)*sizeof(double);
                                long index = iu + iv*xaxis + iw*xaxis*yaxis;
                                fseek(file.pFilereal, global_index, SEEK_SET);
                                fwrite(&gridss_real[index], 1, sizeof(double), file.pFilereal);
                           }
-                      for (int iw=0; iw<num_w_planes; iw++)
+                      for (int iw=0; iw<param.num_w_planes; iw++)
                         for (int iv=0; iv<yaxis; iv++)
                           for (int iu=0; iu<xaxis; iu++)
                           {
-                               long global_index = (iu + (iv+isector*yaxis)*xaxis + iw*grid_size_x*grid_size_y)*sizeof(double);
+                               long global_index = (iu + (iv+isector*yaxis)*xaxis + iw*param.grid_size_x*param.grid_size_y)*sizeof(double);
                                long index = iu + iv*xaxis + iw*xaxis*yaxis;
                                fseek(file.pFileimg, global_index, SEEK_SET);
                                fwrite(&gridss_img[index], 1, sizeof(double), file.pFileimg);
@@ -359,9 +359,9 @@ void write_grided_data()
               }
               else
               {
-                      for (int iw=0; iw<num_w_planes; iw++)
+                      for (int iw=0; iw<param.num_w_planes; iw++)
                       {
-                          long global_index = (xaxis*isector*yaxis + iw*grid_size_x*grid_size_y)*sizeof(double);
+                          long global_index = (xaxis*isector*yaxis + iw*param.grid_size_x*param.grid_size_y)*sizeof(double);
                           long index = iw*xaxis*yaxis;
                           fseek(file.pFilereal, global_index, SEEK_SET);
                           fwrite(&gridss_real[index], xaxis*yaxis, sizeof(double), file.pFilereal);
@@ -371,11 +371,11 @@ void write_grided_data()
               }
           }
        #else
-          for (int iw=0; iw<num_w_planes; iw++)
-             for (int iv=0; iv<grid_size_y; iv++)
-               for (int iu=0; iu<grid_size_x; iu++)
+          for (int iw=0; iw<param.num_w_planes; iw++)
+             for (int iv=0; iv<param.grid_size_y; iv++)
+               for (int iu=0; iu<param.grid_size_x; iu++)
                 {
-                      long index = 2*(iu + iv*grid_size_x + iw*grid_size_x*grid_size_y);
+                      long index = 2*(iu + iv*param.grid_size_x + iw*param.grid_size_x*param.grid_size_y);
                       fwrite(&gridtot[index], 1, sizeof(double), file.pFilereal);
                       fwrite(&gridtot[index+1], 1, sizeof(double), file.pFileimg);
                       //double v_norm = sqrt(gridtot[index]*gridtot[index]+gridtot[index+1]*gridtot[index+1]);
@@ -389,6 +389,7 @@ void write_grided_data()
      #ifdef USE_MPI
         MPI_Win_fence(0,slabwin);
      #endif
+
    #endif //WRITE_DATA 
 
 }