def explore_parameter_space(Svals, Nvals, ncomm, bisec, transect=False): for S in Svals: for N in Nvals: beta = mete.get_beta(S, N) if exp(-beta) < 1: comms = [mete.sim_spatial_whole(S, N, bisec, transect=transect, beta=beta) for i in range(0, ncomm)] output_results(comms, S, N, ncomm, bisec, transect, None, None) print S, N
S = int(sys.argv[1]) N = int(sys.argv[2]) ncomm = int(sys.argv[3]) bisec = int(sys.argv[4]) transect = str2bool(sys.argv[5]) abu = sys.argv[6] shrt_name = sys.argv[7] else: S = 10 N = 100 ncomm = 1 bisec = 9 transect = False abu = 'None' shrt_name = None if abu != 'None': datafile = open(abu,'r') datareader = csv.reader(datafile) data = [] for row in datareader: data.append(row) abu = [int(x) for x in data[0]] else: abu = None comms = [mete.sim_spatial_whole(S, N, bisec, transect, abu) for i in range(0, ncomm)] output_results(comms, S, N, ncomm, bisec, transect, abu, shrt_name) print 'Generating simulated community, complete!'