#!/bin/bash

home=/sas/scratch/hsc/
basedir=$home/rerun/tutorial/deepCoadd/

patch_x=$1
patch_y=$2
patch="$patch_x,$patch_y"
#echo $basedir $patch

g=$basedir/HSC-G/0/$patch/calexp-HSC-G-0-$patch.fits
r=$basedir/HSC-R/0/$patch/calexp-HSC-R-0-$patch.fits
i=$basedir/HSC-I/0/$patch/calexp-HSC-I-0-$patch.fits
z=$basedir/HSC-Z/0/$patch/calexp-HSC-Z-0-$patch.fits
#echo $g

out="patch"$patch_x$patch_y"_grz.png"

ds9 -xpa no \
    -frame delete all \
    \
    -frame "G" $g \
    -scale limits '-0.1' 50 -scale log \
    -catalog load perseus_redshifts.vot \
    \
    -frame "R" $r \
    -scale limits '-0.1' 60 -scale log \
    \
    -rgb \
    -rgb red $z* \
    -scale limits '-0.05' 35 -scale log \
    -rgb green $r* \
    -scale limits '-0.05' 25 -scale log \
    -rgb blue $g* \
    -scale limits '-0.1' 25 -scale log \
    -rgb lock colorbar yes \
    -rgb lock scale yes \
    -catalog load perseus_redshifts.vot \
    -catalog symbol color green \
    -catalog symbol 1 shape "circle" \
    -catalog close \
    \
    -frame G-R \
    $home/analysis/colormaps/g-r/color_gr-$patch_x$patch_y.fits* \
    -scale limits -0.1 1.3 \
    -scale linear \
    \
    -frame G-I \
    $home/analysis/colormaps/g-i/color_gi-$patch_x$patch_y.fits* \
    -scale limits -0.1 2.0 \
    -scale linear \
    \
    -frame G-Z \
    $home/analysis/colormaps/g-z/color_gz-$patch_x$patch_y.fits* \
    -scale limits -0.1 3.0 \
    -scale linear \
    \
    -frame R-Z \
    $home/analysis/colormaps/r-z/color_rz-$patch_x$patch_y.fits* \
    -scale limits -0.1 3.0 \
    -scale linear \
    \
    -frame STRUCT2 \
    $home/analysis/structuremaps/calexp-HSC-G-0-$patch_x,$patch_y.filtered.structuremap_002.0.fits* \
    -scale limits 0.9 1.1 \
    -scale linear \
    \
    -lock frame image \
    -frame "G"