//Parameters for the coalescence simulation program : simcoal.exe
!!!?pops_per_group*3!!! samples
//Population effective sizes (number of genes 2*diploids)
!!!generate_pop_sizes(?pops_per_group, ?pop_size)!!!
//Samples sizes (number of genes 2*diploids)
!!!dupe(str(?sample_size), ?pops_per_group*3)!!!
//Growth rates	: negative growth implies population expansion
!!!dupe('0', ?pops_per_group*3)!!!
//Number of migration matrices : 0 implies no migration between demes
4
//multimig
!!!generate_goatsie_all_mat(?pops_per_group,?pop_size,?minternal,?me,?mimi)!!!
//mig3
!!!generate_goatsie_three_mat(?pops_per_group,?me,?mimi)!!!
//mig2
!!!generate_goatsie_two_mat(?pops_per_group,?me)!!!
//nothing
!!!generate_null_mat(?pops_per_group*3)!!!
//historical event: time, source, sink, migrants, new deme size, new growth rate, migration matrix index
!!!?pops_per_group*3 - 3 + 5!!! historical events
4000 0 0 1 50 0 3
?expand_gen 0 0 1 0.2 0 3
?split_gen 1 0 1 2 0 3
?split_europe_gen 1 1 1 !!!3.0*?neb/(2*?pop_size)!!! 0 2
!!!?split_europe_gen-1!!! 2 1 1 1 0 2
!!!generate_goatsie_sub_pop_events(?pops_per_group,?hier_split_gen)!!!