set ROPT = 135.997451/138.622451/35.253403/37.004803
set TXTDATA = Takayama21.txt
set GRDDATA = Takayama21.grd
set FIG	   =  Takayama21.ps
set CPTFILE = CPT.cpt
set TOPT = 0/3000/500
set JOPT = M10


gmtset BASEMAP_TYPE plain
xyz2grd ${TXTDATA} -G${GRDDATA} -I0.000625/0.000417 -R${ROPT}
makecpt -Crainbow -T${TOPT} > ${CPTFILE}
psbasemap -Ba1f0.2/a1f0.2:."Check for Takayama 21x21":  -J${JOPT} -R${ROPT} -K  > ${FIG}
grdimage ${GRDDATA} -R -C${CPTFILE}  -J${JOPT} -O -K >> ${FIG}
grdcontour ${GRDDATA} -R -C${CPTFILE}  -J${JOPT} -O >> ${FIG}
