rem ***************************************************************************
rem Make input data file inp_Pcolsize.txt
rem   x  y  za  zb  (za: color, zb: size)
rem ***************************************************************************
echo  10 20 20 0.1 > inp_Pcolsize.txt
echo  20 30 30 0.3 >> inp_Pcolsize.txt
echo  30 50 50 0.5 >> inp_Pcolsize.txt
echo  40 80 80 0.7 >> inp_Pcolsize.txt
echo  50 90 90 1.0 >> inp_Pcolsize.txt
echo  60 70 70 0.8 >> inp_Pcolsize.txt
echo  70 60 60 0.6 >> inp_Pcolsize.txt
echo  80 30 30 0.4 >> inp_Pcolsize.txt
echo  90 15 15 0.2 >> inp_Pcolsize.txt
rem ***************************************************************************
rem Make color pallet file inp_PcolCPT.cpt
rem   z1  R1  G1  B1  z2  R1  G2  B2 (z1/z2: range of value)
rem ***************************************************************************
echo   0   0   0   0  20   0   0   0 > inp_PcolCPT.cpt
echo  20 255   0   0  40 255   0   0 >> inp_PcolCPT.cpt
echo  40   0 255   0  60   0 255   0 >> inp_PcolCPT.cpt
echo  60   0   0 255  80   0   0 255 >> inp_PcolCPT.cpt
echo  80 255   0 255 100 255   0 255 >> inp_PcolCPT.cpt
set range=0/100/0/100
set scale=5/5
set xga=a20
set yga=a20
set xlabel="Name of x-axis"
set ylabel="Name of y-axis"
set inp_dat=inp_Pcolsize.txt
set fig=fig_Pcolsize.eps
gmtset ANOT_FONT_SIZE 10
gmtset LABEL_FONT_SIZE 10
gmtset TICK_LENGTH -0.1c
psbasemap -R%range% -JX%scale% -B%xga%:%xlabel%:/%yga%:%ylabel%:WSen -P -K > %fig%
psxy %inp_dat% -R -JX -W2 -K -O >> %fig%
psxy %inp_dat% -R -JX -SC -Cinp_PcolCPT.cpt -N -O >> %fig%