#!/bin/csh -f
source /usr/local/cshrc/gmt
set INFILE1 = asia.lld
set INFILE2 = control.llduv
set PSFILE  = asia.ps
set TITLE = "some variable of water balance"
set USERNAME = "Kengo MIYAOKA (digitizing) and Kooiti MASUDA (mapping)"
set PROJ   = '-R67.5/0/-175/50/r -Js95/90/1:60000000'
set OFF = '-X1.75 -Y5.0'

gmtset MEASURE_UNIT       cm
gmtset WANT_EURO_FONT     false
gmtset BASEMAP_TYPE       plain
gmtset HEADER_FONT_SIZE   18
gmtset DEGREE_FORMAT      3
gmtset ANOT_MIN_ANGLE     40

pscoast $PROJ -K -P -U/0/-3/"Polar stereographic projection" \
  $OFF -V -Dl -W3/0/255/0 -I1/1/0/0/255 -I2/1/0/0/255 \
  -B10g10a10:::."${TITLE}":WESN \
  > $PSFILE

awk '{print $1, $2}' $INFILE1 | \
  psxy $PROJ -O -K -Sc0.1 -G255/255/0 -V >> $PSFILE

awk '{print $1, $2}' $INFILE2 | \
  psxy $PROJ -O -K -Sc0.25 -W1/255/0/0 -V >> $PSFILE

pstext $PROJ -O -N << EOD >> $PSFILE
 95  1.0  16  0  0  MC $USERNAME 
EOD

\rm -f .gmtcommands .gmtdefaults