add bbl fix to atm interp as well
@@ -256,6 +256,7 @@ def atm_interpolation(
output_metadata['band names'] = band_names
output_metadata['description'] = 'Interpolated atmospheric state'
output_metadata['bands'] = len(band_names)
+ output_metadata['bbl'] = '{' + ','.join([str(1) for n in range(len(band_names))]) + '}'
del output_metadata['fwhm']
del output_metadata['wavelength']
