#!/bin/sh

# Galaxy center from DSS opt visual: 20:56:19.07 +31:26:33.9
# http://www.astrouw.edu.pl/~jskowron/ra-dec
ract=314.0794583
decct=31.4927500

# T120 FOV is 13.2'x13.2' = 0.22 by 0.22 deg2 (2048x2048 at ~0.38"/px B2 = *2 scale = 0.77"/px) 
# -> take a ~10% overlap on each side -> offset = 0.19deg
# 0.19 deg. physical offset in E-W -> 0.20/cos(31) = 0.22 deg. RA offset
raoffset=0.22
decoffset=0.19

# We want a 1.78 aspect ratio and an N-S extension of 1deg => 5 vert pointings -> / 1.8 = 4.5 => 4 in RA: iterate on loop and check the area is well covered

# Total = 30 pointings

# Field number from bottom left to top right going left to right per DEC slice
index=0
# DEC loop
i=-3
while [ $i -le 2 ]; do

	dec=`echo $decct $decoffset $i | awk '{printf "%f", $1 + ( $2 * $3) }'`

	# RA loop
	j=-2
	while [ $j -le 2 ]; do

		ra=`echo $ract $raoffset $j | awk '{printf "%f", $1 + ( $2 * $3) }'`

		index=`echo $index | awk '{printf "%02d", $1 + 1 }'`
		# echo "NGC6992_Field_${index}	$ra	$dec	2000.0"


                # For string format
                ralt=`echo $ra | awk '{printf "%d", $1 }' | wc -c | awk '{print $1 }'`

                # Convert J2000d -> J2000
                # RA
                hr=`echo $ra           | awk '{printf "%02d", ( $1 * 24 ) / 360 }'`     # 24hr = 360dg
                reshr=`echo $ra $hr    | awk '{printf "%.6f", ( ($1 - ($2*360/24))) }'`
                mn=`echo $reshr        | awk '{printf "%02d", $1 * 4 }'`                # 1mn  = 0.25dg (1/4)
                resmn=`echo $reshr $mn | awk '{printf "%.6f", $1 - $2/4 }'`
                sc=`echo $resmn        | awk '{printf "%02.2f", $1 * 240  }'`           # 1sec = 0.004166dg (1/240)
                length=`echo $sc       | wc | awk '{printf "%2d", $3 - 1 }'`
                if [ $length = 4 ]; then
                        sc="0${sc}"
                fi
                # DEC
                dg=`echo $dec            | awk '{printf "%02d", $1 }'`                  # Same degree unit
                am=`echo $dec $dg        | awk '{printf "%02d", ( $1 - $2 ) * 60 }'`    # 1deg = 60' (am)
                resam=`echo $dec $dg $am | awk '{printf "%f", ( $1 - $2 - $3/60 )}'`
                ac=`echo $resam          | awk '{printf "%02d", $1 * 3600 }'`           # 1deg = 3600'' (as)

		echo "NGC6992_${index}	${hr}:${mn}:${sc}	+${dg}:${am}:${ac}"

		j=`echo $j | awk '{print $1 + 1}'`
	done

	i=`echo $i | awk '{print $1 + 1}'`
done
# echo " "