Output plots part 2

Weill, S., et al. « Coupling Water Flow and Solute Transport into a Physically-Based Surface–Subsurface Hydrological Model ». Advances in Water Resources, vol. 34, no 1, janvier 2011, p. 128‑36. DOI.org (Crossref), https://doi.org/10.1016/j.advwatres.2010.10.001.

This example shows how to use pyCATHY object to plot the most common ouputs of the hydrological model.

Estimated time to run the notebook = 5min

Here we need to import cathy_tools class that control the CATHY core files preprocessing and processing We also import cathy_plots to render the results

from pyCATHY import cathy_tools
from pyCATHY.plotters import cathy_plots as cplt

Traceback (most recent call last):
  File "/home/runner/work/pycathy_wrapper/pycathy_wrapper/examples/SSHydro/plot_4b_pyCATHY_outputs.py", line 20, in <module>
    from pyCATHY import cathy_tools
  File "/home/runner/work/pycathy_wrapper/pycathy_wrapper/pyCATHY/__init__.py", line 4, in <module>
    from . import cathy_utils, meshtools, sensitivity
  File "/home/runner/work/pycathy_wrapper/pycathy_wrapper/pyCATHY/meshtools.py", line 1317, in <module>
    from sklearn.metrics import r2_score
ModuleNotFoundError: No module named 'sklearn'

if you add True to verbose, the processor log will be printed in the window shell

path2prj = "../SSHydro/"  # add your local path here
simu = cathy_tools.CATHY(dirName=path2prj,
                    prj_name="weil_exemple_outputs_plot"
                    )

simu.run_preprocessor()

simu.read_inputs('atmbc')
# simu.create_mesh_bounds_df(
#     'nansfdirbc',
#     simu.grid3d["mesh3d_nodes"],
#     t_atmbc,
# )

# print(simu.mesh_bound_cond_df.head())
# print("BC columns:", list(simu.mesh_bound_cond_df.columns))

# Apply no-flow to all lateral/bottom boundaries
simu.update_nansfneubc(no_flow=True)
simu.update_nansfdirbc(no_flow=True)
simu.update_sfbc(no_flow=True)

simu.run_processor(IPRT1=2,
                    DTMIN=1e-2,
                    DTMAX=1e2,
                    DELTAT=5,
                   TRAFLAG=0,
                   verbose=False
                   )

simu.show(prop="hgsfdet")

simu.show(prop="dtcoupling", yprop="Atmpot-d")

simu.show(prop="hgraph")

simu.show(prop="cumflowvol")

To select another time step change the value in the function argument

cplt.show_vtk(
    unit="pressure",
    timeStep=1,
    notebook=False,
    path=simu.workdir + "/weil_exemple_outputs_plot/vtk/",
)

cplt.show_vtk(

unit=”saturation”, timeStep=1, notebook=False, path=simu.workdir + “/my_cathy_prj/vtk/”,

)

df_recharge = simu.read_outputs('recharge')
xr_recharge = df_recharge.set_index(['time','X','Y']).to_xarray()

import matplotlib.pyplot as plt

fig, ax = plt.subplots()
xr_recharge['recharge'].isel(time=0).plot.imshow(ax=ax)


fig, ax = plt.subplots()
xr_recharge['recharge'].isel(time=1).plot.imshow(ax=ax)