#!/bin/bash

# Script to compute VIIRS cloud cluster statistics within CrIS FOVS
# and add the cluster data to the CrIS BUFR file
# 
# VIIRS SDR files are read from $PWD/viirs_sdr
# CRIS SDR files are read from $PWD/cris_sdr

# Set up CSPP SDR (to use nagg)
export CSPP_SDR_HOME=$HOME/dbvm/apps/sdr/sdr
export PATH=$PATH:$CSPP_SDR_HOME/common/local/bin

# Create agggregated VIIRS SDR file containing SVM05, SVM07, SVM11, SVM12, SVM15, GMTCO
nagg --onefile -n 99 -t "SVM05,SVM07,SVM12,SVM15,SVM16" -g GMTCO viirs_sdr/SVM{05,07,12,15,16}*h5 viirs_sdr/GMTCO*h5
viirs_agg=$(find . -name GMTCO-SVM05\*h5)

# Create aggregated CrIS SDR file containing SCRIF and GCRSO
nagg --onefile -n 99 -t SCRIF -g GCRSO cris_sdr/SCRIF*h5 cris_sdr/GCRSO*h5
cris_agg=$(find . -name GCRSO-SCRIF\*h5)

# Set up AAPP (containerized version)
export AAPP_INSTALL_DIR=$HOME/dbvm/aapp_build/singularity/v2/v8.8/build
export PATH=$AAPP_INSTALL_DIR/scripts:$PATH

# Create CrIS Level 1C
aapp_run.bash cris_sdr -o cris.l1c -g $cris_agg $cris_agg

# Create VIIRS clusters
aapp_run.bash viirs_to_cris -D -C viirs_clusters.h5 -g $viirs_agg cris.l1c $viirs_agg 

# Create CrIS BUFR (cris.bufr)
aapp_run.bash eccodes_encodebufr_1c -i cris.l1c CRISFSR

# Add VIIRS clusters to CrIS BUFR (cris_clusters.bufr)
aapp_run.bash cris_bufr_add_clusters cris.bufr viirs_clusters.h5 cris_clusters.bufr

# Clean up
rm GMTCO-SVM05*h5 GCRSO-SCRIF*h5 cris.l1c viirs_clusters.h5 viirs_to_cris-M.log cluster*H5 FORMAP*H5 viirs_location.H5