#!/bin/sh
# superimpose Antarctic map over US to show relative sizes
scale='1:45000000'
out='ant_over_us_map.ps'
add_file_loc='/home/braup/lib/ADD_GMT'


# antarctica
region='-180/180/-90/-50'
proj="s0/-90/$scale"
iceshelf_colors='180'
icesheet_colors='150'
outline='0.05p/0/0/0'
pscoast -R$region -J$proj -Dl -G$iceshelf_colors -K -V -P > $out
# The following data are from the Antarctic Digital Database (ADD)
# Draw coastline and border
gzip -dc ${add_file_loc}/ADDant.txt.gz | psxy -J$proj -R$region \
          -M            \
          -G$icesheet_colors -K  -O -P  -V    >> $out

# Add ice shelves
gzip -dc ${add_file_loc}/ADDris.txt.gz |
  psxy -J$proj -R$region -G${iceshelf_colors} -M -K -O -P -V >> $out
gzip -dc ${add_file_loc}/ADDfris.txt.gz |
  psxy -J$proj -R$region -G${iceshelf_colors}  -M -K -O -P -V >> $out
gzip -dc ${add_file_loc}/ADDotheris.txt.gz |
  psxy -J$proj -R$region  -G${iceshelf_colors} -M -K -O -P -V >> $out
gzip -dc  ${add_file_loc}/ADDfrisisl.txt.gz |
  psxy -J$proj -R$region -G${icesheet_colors}  -M -K -O -P -V >> $out
gzip -dc ${add_file_loc}/ADDrisisl.txt.gz |
  psxy -J$proj -R$region -G${icesheet_colors} \
          -M -K -O -P >> $out

echo "-34 -54 10 0 13 RM South Georgia Island" | pstext -R$region -J$proj -G0 -N -V -K -O >> $out

# US
region='-130/-60/20/55.0'
pscoast -R$region -Jb-102/40/35/50/$scale -Dl -A5000 -W1p/0/0/0 -N1/1.5p/0/0/0 -N2/0.5p/0/0/0 -X3.5c -Y5.5c -P -V -O >> $out