Beispiel #1
0
def write_parac(input_c_par, param_name, var_ar, main_directory, d):

    for e in var_ar:
        mkdir_p(main_directory + param_name + "/" + e)
        c_par = open(name_c_par, "r")
        print name_c_par
        path_new = main_directory + param_name + "/" + e + "/"
        c_new = open(path_new + "param.in", "w")
        for linee in c_par.readlines():
            if linee == "   output_path  .\n":
                c_new.write("output_path " + path_new[0:-1] + "\n")
            # 				print "   output_path "+ path_new[0:-1]
            elif (param_name in linee) and (len(linee) < 50):
                c_new.write("   " + param_name + " " + e + "\n")
            # 				print "   "+param_name+" "+e

            else:
                c_new.write(linee)

            if linee[0:12] == "   output_dt":
                timestep = float(linee[12:-1])
                print timestep

            if linee[0:8] == "   x_max":
                length = float(linee[8:-1])
            if linee[0:9] == "   n_cell":
                n_cell = float(linee[9:-1])

        c_par.close()
        c_new.close()
        c = {}
        for key in d:
            c[key] = d[key]
            # print(c['
        c["solver_call"] = d["solver_call"]
        c["solver_input"] = path_new + "param.in"
        c["solver_dir"] = d["solver_dir"]
        c["output_path"] = path_new
        c["timestep"] = timestep
        c["cell_length"] = length / n_cell
        execute_it(c)

        m_analyse.plott_all(c)
        m_global.plot_energy(c)
Beispiel #2
0
# def execute_them(list_input):


script, input = argv

d = read_parameter(input)


if rs(d["var_param"]) == "no":
    if rs(d["execute"]) == "yes":
        execute_it(d)

    add_units(d)

    m_analyse.plott_all(d)
    m_global.plot_energy(d)

else:
    varfile = rs(d["parameter_set_path"]) + rs(d["parameter_set"])
    vf = open(varfile, "r")
    set_var = []
    for line in vf.readlines():
        set_var = set_var + [line[:-1]]

    vf.close()
    for element in set_var:
        varpar = element

        param_file = open(rs(d["parameter_set_path"]) + varpar, "r")
        i = 0
        var_ar = []