GDAL netCDF works in Anaconda, add GDAL commands to file