if os.path.exists(str(path_to_project) + "/" + "head_ogs_" + str(item) + "_" + str(which) + ".txt"): checker.append(True) else: checker.append(False) if all(checker) == True and checker != []: print("All time series have already been extracted. Continuing without checking if content is correct.") else: # extract the time series from the tec files print("Extracting time series...") extract_timeseries(path_to_project, which="mean", process="GROUNDWATER_FLOW") # extract the rfd curve time_time_series, recharge_time_series = extract_rfd( path=path_to_project, rfd=1 ) # plot the time series vs recharge plot_head_timeseries_vs_recharge(path=path_to_project) # write OGS input parameters in DataFrame, but don't return kf because it is ditrubuted Ss, time_step_size, time_steps = get_ogs_parameters(path_to_project, noKf=True) # load from a file with information about the generated field field_info = open(path_to_project + "/field_info" + ".dat", "r") for line in field_info: dim, var, len_scale, mean, seed, geomean, harmean, arimean = line.split() field_info.close() kf_geo = float(geomean) kf_har = float(harmean) kf_ari = float(arimean) S = Ss * aquifer_thickness T_in_geo = kf_geo * aquifer_thickness T_in_har = kf_har * aquifer_thickness T_in_ari = kf_ari * aquifer_thickness # make directory for results
""" script to execute the extract_timeseries() from transect_plot.py. """ from transect_plot import extract_timeseries import sys try: path_to_project = sys.argv[1] except IndexError: path_to_project = input("Insert path to project: ") try: which = sys.argv[2] process = sys.argv[3] rfd = sys.argv[4] plot = sys.argv[5] except IndexError: which = "mean" process = "GROUNDWATER_FLOW" rfd = 1 plot = True rfd = int(rfd) extract_timeseries(path_to_project, which=which, process=process, rfd=rfd) if plot == True: from transect_plot import plot_head_timeseries_vs_recharge plot_head_timeseries_vs_recharge(path_to_project, which=which)