#!/bin/csh

### configuration ###
set OUT=metsite_HKD.ps
set WEST=139
set EAST=146
set NORTH=46
set SOUTH=41
set MARKSIZE=0.2
set FONTSIZE=9
set SCALE=1.8

### site data ####
cat <<EOF >site.txt
140.751 41.802  Hakodate
140.752 42.886  Kucchan
140.969 42.302  Muroran
141.319 43.051  Sapporo
141.784 43.202  Iwamizawa
142.368 43.767  Asahikawa
142.769 42.151  Urakawa
143.202 42.917  Obihiro
143.352 44.335  Monbetsu
EOF

### main routine ###
set SCALEX = `echo "$EAST * 0.8 + $WEST * 0.2" | bc -l`
set SCALEY = `echo "$SOUTH * 0.85 + $NORTH * 0.15" | bc -l`
set SCALELAT = `echo "( $NORTH + $SOUTH ) / 2" | bc -l`

psbasemap -Jm${SCALE} -R${WEST}/${EAST}/${SOUTH}/${NORTH} -Ba2f1 -X5.0 -Y10.0  -P -K -V  >! $OUT
pscoast -Jm -R -Dh -W1 -V -K -Na -S200/230/250 -Lf${SCALEX}/${SCALEY}/${SCALELAT}/100 -O >> $OUT
psxy -Jm -Sc${MARKSIZE} -G250/0/0 -R -V -K -O < site.txt >> $OUT

#(x,y,size,angle,fontno,justify,text
cat site.txt | awk '{print $1-0.27, $2-0.18, FONTSIZE, 0, 0, 0, $3}' FONTSIZE=$FONTSIZE| \
pstext -Jm -R -O -K >> $OUT