+
diff --git a/notebooks/3D_simulation.ipynb b/notebooks/3D_simulation.ipynb
index 5d71c12..78e40c5 100644
--- a/notebooks/3D_simulation.ipynb
+++ b/notebooks/3D_simulation.ipynb
@@ -149,7 +149,9 @@
"outputs": [],
"source": [
"for i in model.builder.iter_node():\n",
- " i.rift=[[182,175]]"
+ " i.rift=[[182,175]]\n",
+ " i.bflux = False\n",
+ " i.adiab = 0.3e-3 "
]
},
{
@@ -158,7 +160,9 @@
"metadata": {},
"outputs": [],
"source": [
- "model.simulator.simulate_every = 1"
+ "model.simulator.simulate_every = 1\n",
+ "model.parameters.HPdcr = 1e11 # set to \"infinite\" decay\n",
+ "model.parameters.bflux = False # set to \"infinite\" decay\n"
]
},
{
@@ -230,7 +234,7 @@
"# interpolate and extrapolate the missing nodes\n",
"# find nearby nodes from the array indexer_full_sim, which is sorted by x-index\n",
"import itertools\n",
- "from subsheat3D.fixed_mesh_model import interpolateNode\n",
+ "from warmth3d.fixed_mesh_model import interpolateNode\n",
"for ni in range(len(model.builder.nodes)):\n",
" for nj in range(len(model.builder.nodes[ni])):\n",
" if model.builder.nodes[ni][nj] is False:\n",
@@ -269,8 +273,7 @@
" model.builder.nodes[ni][nj] = node\n",
" # if (model.builder.nodes[ni][nj].Y>12000):\n",
" # model.builder.nodes[ni][nj].crustRHP = (2e-6) * 4\n",
- " # model.builder.nodes[ni][nj].rhp = (2e-6) * 4\n",
- " model.builder.nodes[ni][nj].crustRHP = (2e0) * 1\n"
+ " # model.builder.nodes[ni][nj].rhp = (2e-6) * 4\n"
]
},
{
@@ -308,7 +311,7 @@
" os.mkdir('temp')\n",
"except FileExistsError:\n",
" pass\n",
- "run(model,start_time=model.parameters.time_start,end_time=0)\n",
+ "mm2 = run(model,start_time=model.parameters.time_start,end_time=0)\n",
"\n"
]
},
@@ -317,12 +320,51 @@
"execution_count": null,
"metadata": {},
"outputs": [],
- "source": []
+ "source": [
+ "import numpy as np\n",
+ "hx = model.builder.grid.num_nodes_x // 2\n",
+ "hy = model.builder.grid.num_nodes_y // 2\n",
+ "# hx = model.builder.grid.num_nodes_x - 1 - pad\n",
+ "# hy = model.builder.grid.num_nodes_y - 1 - pad\n",
+ "\n",
+ "nn = model.builder.nodes[hy][hx]\n",
+ "dd = nn._depth_out[:,0]\n",
+ "\n",
+ "node_ind = hy*model.builder.grid.num_nodes_x + hx\n",
+ "v_per_n = int(mm2.mesh_vertices.shape[0]/(model.builder.grid.num_nodes_y*model.builder.grid.num_nodes_x))\n",
+ "\n",
+ "temp_1d = np.nan_to_num(nn.temperature_out[:,0], nan=5.0)\n",
+ "temp_3d_ind = np.array([ np.where([mm2.mesh_reindex==i])[1][0] for i in range(node_ind*v_per_n, (node_ind+1)*v_per_n) ] )\n",
+ "dd_mesh = mm2.mesh.geometry.x[temp_3d_ind,2]\n",
+ "temp_3d_mesh = mm2.u_n.x.array[temp_3d_ind]\n",
+ "\n",
+ "temp_1d_at_mesh_pos = np.interp(dd_mesh, dd, temp_1d)\n",
+ "dd_subset = np.where(dd_mesh<5000)\n",
+ "print(f'Max. absolute error in temperature at 3D mesh vertex positions: {np.amax(np.abs(temp_1d_at_mesh_pos-temp_3d_mesh))}')\n",
+ "print(f'Max. absolute error at depths < 5000m: {np.amax(np.abs(temp_1d_at_mesh_pos[dd_subset]-temp_3d_mesh[dd_subset]))}')\n"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "%matplotlib\n",
+ "import matplotlib as plt\n",
+ "plt.use('TkAgg')\n",
+ "\n",
+ "plt.pyplot.plot( dd, temp_1d, label=f'1D simulation'); \n",
+ "# plt.pyplot.plot( dd, temp_3d, label=f'3D simulation'); \n",
+ "plt.pyplot.plot( dd_mesh, temp_3d_mesh, 'o', label=f'3D simulation (nodes)'); \n",
+ "plt.pyplot.legend(loc=\"lower right\", prop={'size': 7})\n",
+ "plt.pyplot.show()\n"
+ ]
}
],
"metadata": {
"kernelspec": {
- "display_name": ".venv",
+ "display_name": "Python 3.9.10 64-bit",
"language": "python",
"name": "python3"
},
@@ -336,7 +378,12 @@
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
- "version": "3.10.13"
+ "version": "3.9.10"
+ },
+ "vscode": {
+ "interpreter": {
+ "hash": "e51a2b4e99f2ffc4c1b139fefabaf83a9a6f27184f9cac248baef6121453d47a"
+ }
}
},
"nbformat": 4,
diff --git a/searchindex.js b/searchindex.js
index fdbfce5..2dc47e2 100644
--- a/searchindex.js
+++ b/searchindex.js
@@ -1 +1 @@
-Search.setIndex({"docnames": ["apireference", "example", "index", "notebooks/3D_simulation", "notebooks/Build_within_Python", "theory"], "filenames": ["apireference.rst", "example.rst", "index.rst", "notebooks/3D_simulation.ipynb", "notebooks/Build_within_Python.ipynb", "theory.md"], "titles": ["API Reference", "Example", "Welcome to warmth\u2019s documentation!", "<no title>", "Build 1D model", "Background theory"], "terms": {"class": [0, 5], "warmth": [0, 3, 4, 5], "main": 0, "properti": [0, 5], "return": 0, "type": 0, "load": 0, "path": [0, 3, 4], "str": 0, "overrid": 0, "current": [0, 3, 5], "data": [0, 3, 4, 5], "build": [0, 1, 2, 5], "single_nod": [0, 4], "1d": [0, 1, 2, 3], "locat": 0, "hc": [0, 3], "initi": 0, "crustal": [0, 5], "thick": [0, 4, 5], "m": [0, 3, 4], "float": 0, "hlith": [0, 3], "depth": [0, 4], "base": [0, 4, 5], "lithospher": [0, 5], "thermal": 0, "asthenospher": 0, "boundari": 0, "kcrust": [0, 3], "conduct": [0, 4], "crust": [0, 5], "20c": 0, "w": 0, "k": 0, "2": [0, 3, 4], "klith": [0, 3], "mantl": 0, "kasth": [0, 3], "rhp": [0, 3, 4], "radiogen": 0, "heat": [0, 2, 4], "product": [0, 3], "3": [0, 3, 4], "crustsolid": [0, 3], "densiti": 0, "lithsolid": [0, 3], "asthsolid": [0, 3], "th": 0, "t0": [0, 3], "seab": 0, "temperatur": [0, 4], "c": [0, 4], "tm": [0, 3], "lab": 0, "qbase": [0, 3, 4], "flow": [0, 4], "moho": 0, "sediments_input": [0, 4], "present": [0, 5], "dai": [0, 5], "sediment": [0, 4, 5], "see": [0, 5], "single_node_sediments_inputs_templ": [0, 4], "pd": [0, 4], "datafram": [0, 4], "x": [0, 3, 4], "node": [0, 2, 3, 4], "y": [0, 3], "paleowd": [0, 4], "paleo": 0, "water": 0, "multi": 0, "rift": [0, 3, 4], "np": [0, 4], "ndarrai": 0, "float64": 0, "episod": 0, "list": [0, 3], "int": [0, 4], "water_depth_differ": 0, "differ": 0, "between": 0, "observ": 0, "sediment_fill_margin": 0, "maximum": 0, "when": [0, 4], "fit": 0, "i": [0, 3, 4, 5], "consid": 0, "achiev": 0, "total_beta_test": 0, "total": [0, 4], "number": 0, "beta": [0, 4, 5], "factor": [0, 5], "test": [0, 5], "error": 0, "from": [0, 3, 4, 5], "none": [0, 3], "simulated_at": 0, "timestamp": 0, "finish": 0, "bool": 0, "whether": 0, "found": 0, "true": [0, 3, 4], "smaller": 0, "than": 0, "accept": 0, "clean": 0, "up": 0, "input": [0, 3, 5], "grid": [0, 2, 3], "origin_x": 0, "origin_i": 0, "num_nodes_x": 0, "num_nodes_i": 0, "step_x": 0, "step_i": 0, "defin": [0, 3, 5], "geometri": 0, "3d": [0, 2, 3], "location_grid": [0, 3], "all": 0, "A": 0, "2d": 0, "arrai": [0, 3, 4], "make_grid_arr": 0, "templat": 0, "store": 0, "object": [0, 4], "indexing_arr": 0, "indic": 0, "arri": 0, "dump": 0, "filepath": 0, "save": [0, 3], "file": [0, 5], "util": 0, "creat": [0, 4], "singl": [0, 4], "input_horizons_templ": [0, 3], "us": [0, 2, 3], "map": 0, "emtpi": 0, "append": [0, 4], "extract_nod": [0, 3], "thread": 0, "formatfil": 0, "irap_binari": [0, 3], "facies_dict": 0, "dict": 0, "extract": 0, "concurr": 0, "process": 0, "directori": 0, "option": 0, "format": [0, 5], "support": 0, "xtgeo": [0, 5], "default": 0, "lithologi": 0, "valu": [0, 4, 5], "faci": 0, "rais": 0, "valueerror": 0, "invalid": [0, 4], "tabl": [0, 3], "check": 0, "self": 0, "input_horizon": [0, 3], "sort": [0, 3], "ascend": 0, "ag": [0, 4], "define_geometri": [0, 3], "xinc": [0, 3], "yinc": [0, 3], "fformat": [0, 3], "overwrit": 0, "distanc": 0, "direct": 0, "iter_nod": [0, 3], "iter": 0, "yield": 0, "set_eustatic_sea_level": [0, 3, 4], "sealevel": 0, "set": [0, 3], "eustat": 0, "sea": 0, "level": 0, "correct": 0, "subsid": [0, 5], "initial_hlith_max": 0, "allow": 0, "calibr": 0, "initial_hlith_min": 0, "minimum": 0, "time_start": [0, 3, 4], "start": 0, "ma": [0, 4], "time_end": [0, 4], "end": [0, 3], "positive_down": 0, "ar": [0, 5], "posit": [0, 5], "downward": 0, "positv": 0, "time_step_ma": 0, "time": [0, 3, 4], "step": [0, 4], "solv": [0, 5], "1": [0, 3, 4], "forward_model": [0, 4], "current_nod": 0, "simulate_single_nod": 0, "simulate_continent": 0, "run": [0, 3, 4, 5], "static": 0, "compact": 0, "top_m": 0, "base_m": 0, "phi0": 0, "phi_decai": 0, "base_maximum_burial_depth_m": 0, "0": [0, 3, 4], "top": [0, 4], "poros": 0, "surfac": [0, 3], "fraction": 0, "exponenti": 0, "decai": [0, 4], "burial": 0, "To": [0, 5], "make": [0, 3], "sure": [0, 3], "irrevers": 0, "base_result": 0, "decompact": 0, "without": 0, "pore": 0, "space": 0, "For": [0, 5], "calcul": 0, "deposit": 0, "rate": 0, "grain_thick": [0, 4], "add_sedi": 0, "sedrat": 0, "sed": 0, "xsed_old": 0, "tsed_old": 0, "hpsed_old": 0, "idsed_old": 0, "int32": 0, "tupl": 0, "take": 0, "care": 0, "thi": [0, 5], "unit": 0, "referenc": 0, "sedimentari": 0, "column": 0, "xsed": 0, "m3": 0, "id": 0, "sedflag": 0, "exist": [0, 3], "tsed": 0, "hpsed": 0, "ids": 0, "calculate_new_temperatur": 0, "coord_crust_lith": 0, "t_old": 0, "hp": 0, "new": [0, 4], "profil": 0, "If": 0, "coord": 0, "vertic": 0, "mesh": [0, 3, 5], "previou": [0, 3], "_description_": 0, "implicit_euler_solv": 0, "t_start": 0, "coord_al": 0, "density_al": 0, "conductivity_pack": 0, "sourc": [0, 5], "k_last_nod": 0, "hp_last_nod": 0, "equat": [0, 2], "backward": 0, "euler": 0, "scheme": 0, "condit": 0, "whole": 0, "term": 0, "postprocess": [0, 4], "sediments_id": 0, "sediment_input": 0, "k_crust": 0, "k_lith": 0, "k_asth": 0, "typeddict": 0, "resultvalu": 0, "requir": [0, 5], "kei": 0, "layerid": 0, "top_crust": 0, "geolog": 0, "top_lithospher": 0, "top_asthenospher": 0, "crust_thick": 0, "lithosphere_thick": 0, "cell": [0, 3], "sediment_id": 0, "filter": 0, "layer": [0, 5], "centr": 0, "center": 0, "sediment_poros": 0, "effective_conduct": [0, 4], "effect": 0, "heatflow": 0, "basement_heatflow": [0, 4], "basement": [0, 4], "model": [1, 2, 3, 5], "background": 2, "theori": 2, "pre": 2, "requisit": 2, "simul": [2, 3, 4], "dolfinx": 2, "resqml": 2, "output": 2, "depend": 2, "exampl": 2, "api": 2, "refer": 2, "builder": [2, 3, 4], "paramet": [2, 3, 4], "forward": 2, "result": [2, 4, 5], "index": [2, 3, 4], "modul": 2, "search": 2, "page": 2, "download": [3, 4], "notebook": [3, 4], "import": [3, 4], "pathlib": 3, "1e": 3, "6": [3, 4], "30e3": 3, "03": 3, "30e": [3, 4], "maps_dir": 3, "add": [3, 4], "you": 3, "can": [3, 5], "other": [3, 5], "method": 3, "well": 3, "loc": 3, "gri": [3, 5], "onlap": 3, "66": [3, 4], "100": [3, 4], "163": 3, "eros": 3, "4": [3, 4], "168": 3, "5": [3, 4], "170": 3, "182": 3, "inc": 3, "2000": 3, "__dict__": 3, "haq87": [3, 4], "175": 3, "simulate_everi": 3, "n_valid_nod": 3, "print": 3, "nameerror": 3, "traceback": 3, "most": 3, "recent": 3, "call": 3, "last": 3, "workspac": [3, 4], "doc": 3, "3d_simul": 3, "ipynb": 3, "14": 3, "line": 3, "gt": [3, 4], "lt": [3, 4], "href": 3, "39": 3, "vscode": 3, "codespac": 3, "2bobscur": 3, "meme": 3, "qpw4jj4pgvp34j7v": 3, "x45sdnnjb2rllxjlbw90zq": 3, "name": 3, "purg": 3, "interpol": [3, 5], "extrapol": 3, "miss": 3, "find": 3, "nearbi": 3, "indexer_full_sim": 3, "which": [3, 5], "itertool": 3, "subsheat3d": 3, "fixed_mesh_model": 3, "interpolatenod": 3, "ni": 3, "rang": 3, "len": 3, "nj": 3, "fals": [3, 4], "closest_x_up": 3, "j": 3, "matching_x": 3, "break": 3, "closest_x_down": 3, "closest_y_up": 3, "matching_i": 3, "closest_y_down": 3, "interpolationnod": 3, "nn": 3, "els": 3, "itself": 3, "same": 3, "member": 3, "variabl": 3, "12000": 3, "crustrhp": [3, 4], "2e": [3, 4], "2e0": 3, "15": 3, "mesh_model": 3, "o": [3, 5], "try": 3, "mkdir": 3, "except": 3, "fileexistserror": 3, "pass": [3, 4], "temp": 3, "start_tim": 3, "end_tim": 3, "rebuild": 3, "reload": 3, "tti": 3, "nodeparameters1d": 3, "shf": 3, "30000": 3, "hw": 3, "3600": 3, "130000": 3, "109": 3, "crustliquid": 3, "2500": 3, "2800": 3, "lithliquid": 3, "2700": 3, "3300": 3, "asthliquid": 3, "3200": 3, "1330": 3, "buil": 3, "buildvertic": 3, "The": [3, 5], "kernel": 3, "crash": 3, "while": 3, "execut": 3, "code": 3, "pleas": 3, "review": 3, "": 3, "identifi": 3, "possibl": 3, "caus": 3, "failur": 3, "click": 3, "http": 3, "aka": 3, "vscodejupyterkernelcrash": 3, "here": 3, "more": 3, "info": 3, "view": 3, "jupyt": 3, "command": 3, "viewoutput": 3, "log": 3, "further": [3, 5], "detail": 3, "sy": 4, "numpi": 4, "panda": 4, "node_templ": 4, "h1": 4, "152": 4, "500000": 4, "301755e": 4, "09": 4, "620000": 4, "500": 4, "2720": 4, "2448": 4, "h2": 4, "810": 4, "20": 4, "538462": 4, "079433e": 4, "599730": 4, "490": 4, "2708": 4, "2437": 4, "h3": 4, "1608": 4, "h4": 4, "1973": 4, "h5": [4, 5], "2262": 4, "145": 4, "h6": 4, "2362": 4, "904762": 4, "719506e": 4, "10": 4, "447705": 4, "415": 4, "2618": 4, "2356": 4, "h7": 4, "2427": 4, "160": 4, "horizon": [4, 5], "from_dict": 4, "topag": 4, "k_cond": 4, "phi": 4, "solidu": 4, "liquidu": 4, "7": 4, "concat": 4, "ignore_index": 4, "py": [4, 5], "199": 4, "futurewarn": 4, "In": 4, "futur": 4, "version": [4, 5], "constructor": 4, "infer": 4, "numer": 4, "dtype": 4, "sequenc": 4, "match": [4, 5], "seri": [4, 5], "behavior": 4, "check_ascend": 4, "df": 4, "appli": 4, "lambda": 4, "is_monotonic_increas": 4, "baseag": 4, "phi_mean": 4, "658": 4, "310357": 4, "528332": 4, "798": 4, "511062": 4, "359571": 4, "200000": 4, "365": 4, "332780": 4, "088275": 4, "289": 4, "221878": 4, "232256": 4, "078943": 4, "210571": 4, "65": 4, "053525": 4, "176536": 4, "iloc": 4, "200": 4, "onli": [4, 5], "mulit": 4, "8": 4, "9": 4, "parallel": [4, 5], "cpu": 4, "user": 4, "612": 4, "n": 4, "wall": 4, "624": 4, "24": 4, "11": 4, "matplotlib": 4, "inlin": 4, "pyplot": 4, "plt": 4, "plot": 4, "1000": 4, "gca": 4, "invert_xaxi": 4, "xlabel": 4, "ylabel": 4, "mw": 4, "m2": 4, "show": 4, "runtimewarn": 4, "encount": 4, "divid": 4, "v": 4, "initial_poro": 4, "initial_decai": 4, "phi1": 4, "12": 4, "ipywidget": 4, "widget": 4, "def": 4, "plot_cell_bound": 4, "t": 4, "invert_yaxi": 4, "ylim": 4, "bottom": 4, "interact": 4, "intslid": 4, "max": 4, "min": 4, "function": [4, 5], "__main__": 4, "17": 4, "plot_cell_properti": 4, "under": 5, "activ": 5, "develop": 5, "yhei": 5, "uniform": 5, "reactangular": 5, "nodegrid": 5, "structur": 5, "provid": 5, "read": 5, "sedimentstack": 5, "complet": 5, "some": 5, "script": 5, "1dsed": 5, "each": 5, "pickl": 5, "separ": 5, "improv": 5, "perform": 5, "regularli": 5, "updat": 5, "petsc": 5, "solver": 5, "packag": 5, "part": 5, "fenic": 5, "project": 5, "comput": 5, "built": 5, "hexahedra": 5, "everi": 5, "rectangl": 5, "e": 5, "aesthenospher": 5, "subdivid": 5, "tetrahedra": 5, "manag": 5, "uniformnodegridfixedsizemeshmodel": 5, "demonstr": 5, "warmth3d_mapa_exampl": 5, "note": 5, "definit": 5, "should": 5, "solut": 5, "write": 5, "1m": 5, "year": 5, "xdmf": 5, "visual": 5, "paraview": 5, "unstructur": 5, "pair": 5, "epc": 5, "resqpy_help": 5, "modifi": 5, "resqpi": 5, "librari": 5, "visualis": 5, "3rd": 5, "parti": 5, "plug": 5, "instal": 5, "fespp": 5, "linux": 5, "ha": 5, "compil": 5, "apt": 5, "get": 5, "pip": 5, "now": 5, "chang": 5, "yet": 5, "merg": 5, "meshio": 5, "ad": 5}, "objects": {"warmth": [[0, 0, 0, "-", "build"], [0, 0, 0, "-", "forward_modelling"], [0, 0, 0, "-", "model"], [0, 0, 0, "-", "parameters"], [0, 0, 0, "-", "postprocessing"], [0, 0, 0, "-", "simulator"]], "warmth.build": [[0, 1, 1, "", "Builder"], [0, 1, 1, "", "Grid"], [0, 1, 1, "", "single_node"]], "warmth.build.Builder": [[0, 2, 1, "", "define_geometry"], [0, 2, 1, "", "extract_nodes"], [0, 3, 1, "", "input_horizons_template"], [0, 2, 1, "", "iter_node"], [0, 3, 1, "", "locations"], [0, 2, 1, "", "set_eustatic_sea_level"], [0, 3, 1, "", "single_node_sediments_inputs_template"]], "warmth.build.Grid": [[0, 2, 1, "", "dump"], [0, 3, 1, "", "indexing_arr"], [0, 3, 1, "", "location_grid"], [0, 2, 1, "", "make_grid_arr"]], "warmth.build.single_node": [[0, 4, 1, "", "T0"], [0, 4, 1, "", "Tm"], [0, 4, 1, "", "X"], [0, 4, 1, "", "Y"], [0, 4, 1, "", "asthsolid"], [0, 4, 1, "", "crustsolid"], [0, 4, 1, "", "error"], [0, 3, 1, "", "fitting"], [0, 4, 1, "", "hLith"], [0, 4, 1, "", "hc"], [0, 4, 1, "", "kAsth"], [0, 4, 1, "", "kCrust"], [0, 4, 1, "", "kLith"], [0, 4, 1, "", "lithsolid"], [0, 4, 1, "", "paleoWD"], [0, 4, 1, "", "qbase"], [0, 3, 1, "", "result"], [0, 4, 1, "", "rhp"], [0, 4, 1, "", "rift"], [0, 4, 1, "", "sediment_fill_margin"], [0, 3, 1, "", "sediments"], [0, 4, 1, "", "sediments_inputs"], [0, 4, 1, "", "simulated_at"], [0, 4, 1, "", "total_beta_tested"], [0, 4, 1, "", "water_depth_difference"]], "warmth.forward_modelling": [[0, 1, 1, "", "Forward_model"]], "warmth.forward_modelling.Forward_model": [[0, 2, 1, "", "add_sediments"], [0, 2, 1, "", "calculate_new_temperature"], [0, 2, 1, "", "compaction"], [0, 2, 1, "", "decompaction"], [0, 2, 1, "", "implicit_euler_solve"], [0, 2, 1, "", "simulate_continental"], [0, 2, 1, "", "simulate_single_node"]], "warmth.model": [[0, 1, 1, "", "Model"]], "warmth.model.Model": [[0, 3, 1, "", "builder"], [0, 2, 1, "", "load"], [0, 3, 1, "", "parameters"], [0, 3, 1, "", "simulator"]], "warmth.parameters": [[0, 1, 1, "", "Parameters"]], "warmth.parameters.Parameters": [[0, 3, 1, "", "initial_hLith_max"], [0, 3, 1, "", "initial_hLith_min"], [0, 3, 1, "", "positive_down"], [0, 3, 1, "", "time_end"], [0, 3, 1, "", "time_start"], [0, 3, 1, "", "time_step_Ma"]], "warmth.postprocessing": [[0, 1, 1, "", "Results"]], "warmth.postprocessing.Results": [[0, 3, 1, "", "ages"], [0, 2, 1, "", "basement_heatflow"], [0, 2, 1, "", "crust_thickness"], [0, 2, 1, "", "depth"], [0, 2, 1, "", "effective_conductivity"], [0, 2, 1, "", "heatflow"], [0, 2, 1, "", "lithosphere_thickness"], [0, 5, 1, "", "resultValues"], [0, 2, 1, "", "sediment_ids"], [0, 2, 1, "", "sediment_porosity"], [0, 2, 1, "", "temperature"], [0, 2, 1, "", "top_asthenosphere"], [0, 2, 1, "", "top_crust"], [0, 2, 1, "", "top_lithosphere"]], "warmth.simulator": [[0, 1, 1, "", "Simulator"]]}, "objtypes": {"0": "py:module", "1": "py:class", "2": "py:method", "3": "py:property", "4": "py:attribute", "5": "py:typeddict"}, "objnames": {"0": ["py", "module", "Python module"], "1": ["py", "class", "Python class"], "2": ["py", "method", "Python method"], "3": ["py", "property", "Python property"], "4": ["py", "attribute", "Python attribute"], "5": ["py", "typeddict", "Python typeddict"]}, "titleterms": {"api": 0, "refer": 0, "model": [0, 4], "builder": 0, "paramet": 0, "forward": 0, "simul": [0, 5], "result": 0, "exampl": 1, "content": [1, 2], "welcom": 2, "warmth": 2, "": 2, "document": 2, "indic": 2, "tabl": 2, "build": 4, "1d": [4, 5], "background": 5, "theori": 5, "pre": 5, "requisit": 5, "grid": 5, "node": 5, "3d": 5, "heat": 5, "equat": 5, "us": 5, "dolfinx": 5, "resqml": 5, "output": 5, "depend": 5}, "envversion": {"sphinx.domains.c": 3, "sphinx.domains.changeset": 1, "sphinx.domains.citation": 1, "sphinx.domains.cpp": 9, "sphinx.domains.index": 1, "sphinx.domains.javascript": 3, "sphinx.domains.math": 2, "sphinx.domains.python": 4, "sphinx.domains.rst": 2, "sphinx.domains.std": 2, "nbsphinx": 4, "sphinx": 60}, "alltitles": {"API Reference": [[0, "api-reference"]], "Model": [[0, "module-warmth.model"]], "Model builder": [[0, "module-warmth.build"]], "Parameters": [[0, "module-warmth.parameters"]], "Forward model": [[0, "module-warmth.forward_modelling"]], "Simulator": [[0, "module-warmth.simulator"]], "Results": [[0, "module-warmth.postprocessing"]], "Example": [[1, "example"]], "Contents:": [[1, null], [2, null]], "Welcome to warmth\u2019s documentation!": [[2, "welcome-to-warmth-s-documentation"]], "Indices and tables": [[2, "indices-and-tables"]], "Build 1D model": [[4, "Build-1D-model"]], "Background theory": [[5, "background-theory"]], "Pre-requisite: grid of 1D node simulation": [[5, "pre-requisite-grid-of-1d-node-simulation"]], "3D heat equation simulation using dolfinx": [[5, "d-heat-equation-simulation-using-dolfinx"]], "RESQML output": [[5, "resqml-output"]], "3D dependencies": [[5, "d-dependencies"]]}, "indexentries": {"builder (class in warmth.build)": [[0, "warmth.build.Builder"]], "forward_model (class in warmth.forward_modelling)": [[0, "warmth.forward_modelling.Forward_model"]], "grid (class in warmth.build)": [[0, "warmth.build.Grid"]], "model (class in warmth.model)": [[0, "warmth.model.Model"]], "parameters (class in warmth.parameters)": [[0, "warmth.parameters.Parameters"]], "results (class in warmth.postprocessing)": [[0, "warmth.postprocessing.Results"]], "results.resultvalues (typeddict in warmth.postprocessing)": [[0, "warmth.postprocessing.Results.resultValues"]], "simulator (class in warmth.simulator)": [[0, "warmth.simulator.Simulator"]], "t0 (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.T0"]], "tm (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.Tm"]], "x (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.X"]], "y (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.Y"]], "add_sediments() (warmth.forward_modelling.forward_model method)": [[0, "warmth.forward_modelling.Forward_model.add_sediments"]], "ages (warmth.postprocessing.results property)": [[0, "warmth.postprocessing.Results.ages"]], "asthsolid (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.asthsolid"]], "basement_heatflow() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.basement_heatflow"]], "builder (warmth.model.model property)": [[0, "warmth.model.Model.builder"]], "calculate_new_temperature() (warmth.forward_modelling.forward_model method)": [[0, "warmth.forward_modelling.Forward_model.calculate_new_temperature"]], "compaction() (warmth.forward_modelling.forward_model static method)": [[0, "warmth.forward_modelling.Forward_model.compaction"]], "crust_thickness() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.crust_thickness"]], "crustsolid (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.crustsolid"]], "decompaction() (warmth.forward_modelling.forward_model static method)": [[0, "warmth.forward_modelling.Forward_model.decompaction"]], "define_geometry() (warmth.build.builder method)": [[0, "warmth.build.Builder.define_geometry"]], "depth() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.depth"]], "dump() (warmth.build.grid method)": [[0, "warmth.build.Grid.dump"]], "effective_conductivity() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.effective_conductivity"]], "error (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.error"]], "extract_nodes() (warmth.build.builder method)": [[0, "warmth.build.Builder.extract_nodes"]], "fitting (warmth.build.single_node property)": [[0, "warmth.build.single_node.fitting"]], "hlith (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.hLith"]], "hc (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.hc"]], "heatflow() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.heatflow"]], "implicit_euler_solve() (warmth.forward_modelling.forward_model method)": [[0, "warmth.forward_modelling.Forward_model.implicit_euler_solve"]], "indexing_arr (warmth.build.grid property)": [[0, "warmth.build.Grid.indexing_arr"]], "initial_hlith_max (warmth.parameters.parameters property)": [[0, "warmth.parameters.Parameters.initial_hLith_max"]], "initial_hlith_min (warmth.parameters.parameters property)": [[0, "warmth.parameters.Parameters.initial_hLith_min"]], "input_horizons_template (warmth.build.builder property)": [[0, "warmth.build.Builder.input_horizons_template"]], "iter_node() (warmth.build.builder method)": [[0, "warmth.build.Builder.iter_node"]], "kasth (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.kAsth"]], "kcrust (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.kCrust"]], "klith (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.kLith"]], "lithosphere_thickness() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.lithosphere_thickness"]], "lithsolid (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.lithsolid"]], "load() (warmth.model.model method)": [[0, "warmth.model.Model.load"]], "location_grid (warmth.build.grid property)": [[0, "warmth.build.Grid.location_grid"]], "locations (warmth.build.builder property)": [[0, "warmth.build.Builder.locations"]], "make_grid_arr() (warmth.build.grid method)": [[0, "warmth.build.Grid.make_grid_arr"]], "module": [[0, "module-warmth.build"], [0, "module-warmth.forward_modelling"], [0, "module-warmth.model"], [0, "module-warmth.parameters"], [0, "module-warmth.postprocessing"], [0, "module-warmth.simulator"]], "paleowd (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.paleoWD"]], "parameters (warmth.model.model property)": [[0, "warmth.model.Model.parameters"]], "positive_down (warmth.parameters.parameters property)": [[0, "warmth.parameters.Parameters.positive_down"]], "qbase (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.qbase"]], "result (warmth.build.single_node property)": [[0, "warmth.build.single_node.result"]], "rhp (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.rhp"]], "rift (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.rift"]], "sediment_fill_margin (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.sediment_fill_margin"]], "sediment_ids() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.sediment_ids"]], "sediment_porosity() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.sediment_porosity"]], "sediments (warmth.build.single_node property)": [[0, "warmth.build.single_node.sediments"]], "sediments_inputs (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.sediments_inputs"]], "set_eustatic_sea_level() (warmth.build.builder method)": [[0, "warmth.build.Builder.set_eustatic_sea_level"]], "simulate_continental() (warmth.forward_modelling.forward_model method)": [[0, "warmth.forward_modelling.Forward_model.simulate_continental"]], "simulate_single_node() (warmth.forward_modelling.forward_model method)": [[0, "warmth.forward_modelling.Forward_model.simulate_single_node"]], "simulated_at (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.simulated_at"]], "simulator (warmth.model.model property)": [[0, "warmth.model.Model.simulator"]], "single_node (class in warmth.build)": [[0, "warmth.build.single_node"]], "single_node_sediments_inputs_template (warmth.build.builder property)": [[0, "warmth.build.Builder.single_node_sediments_inputs_template"]], "temperature() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.temperature"]], "time_end (warmth.parameters.parameters property)": [[0, "warmth.parameters.Parameters.time_end"]], "time_start (warmth.parameters.parameters property)": [[0, "warmth.parameters.Parameters.time_start"]], "time_step_ma (warmth.parameters.parameters property)": [[0, "warmth.parameters.Parameters.time_step_Ma"]], "top_asthenosphere() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.top_asthenosphere"]], "top_crust() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.top_crust"]], "top_lithosphere() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.top_lithosphere"]], "total_beta_tested (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.total_beta_tested"]], "warmth.build": [[0, "module-warmth.build"]], "warmth.forward_modelling": [[0, "module-warmth.forward_modelling"]], "warmth.model": [[0, "module-warmth.model"]], "warmth.parameters": [[0, "module-warmth.parameters"]], "warmth.postprocessing": [[0, "module-warmth.postprocessing"]], "warmth.simulator": [[0, "module-warmth.simulator"]], "water_depth_difference (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.water_depth_difference"]]}})
\ No newline at end of file
+Search.setIndex({"docnames": ["apireference", "example", "index", "notebooks/3D_simulation", "notebooks/Build_within_Python", "theory"], "filenames": ["apireference.rst", "example.rst", "index.rst", "notebooks/3D_simulation.ipynb", "notebooks/Build_within_Python.ipynb", "theory.md"], "titles": ["API Reference", "Example", "Welcome to warmth\u2019s documentation!", "<no title>", "Build 1D model", "Background theory"], "terms": {"class": [0, 5], "warmth": [0, 3, 4, 5], "main": 0, "properti": [0, 5], "return": 0, "type": 0, "load": 0, "path": [0, 3, 4], "str": 0, "overrid": 0, "current": [0, 3, 5], "data": [0, 3, 4, 5], "build": [0, 1, 2, 5], "single_nod": [0, 4], "1d": [0, 1, 2, 3], "locat": 0, "hc": [0, 3], "initi": 0, "crustal": [0, 5], "thick": [0, 4, 5], "m": [0, 3, 4], "float": 0, "hlith": [0, 3], "depth": [0, 3, 4], "base": [0, 4, 5], "lithospher": [0, 5], "thermal": 0, "asthenospher": 0, "boundari": 0, "kcrust": [0, 3], "conduct": [0, 4], "crust": [0, 5], "20c": 0, "w": 0, "k": 0, "2": [0, 3, 4], "klith": [0, 3], "mantl": 0, "kasth": [0, 3], "rhp": [0, 3, 4], "radiogen": 0, "heat": [0, 2, 4], "product": [0, 3], "3": [0, 3, 4], "crustsolid": [0, 3], "densiti": 0, "lithsolid": [0, 3], "asthsolid": [0, 3], "th": 0, "t0": [0, 3], "seab": 0, "temperatur": [0, 3, 4], "c": [0, 4], "tm": [0, 3], "lab": 0, "qbase": [0, 3, 4], "flow": [0, 4], "moho": 0, "sediments_input": [0, 4], "present": [0, 5], "dai": [0, 5], "sediment": [0, 4, 5], "see": [0, 5], "single_node_sediments_inputs_templ": [0, 4], "pd": [0, 4], "datafram": [0, 4], "x": [0, 3, 4], "node": [0, 2, 3, 4], "y": [0, 3], "paleowd": [0, 4], "paleo": 0, "water": 0, "multi": 0, "rift": [0, 3, 4], "np": [0, 3, 4], "ndarrai": 0, "float64": 0, "episod": 0, "list": [0, 3], "int": [0, 3, 4], "water_depth_differ": 0, "differ": 0, "between": 0, "observ": 0, "sediment_fill_margin": 0, "maximum": 0, "when": [0, 4], "fit": 0, "i": [0, 3, 4, 5], "consid": 0, "achiev": 0, "total_beta_test": 0, "total": [0, 4], "number": 0, "beta": [0, 4, 5], "factor": [0, 5], "test": [0, 5], "error": [0, 3], "from": [0, 3, 4, 5], "none": [0, 3], "simulated_at": 0, "timestamp": 0, "finish": 0, "bool": 0, "whether": 0, "found": 0, "true": [0, 3, 4], "smaller": 0, "than": 0, "accept": 0, "clean": 0, "up": 0, "input": [0, 3, 5], "grid": [0, 2, 3], "origin_x": 0, "origin_i": 0, "num_nodes_x": [0, 3], "num_nodes_i": [0, 3], "step_x": 0, "step_i": 0, "defin": [0, 3, 5], "geometri": [0, 3], "3d": [0, 2, 3], "location_grid": [0, 3], "all": 0, "A": 0, "2d": 0, "arrai": [0, 3, 4], "make_grid_arr": 0, "templat": 0, "store": 0, "object": [0, 4], "indexing_arr": 0, "indic": 0, "arri": 0, "dump": 0, "filepath": 0, "save": [0, 3], "file": [0, 5], "util": 0, "creat": [0, 4], "singl": [0, 4], "input_horizons_templ": [0, 3], "us": [0, 2, 3], "map": 0, "emtpi": 0, "append": [0, 4], "extract_nod": [0, 3], "thread": 0, "formatfil": 0, "irap_binari": [0, 3], "facies_dict": 0, "dict": 0, "extract": 0, "concurr": 0, "process": 0, "directori": 0, "option": 0, "format": [0, 5], "support": 0, "xtgeo": [0, 5], "default": 0, "lithologi": 0, "valu": [0, 4, 5], "faci": 0, "rais": 0, "valueerror": 0, "invalid": [0, 4], "tabl": [0, 3], "check": 0, "self": 0, "input_horizon": [0, 3], "sort": [0, 3], "ascend": 0, "ag": [0, 4], "define_geometri": [0, 3], "xinc": [0, 3], "yinc": [0, 3], "fformat": [0, 3], "overwrit": 0, "distanc": 0, "direct": 0, "iter_nod": [0, 3], "iter": 0, "yield": 0, "set_eustatic_sea_level": [0, 3, 4], "sealevel": 0, "set": [0, 3], "eustat": 0, "sea": 0, "level": 0, "correct": 0, "subsid": [0, 5], "initial_hlith_max": 0, "allow": 0, "calibr": 0, "initial_hlith_min": 0, "minimum": 0, "time_start": [0, 3, 4], "start": 0, "ma": [0, 4], "time_end": [0, 4], "end": [0, 3], "positive_down": 0, "ar": [0, 5], "posit": [0, 3, 5], "downward": 0, "positv": 0, "time_step_ma": 0, "time": [0, 3, 4], "step": [0, 4], "solv": [0, 5], "1": [0, 3, 4], "forward_model": [0, 4], "current_nod": 0, "simulate_single_nod": 0, "simulate_continent": 0, "run": [0, 3, 4, 5], "static": 0, "compact": 0, "top_m": 0, "base_m": 0, "phi0": 0, "phi_decai": 0, "base_maximum_burial_depth_m": 0, "0": [0, 3, 4], "top": [0, 4], "poros": 0, "surfac": [0, 3], "fraction": 0, "exponenti": 0, "decai": [0, 3, 4], "burial": 0, "To": [0, 5], "make": [0, 3], "sure": [0, 3], "irrevers": 0, "base_result": 0, "decompact": 0, "without": 0, "pore": 0, "space": 0, "For": [0, 5], "calcul": 0, "deposit": 0, "rate": 0, "grain_thick": [0, 4], "add_sedi": 0, "sedrat": 0, "sed": 0, "xsed_old": 0, "tsed_old": 0, "hpsed_old": 0, "idsed_old": 0, "int32": 0, "tupl": 0, "take": 0, "care": 0, "thi": [0, 5], "unit": 0, "referenc": 0, "sedimentari": 0, "column": 0, "xsed": 0, "m3": 0, "id": 0, "sedflag": 0, "exist": [0, 3], "tsed": 0, "hpsed": 0, "ids": 0, "calculate_new_temperatur": 0, "coord_crust_lith": 0, "t_old": 0, "hp": 0, "new": [0, 4], "profil": 0, "If": 0, "coord": 0, "vertic": 0, "mesh": [0, 3, 5], "previou": [0, 3], "_description_": 0, "implicit_euler_solv": 0, "t_start": 0, "coord_al": 0, "density_al": 0, "conductivity_pack": 0, "sourc": [0, 5], "k_last_nod": 0, "hp_last_nod": 0, "t_base": 0, "equat": [0, 2], "backward": 0, "euler": 0, "scheme": 0, "condit": 0, "whole": 0, "term": 0, "postprocess": [0, 4], "sediments_id": 0, "sediment_input": 0, "k_crust": 0, "k_lith": 0, "k_asth": 0, "typeddict": 0, "resultvalu": 0, "requir": [0, 5], "kei": 0, "layerid": 0, "top_crust": 0, "geolog": 0, "top_lithospher": 0, "top_asthenospher": 0, "crust_thick": 0, "lithosphere_thick": 0, "cell": [0, 3], "sediment_id": 0, "filter": 0, "layer": [0, 5], "centr": 0, "center": 0, "sediment_poros": 0, "effective_conduct": [0, 4], "effect": 0, "heatflow": 0, "basement_heatflow": [0, 4], "basement": [0, 4], "model": [1, 2, 3, 5], "background": 2, "theori": 2, "pre": 2, "requisit": 2, "simul": [2, 3, 4], "dolfinx": 2, "resqml": 2, "output": 2, "depend": 2, "exampl": 2, "api": 2, "refer": 2, "builder": [2, 3, 4], "paramet": [2, 3, 4], "forward": 2, "result": [2, 4, 5], "index": [2, 3, 4], "modul": 2, "search": 2, "page": 2, "download": [3, 4], "notebook": [3, 4], "import": [3, 4], "pathlib": 3, "1e": 3, "6": [3, 4], "30e3": 3, "03": 3, "30e": [3, 4], "maps_dir": 3, "add": [3, 4], "you": 3, "can": [3, 5], "other": [3, 5], "method": 3, "well": 3, "loc": 3, "gri": [3, 5], "onlap": 3, "66": [3, 4], "100": [3, 4], "163": 3, "eros": 3, "4": [3, 4], "168": 3, "5": [3, 4], "170": 3, "182": 3, "inc": 3, "2000": 3, "__dict__": 3, "haq87": [3, 4], "175": 3, "bflux": 3, "fals": [3, 4], "adiab": 3, "3e": 3, "simulate_everi": 3, "hpdcr": 3, "1e11": 3, "infinit": 3, "n_valid_nod": 3, "print": 3, "nameerror": 3, "traceback": 3, "most": 3, "recent": 3, "call": 3, "last": 3, "workspac": [3, 4], "doc": 3, "3d_simul": 3, "ipynb": 3, "14": 3, "line": 3, "gt": [3, 4], "lt": [3, 4], "href": 3, "39": 3, "vscode": 3, "codespac": 3, "2bobscur": 3, "meme": 3, "qpw4jj4pgvp34j7v": 3, "x45sdnnjb2rllxjlbw90zq": 3, "name": 3, "purg": 3, "interpol": [3, 5], "extrapol": 3, "miss": 3, "find": 3, "nearbi": 3, "indexer_full_sim": 3, "which": [3, 5], "itertool": 3, "warmth3d": 3, "fixed_mesh_model": 3, "interpolatenod": 3, "ni": 3, "rang": 3, "len": 3, "nj": 3, "closest_x_up": 3, "j": 3, "matching_x": 3, "break": 3, "closest_x_down": 3, "closest_y_up": 3, "matching_i": 3, "closest_y_down": 3, "interpolationnod": 3, "nn": 3, "els": 3, "itself": 3, "same": 3, "member": 3, "variabl": 3, "12000": 3, "crustrhp": [3, 4], "2e": [3, 4], "15": 3, "mesh_model": 3, "o": [3, 5], "try": 3, "mkdir": 3, "except": 3, "fileexistserror": 3, "pass": [3, 4], "temp": 3, "mm2": 3, "start_tim": 3, "end_tim": 3, "rebuild": 3, "reload": 3, "tti": 3, "nodeparameters1d": 3, "shf": 3, "30000": 3, "hw": 3, "3600": 3, "130000": 3, "109": 3, "crustliquid": 3, "2500": 3, "2800": 3, "lithliquid": 3, "2700": 3, "3300": 3, "asthliquid": 3, "3200": 3, "1330": 3, "buil": 3, "buildvertic": 3, "The": [3, 5], "kernel": 3, "crash": 3, "while": 3, "execut": 3, "code": 3, "pleas": 3, "review": 3, "": 3, "identifi": 3, "possibl": 3, "caus": 3, "failur": 3, "click": 3, "http": 3, "aka": 3, "vscodejupyterkernelcrash": 3, "here": 3, "more": 3, "info": 3, "view": 3, "jupyt": 3, "command": 3, "viewoutput": 3, "log": 3, "further": [3, 5], "detail": 3, "numpi": [3, 4], "hx": 3, "hy": 3, "pad": 3, "dd": 3, "_depth_out": 3, "node_ind": 3, "v_per_n": 3, "mesh_vertic": 3, "shape": 3, "temp_1d": 3, "nan_to_num": 3, "temperature_out": 3, "nan": 3, "temp_3d_ind": 3, "where": 3, "mesh_reindex": 3, "dd_mesh": 3, "temp_3d_mesh": 3, "u_n": 3, "temp_1d_at_mesh_po": 3, "interp": 3, "dd_subset": 3, "5000": 3, "f": 3, "max": [3, 4], "absolut": 3, "vertex": 3, "amax": 3, "ab": 3, "5000m": 3, "matplotlib": [3, 4], "plt": [3, 4], "tkagg": 3, "pyplot": [3, 4], "plot": [3, 4], "label": 3, "temp_3d": 3, "legend": 3, "lower": 3, "right": 3, "prop": 3, "size": 3, "7": [3, 4], "show": [3, 4], "sy": 4, "panda": 4, "node_templ": 4, "h1": 4, "152": 4, "500000": 4, "301755e": 4, "09": 4, "620000": 4, "500": 4, "2720": 4, "2448": 4, "h2": 4, "810": 4, "20": 4, "538462": 4, "079433e": 4, "599730": 4, "490": 4, "2708": 4, "2437": 4, "h3": 4, "1608": 4, "h4": 4, "1973": 4, "h5": [4, 5], "2262": 4, "145": 4, "h6": 4, "2362": 4, "904762": 4, "719506e": 4, "10": 4, "447705": 4, "415": 4, "2618": 4, "2356": 4, "h7": 4, "2427": 4, "160": 4, "horizon": [4, 5], "from_dict": 4, "topag": 4, "k_cond": 4, "phi": 4, "solidu": 4, "liquidu": 4, "concat": 4, "ignore_index": 4, "py": [4, 5], "199": 4, "futurewarn": 4, "In": 4, "futur": 4, "version": [4, 5], "constructor": 4, "infer": 4, "numer": 4, "dtype": 4, "sequenc": 4, "match": [4, 5], "seri": [4, 5], "behavior": 4, "check_ascend": 4, "df": 4, "appli": 4, "lambda": 4, "is_monotonic_increas": 4, "baseag": 4, "phi_mean": 4, "658": 4, "310357": 4, "528332": 4, "798": 4, "511062": 4, "359571": 4, "200000": 4, "365": 4, "332780": 4, "088275": 4, "289": 4, "221878": 4, "232256": 4, "078943": 4, "210571": 4, "65": 4, "053525": 4, "176536": 4, "iloc": 4, "200": 4, "onli": [4, 5], "mulit": 4, "8": 4, "9": 4, "parallel": [4, 5], "cpu": 4, "user": 4, "612": 4, "n": 4, "wall": 4, "624": 4, "24": 4, "11": 4, "inlin": 4, "1000": 4, "gca": 4, "invert_xaxi": 4, "xlabel": 4, "ylabel": 4, "mw": 4, "m2": 4, "runtimewarn": 4, "encount": 4, "divid": 4, "v": 4, "initial_poro": 4, "initial_decai": 4, "phi1": 4, "12": 4, "ipywidget": 4, "widget": 4, "def": 4, "plot_cell_bound": 4, "t": 4, "invert_yaxi": 4, "ylim": 4, "bottom": 4, "interact": 4, "intslid": 4, "min": 4, "function": [4, 5], "__main__": 4, "17": 4, "plot_cell_properti": 4, "under": 5, "activ": 5, "develop": 5, "yhei": 5, "uniform": 5, "reactangular": 5, "nodegrid": 5, "structur": 5, "provid": 5, "read": 5, "sedimentstack": 5, "complet": 5, "some": 5, "script": 5, "1dsed": 5, "each": 5, "pickl": 5, "separ": 5, "improv": 5, "perform": 5, "regularli": 5, "updat": 5, "petsc": 5, "solver": 5, "packag": 5, "part": 5, "fenic": 5, "project": 5, "comput": 5, "built": 5, "hexahedra": 5, "everi": 5, "rectangl": 5, "e": 5, "aesthenospher": 5, "subdivid": 5, "tetrahedra": 5, "manag": 5, "uniformnodegridfixedsizemeshmodel": 5, "demonstr": 5, "warmth3d_mapa_exampl": 5, "note": 5, "definit": 5, "should": 5, "solut": 5, "write": 5, "1m": 5, "year": 5, "xdmf": 5, "visual": 5, "paraview": 5, "unstructur": 5, "pair": 5, "epc": 5, "resqpy_help": 5, "modifi": 5, "resqpi": 5, "librari": 5, "visualis": 5, "3rd": 5, "parti": 5, "plug": 5, "instal": 5, "fespp": 5, "linux": 5, "ha": 5, "compil": 5, "apt": 5, "get": 5, "pip": 5, "now": 5, "chang": 5, "yet": 5, "merg": 5, "meshio": 5, "ad": 5}, "objects": {"warmth": [[0, 0, 0, "-", "build"], [0, 0, 0, "-", "forward_modelling"], [0, 0, 0, "-", "model"], [0, 0, 0, "-", "parameters"], [0, 0, 0, "-", "postprocessing"], [0, 0, 0, "-", "simulator"]], "warmth.build": [[0, 1, 1, "", "Builder"], [0, 1, 1, "", "Grid"], [0, 1, 1, "", "single_node"]], "warmth.build.Builder": [[0, 2, 1, "", "define_geometry"], [0, 2, 1, "", "extract_nodes"], [0, 3, 1, "", "input_horizons_template"], [0, 2, 1, "", "iter_node"], [0, 3, 1, "", "locations"], [0, 2, 1, "", "set_eustatic_sea_level"], [0, 3, 1, "", "single_node_sediments_inputs_template"]], "warmth.build.Grid": [[0, 2, 1, "", "dump"], [0, 3, 1, "", "indexing_arr"], [0, 3, 1, "", "location_grid"], [0, 2, 1, "", "make_grid_arr"]], "warmth.build.single_node": [[0, 4, 1, "", "T0"], [0, 4, 1, "", "Tm"], [0, 4, 1, "", "X"], [0, 4, 1, "", "Y"], [0, 4, 1, "", "asthsolid"], [0, 4, 1, "", "crustsolid"], [0, 4, 1, "", "error"], [0, 3, 1, "", "fitting"], [0, 4, 1, "", "hLith"], [0, 4, 1, "", "hc"], [0, 4, 1, "", "kAsth"], [0, 4, 1, "", "kCrust"], [0, 4, 1, "", "kLith"], [0, 4, 1, "", "lithsolid"], [0, 4, 1, "", "paleoWD"], [0, 4, 1, "", "qbase"], [0, 3, 1, "", "result"], [0, 4, 1, "", "rhp"], [0, 4, 1, "", "rift"], [0, 4, 1, "", "sediment_fill_margin"], [0, 3, 1, "", "sediments"], [0, 4, 1, "", "sediments_inputs"], [0, 4, 1, "", "simulated_at"], [0, 4, 1, "", "total_beta_tested"], [0, 4, 1, "", "water_depth_difference"]], "warmth.forward_modelling": [[0, 1, 1, "", "Forward_model"]], "warmth.forward_modelling.Forward_model": [[0, 2, 1, "", "add_sediments"], [0, 2, 1, "", "calculate_new_temperature"], [0, 2, 1, "", "compaction"], [0, 2, 1, "", "decompaction"], [0, 2, 1, "", "implicit_euler_solve"], [0, 2, 1, "", "simulate_continental"], [0, 2, 1, "", "simulate_single_node"]], "warmth.model": [[0, 1, 1, "", "Model"]], "warmth.model.Model": [[0, 3, 1, "", "builder"], [0, 2, 1, "", "load"], [0, 3, 1, "", "parameters"], [0, 3, 1, "", "simulator"]], "warmth.parameters": [[0, 1, 1, "", "Parameters"]], "warmth.parameters.Parameters": [[0, 3, 1, "", "initial_hLith_max"], [0, 3, 1, "", "initial_hLith_min"], [0, 3, 1, "", "positive_down"], [0, 3, 1, "", "time_end"], [0, 3, 1, "", "time_start"], [0, 3, 1, "", "time_step_Ma"]], "warmth.postprocessing": [[0, 1, 1, "", "Results"]], "warmth.postprocessing.Results": [[0, 3, 1, "", "ages"], [0, 2, 1, "", "basement_heatflow"], [0, 2, 1, "", "crust_thickness"], [0, 2, 1, "", "depth"], [0, 2, 1, "", "effective_conductivity"], [0, 2, 1, "", "heatflow"], [0, 2, 1, "", "lithosphere_thickness"], [0, 5, 1, "", "resultValues"], [0, 2, 1, "", "sediment_ids"], [0, 2, 1, "", "sediment_porosity"], [0, 2, 1, "", "temperature"], [0, 2, 1, "", "top_asthenosphere"], [0, 2, 1, "", "top_crust"], [0, 2, 1, "", "top_lithosphere"]], "warmth.simulator": [[0, 1, 1, "", "Simulator"]]}, "objtypes": {"0": "py:module", "1": "py:class", "2": "py:method", "3": "py:property", "4": "py:attribute", "5": "py:typeddict"}, "objnames": {"0": ["py", "module", "Python module"], "1": ["py", "class", "Python class"], "2": ["py", "method", "Python method"], "3": ["py", "property", "Python property"], "4": ["py", "attribute", "Python attribute"], "5": ["py", "typeddict", "Python typeddict"]}, "titleterms": {"api": 0, "refer": 0, "model": [0, 4], "builder": 0, "paramet": 0, "forward": 0, "simul": [0, 5], "result": 0, "exampl": 1, "content": [1, 2], "welcom": 2, "warmth": 2, "": 2, "document": 2, "indic": 2, "tabl": 2, "build": 4, "1d": [4, 5], "background": 5, "theori": 5, "pre": 5, "requisit": 5, "grid": 5, "node": 5, "3d": 5, "heat": 5, "equat": 5, "us": 5, "dolfinx": 5, "resqml": 5, "output": 5, "depend": 5}, "envversion": {"sphinx.domains.c": 3, "sphinx.domains.changeset": 1, "sphinx.domains.citation": 1, "sphinx.domains.cpp": 9, "sphinx.domains.index": 1, "sphinx.domains.javascript": 3, "sphinx.domains.math": 2, "sphinx.domains.python": 4, "sphinx.domains.rst": 2, "sphinx.domains.std": 2, "nbsphinx": 4, "sphinx": 60}, "alltitles": {"API Reference": [[0, "api-reference"]], "Model": [[0, "module-warmth.model"]], "Model builder": [[0, "module-warmth.build"]], "Parameters": [[0, "module-warmth.parameters"]], "Forward model": [[0, "module-warmth.forward_modelling"]], "Simulator": [[0, "module-warmth.simulator"]], "Results": [[0, "module-warmth.postprocessing"]], "Example": [[1, "example"]], "Contents:": [[1, null], [2, null]], "Welcome to warmth\u2019s documentation!": [[2, "welcome-to-warmth-s-documentation"]], "Indices and tables": [[2, "indices-and-tables"]], "Build 1D model": [[4, "Build-1D-model"]], "Background theory": [[5, "background-theory"]], "Pre-requisite: grid of 1D node simulation": [[5, "pre-requisite-grid-of-1d-node-simulation"]], "3D heat equation simulation using dolfinx": [[5, "d-heat-equation-simulation-using-dolfinx"]], "RESQML output": [[5, "resqml-output"]], "3D dependencies": [[5, "d-dependencies"]]}, "indexentries": {"builder (class in warmth.build)": [[0, "warmth.build.Builder"]], "forward_model (class in warmth.forward_modelling)": [[0, "warmth.forward_modelling.Forward_model"]], "grid (class in warmth.build)": [[0, "warmth.build.Grid"]], "model (class in warmth.model)": [[0, "warmth.model.Model"]], "parameters (class in warmth.parameters)": [[0, "warmth.parameters.Parameters"]], "results (class in warmth.postprocessing)": [[0, "warmth.postprocessing.Results"]], "results.resultvalues (typeddict in warmth.postprocessing)": [[0, "warmth.postprocessing.Results.resultValues"]], "simulator (class in warmth.simulator)": [[0, "warmth.simulator.Simulator"]], "t0 (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.T0"]], "tm (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.Tm"]], "x (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.X"]], "y (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.Y"]], "add_sediments() (warmth.forward_modelling.forward_model method)": [[0, "warmth.forward_modelling.Forward_model.add_sediments"]], "ages (warmth.postprocessing.results property)": [[0, "warmth.postprocessing.Results.ages"]], "asthsolid (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.asthsolid"]], "basement_heatflow() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.basement_heatflow"]], "builder (warmth.model.model property)": [[0, "warmth.model.Model.builder"]], "calculate_new_temperature() (warmth.forward_modelling.forward_model method)": [[0, "warmth.forward_modelling.Forward_model.calculate_new_temperature"]], "compaction() (warmth.forward_modelling.forward_model static method)": [[0, "warmth.forward_modelling.Forward_model.compaction"]], "crust_thickness() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.crust_thickness"]], "crustsolid (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.crustsolid"]], "decompaction() (warmth.forward_modelling.forward_model static method)": [[0, "warmth.forward_modelling.Forward_model.decompaction"]], "define_geometry() (warmth.build.builder method)": [[0, "warmth.build.Builder.define_geometry"]], "depth() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.depth"]], "dump() (warmth.build.grid method)": [[0, "warmth.build.Grid.dump"]], "effective_conductivity() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.effective_conductivity"]], "error (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.error"]], "extract_nodes() (warmth.build.builder method)": [[0, "warmth.build.Builder.extract_nodes"]], "fitting (warmth.build.single_node property)": [[0, "warmth.build.single_node.fitting"]], "hlith (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.hLith"]], "hc (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.hc"]], "heatflow() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.heatflow"]], "implicit_euler_solve() (warmth.forward_modelling.forward_model method)": [[0, "warmth.forward_modelling.Forward_model.implicit_euler_solve"]], "indexing_arr (warmth.build.grid property)": [[0, "warmth.build.Grid.indexing_arr"]], "initial_hlith_max (warmth.parameters.parameters property)": [[0, "warmth.parameters.Parameters.initial_hLith_max"]], "initial_hlith_min (warmth.parameters.parameters property)": [[0, "warmth.parameters.Parameters.initial_hLith_min"]], "input_horizons_template (warmth.build.builder property)": [[0, "warmth.build.Builder.input_horizons_template"]], "iter_node() (warmth.build.builder method)": [[0, "warmth.build.Builder.iter_node"]], "kasth (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.kAsth"]], "kcrust (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.kCrust"]], "klith (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.kLith"]], "lithosphere_thickness() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.lithosphere_thickness"]], "lithsolid (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.lithsolid"]], "load() (warmth.model.model method)": [[0, "warmth.model.Model.load"]], "location_grid (warmth.build.grid property)": [[0, "warmth.build.Grid.location_grid"]], "locations (warmth.build.builder property)": [[0, "warmth.build.Builder.locations"]], "make_grid_arr() (warmth.build.grid method)": [[0, "warmth.build.Grid.make_grid_arr"]], "module": [[0, "module-warmth.build"], [0, "module-warmth.forward_modelling"], [0, "module-warmth.model"], [0, "module-warmth.parameters"], [0, "module-warmth.postprocessing"], [0, "module-warmth.simulator"]], "paleowd (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.paleoWD"]], "parameters (warmth.model.model property)": [[0, "warmth.model.Model.parameters"]], "positive_down (warmth.parameters.parameters property)": [[0, "warmth.parameters.Parameters.positive_down"]], "qbase (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.qbase"]], "result (warmth.build.single_node property)": [[0, "warmth.build.single_node.result"]], "rhp (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.rhp"]], "rift (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.rift"]], "sediment_fill_margin (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.sediment_fill_margin"]], "sediment_ids() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.sediment_ids"]], "sediment_porosity() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.sediment_porosity"]], "sediments (warmth.build.single_node property)": [[0, "warmth.build.single_node.sediments"]], "sediments_inputs (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.sediments_inputs"]], "set_eustatic_sea_level() (warmth.build.builder method)": [[0, "warmth.build.Builder.set_eustatic_sea_level"]], "simulate_continental() (warmth.forward_modelling.forward_model method)": [[0, "warmth.forward_modelling.Forward_model.simulate_continental"]], "simulate_single_node() (warmth.forward_modelling.forward_model method)": [[0, "warmth.forward_modelling.Forward_model.simulate_single_node"]], "simulated_at (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.simulated_at"]], "simulator (warmth.model.model property)": [[0, "warmth.model.Model.simulator"]], "single_node (class in warmth.build)": [[0, "warmth.build.single_node"]], "single_node_sediments_inputs_template (warmth.build.builder property)": [[0, "warmth.build.Builder.single_node_sediments_inputs_template"]], "temperature() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.temperature"]], "time_end (warmth.parameters.parameters property)": [[0, "warmth.parameters.Parameters.time_end"]], "time_start (warmth.parameters.parameters property)": [[0, "warmth.parameters.Parameters.time_start"]], "time_step_ma (warmth.parameters.parameters property)": [[0, "warmth.parameters.Parameters.time_step_Ma"]], "top_asthenosphere() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.top_asthenosphere"]], "top_crust() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.top_crust"]], "top_lithosphere() (warmth.postprocessing.results method)": [[0, "warmth.postprocessing.Results.top_lithosphere"]], "total_beta_tested (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.total_beta_tested"]], "warmth.build": [[0, "module-warmth.build"]], "warmth.forward_modelling": [[0, "module-warmth.forward_modelling"]], "warmth.model": [[0, "module-warmth.model"]], "warmth.parameters": [[0, "module-warmth.parameters"]], "warmth.postprocessing": [[0, "module-warmth.postprocessing"]], "warmth.simulator": [[0, "module-warmth.simulator"]], "water_depth_difference (warmth.build.single_node attribute)": [[0, "warmth.build.single_node.water_depth_difference"]]}})
\ No newline at end of file