Raster created by python GDAL interpreted upside down

I am using CRS which has axis in left bottom corner. Code below produce Geotiff which is shown allright in QGIS but when I open it in other software it is upside down. Also gdal_calc takes it as upside down when calculating. Is there some easy way how to make sure this grid will be interpreted right as it is if I create it using gdal_translate from xyz?

geotransform=(min_x,min_cell_size,0,min_y,0,min_cell_size)
nrows,ncols = np.shape(zi)
output_raster = gdal.GetDriverByName('GTiff').Create(output_reaster,ncols, nrows, 1 ,gdal.GDT_Float32)
output_raster.SetGeoTransform(geotransform)
srs = osr.SpatialReference()
srs.ImportFromEPSG(28355)
output_raster.SetProjection(srs.ExportToWkt())
output_raster.GetRasterBand(1).WriteArray(zi)