%% Plot USGS turbidity
load USGSturbid1.mat

SF1=SFBSSalcmglalldata;    SF8=SFBSSnscarmglalldata;  
SF2=SFBSSm17mglalldata;    SF9=SFBSSnscm1mglalldata;  
SF3=SFBSSmdmbmglalldata;   SF10=SFBSSnsmalmglalldata;  
SF4=SFBSSmdsmbmglalldata;  SF11=SFBSSnspspmglalldata;  
SF5=SFBSSmmicmglalldata;   SF12=SFBSSnsrchmglalldata;  
SF6=SFBSSmp24mglalldata;   
SF7=SFBSSnsbenmglalldata;

figure
for i=1:12
    eval(['y=SF' num2str(i) '(:,2)+(' num2str(i) '*1000);']);
eval(['plot(SF' num2str(i) '(:,1), y)']);
hold on
end