Reading in custom-written raster bands causes seg fault

download.sh

wget https://www.unidata.ucar.edu/software/netcdf/examples/tos_O1_2001-2002.nc

test.py

import rasterio as rio

with rio.drivers():
    with rio.open('tos_O1_2001-2002.nc', 'r') as src:
        raster = src.read()[0]
        profile = src.profile
        
    profile['count'] = 2
    
    with rio.open('test.nc', 'w', **profile) as dst:
        dst.write_band(1, raster)
        dst.write_band(2, raster)
        
    with rio.open('test.nc') as src:
        rasters = src.read()

Two quick comments:

• you can read just the first band using src.read(1) rather than reading the entire dataset into a 3D array
• you can write bands using dst.write(raster, 1) which is preferred over the write_band method.