close all
clear all
dirnem='nemed/';
nemlist='liste_perso';
format long
fid=fopen([dirnem nemlist],'r');
count=1;
float_obs=[];
index=[];
while 1
            tline = fgetl(fid);
                 disp(tline)
            if ~ischar(tline), break, end
            lon=ncread([dirnem tline],'longitudes');
            lat=ncread([dirnem tline],'latitudes');
            pdt=ncread([dirnem tline],'time');
            
            M=zeros(size(lat,1),4);
            M(:,1)=count;
            M(:,2)=pdt/365.25;
            M(:,3)=lon;
            M(:,4)=lat;
            [m1,n1]=size(M);
            index=[index;m1]
            float_obs=[float_obs;M];
            count=count+1;
end

[m,n]=size(float_obs)
fclose(fid);


% figure (1)
%     I=find(float_obs(:,1)==1);
%     plot(float_obs(I,3),float_obs(I,4),'m');
%     hold on;
%   
%     I2=find(float_obs(:,1)==2);
%     I2
%     size2=size(I2)
%     plot(float_obs(I2,3),float_obs(I2,4),'r');
%     hold on;
%     
%     I=find(float_obs(:,1)==3);
%     plot(float_obs(I,3),float_obs(I,4),'b');
%     hold on;
%     
%     I=find(float_obs(:,1)==4);
%     plot(float_obs(I,3),float_obs(I,4),'g');
%     hold on;
%     
%     I=find(float_obs(:,1)==5);
%     plot(float_obs(I,3),float_obs(I,4),'g');
%     hold on;
%     
%     I6=find(float_obs(:,1)==6);
%     I6
%     size6=size(I6)
%     plot(float_obs(I6,3),float_obs(I6,4),'c');
%     hold on;
%     
%     I=find(float_obs(:,1)==7);
%     plot(float_obs(I,3),float_obs(I,4),'k');
%     hold on;
%     
%     I=find(float_obs(:,1)==8);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
    
    
%     figure(2)
%     axis([0 37 2008 2012])
% hold on
% grid on 
% hold on
%     
%     I=find(float_obs(:,1)==1);
%     plot(float_obs(I,1),float_obs(I,2),'m');
%     hold on;
%   
%     I=find(float_obs(:,1)==2);
%     plot(float_obs(I,1),float_obs(I,2),'r');
%     hold on;
%     
%     I=find(float_obs(:,1)==3);
%     plot(float_obs(I,1),float_obs(I,2),'b');
%     hold on;
%     
%     I=find(float_obs(:,1)==4);
%     plot(float_obs(I,1),float_obs(I,2),'m');
%     hold on;
%     
%     I=find(float_obs(:,1)==5);
%     plot(float_obs(I,1),float_obs(I,2),'g');
%     hold on;
%     
%     I=find(float_obs(:,1)==6);
%     plot(float_obs(I,1),float_obs(I,2),'c');
%     hold on;
%     
%     I=find(float_obs(:,1)==7);
%     plot(float_obs(I,1),float_obs(I,2),'k');
%     hold on;
%     
%     
%     I=find(float_obs(:,1)==8);
%     plot(float_obs(I,1),float_obs(I,2),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==9);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==10);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==11);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==12);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==13);
%     plot(float_obs(I,3),float_obs(I,4),'--r');
%     hold on;
%     
%      I=find(float_obs(:,1)==14);
%     plot(float_obs(I,3),float_obs(I,4),'--k');
%     hold on;
%     
%      I=find(float_obs(:,1)==15);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==16);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==17);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==18);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==19);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==20);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==21);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==22);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==23);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==24);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==25);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==26);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==27);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==28);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==29);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==30);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==31);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==32);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==33);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==34);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==35);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%     
%      I=find(float_obs(:,1)==36);
%     plot(float_obs(I,3),float_obs(I,4),'--b');
%     hold on;
%  
    
    
   firstcol=ones(m,1);
    
    float_obs(:,5)=float_obs(:,1);
    float_obs(:,1)=ones(m,1);
    float_obs(:,6)=float_obs(:,3);
    float_obs(:,3)=float_obs(:,5);
    float_obs(:,7)=float_obs(:,4);
    float_obs(:,4)=float_obs(:,6);
    
    

    
    
    outfile='obs_float/nemed_observations.dat';
format long
 files=fopen(outfile,'w');
 fprintf(files,'%d %d %d %f\n',float_obs(:,1:4)');
 
     float_obs(:,1)=firstcol+1;

 fprintf(files,'%d %d %d %f\n',float_obs(:,[1,2,3,7])');

 fclose(files);
 
 
load obs_float/meshgrid2.dat

Nlat=58; %Phi
Nlong=87; %Lambda
GR_long=reshape(meshgrid2(:,1),Nlong,Nlat); 
GR_lat=reshape(meshgrid2(:,2),Nlong,Nlat); 
load obs_float/nemed_observations.dat
I=find(nemed_observations(:,1)==1);
hold on
xi=interp1(GR_long(:,1),1:Nlong,nemed_observations(I,end));
xi(1:10)
I=find(nemed_observations(:,1)==2);
yi=interp1(GR_lat(1,:),1:Nlat,nemed_observations(I,end));

    float_obs(:,1)=firstcol;

    float_obs(:,4)=xi;
    float_obs(:,7)=yi;

    
 outfile='obs_float/nemed_observations_2.dat';
format long
 files=fopen(outfile,'w');
 fprintf(files,'%d %d %d %f\n',float_obs(:,1:4)');
 
     float_obs(:,1)=firstcol+1;

 fprintf(files,'%d %d %d %f\n',float_obs(:,[1,2,3,7])');

 fclose(files);

%  
%  outfile='scripts/matlab_toolbox/nemed_observations_final.dat';
% load scripts/matlab_toolbox/nemed_observations_2.dat;
% [d,f]=size(nemed_observations_2)
% floater1=zeros(d,f);
% count1=1;
% count2=1;
% nboff=8;
% for i=1:d/2
%         floater(count1,:)=nemed_observations_2(i,:);
%         floater(count1+1,:)=nemed_observations_2(i+d/2,:);
% 
%         count1=count1+2;
%    
% end
% 
% fid=fopen(outfile,'w');
% fprintf(fid,'%d %d %d %f\n',floater');
% fclose(fid);
%  
% 




