+main3d :: IO ()
+main3d = do
+ (s:_) <- getArgs
+ let scale = read s :: Int
+ let zoom_factor = (scale, scale, scale)
+ let out_file = "output.bin"
+ arr <- read_word16s in_file
+ let arr' = swap_bytes arr
+ let arrMRI = R.reshape mri_shape arr'
+ let dbl_data = R.computeS $ R.map fromIntegral arrMRI
+ let output = zoom dbl_data zoom_factor
+ let word16_output = R.computeS $ round_array output
+ write_word16s out_file word16_output