elevation=lambda(1).zd2;
for i=1%4:2:13
    %elevation(j)
    count=1;
%     h=figure;
    
    for j=4
        elevation(j)
        lambda(i).trans
        psd=lambda(i).psd;
        f=lambda(i).fqcy;
%          subplot(4,2,count)
         
         semilogx(f(:,j),psd(:,j))
%          count=count+1;
%          
%          zdata=lambda(i).z_data;
%          xdd=lambda(i).xdd;
%          
%          
%          subplot(4,2,count)
%          plot(xdd(j,:),zdata(:,j))
%                   count=count+1;

        1./ginput
    end
    %text(-15573.34,60.41,['vol0034321 rect at el ',num2str(elevation(j)),'m']);
    %saveas(h,['C:\Documents and Settings\sylvie lorsolo\My Documents\research\frances\images\wavelength analysis\rect031007',num2str(elevation(j)),'.fig'])

end