diff --git a/flash2d/src/options-sphere/ADV_wind_stacked.f90 b/flash2d/src/options-sphere/ADV_wind_stacked.f90
index af2458c2e67112f8a1700c2438bcd0f6d0919e6a..8d7efcf3670e0cc5c130abf40b51676e8818c769 100644
--- a/flash2d/src/options-sphere/ADV_wind_stacked.f90
+++ b/flash2d/src/options-sphere/ADV_wind_stacked.f90
@@ -204,14 +204,14 @@
       TYPE (control_struct), INTENT(in)      :: p_control
       
 !---------- local variables
-      CHARACTER (len=io_fillen)              :: c_uname, c_vname, c_wname
-      INTEGER (KIND = GRID_SI)               :: i_alct, i_ncstat
-      INTEGER (KIND = GRID_SI)               :: i_filexid, i_fileyid, i_filezid
-      INTEGER (KIND = GRID_SI)               :: i_dimid, i_varid
-      INTEGER (KIND = GRID_SI)               :: i_tm, i_ln, i_lt, i_ly
-      REAL (KIND = GRID_SR), DIMENSION(:,:,:,:), ALLOCATABLE :: r_auxx
-      REAL (KIND = GRID_SR), DIMENSION(:,:,:,:), ALLOCATABLE :: r_auxy
-      REAL (KIND = GRID_SR), DIMENSION(:,:,:,:), ALLOCATABLE :: r_auxz
+      CHARACTER (len=io_fillen)               :: c_uname, c_vname, c_wname
+      INTEGER (KIND = GRID_SI)                :: i_alct, i_ncstat
+      INTEGER (KIND = GRID_SI)                :: i_filexid, i_fileyid, i_filezid
+      INTEGER (KIND = GRID_SI)                :: i_dimid, i_varid
+      INTEGER (KIND = GRID_SI)                :: i_tm, i_ln, i_lt, i_ly
+      REAL*4, DIMENSION(:,:,:,:), ALLOCATABLE :: r_auxx
+      REAL*4, DIMENSION(:,:,:,:), ALLOCATABLE :: r_auxy
+      REAL*4, DIMENSION(:,:,:,:), ALLOCATABLE :: r_auxz
 
 !---------- initialize file names
       WRITE(GRID_parameters%ioout,*) &
@@ -572,11 +572,11 @@
 
       IF(r_dy /= 0.0_GRID_SR) THEN
         r_dyi= 1._GRID_SR/r_dy
-        r_inter(1)= (r_hy* r_l1+ r_ly* r_h1)* r_dyi* r_scalx
-        r_inter(2)= (r_hy* r_l2+ r_ly* r_h2)* r_dyi* r_scaly
+        r_inter(1)= (r_hy* r_l1+ r_ly* r_h1)* r_dyi* r_scaly
+        r_inter(2)= (r_hy* r_l2+ r_ly* r_h2)* r_dyi* r_scalx
       ELSE
-        r_inter(1)= r_l1* r_scalx
-        r_inter(2)= r_l2* r_scaly
+        r_inter(1)= r_l1* r_scaly
+        r_inter(2)= r_l2* r_scalx
       END IF
 
 !---------- calculate wind vector in kartesian coordinates