#!/bin/csh -f
# Create Orthographic map of Satellite derived bathymetry
set GMT_GRIDDIR = $DAR/geosat/grdraster.info
grdraster 8 -R0/360/-90/90 -Gworld_relief_30s.grd -V

pscoast -R0/360/-90/90 -JG280/30/6i -Bg30/g15 -Dc -A5000 -Gwhite -SDarkTurquoise -P >! map.ps

set grd = GlobalTopography_30s.grd
makecpt -Cglobe  > globe.cpt

set grd = world_relief_2m.grd
#grdgradient $grd -G$grd.grad -A315 -Nt1 -V
set cpt = -C/usr/local/gmt/share/cpt/GMT_globe.cpt
set R = -R0/360/-90/90 
grdimage $grd $cpt $R -JG200/22/6i -Bg30/g30 -K -V >! globe.ps
psxy inset.xy $R -JG200/22/6i -O -V >> globe.ps

pscoast -R0/360/-90/90  -JG200/22/6i -Bg30/g30 -Dc -A5000 -Gblack -Swhite -P -K -V >! map.ps
psxy inset.xy -R0/360/-90/90  -JG200/22/6i -M -L -O -K -V >> map.ps
awk '{print $1,$2}' atoll.xy | psxy $R -JG200/22/6i -Sa0.1i -G255/0/0 -O -V >> map.ps

grdraster 22 -Rd-180/-155/-10/25 -Gmuller_age.grd -V

set R = -R0/360/-90/90 
pscoast $R -JG180/10/6i -Bg30/g30 -Dc -A5000 -Gblack -Swhite -P -K -V >! map.ps


pscoast -R0/360/-90/90  -JG180/10/6i -Bg30/g30 -Dc -A5000 -Gblack -Swhite -P -K -V >! map.ps
psxy inset.xy -R0/360/-90/90  -JG180/10/6i -M -L -O -K -V >> map.ps
awk '{print $1,$2}' atoll.xy | psxy $R -JG180/10/6i -Sa0.1i -G255/0/0 -O -V >> map.ps

#grdgradient $grd -G$grd.grad -A315 -Nt1 -V
set grd = world_relief_2m.grd
set cpt = -C/usr/local/gmt/share/cpt/GMT_globe.cpt
set R = -R0/360/-90/90
set J = -JG180/10/6i
set B = -Bg30/g30
grdimage $grd $cpt $R $J $B -K -V >! gradsemglobe.ps
psxy inset.xy $R  $J -M -L -O -K -V >> gradsemglobe.ps
awk '{print $1,$2}' atoll.xy | psxy $R $J -Sa0.1i -G255/0/0 -O -V >> gradsemglobe.ps

