#!/bin/csh 
# plot index map for Arctic Ocean

gmtset `cat gmt.plot.parameters` D_FORMAT %.4lf
set scale = "-JS0/90/6.0i -P"
set box = -R230/79/180/70r

mapproject -R-45/53:49:1.4687/135/53:49:1.4687r -Js0/90/75/1:5805000 -I -Fk -V -C << END
-1851.637 2251.637
859.0634 -459.0634
END

makecpt -Cglobe > ../figures/data/ibcao.cpt

grdimage ../figures/data/IBCAO.new.grd -JX6i -R-1851.637/859.0634/-459.0634/2251.637 -C../figures/data/ibcao.cpt -P -K > IBCAO.ps

psxy -JS0/90/6i -R283.9242/72.7438/159.1167/68.29398r -W8/255/255/0 -K -O << END >> IBCAO.ps
230 79
180 83
180 70
202.5 68.5
230 79
END

cat ../figures/data/"Healy 05-03"/track.pts | psxy -JS0/90/6i  -R283.9242/72.7438/159.1167/68.29398r  -W2 -K -O >> IBCAO.ps

foreach file (../figures/data/"Healy 05-03"/05*.pts)
	cat "$file" | awk '{printf"%s %s\n",$2,$3}' | psxy -JS0/90/6i  -R283.9242/72.7438/159.1167/68.29398r  -W8/255/125/255 -K -O >> IBCAO.ps
end

cat ../figures/data/"Healy 05-03"/0520a.pts | awk '{printf"%s %s\n",$2,$3}' | psxy -JS0/90/6i  -R283.9242/72.7438/159.1167/68.29398r  -W8/55/200/55 -K -O >> IBCAO.ps
cat ../figures/data/"Healy 05-03"/0520b.pts | awk '{printf"%s %s\n",$2,$3}' | psxy -JS0/90/6i  -R283.9242/72.7438/159.1167/68.29398r  -W8/55/200/55 -K -O >> IBCAO.ps

pstext -JX6i -R0/6/0/6 -G255/255/0 -K -O << END >> IBCAO.ps
2.1 2.8 12 0 1 1 Figure 6
END

pstext -JX6i -R0/6/0/6 -G128 -K -O << END >> IBCAO.ps
2.85  4.1 10 45 1 1 Chukchi 
2.9 3.9 10 45 1 1 Borderland
END

pstext -JX6i -R0/6/0/6 -K -O << END >> IBCAO.ps
4.3 5.15 10 0 1 1 USCGC Healy 05-03
4.5 5.0 10 0 1 1 MCS Profiles
4.5 4.85 10 0 1 1 Track
4.5 4.7 10 0 1 1 Figure 1 (line 05-20)
2.8 3.70 8 90 1 1 Northwind Ridge
3.1 3.9 10 0 1 1 Chukchi
3.1 3.75 10 0 1 1 Cap
3.1 5.2 12 0 1 1 Chukchi 
3.1 5.0 12 0 1 1 Shelf
1.7 3.8 14 0 1 1 Canada 
1.7 3.60 14 0 1 1 Basin
4.0 3.95 6 90 1 1 Mendeleev Ridge
END

psxy -JX6i -R0/6/0/6  -W8/255/125/255  -K -O << END >> IBCAO.ps
4.2 4.9 
4.45 5.1
END

psxy -JX6i -R0/6/0/6  -W2  -K -O << END >> IBCAO.ps
4.2 4.75 
4.45 4.95
END

psxy -JX6i -R0/6/0/6  -W8/55/200/55  -K -O << END >> IBCAO.ps
4.2 4.6
4.45 4.8
END

psbasemap -JS0/90/6i  -R283.9242/72.7438/159.1167/68.29398r  -Ba45g90/a10g5WEsn -O >> IBCAO.ps