real    array(i0dum:ildum, j0dum:jldum, k0dum:kldum)
real, dimension(:,:,:), allocatable :: arrayNewOrder
allocate( arrayNewOrder(k0dum:kldum, j0dum:jldum, i0dum:ildum) )

arrayNewOrder= reshape(array, shape(arrayNewOrder), order=[1,3,2])

...
do some shit
...

array= reshape(arrayNewOrder, shape(array), order=[1,3,2,4])