def main():
gb = import_gb("cut", 2)
param = get_param()
scheme = Scheme(gb)
scheme.set_data(param)
# exporter
save = pp.Exporter(gb, "case5", folder_name="solution")
vars_to_save = scheme.vars_to_save()
# post process
scheme.extract()
save.write_vtk(vars_to_save, time_step=0)
for i in np.arange(param["time"]["num_steps"]):
print("processing", i, "step at time", i * param["time"]["step"])
# do one step of the splitting scheme
scheme.one_step_splitting_scheme()
# post process
time = param["time"]["step"] * (i + 1)
save.write_vtk(vars_to_save, time_step=time)
time = np.arange(param["time"]["num_steps"] + 1) * param["time"]["step"]
save.write_pvd(time)
def main():
mesh_size = np.power(2., -7)
gb = create_gb(mesh_size)
param = get_param()
scheme = Scheme(gb)
scheme.set_data(param)
# exporter
folder = "./case2/"
save = pp.Exporter(gb, "case2", folder_name=folder+"solution_ml")
vars_to_save = scheme.vars_to_save()
# post process
save.write_vtk(vars_to_save, time_step=0)
for i in np.arange(param["time"]["num_steps"]):
time = (i+1)*param["time"]["step"]
print("processing", i, "step at time", time)
# do one step of the splitting scheme
scheme.one_step_splitting_scheme(time)
# post process
save.write_vtk(vars_to_save, time_step=time)
time = np.arange(param["time"]["num_steps"]+1)*param["time"]["step"]
save.write_pvd(time)