real conversion in wind interpolation

flash2d/src/options-sphere/ADV_wind_stacked.f90

@@ -428,12 +428,9 @@
       WRITE(GRID_parameters%ioout,*) &
         '***** reshaping core data arrays ...                  *****'
 
-!       r_flowx = reshape(r_auxx, shape(r_flowx), order=(/4,3,2,1/))
-!       r_flowy = reshape(r_auxy, shape(r_flowy), order=(/4,3,2,1/))
-!       r_flowz = reshape(r_auxz, shape(r_flowz), order=(/4,3,2,1/))
-      r_flowx = r_auxx
-      r_flowy = r_auxy
-      r_flowz = r_auxz
+      r_flowx = reshape(r_auxx, shape(r_flowx), order=(/4,3,2,1/))
+      r_flowy = reshape(r_auxy, shape(r_flowy), order=(/4,3,2,1/))
+      r_flowz = reshape(r_auxz, shape(r_flowz), order=(/4,3,2,1/))
       DEALLOCATE(r_auxx, r_auxy, r_auxz)
 
 !---------- Fix mask values