%plot points on track line indexed into colorbar
% modified 8/29/01 on NBP0104
% want to plot over Etopo2.
%;
% x = lon, y=lat, f=data
cruiseids=['9808';]
[rcid, ccid]=size(cruiseids);

for idex=1:rcid
   cruise=cruiseids(idex,:)
   load(['ESS' cruise '.mat']) % day fluor sal temp xlon ylat GB

	%eval(['GODAT=DAT' num2str(numlist(k)) ';'])
GODAT=f';
gbmapls2000
hold on

%% PLOT SCALED DATA (64 colors) ON TOP OF SOME MAP
factorby=(mean(GODAT)-min(GODAT))./7
pieceplot1=plot(x, y, 'k');
grayt=gray;
set(pieceplot1, 'linew', 8)

for j=1:64
hiGO=ceil((1+GODAT)*10);
ec=grayt(hiGO, :);
i=find(abs(hiGO)>=j & abs(hiGO)<(j+1));

%pieceplot=plot(xreplot(i), yreplot(i), coldots(j,:));
pieceplot=plot(x(i), y(i), 'w.');

set(pieceplot, 'markersize', 5, 'markerfacecol', grayt(j,:), ...
'markeredgecol', grayt(j,:), 'marker', 'o')
ax1=gca;
end
%title(['Variation in  Fluorescence on ' cruise])
h=colorbar;colormap(gray);set(h, 'yticklabel',...
'0|0.5|1.0|1.5|2.0|2.5|3.0|3.5|4.0|4.5|5.0 VOLTS')

		xf=xlabel('Degrees Longitude ')
		yf=ylabel('Degrees Latitude ')
		set(xf, 'fontsize', 16, 'fontweight', 'bold');
		set(yf, 'fontsize', 16, 'fontweight', 'bold');
		set(ax1, 'fontsize', 16, 'fontweight', 'bold');
		orient landscape
        print -djpeg99 dataorth.jpg
		%axis tight
	%	[m,map]=capture;
%		imwrite(m,map, [cruise 'fluorortho.tif']);
%		print -dtiffnocompression fluortho.tif
%		print -dpict fluortho.pict
%%print -dps2
%%end
end
