#!/bin/sh
############################################################
#to    draw global map
#on    31/May/2002
#by    nhanasaki
#at    IIS,UT
############################################################
GRD=temp.grd
XYZ=temp.xyz
############################################################
#---Get File name
############################################################
ONEFILE=$1
CPTFILE=$2
PSFILE=$3
TITLE1=$4
TITLE2=$5
PFLAG=$6
OPT=$7
if [ "$3" = "" ]; then
  echo "Usage: cmap2ps CMAP CPT PS [Title1] [Title2] [-P] [OPT]"
  echo "OPT:   whitesea  blacksea  bluesea  graysea"
  echo "OPT:   whiteland blackland"
  exit
fi
if [ "$TITLE1" = "." ]; then
  TITLE1=
fi
if [ "$TITLE2" = "." ]; then
  TITLE2=
fi
if [ "$PFLAG" = "." ]; then
  PFLAG=
fi
############################################################
#---Convert one file to xyz file
############################################################
cmap2xyz $ONEFILE $XYZ
############################################################
#---Map Projection
############################################################
D=2.5
XMIN=-180.00
XMAX=180.00
XMINL=-180
XMAXL=0
XMINR=0
XMAXR=180
YMIN=-90.00
YMAX=90.00
XWID=21.0
YWID=10.5
XWIDHLF=10.5
DXa=30.0		#a:Anotation
DXf=30.0		#f:Frame
DXg=10.0		#G:Grid
DYa=30.0		#a:Anotation	(JPN:Memori)
DYf=30.0		#f:Frame
DYg=10.0		#g:Grid
############################################################
#---Short Cuts
############################################################
RFLAG="-R${XMIN}/${XMAX}/${YMIN}/${YMAX}"
RFLAGL="-R${XMINL}/${XMAXL}/${YMIN}/${YMAX}"
RFLAGR="-R${XMINR}/${XMAXR}/${YMIN}/${YMAX}"
JFLAG="-JX${XWID}d/${YWID}d"
JFLAGHLF="-JX${XWIDHLF}d/${YWID}d"
BFLAG="-Ba${DXa}f${DXf}::/a${DYa}f${DYf}::WSne"
############################################################
#---Jobs---
############################################################
psbasemap   $RFLAG $JFLAG $BFLAG -X5.0 -Y5.0   $PFLAG  -K  > $PSFILE 
awk '{print $1, $2, $3}' $XYZ | \
xyz2grd $RFLAG -G$GRD -I${D}/${D} -F
grdimage -O $RFLAG $JFLAG $GRD -C$CPTFILE              -K >> $PSFILE 
psscale  -O -C$CPTFILE -D10.5/-1.5/21.0/0.5h -L          -K >> $PSFILE
if [ "$TITLE1" = "" ]; then
  echo No main title
else
  echo With main title ${TITLE1}
  pstext -O $RFLAG $JFLAG -N                  -K   << EOF >> $PSFILE
  0 110 20 0.0 1 6 ${TITLE1}
EOF
fi
if [ "$TITLE2" = "" ]; then
  echo No sub title
else
  echo With sub title ${TITLE2}
  pstext -O $RFLAG $JFLAG -N                   -K  << EOF >> $PSFILE
  150 110 20 0.0 1 6 ${TITLE2}		
EOF
fi
if [ "$OPT" = "" ]; then
  pscoast -O $RFLAGL $JFLAGHLF -Dc -W5                          -K >> $PSFILE
  pscoast -O $RFLAGR $JFLAGHLF -Dc -W5 -X$XWIDHLF                  >> $PSFILE
elif [ "$OPT" = "graysea" ]; then
  pscoast -O $RFLAGL $JFLAGHLF -Dc -W5            -S51/51/51    -K >> $PSFILE
  pscoast -O $RFLAGR $JFLAGHLF -Dc -W5 -X$XWIDHLF -S51/51/51       >> $PSFILE
elif [ "$OPT" = "bluesea" ]; then
  pscoast -O $RFLAGL $JFLAGHLF -Dc -W5            -S128/128/255 -K >> $PSFILE
  pscoast -O $RFLAGR $JFLAGHLF -Dc -W5 -X$XWIDHLF -S128/128/255    >> $PSFILE
elif [ "$OPT" = "blacksea" ]; then
  pscoast -O $RFLAGL $JFLAGHLF -Dc -W5            -S0/0/0       -K >> $PSFILE
  pscoast -O $RFLAGR $JFLAGHLF -Dc -W5 -X$XWIDHLF -S0/0/0          >> $PSFILE
elif [ "$OPT" = "whitesea" ]; then
  pscoast -O $RFLAGL $JFLAGHLF -Dc -W5            -S255/255/255 -K >> $PSFILE
  pscoast -O $RFLAGR $JFLAGHLF -Dc -W5 -X$XWIDHLF -S255/255/255    >> $PSFILE
elif [ "$OPT" = "blackland" ]; then
  pscoast -O $RFLAGL $JFLAGHLF -Dc -W5            -G0/0/0       -K >> $PSFILE
  pscoast -O $RFLAGR $JFLAGHLF -Dc -W5 -X$XWIDHLF -G0/0/0          >> $PSFILE
elif [ "$OPT" = "whiteland" ]; then
  pscoast -O $RFLAGL $JFLAGHLF -Dc -W5            -G255/255/255 -K >> $PSFILE
  pscoast -O $RFLAGR $JFLAGHLF -Dc -W5 -X$XWIDHLF -G255/255/255    >> $PSFILE
fi
############################################################
#
############################################################
/bin/rm -f .gmtcommands
exit 0





