#!/usr/bin/env python3
# -*- coding: utf-8 -*-
# Composite 'snow_age' does not work in SatPy 0.37.1 (it used to work in SatPy 0.36.0)

# We certainly need this
from satpy.scene import Scene
# Why to hell is it not working?
from satpy.utils import debug_on
debug_on()

# Available DAY composites: See files ../satpy/etc/composites/visir.yaml and viirs.yaml
# *************************************************************************************
# ['ash', 'cloudtop_daytime', 'dust', 'false_color', 'fire_temperature', 'fire_temperature_39refl',
#  'fire_temperature_awips', 'fire_temperature_eumetsat', 'fog', 'ir108_3d', 'ir_cloud_day',
#  'natural_color', 'natural_color_sun_lowres', 'natural_with_night_fog', 'night_fog',
#  'ocean_color', 'overview', 'snow_age', 'snow_lowres', 'true_color', 'true_color_crefl',
#  'true_color_lowres', 'true_color_lowres_crefl', 'true_color_lowres_land',
#  'true_color_lowres_marine_tropical', 'true_color_raw']

composite = 'snow_age'

# These two segments should show Italy left and Spain with Gibraltar right (upside down)
files = ['DATA/SVMC_npp_d20221105_t1224383_e1226025_b57121_c20221105123534000154_eum_ops.h5',
         'DATA/SVMC_npp_d20221105_t1226037_e1227279_b57121_c20221105123541000704_eum_ops.h5']

# Swath image size is 3200 X (768) * NumberOfSegments (MX-Bands)
scn = Scene(filenames = files, reader = 'viirs_compact')
scn.load([composite])
scn.save_dataset(composite, 'IMAGES/SNPP-'+composite+'.png')