data=temp.out
out=contour1

1:
# 2D contour map using grdcontour
# Plot ocean temperature at depth=0 for January using WOA data
#
	makecpt  -T20/30/1 -Crainbow > color.cpt
	pscoast -P -JM12 -Dl -R40/110/-30/30 -B20WeSn -W3 -X5 -Y10 -K > $(out).ps
	awk '$$1 == 1 && $$2 == 0 && $$5 > 0 {print $$4, $$3, $$5}' $(data) > junk
	surface junk -R -I0.1 -Gjunk.grd
	grdcontour junk.grd -J -G2 -C1 -A1f8 -W3 -O -K >> $(out).ps
	pscoast -J -G200 -R -B -W2 -O >> $(out).ps
	convert $(out).ps $(out).jpg