#
# magind_cor0.sh
#
 gmtset MEASURE_UNIT inch
#
 cat <<END > mag1.cpt
# Standard Rainbow colortable
# COLOR_MODEL = HSV
-800	220	1	1	-400	200	1	1
-400	200	1	1	-200	180	1	1
-200	180	1	1	-100	160	1	1
-100	160	1	1	-50	140	1	1
-50	140	1	1	-20	120	1	1
-20	120	1	1	0	100	1	1
0	100	1	1	20	80	1	1
20	80	1	1	50	60	1	1
50	60	1	1	100	40	1	1
100	40	1	1	200	20	1	1
200	20	1	1	400	0	1	1
400	360	1	1	800	340	1	1
F	340	1	1
B	220	1	1
END
#
 tail +2 indcor.res | blockmean -R75/115/-12/13 -I5m -V > indcor.dat
#
# gridding
 surface indcor.dat -R75/115/-12/13 -I5m -Gindcor.grd -T0.25 -V
#
# area
 psbasemap -R -Jm1:20000000.0 -X1.5 -Y1.8 -Bf5a10WSne -K -V > indcor.ps
#
# gradient calculation
 grdgradient indcor.grd -Gindgrad.grd -A135 -Nt -M -V
#
# shaded color image
 grdimage indcor.grd -Cmag1.cpt -Iindgrad.grd -Jm -R -O -K -V >> indcor.ps
#
# mask area
 psmask indcor.dat -I5m -Jm -G255 -S30k -O -K -R -B -N -T -V >> indcor.ps
#
# draw a color palette
 psscale -Cmag1.cpt -D4.0/-0.5/6.0/0.40h -B::/:nT: -I -L -O -V -K >> indcor.ps
#
# coast lines
 pscoast -R -Jm -Dl -G200/200/200 -W2 -O -V >> indcor.ps
#
 rm indcor.dat indcor.grd indgrad.grd mag1.cpt
#