#!/bin/csh
# Name:    c.hw08-3
# Author:  Marvin Simkin
# Date:    2002-11-03
# Purpose: Generate plots for homework 8 part 3

# ARGUMENTS (required)
set WESTDEG = $1
set EASTDEG = $2
set SOUTHDEG = $3
set NORTHDEG = $4

# for geophiz...
source /home/guest/.login
source /home/guest/.cshrc


set cptfile = etopo5.cpt
set ingrd = etopo5.grd

set output = hw08-3.ps
/usr/bin/rm $output
touch $output


# ------------------------------------------------------------------------------
# start the postscript file with a title at absolute page coordinates 1, 9.5
##psbasemap -K -B1 -JX6.5i/0.5i -P -R0/1/0/1 -Xa1 -Ya9.5 >> $output

pstext -K -JX6.5i/0.5i -N -P -R0/1/0/1 -Xa1 -Ya9.5 << TITLE >> $output
0.5 0.7 14 0 1 CM Final Project - Area of Study
TITLE


# ------------------------------------------------------------------------------
# panel 1 at absolute page coordinates 1, 1 size 6.5, 8
##set region = -120/-105/25/40
##set region = -118/-108/30/40
##set region = -116/-109/32/38
set region = -115.5/-108.5/31.5/37.5

psbasemap -O -K -B1 -JX6.5i/8i -P -R$region -Xa1 -Ya1 >> $output

grdimage $ingrd -C$cptfile -O -K -R$region -JX6.5i/8id -E100 \
 -Xa1 -Ya1 -P >> $output

pscoast -O -K -R$region -JX6.5i/8id -Df -Na -W3/0/0/0 -I1/1p/0/0/255 -I2/0.25p/0/0/255 -P -Xa1 -Ya1 >> $output


psxy -O -K -R$region -JX6.5i/8id -W5/255/0/0 -P -Xa1 -Ya1 << LINES >> $output
$WESTDEG $SOUTHDEG
$WESTDEG $NORTHDEG
$EASTDEG $NORTHDEG
$EASTDEG $SOUTHDEG
$WESTDEG $SOUTHDEG
LINES




# ------------------------------------------------------------------------------
# to close the plot
pstext -O -JX6.5i/0.5i -N -P -R0/1/0/1 -Xa1 -Ya9.5 << TITLE >> $output
0.5 0.30 6 0 1 CM Phoenix, Arizona, USA
TITLE