% plot fl vs Press for inject casts 
% see various regimes of fl

load INJ_003_fl
load INJ_003_spd

figure(1)
plot(fl,press,'.')
axis ij
hold on
aa=find(fl>1.5 & press>95);
bb=find(fl<1.5 & fl>.9& press > 100);

plot(fl(aa),press(aa),'g.',fl(bb),press(bb),'r.')
xlabel('Chlorophyll (\mug/l)')
ylabel('Pressure (dbar)')
title('INJ\_003')

figure(2)
plot(yearday,fl,'.')
hold on
plot(yearday(aa),fl(aa),'g.',yearday(bb),fl(bb),'r.')
ylabel('Chlorophyll (\mug/l)')
xlabel('Yearday (2005)')
title('INJ\_003')

figure(5)
cc=find(lon > -60);lon(cc)=nan;lat(cc)=nan;
plot(lon,lat,'.',lon(aa),lat(aa),'.',lon(bb),lat(bb),'.')
title('INJ\_003')

%%%%%%%%%%%%%%%%%%%%%%%%%%%
load INJ_004_fl
figure(3)
plot(fl,press,'.')
axis ij
hold on
aa=find(fl < 1.5 & press < 90);
bb=find(yearday < 203.7103 & fl > 1.5);
plot(fl(aa),press(aa),'g.',fl(bb),press(bb),'r.')
xlabel('Chlorophyll (\mug/l)')
ylabel('Pressure (dbar)')
title('INJ\_004')

figure(4)
plot(yearday,fl,'.')
hold on
plot(yearday(aa),fl(aa),'g.',yearday(bb),fl(bb),'r.')
ylabel('Chlorophyll (\mug/l)')
xlabel('Yearday (2005)')
title('INJ\_004')

figure(6)
plot(lon,lat,'.',lon(aa),lat(aa),'.',lon(bb),lat(bb),'.')
title('INJ\_004')
