#!/bin/sh

# Set common options
PROJECTION=E0.0/-90.0/20.0c
REGION=-180/180/-90/-50
FILENAME=Antarctic_Ocean_location_map.ps

# Get Natural Earth datafiles
echo "Downloading datafiles from naturalearthdata.com"
wget -nv -nc \
 http://www.naturalearthdata.com/http//www.naturalearthdata.com/download/50m/physical/50m-coastline.zip\
 http://www.naturalearthdata.com/http//www.naturalearthdata.com/download/50m/physical/50m-land.zip\
 http://www.naturalearthdata.com/http//www.naturalearthdata.com/download/50m/cultural/50m-admin-0-boundary-lines.zip\
 http://www.naturalearthdata.com/http//www.naturalearthdata.com/download/50m/physical/50m-antarctic-ice-shelves-polys.zip\
 http://www.naturalearthdata.com/http//www.naturalearthdata.com/download/50m/physical/50m-antarctic-ice-shelves-lines.zip
echo "...done\n"

echo "Inflating archives"
unzip 50m\*.zip
echo "...done\n"

echo "Converting shapefiles to GMT"
ogr2ogr -F GMT 50m_land.gmt 50m_land.shp
ogr2ogr -F GMT 50m_coastline.gmt 50m_coastline.shp
ogr2ogr -F GMT 50m_borders.gmt ne_50m_admin_0_boundary_lines_land.shp
ogr2ogr -F GMT 50m_ice.gmt 50m_antarctic_ice_shelves_polys.shp
ogr2ogr -F GMT 50m_ice_borders.gmt 50m_antarctic_ice_shelves_lines.shp
echo "...done\n"

echo "Creating maps"
GMT psxy -m -J$PROJECTION -R$REGION -Bg15 --POLAR_CAP=none --BASEMAP_TYPE=plain --BASEMAP_AXES=N --FRAME_PEN=0.15p,black 50m_land.gmt -K -G254/254/233 > $FILENAME
GMT psxy -m -J$PROJECTION -R$REGION 50m_ice.gmt -O -K -G236/249/255 >> $FILENAME
GMT psxy -m -J$PROJECTION -R$REGION 50m_coastline.gmt -O -K -W0.535p,9/120/171  >> $FILENAME
GMT psxy -m -J$PROJECTION -R$REGION 50m_ice_borders.gmt -O -K -W0.27p,9/120/171  >> $FILENAME
GMT psxy -m -J$PROJECTION -R$REGION 50m_borders.gmt -O -W100/100/100  >> $FILENAME
echo "...done\n"

echo "You can now open $FILENAME in Inkscape to clean up the file and save it as SVG. You should add an ocean coloured background and simplify the graticules using CTRL+L."