out_species.writerow(out_row_species)
file_species.close()
stop_sim = datetime.now()
elapsed_sim = stop_sim - start_sim
print 'time for simulation', elapsed_sim
#ecosystem.draw_species_distribution()
#layout = nx.graphviz_layout(net, prog='dot', args='-Gnodesep=.07, -Granksep=.1, -Grankdir=BT')
layout = nx.circular_layout(net)
#fig = figure()
#network_plot = fig.add_subplot(111)
#nx.draw_networkx(net, layout, ax=network_plot)
print 'S original =', net.order(), 'L original =', net.size(
), 'C original =', net.connectance()
#
# net_final = ecosystem.realised_net.copy()
# layout2 = nx.circular_layout(net_final)
#
# fig2 = figure()
# network_plot2 = fig2.add_subplot(111)
# nx.draw_networkx(net_final, layout2, ax=network_plot2)
#
# print 'S realised =', net_final.order(), 'L realised =', net_final.size(), 'C realised =', net_final.connectance()
#show()
