Skip to content
Snippets Groups Projects
Commit 33f74ea0 authored by Julian Heinze's avatar Julian Heinze
Browse files

tiny fixes

parent 5348fd9c
No related branches found
No related tags found
1 merge request!240feflowlib: rename boundary conditions to subdomains
Pipeline #27340 passed with warnings
...@@ -48,7 +48,7 @@ ogs_sim_res = ms.mesh(ms.timesteps[-1]) ...@@ -48,7 +48,7 @@ ogs_sim_res = ms.mesh(ms.timesteps[-1])
ot.plot.contourf(ogs_sim_res, ot.variables.temperature) ot.plot.contourf(ogs_sim_res, ot.variables.temperature)
# %% # %%
# 6. The boundary meshes are manipulated to alter the model. # 6. The boundary meshes are manipulated to alter the model.
# The original boundary conditions are shown in this example: :ref:`sphx_glr_auto_examples_howto_conversions_plot_F_HT_simulation.py` # The original boundary conditions are shown in this example: :ref:`sphx_glr_auto_examples_howto_conversions_plot_F_feflowlib_HT_simulation.py`
# 6.1 The Dirichlet boundary conditions for the hydraulic head are set to 0. Therefore, no water flows from the left to the right edge. # 6.1 The Dirichlet boundary conditions for the hydraulic head are set to 0. Therefore, no water flows from the left to the right edge.
bc_flow = feflow_model.subdomains["P_BC_FLOW"]["P_BC_FLOW"] bc_flow = feflow_model.subdomains["P_BC_FLOW"]["P_BC_FLOW"]
feflow_model.subdomains["P_BC_FLOW"]["P_BC_FLOW"][bc_flow == 10] = 0 feflow_model.subdomains["P_BC_FLOW"]["P_BC_FLOW"][bc_flow == 10] = 0
......
...@@ -75,6 +75,7 @@ hydraulic_head_diff = ot.variables.Scalar( ...@@ -75,6 +75,7 @@ hydraulic_head_diff = ot.variables.Scalar(
data_name="HEAD_difference", data_unit="m", output_unit="m" data_name="HEAD_difference", data_unit="m", output_unit="m"
) )
ot.plot.contourf(diff_mesh, hydraulic_head_diff, vmin=-1.5e-9, vmax=1.5e-9) ot.plot.contourf(diff_mesh, hydraulic_head_diff, vmin=-1.5e-9, vmax=1.5e-9)
# %%
# Plot differences in temperature. # Plot differences in temperature.
feflow_model.mesh["temperature"] = feflow_model.mesh["P_TEMP"] feflow_model.mesh["temperature"] = feflow_model.mesh["P_TEMP"]
# Plot differences in temperature. # Plot differences in temperature.
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment