D=zeros(17520,1);H=zeros(17520,1);M=zeros(17520,1);

x=reshape(D,48,365);

for i=1:365

   x(:,i)=i;

end

D=reshape(x,17520,1);

h=[0 0 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 11 11 12 12 13 13 14 14 15 15 16 16 17 17 18 18 19 19 20 20 21 21 22 22 23 23]';

y=reshape(H,48,365);

for i=1:365

   y(:,i)=h;

end

H=reshape(y,17520,1);

z=reshape(M,2,17520/2);

z(1,:)=0;z(2,:)=30;M=reshape(z,17520,1);

date=[D H M];

%Data=load('E:模型模式FASST唐古拉测试模拟结果输出tgl_fluxes.out');

flux=tglfluxes;

Flux=[];

for j=1:17520

   a=find(flux(:,2)==date(j,1) & flux(:,3)==date(j,2) & flux(:,4)==date(j,3));

   a=min(a);

   a

   b=flux(a,:);    

   Flux=[Flux;b];

end