# mkdir -p 2mass-cut
# python cut_2mass.py
ln -s ../4000/the-making-of/2mass-cut .


mkdir -p 2mass-split
hpsplit -o 2mass-split/2mass-hp%02i.fits -n 2 -m 1 2mass-cut/2mass-hp* > hpsplit.log 2>hpsplit.err &

# mkdir -p 2mass-cut-J15
# for ((x=0; x<972; x++)); do
#   echo $x
#   echo fitscopy $(printf 2mass-cut/2mass-hp%03i.fits $x)"[j_mag < 15]" $(printf 2mass-cut-J15/2mass-hp%03i.fits $x)
#   fitscopy $(printf 2mass-cut/2mass-hp%03i.fits $x)"[j_mag < 15]" $(printf 2mass-cut-J15/2mass-hp%03i.fits $x)
# done
# 
# cd 2mass-cut-J15
# cp 2mass-hp000.fits 2mass-all.fits
# for ((x=1; x<972; x++)); do
#   echo $x
#   tabmerge $(printf 2mass-hp%03i.fits+1 $x) 2mass-all.fits+1
# done
ln -s ../4000/the-making-of/2mass-cut-J15 .


mkdir -p 2mass-split-1
hpsplit -o 2mass-split-1/2mass-hp%02i.fits -n 1 -m 10 2mass-cut/2mass-hp* > hpsplit1.log 2>hpsplit1.err &

python index.py

## I had a bug where the index files were still named index-40* and had headers INDEXID = 40*;
## convert  40 -> 42.

for x in index-4000-*.fits; do modhead $x INDEXID 4200; done
for x in index-4001-*.fits; do modhead $x INDEXID 4201; done
for x in index-4002-*.fits; do modhead $x INDEXID 4202; done
for x in index-4003-*.fits; do modhead $x INDEXID 4203; done
for x in index-4004-*.fits; do modhead $x INDEXID 4204; done
for x in index-4005-*.fits; do modhead $x INDEXID 4205; done
for x in index-4006-*.fits; do modhead $x INDEXID 4206; done
for x in index-4007-*.fits; do modhead $x INDEXID 4207; done

modhead index-4008.fits INDEXID 4208
modhead index-4009.fits INDEXID 4209
for ((i=10; i<20; i++)); do modhead index-40${i}.fits INDEXID 42${i}; done

for x in index*.fits; do mv -v $x $(echo $x | sed s/40/42/); done


# Tests
# solve-field --backend-config cfg2 -D tests --continue \
# 			-L 0.391 -H 0.401 -u app --odds-to-tune-up 1e100 --odds-to-reject 1e-9 \
# 			-d 1-30 -w 2048 -e 1361 \
# 			--fields 1-10000 -s mag -a \
# 			--no-remove-lines --no-fits2fits --no-tweak --no-plots \
# 			-N none -U none -R none -B none -W none \
# 			~/binary/sdss-tests-astrometry.net/sdss-e-r-00.fits > e-r-00.log
# 
# 
# 
# $ printsolved -u -p tests/sdss-e-r-00.solved 
# File tests/sdss-e-r-00.solved
# tests/sdss-e-r-00.solved: 39/10000 (0.390000 %) unsolved
# 92 251 285 806 1040 1402 1661 2310 2659 2676 2835 3298 3470 4010 4033 4171 4261 5241 5400 5563 5618 5703 5882 5915 6410 6606 6704 6849 7246 7302 7385 7394 7439 7461 7470 8514 8929 9257 9732 
# Total: 39/10000 (0.390000 %) unsolved
# 
# for ((i=1; i<=18; i++)); do
# I=$(printf %02i $i);
# solve-field --backend-config cfg2 -D tests --continue \
# 		   	-L 0.391 -H 0.401 -u app --odds-to-tune-up 1e100 --odds-to-reject 1e-9 \
# 			-d 1-30 -w 2048 -e 1361 \
# 			--fields 1-10000 -s mag -a \
# 			--no-remove-lines --no-fits2fits --no-tweak --no-plots \
# 			-N none -U none -R none -B none -W none \
# 			~/binary/sdss-tests-astrometry.net/sdss-e-r-${I}.fits > e-r-${I}.log &
# done
# wait
# 
# 
# 
# $ printsolved -u -p tests/*.solved -s
# Total: 584/182221 (0.320490 %) unsolved
# 
# 
# for ((i=0; i<=18; i++)); do
# I=$(printf %02i $i);
# solve-field --backend-config cfg2 -D tests2 --continue \
# 		   	-L 0.391 -H 0.401 -u app --odds-to-tune-up 1e100 --odds-to-reject 1e-9 \
# 			-d 1-100 -w 2048 -e 1361 -s mag -a \
# 			--fields $(printsolved -j -u tests/sdss-e-r-${I}.solved | tr " " ,) \
# 			--no-remove-lines --no-fits2fits --no-tweak --no-plots \
# 			-N none -U none -R none -B none -W none \
# 			~/binary/sdss-tests-astrometry.net/sdss-e-r-${I}.fits > e-r-${I}-stage2.log &
# done
# 
# $ printsolved -j tests2/*.solved | wc -w
# 303
# 
# 584 - 303 = 281
#