llc_hires/llc_4320/timing.m

% ls -l --full-time pickup_000*data > ~/pickup.txt
fid=fopen('pickup.txt');
tme=[];
while(~feof(fid))
    tmp=fgetl(fid);
    tme=[tme datenum(tmp([44:62]))];
end
fid=fclose(fid);
difftme=diff(tme(3:end))*24*60;
tme(1:3)=[];

clf reset
orient landscape
wysiwyg
plot(tme,difftme,'b.')
axis([datenum(tme(3))-1/24 datenum(tme(end))+1/24 50 85])
grid
set(gca,'xtick',736437:736700)
datetick('x',7,'keeplimits','keepticks')
title('Wallclock minutes per 12 hours of simulation')
ylabel('Wallclock minutes')
xlabel('April / May 2016')
hold on
it=find(difftme>70);
plot(tme(it),difftme(it),'b*')
it=find(difftme>85);
plot(tme(it),84.8,'r*')
print -depsc timing
1	dimitri	1.1	% ls -l --full-time pickup_000*data > ~/pickup.txt
2			fid=fopen('pickup.txt');
3			tme=[];
4			while(~feof(fid))
5			tmp=fgetl(fid);
6			tme=[tme datenum(tmp([44:62]))];
7			end
8			fid=fclose(fid);
9			difftme=diff(tme(3:end))2460;
10			tme(1:3)=[];
11
12			clf reset
13			orient landscape
14			wysiwyg
15			plot(tme,difftme,'b.')
16			axis([datenum(tme(3))-1/24 datenum(tme(end))+1/24 50 85])
17			grid
18			set(gca,'xtick',736437:736700)
19			datetick('x',7,'keeplimits','keepticks')
20			title('Wallclock minutes per 12 hours of simulation')
21			ylabel('Wallclock minutes')
22			xlabel('April / May 2016')
23			hold on
24			it=find(difftme>70);
25			plot(tme(it),difftme(it),'b*')
26			it=find(difftme>85);
27			plot(tme(it),84.8,'r*')
28			print -depsc timing