I have produced new EROD grids from soil surveys via GEOTIFF. I have used gdal_translate, but since erod as three bands, I run into issues. I can either get one band in and the file is okay (but still unacceptable to geogrid, maybe because it lacks a Z dimension?), or using -b 1 -b 2 -b 3, I get a file that is not anything close to what it should be, usually checkered (I used all the banding options [BSQ/BIP/BIL]). What is the best way to write new binary files for geogrid that have a z dimension?