# Restrict contours using -L option
# Plot ocean temperature at depth=0 for January using WOA data
#
###############################################################

gmt pscoast -P -JM12 -Dl -R40/110/-30/30 -B20WeSn -W1 -X5 -Y10 -K > fig.ps
awk '$3 > 0 {print $1, $2, $3}' temp-jan.out > junk
gmt surface junk -R -I0.1 -Gjunk.grd
gmt grdcontour junk.grd -J -C0.2 -L28/30 -A0.2f9 -W1 -O -K >> fig.ps
gmt pscoast -J -Dl -G200 -R -B -W1 -O >> fig.ps