def export_tables(silent=True): if config.plotSave: caption = "Infiltration and pumping rates for example {}.".format( sim_name ) headings = ( "Stress period", "Infiltration rate", "Pumping rate", ) fpth = os.path.join("..", "tables", "{}-02.tex".format(sim_name)) dtype = [ ("nper", "U30"), ("infilt", "U30"), ("rate", "U30"), ] arr = np.zeros(nper, dtype=dtype) for n in range(nper): arr["nper"][n] = bt.int_format(n + 1) arr["infilt"][n] = bt.exp_format(infiltration[n]) arr["rate"][n] = bt.float_format(well_rates[n]) if not silent: print("creating...'{}'".format(fpth)) col_widths = (0.1, 0.30, 0.30) bt.build_table( caption, fpth, arr, headings=headings, col_widths=col_widths )
def export_tables(silent=True): if config.plotSave: name = list(parameters.keys())[1] caption = "Aquifer properties for example {}.".format(sim_name) headings = ( "Layer", "Thickness", "Hydraulic conductivity", "Initial head", ) fpth = os.path.join("..", "tables", "{}-01.tex".format(sim_name)) dtype = [ ("k", "U30"), ("thickness", "U30"), ("k33", "U30"), ("h0", "U30"), ] arr = np.zeros(nlay, dtype=dtype) for k in range(nlay): arr["k"][k] = bt.int_format(k + 1) arr["thickness"][k] = bt.float_format(zelevs[k] - zelevs[k + 1]) arr["k33"][k] = bt.exp_format(k33[k]) arr["h0"][k] = bt.float_format(strt[k]) if not silent: print("creating...'{}'".format(fpth)) col_widths = (0.1, 0.15, 0.30, 0.25) bt.build_table( caption, fpth, arr, headings=headings, col_widths=col_widths ) caption = "Interbed properties for example {}.".format(sim_name) headings = ( "Interbed", "Layer", "Thickness", "Initial stress", ) fpth = os.path.join("..", "tables", "{}-02.tex".format(sim_name)) dtype = [ ("ib", "U30"), ("k", "U30"), ("thickness", "U30"), ("pcs0", "U30"), ] arr = np.zeros(len(ib_ctype), dtype=dtype) for idx, ctype in enumerate(ib_ctype): arr["ib"][idx] = bt.int_format(idx + 1) arr["k"][idx] = bt.int_format(ib_cellid[idx][0] + 1) if ctype == "nodelay": arr["thickness"][idx] = bt.float_format(ib_thickness[idx]) else: b = ib_thickness[idx] * ib_rnb[idx] arr["thickness"][idx] = bt.float_format(b) arr["pcs0"][idx] = bt.float_format(parameters[name]["pcs0"][idx]) if not silent: print("creating...'{}'".format(fpth)) bt.build_table(caption, fpth, arr, headings=headings) caption = "Aquifer storage properties for example {}.".format(sim_name) headings = ( "Layer", "Specific Storage", ) fpth = os.path.join("..", "tables", "{}-03.tex".format(sim_name)) dtype = [("k", "U30"), ("ss", "U30")] arr = np.zeros(4, dtype=dtype) for idx, k in enumerate((4, 6, 11, 13)): arr["k"][idx] = bt.int_format(k + 1) arr["ss"][idx] = bt.exp_format(parameters[name]["cg_ske"][k]) if not silent: print("creating...'{}'".format(fpth)) col_widths = (0.1, 0.25) bt.build_table( caption, fpth, arr, headings=headings, col_widths=col_widths ) caption = "Interbed storage properties for example {}.".format( sim_name ) headings = ( "Interbed", "Layer", "Inelastic \\newline Specific \\newline Storage", "Elastic \\newline Specific \\newline Storage", ) fpth = os.path.join("..", "tables", "{}-04.tex".format(sim_name)) dtype = [ ("ib", "U30"), ("k", "U30"), ("ssv", "U30"), ("sse", "U30"), ] arr = np.zeros(len(ib_ctype), dtype=dtype) for idx, ctype in enumerate(ib_ctype): arr["ib"][idx] = bt.int_format(idx + 1) arr["k"][idx] = bt.int_format(ib_cellid[idx][0] + 1) arr["ssv"][idx] = bt.exp_format(parameters[name]["ssv"][idx]) arr["sse"][idx] = bt.exp_format(parameters[name]["sse"][idx]) if not silent: print("creating...'{}'".format(fpth)) col_widths = (0.2, 0.2, 0.2, 0.2) bt.build_table( caption, fpth, arr, headings=headings, col_widths=col_widths, ) return