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!'