close all
clear all

load obs_float/nemed_observations_2.dat

figure (1)
    I1=find( nemed_observations_2(:,3)==1);
    II1=nemed_observations_2(I1,:)
    Iu1=find(II1(:,1)==1);
    Iv1=find(II1(:,1)==2);
    size1u=size(Iu1)
    size1v=size(Iv1)    
    plot(II1(Iu1,4),II1(Iv1,4),'--r');
    hold on;

    I2=find( nemed_observations_2(:,3)==2);
    II2=nemed_observations_2(I2,:)
    Iu2=find(II2(:,1)==1);
    Iv2=find(II2(:,1)==2);
    size2u=size(Iu2)
    size2v=size(Iv2)    
    plot(II2(Iu2,4),II2(Iv2,4),'r');
    hold on;
    
    I3=find( nemed_observations_2(:,3)==3);
    II3=nemed_observations_2(I3,:)
    Iu3=find(II3(:,1)==1);
    Iv3=find(II3(:,1)==2);
    size3u=size(Iu3)
    size3v=size(Iv3)    
    plot(II3(Iu3,4),II3(Iv3,4),'b');
    hold on;
    
    I4=find( nemed_observations_2(:,3)==4);
    II4=nemed_observations_2(I4,:)
    Iu4=find(II4(:,1)==1);
    Iv4=find(II4(:,1)==2);
    size4u=size(Iu4)
    size4v=size(Iv4)    
    plot(II4(Iu4,4),II4(Iv4,4),'m');
    hold on;
    
    I5=find( nemed_observations_2(:,3)==5);
    II5=nemed_observations_2(I5,:)
    Iu5=find(II5(:,1)==1);
    Iv5=find(II5(:,1)==2);
    size5u=size(Iu5)
    size5v=size(Iv5)    
    plot(II5(Iu5,4),II5(Iv5,4),'g');
    hold on;
    
    I6=find( nemed_observations_2(:,3)==6);
    II6=nemed_observations_2(I6,:)
    Iu6=find(II6(:,1)==1);
    Iv6=find(II6(:,1)==2);
    size6u=size(Iu6)
    size6v=size(Iv6)    
    plot(II6(Iu6,4),II6(Iv6,4),'c');
    hold on;
    
    I7=find( nemed_observations_2(:,3)==7);
    II7=nemed_observations_2(I7,:)
    Iu7=find(II7(:,1)==1);
    Iv7=find(II7(:,1)==2);
    size7u=size(Iu7)
    size7v=size(Iv7)    
    plot(II7(Iu7,4),II7(Iv7,4),'k');
    hold on;

% % % % % %     I8=find( nemed_observations_2(:,3)==8);
% % % % % %     II8=nemed_observations_2(I8,:)
% % % % % %     Iu8=find(II8(:,1)==1);
% % % % % %     Iv8=find(II8(:,1)==2);
% % % % % %     size8u=size(Iu8)
% % % % % %     size8v=size(Iv8)    
% % % % % %     plot(II8(Iu8,4),II8(Iv8,4),'--b');
% % % % % %     hold on;