Replace gdal_calc.py with BandMath