namimuki=240
yyyy=2012;mm=11;dd=30;hh=18;
date_next=(yyyy)*1000000+(mm)*10000+(dd)*100+(hh);
for j=2:83
Cg=dlmread(['Cg_',num2str(date_next),'.txt'],'',0,0);
Cg(Cg<0)=0;
matrixlonlat=dlmread(['matrixlonlat',num2str(namimuki),'.txt'],'',0,0);
distance=dlmread([num2str(namimuki),'distance.txt'],'',0,0);
matrixlat=matrixlonlat(:,1);
matrixlon=matrixlonlat(:,2);
Cg_next=Cg(matrixlat(j,:),matrixlon(j,:));
time=round(distance(j)*1000./Cg_next/3600);
jikoku_now=datetime([yyyy,mm,dd]) + hours(hh);
jikoku_next=jikoku_now-hours(time);
jikoku_next=datevec(jikoku_next');
year=jikoku_next(:,1);
month=jikoku_next(:,2);
day=jikoku_next(:,3);
hour=jikoku_next(:,4);
format long
date_next=year*1000000+month*10000+day*100+hour
save(['time',num2str(date_next),'.txt']);
goto 11
return
end