Skip to content

Commit

Permalink
Add script demo v0.3 to scripts
Browse files Browse the repository at this point in the history
  • Loading branch information
jnsebgosselin committed Apr 27, 2022
1 parent 1b4b1e9 commit 1c005f8
Showing 1 changed file with 53 additions and 0 deletions.
53 changes: 53 additions & 0 deletions pyhelp/scripts/example_pyhelp_outputs_v0.3.py
Original file line number Diff line number Diff line change
@@ -0,0 +1,53 @@
# -*- coding: utf-8 -*-
"""
Created on Tue Apr 26 09:21:49 2022
@author: Jean-Sébastien Gosselin
"""

from pyhelp.output import HelpOutput
import matplotlib.pyplot as plt

plt.close('all')
dirname = "D:/Projets/pyhelp/example/"

# On charge en mémoire les résultats PyHELP d'une simulation antérieure
# dont nous avions sauvegarder les résultats sur le disque.

# À noter que cela revient au même que d'utiliser l'output de la commande
# "output_help = helpm.calc_help_cells()" directement.
output_help = HelpOutput(dirname + "help_example.out")

# On peut ensuite exporter, dans un fichier CSV, les moyennes annuelles du
# bilan hydrologique.

# On peut également spécifier (disponible dans pyhelp 0.3) les années pour
# lesquelles ont souhaite calculer les moyennes annuelles du bilan. Par
# exemple, pour calculer les moyennes entre 2003 et 2009 inclusivement, on
# écrirait :
output_help.save_to_csv(
dirname + "help_example_yearly_2003-2009.csv",
year_from=2003,
year_to=2007)

# On peut également produires des graphiques des moyennes annuelles et
# ensuelles du bilan hydrologique de même que les valeurs annuelles du
# bilan hydrologique à l'aide des 3 fonctions ci-dessous. Ici encore,
# Il est possible de spécifier la période pour laquelle ont veut calculer
# les moyennes ou afficher les valeurs annuelles.
output_help.plot_area_monthly_avg(year_from=2003, year_to=2007)
output_help.plot_area_yearly_avg(year_from=2005, year_to=2005)
output_help.plot_area_yearly_series(year_from=2004, year_to=2009)

# Les valeurs numériques derrière les graphiques produits ci-dessus
# peuvent être accédées avec les fonctions suivantes:
yearly_avg = output_help.calc_area_yearly_avg()
monthly_avg = output_help.calc_area_monthly_avg()

# Il est très facile d'exporter ces résultats dans des csv. Par exemple,
# pour exporter les valeurs mensuelles moyennes de recharge, nous
# ferions tout simplement:
monthly_avg['rechg'].to_csv(dirname + "monthly_rechg.csv")

# Similairement, pour exporter les valeurs moyennes annuelles
# d'évapotranspiration, on pourrait faire:
yearly_avg['evapo'].to_csv(dirname + "yearly_evapo.csv")

0 comments on commit 1c005f8

Please sign in to comment.