#!/bin/bash
# A simple top-view of the inner Solar system
# (c) SsODNet, B.Carry & J.Berthier, Lagrange/OCA & LTE/Obs. Paris

if [[ ! -f ssoBFT-latest_Asteroid.parquet ]]; then
  read -p "ssoBFT not found. Download it [~800MB]? [y/n]: " resp
  if [[ ${resp,,} = "y" ]]; then
    echo "Downloading ssoBFT ..."
    wget -q https://ssp.imcce.fr/data/ssoBFT-latest_Asteroid.parquet
  else
    exit 1
  fi
fi

topcat -stilts plot2plane \
   xpix=644 ypix=568 \
   xlabel='x\ /\ au' ylabel='y\ /\ au' aspect=1 texttype=latex \
    fontsize=16 fontstyle=serif \
   xmin=-6.75 xmax=6.75 ymin=-6 ymax=6 \
   legend=true legpos=0.03,0.962 \
   auxvisible=false \
   in=ssoBFT-latest_Asteroid.parquet \
    x='position[0]' y='position[1]' shading=translucent \
   layer_1=Mark \
      icmd_1='select "contains(class,\"MB\")"' \
      color_1=light_grey \
      leglabel_1='Main\ belt' \
   layer_2=Mark \
      icmd_2='select "contains(class,\"Hilda\")"' \
      color_2=black \
      leglabel_2=Hilda \
   layer_3=Mark \
      icmd_3='select "contains(class,\"L4\")"' \
      color_3=green \
      leglabel_3=Trojan/L4 \
   layer_4=Mark \
      icmd_4='select "contains(class,\"L5\")"' \
      color_4=blue \
      leglabel_4=Trojan/L5 \
   layer_5=Mark \
      icmd_5='select "contains(class,\"NEA\")"' \
      leglabel_5=NEA