#!/bin/bash
#
# wrapper for union_bigwig
#

DRYRUN=${DRYRUN:-"1"};

for A in CA1 CA2 CA3 DG; do
   for B in CB DE; do
      for STRAND in pos neg; do
         OUTBW="${A}_${B}.union.${STRAND}.bw";
         INBW=$(ls ../NS*${A}${B}*STAR_mm10*${STRAND}.bw);
         if [ "1" == "${DRYRUN}" ]; then
            echo -e "Skipping due to DRYRUN=1:${BOLD}GENOME=mm10 ./union_bigwig_bwtools.sh ${OUTBW} ${INBW[*]}${RESET}";
         else
            if [ -f "${OUTBW}" ]; then
               echo -e "${FYELLOW}Output BW already exists:${BOLD}${OUTBW}${RESET}";
            else
               #DRYRUN=${DRYRUN} GENOME=mm10 ./union_bigwig_bwtools.sh CA1_CB.neg.bw ../NS50211CA1CB492_1.sickle.merged.cutadapt.STAR_mm10.neg.bw ../NS50217CA1CB502_1.sickle.merged.cutadapt.STAR_mm10.neg.bw ../NS50218CA1CB496_1.sickle.merged.cutadapt.STAR_mm10.neg.bw
#            echo -e "${FCYAN}DRYRUN=${DRYRUN} GENOME=mm10 ./union_bigwig_bwtools.sh ${OUTBW} ${INBW[*]}${RESET}";
               time DRYRUN=${DRYRUN} GENOME=mm10 ./union_bigwig_bwtools.sh ${OUTBW} ${INBW[*]};
            fi;
         fi;
      done;
   done;
done;