function [ TE ] = calcExtinctions ( D ) %calculate the mean times to extinction; TE==0 is not extinct
[ nReps maxTime nPops ] = size ( D );
ME = zeros ( nReps, 1 ); %metapopulation extinction - 1=yes; 0=no
TE = zeros ( nReps, 1 ); %time to metapopulation extinction
for irep = 1 : nReps
Dr = squeeze ( D ( irep, :, : ) );
extpops = find ( Dr ( end, : ) == 0 ); %indicies to the populations that went extinct
numextpops = length ( extpops ); %number of populations that went extinct
ME ( irep ) = numextpops==nPops;
if ME ( irep ) %if the metapopulation went extinct - calculate the time at which it happened
mpeind = find ( all ( ( Dr == 0 )' ) );
TE ( irep ) = mpeind ( 1 )+1;
end
end