#!/usr/bin/env python

#=============================================================================
# Imports
#=============================================================================
import os
import os.path as osp
import xarray as xr
import numpy as np

#=============================================================================
# Compute sums
#=============================================================================

# Open the bpch dataset
ds = xr.open_dataset('modis_emissivity_tmp.4x5.nc')
print(ds)

for i in range(0,16):
    ds_temp = ds.isel(lev=i)
    ds_temp.to_netcdf('modis_emissivity.band{}.4x5.nc'.format(str(i+1).zfill(2)))