{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"scrolled": true
},
"outputs": [],
"source": [
"from idf_analysis.idf_class import IntensityDurationFrequencyAnalyse\n",
"from idf_analysis.definitions import *\n",
"import pandas as pd\n",
"from os import path\n",
"%matplotlib inline\n",
"import matplotlib.pyplot as plt\n",
"plt.style.use('bmh')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Intensity Duration Frequency Analyse"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Parameter\n",
"\n",
"**series_kind**:\n",
"\n",
"`SERIES.PARTIAL` = Partielle Serie (partial duration series, PDS) (peak over threshold, POT)\n",
"\n",
"`SERIES.ANNUAL` = Jährliche Serie (annual maximum series, AMS)\n",
"\n",
"**worksheet**:\n",
"\n",
"`METHOD.KOSTRA`:\n",
"- DWA-A 531\n",
"- KOSTRA - empfohlen\n",
"- Stützstellen: 60 min und 12 h\n",
"\n",
"`METHOD.CONVECTIVE_ADVECTIVE`:\n",
"- DWA-A 531\n",
"- Unterscheidung in überwiegend konvektiv und advektiv verursachte Starkregen\n",
"- Stützstellen: 3 h und 24 h\n",
"\n",
"`METHOD.ATV`:\n",
"- ATV-A 121\n",
"- Stützstellen: 3 h und 48 h\n",
"\n",
"**extended_durations** = Inkludiert die Dauerstufen `[0.75d, 1d, 2d, 3d, 4d, 5d, 6d]` in der Analyse (in d=Tage)\n",
"\n",
"Standardmäßig berechnete Dauerstufen `[5m, 10m, 15m, 20m, 30m, 45m, 60m, 1.5h, 3h, 4.5h, 6h, 7.5h, 10h, 12h]`"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"scrolled": true
},
"outputs": [],
"source": [
"idf = IntensityDurationFrequencyAnalyse(series_kind=SERIES.PARTIAL, worksheet=METHOD.KOSTRA, extended_durations=True)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"I used the rain-time-series from ehyd.gv.at with the ID 112086 (Graz-Andritz) [created with the ehyd-tools package](https://github.com/MarkusPic/ehyd_tools/blob/main/example/example_python_api.ipynb). "
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"scrolled": true
},
"outputs": [],
"source": [
"data = pd.read_parquet('ehyd_112086.parquet')"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"scrolled": true
},
"outputs": [],
"source": [
"output_directory = 'ehyd_112086_idf_data'"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/plain": " precipitation\ndatetime \n2007-09-18 11:09:00 0.1\n2007-09-18 11:12:00 0.1\n2007-09-18 11:13:00 0.1\n2007-09-18 11:14:00 0.1\n2007-09-18 11:20:00 0.1",
"text/html": "
\n\n
\n \n \n | \n precipitation | \n
\n \n datetime | \n | \n
\n \n \n \n 2007-09-18 11:09:00 | \n 0.1 | \n
\n \n 2007-09-18 11:12:00 | \n 0.1 | \n
\n \n 2007-09-18 11:13:00 | \n 0.1 | \n
\n \n 2007-09-18 11:14:00 | \n 0.1 | \n
\n \n 2007-09-18 11:20:00 | \n 0.1 | \n
\n \n
\n
"
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"data.head()"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/plain": " precipitation\ndatetime \n2016-12-28 22:24:00 0.006\n2016-12-28 22:25:00 0.092\n2016-12-28 22:26:00 0.006\n2016-12-28 22:45:00 0.090\n2016-12-28 22:46:00 0.006",
"text/html": "\n\n
\n \n \n | \n precipitation | \n
\n \n datetime | \n | \n
\n \n \n \n 2016-12-28 22:24:00 | \n 0.006 | \n
\n \n 2016-12-28 22:25:00 | \n 0.092 | \n
\n \n 2016-12-28 22:26:00 | \n 0.006 | \n
\n \n 2016-12-28 22:45:00 | \n 0.090 | \n
\n \n 2016-12-28 22:46:00 | \n 0.006 | \n
\n \n
\n
"
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"data.tail()"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"scrolled": true
},
"outputs": [],
"source": [
"idf.set_series(data['precipitation'])"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Bei jeder neuen Berechnung werden Zwischenergebnisse erstellt, welche nur abhängig von der gewählten Serie `series_kind` und der angegebenen/benötigten Dauerstufen sind. Dieser Vorgang dauert einige Sekunden.\n",
"Auserdem enthalten diese Zwischenergebnisse die Parameter, die zur Berechnung der Regenhöhe und Regenspende benötigt werden.\n",
"Hier sind bereist die Berechnungsverfahren und Stückpunkte laut dem gewählten `worksheet` berücksichtigt."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Um Zeit zu sparen, gibt es die Möglichkeit, die Parameter zwischenzuspeichern und bei erneutem Aufrufen des Skripts werden diese Parameter nicht mehr berechnet, sondern aus der Datei gelesen."
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"scrolled": true
},
"outputs": [],
"source": [
"idf.auto_save_parameters(path.join(output_directory, 'idf_parameters.yaml'))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Abgerufen können diese Zwischenergebnisse mit:"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/plain": ""
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"idf.parameters"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Berechnungen"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/plain": "19.031596336052708"
},
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"idf.depth_of_rainfall(duration=15, return_period=1)"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {
"scrolled": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Resultierende Regenhöhe h_N(T_n=1.0a, D=15.0min) = 19.03 mm\n"
]
}
],
"source": [
"print('Resultierende Regenhöhe h_N(T_n={t:0.1f}a, D={d:0.1f}min) = {h:0.2f} mm'\n",
" ''.format(t=1, d=15, h=idf.depth_of_rainfall(15, 1)))"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/plain": "211.46218151169674"
},
"execution_count": 12,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"idf.rain_flow_rate(duration=15, return_period=1)"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {
"scrolled": true
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Resultierende Regenspende r_N(T_n=1.0a, D=15.0min) = 211.46 L/(s*ha)\n"
]
}
],
"source": [
"print('Resultierende Regenspende r_N(T_n={t:0.1f}a, D={d:0.1f}min) = {r:0.2f} L/(s*ha)'\n",
" ''.format(t=1, d=15, r=idf.rain_flow_rate(15, 1)))"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/plain": "11.410836729727"
},
"execution_count": 14,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"idf.r_720_1()"
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/plain": "0.1430180144131331"
},
"execution_count": 15,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"idf.get_return_period(height_of_rainfall=10, duration=15)"
]
},
{
"cell_type": "code",
"execution_count": 16,
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"text/plain": "5.433080747189968"
},
"execution_count": 16,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"idf.get_duration(height_of_rainfall=10, return_period=1)"
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {
"scrolled": false
},
"outputs": [
{
"data": {
"text/plain": " 1 2 3 5 10 20 25 30 50 \\\n5 9.39 10.97 11.89 13.04 14.61 16.19 16.69 17.11 18.26 \n10 15.15 17.62 19.06 20.88 23.35 25.82 26.62 27.27 29.09 \n15 19.03 22.25 24.13 26.51 29.72 32.94 33.98 34.83 37.20 \n20 21.83 25.71 27.99 30.85 34.73 38.62 39.87 40.89 43.75 \n30 25.60 30.66 33.62 37.35 42.41 47.47 49.10 50.43 54.16 \n45 28.92 35.51 39.37 44.23 50.83 57.42 59.54 61.28 66.14 \n60 30.93 38.89 43.54 49.40 57.36 65.31 67.88 69.97 75.83 \n90 33.37 41.74 46.64 52.80 61.17 69.54 72.23 74.43 80.60 \n120 35.22 43.90 48.97 55.36 64.03 72.70 75.49 77.78 84.17 \n180 38.01 47.13 52.46 59.18 68.30 77.42 80.36 82.76 89.48 \n240 40.12 49.57 55.10 62.06 71.51 80.97 84.01 86.49 93.46 \n360 43.29 53.23 59.04 66.37 76.31 86.25 89.45 92.06 99.39 \n540 46.71 57.16 63.27 70.98 81.43 91.89 95.25 98.00 105.71 \n720 49.29 60.13 66.47 74.45 85.29 96.12 99.61 102.46 110.44 \n1080 54.41 64.97 71.15 78.94 89.50 100.06 103.46 106.24 114.02 \n1440 58.02 67.72 73.39 80.54 90.24 99.93 103.05 105.61 112.75 \n2880 66.70 77.41 83.68 91.57 102.29 113.00 116.45 119.26 127.16 \n4320 71.93 85.72 93.78 103.95 117.73 131.52 135.96 139.58 149.75 \n5760 78.95 95.65 105.42 117.72 134.43 151.13 156.50 160.89 173.20 \n7200 83.53 101.38 111.82 124.98 142.83 160.68 166.43 171.12 184.28 \n8640 85.38 104.95 116.40 130.82 150.38 169.95 176.25 181.40 195.82 \n\n 75 100 \n5 19.18 19.83 \n10 30.54 31.56 \n15 39.08 40.42 \n20 46.02 47.63 \n30 57.12 59.22 \n45 69.99 72.73 \n60 80.49 83.79 \n90 85.49 88.96 \n120 89.24 92.84 \n180 94.81 98.60 \n240 98.99 102.91 \n360 105.20 109.33 \n540 111.82 116.16 \n720 116.78 121.28 \n1080 120.20 124.58 \n1440 118.42 122.45 \n2880 133.42 137.87 \n4320 157.81 163.53 \n5760 182.97 189.90 \n7200 194.72 202.13 \n8640 207.27 215.39 ",
"text/html": "\n\n
\n \n \n | \n 1 | \n 2 | \n 3 | \n 5 | \n 10 | \n 20 | \n 25 | \n 30 | \n 50 | \n 75 | \n 100 | \n
\n \n \n \n 5 | \n 9.39 | \n 10.97 | \n 11.89 | \n 13.04 | \n 14.61 | \n 16.19 | \n 16.69 | \n 17.11 | \n 18.26 | \n 19.18 | \n 19.83 | \n
\n \n 10 | \n 15.15 | \n 17.62 | \n 19.06 | \n 20.88 | \n 23.35 | \n 25.82 | \n 26.62 | \n 27.27 | \n 29.09 | \n 30.54 | \n 31.56 | \n
\n \n 15 | \n 19.03 | \n 22.25 | \n 24.13 | \n 26.51 | \n 29.72 | \n 32.94 | \n 33.98 | \n 34.83 | \n 37.20 | \n 39.08 | \n 40.42 | \n
\n \n 20 | \n 21.83 | \n 25.71 | \n 27.99 | \n 30.85 | \n 34.73 | \n 38.62 | \n 39.87 | \n 40.89 | \n 43.75 | \n 46.02 | \n 47.63 | \n
\n \n 30 | \n 25.60 | \n 30.66 | \n 33.62 | \n 37.35 | \n 42.41 | \n 47.47 | \n 49.10 | \n 50.43 | \n 54.16 | \n 57.12 | \n 59.22 | \n
\n \n 45 | \n 28.92 | \n 35.51 | \n 39.37 | \n 44.23 | \n 50.83 | \n 57.42 | \n 59.54 | \n 61.28 | \n 66.14 | \n 69.99 | \n 72.73 | \n
\n \n 60 | \n 30.93 | \n 38.89 | \n 43.54 | \n 49.40 | \n 57.36 | \n 65.31 | \n 67.88 | \n 69.97 | \n 75.83 | \n 80.49 | \n 83.79 | \n
\n \n 90 | \n 33.37 | \n 41.74 | \n 46.64 | \n 52.80 | \n 61.17 | \n 69.54 | \n 72.23 | \n 74.43 | \n 80.60 | \n 85.49 | \n 88.96 | \n
\n \n 120 | \n 35.22 | \n 43.90 | \n 48.97 | \n 55.36 | \n 64.03 | \n 72.70 | \n 75.49 | \n 77.78 | \n 84.17 | \n 89.24 | \n 92.84 | \n
\n \n 180 | \n 38.01 | \n 47.13 | \n 52.46 | \n 59.18 | \n 68.30 | \n 77.42 | \n 80.36 | \n 82.76 | \n 89.48 | \n 94.81 | \n 98.60 | \n
\n \n 240 | \n 40.12 | \n 49.57 | \n 55.10 | \n 62.06 | \n 71.51 | \n 80.97 | \n 84.01 | \n 86.49 | \n 93.46 | \n 98.99 | \n 102.91 | \n
\n \n 360 | \n 43.29 | \n 53.23 | \n 59.04 | \n 66.37 | \n 76.31 | \n 86.25 | \n 89.45 | \n 92.06 | \n 99.39 | \n 105.20 | \n 109.33 | \n
\n \n 540 | \n 46.71 | \n 57.16 | \n 63.27 | \n 70.98 | \n 81.43 | \n 91.89 | \n 95.25 | \n 98.00 | \n 105.71 | \n 111.82 | \n 116.16 | \n
\n \n 720 | \n 49.29 | \n 60.13 | \n 66.47 | \n 74.45 | \n 85.29 | \n 96.12 | \n 99.61 | \n 102.46 | \n 110.44 | \n 116.78 | \n 121.28 | \n
\n \n 1080 | \n 54.41 | \n 64.97 | \n 71.15 | \n 78.94 | \n 89.50 | \n 100.06 | \n 103.46 | \n 106.24 | \n 114.02 | \n 120.20 | \n 124.58 | \n
\n \n 1440 | \n 58.02 | \n 67.72 | \n 73.39 | \n 80.54 | \n 90.24 | \n 99.93 | \n 103.05 | \n 105.61 | \n 112.75 | \n 118.42 | \n 122.45 | \n
\n \n 2880 | \n 66.70 | \n 77.41 | \n 83.68 | \n 91.57 | \n 102.29 | \n 113.00 | \n 116.45 | \n 119.26 | \n 127.16 | \n 133.42 | \n 137.87 | \n
\n \n 4320 | \n 71.93 | \n 85.72 | \n 93.78 | \n 103.95 | \n 117.73 | \n 131.52 | \n 135.96 | \n 139.58 | \n 149.75 | \n 157.81 | \n 163.53 | \n
\n \n 5760 | \n 78.95 | \n 95.65 | \n 105.42 | \n 117.72 | \n 134.43 | \n 151.13 | \n 156.50 | \n 160.89 | \n 173.20 | \n 182.97 | \n 189.90 | \n
\n \n 7200 | \n 83.53 | \n 101.38 | \n 111.82 | \n 124.98 | \n 142.83 | \n 160.68 | \n 166.43 | \n 171.12 | \n 184.28 | \n 194.72 | \n 202.13 | \n
\n \n 8640 | \n 85.38 | \n 104.95 | \n 116.40 | \n 130.82 | \n 150.38 | \n 169.95 | \n 176.25 | \n 181.40 | \n 195.82 | \n 207.27 | \n 215.39 | \n
\n \n
\n
"
},
"execution_count": 17,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"idf.result_table().round(2)"
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {
"scrolled": false
},
"outputs": [
{
"data": {
"text/plain": "return period (a) 1 2 3 5 10 20 25 \\\nfrequency (1/a) 1.000 0.500 0.333 0.200 0.100 0.050 0.040 \nduration (min) \n5 9.39 10.97 11.89 13.04 14.61 16.19 16.69 \n10 15.15 17.62 19.06 20.88 23.35 25.82 26.62 \n15 19.03 22.25 24.13 26.51 29.72 32.94 33.98 \n20 21.83 25.71 27.99 30.85 34.73 38.62 39.87 \n30 25.60 30.66 33.62 37.35 42.41 47.47 49.10 \n45 28.92 35.51 39.37 44.23 50.83 57.42 59.54 \n60 30.93 38.89 43.54 49.40 57.36 65.31 67.88 \n90 33.37 41.74 46.64 52.80 61.17 69.54 72.23 \n120 35.22 43.90 48.97 55.36 64.03 72.70 75.49 \n180 38.01 47.13 52.46 59.18 68.30 77.42 80.36 \n240 40.12 49.57 55.10 62.06 71.51 80.97 84.01 \n360 43.29 53.23 59.04 66.37 76.31 86.25 89.45 \n540 46.71 57.16 63.27 70.98 81.43 91.89 95.25 \n720 49.29 60.13 66.47 74.45 85.29 96.12 99.61 \n1080 54.41 64.97 71.15 78.94 89.50 100.06 103.46 \n1440 58.02 67.72 73.39 80.54 90.24 99.93 103.05 \n2880 66.70 77.41 83.68 91.57 102.29 113.00 116.45 \n4320 71.93 85.72 93.78 103.95 117.73 131.52 135.96 \n5760 78.95 95.65 105.42 117.72 134.43 151.13 156.50 \n7200 83.53 101.38 111.82 124.98 142.83 160.68 166.43 \n8640 85.38 104.95 116.40 130.82 150.38 169.95 176.25 \n\nreturn period (a) 30 50 75 100 \nfrequency (1/a) 0.033 0.020 0.013 0.010 \nduration (min) \n5 17.11 18.26 19.18 19.83 \n10 27.27 29.09 30.54 31.56 \n15 34.83 37.20 39.08 40.42 \n20 40.89 43.75 46.02 47.63 \n30 50.43 54.16 57.12 59.22 \n45 61.28 66.14 69.99 72.73 \n60 69.97 75.83 80.49 83.79 \n90 74.43 80.60 85.49 88.96 \n120 77.78 84.17 89.24 92.84 \n180 82.76 89.48 94.81 98.60 \n240 86.49 93.46 98.99 102.91 \n360 92.06 99.39 105.20 109.33 \n540 98.00 105.71 111.82 116.16 \n720 102.46 110.44 116.78 121.28 \n1080 106.24 114.02 120.20 124.58 \n1440 105.61 112.75 118.42 122.45 \n2880 119.26 127.16 133.42 137.87 \n4320 139.58 149.75 157.81 163.53 \n5760 160.89 173.20 182.97 189.90 \n7200 171.12 184.28 194.72 202.13 \n8640 181.40 195.82 207.27 215.39 ",
"text/html": "\n\n
\n \n \n return period (a) | \n 1 | \n 2 | \n 3 | \n 5 | \n 10 | \n 20 | \n 25 | \n 30 | \n 50 | \n 75 | \n 100 | \n
\n \n frequency (1/a) | \n 1.000 | \n 0.500 | \n 0.333 | \n 0.200 | \n 0.100 | \n 0.050 | \n 0.040 | \n 0.033 | \n 0.020 | \n 0.013 | \n 0.010 | \n
\n \n duration (min) | \n | \n | \n | \n | \n | \n | \n | \n | \n | \n | \n | \n
\n \n \n \n 5 | \n 9.39 | \n 10.97 | \n 11.89 | \n 13.04 | \n 14.61 | \n 16.19 | \n 16.69 | \n 17.11 | \n 18.26 | \n 19.18 | \n 19.83 | \n
\n \n 10 | \n 15.15 | \n 17.62 | \n 19.06 | \n 20.88 | \n 23.35 | \n 25.82 | \n 26.62 | \n 27.27 | \n 29.09 | \n 30.54 | \n 31.56 | \n
\n \n 15 | \n 19.03 | \n 22.25 | \n 24.13 | \n 26.51 | \n 29.72 | \n 32.94 | \n 33.98 | \n 34.83 | \n 37.20 | \n 39.08 | \n 40.42 | \n
\n \n 20 | \n 21.83 | \n 25.71 | \n 27.99 | \n 30.85 | \n 34.73 | \n 38.62 | \n 39.87 | \n 40.89 | \n 43.75 | \n 46.02 | \n 47.63 | \n
\n \n 30 | \n 25.60 | \n 30.66 | \n 33.62 | \n 37.35 | \n 42.41 | \n 47.47 | \n 49.10 | \n 50.43 | \n 54.16 | \n 57.12 | \n 59.22 | \n
\n \n 45 | \n 28.92 | \n 35.51 | \n 39.37 | \n 44.23 | \n 50.83 | \n 57.42 | \n 59.54 | \n 61.28 | \n 66.14 | \n 69.99 | \n 72.73 | \n
\n \n 60 | \n 30.93 | \n 38.89 | \n 43.54 | \n 49.40 | \n 57.36 | \n 65.31 | \n 67.88 | \n 69.97 | \n 75.83 | \n 80.49 | \n 83.79 | \n
\n \n 90 | \n 33.37 | \n 41.74 | \n 46.64 | \n 52.80 | \n 61.17 | \n 69.54 | \n 72.23 | \n 74.43 | \n 80.60 | \n 85.49 | \n 88.96 | \n
\n \n 120 | \n 35.22 | \n 43.90 | \n 48.97 | \n 55.36 | \n 64.03 | \n 72.70 | \n 75.49 | \n 77.78 | \n 84.17 | \n 89.24 | \n 92.84 | \n
\n \n 180 | \n 38.01 | \n 47.13 | \n 52.46 | \n 59.18 | \n 68.30 | \n 77.42 | \n 80.36 | \n 82.76 | \n 89.48 | \n 94.81 | \n 98.60 | \n
\n \n 240 | \n 40.12 | \n 49.57 | \n 55.10 | \n 62.06 | \n 71.51 | \n 80.97 | \n 84.01 | \n 86.49 | \n 93.46 | \n 98.99 | \n 102.91 | \n
\n \n 360 | \n 43.29 | \n 53.23 | \n 59.04 | \n 66.37 | \n 76.31 | \n 86.25 | \n 89.45 | \n 92.06 | \n 99.39 | \n 105.20 | \n 109.33 | \n
\n \n 540 | \n 46.71 | \n 57.16 | \n 63.27 | \n 70.98 | \n 81.43 | \n 91.89 | \n 95.25 | \n 98.00 | \n 105.71 | \n 111.82 | \n 116.16 | \n
\n \n 720 | \n 49.29 | \n 60.13 | \n 66.47 | \n 74.45 | \n 85.29 | \n 96.12 | \n 99.61 | \n 102.46 | \n 110.44 | \n 116.78 | \n 121.28 | \n
\n \n 1080 | \n 54.41 | \n 64.97 | \n 71.15 | \n 78.94 | \n 89.50 | \n 100.06 | \n 103.46 | \n 106.24 | \n 114.02 | \n 120.20 | \n 124.58 | \n
\n \n 1440 | \n 58.02 | \n 67.72 | \n 73.39 | \n 80.54 | \n 90.24 | \n 99.93 | \n 103.05 | \n 105.61 | \n 112.75 | \n 118.42 | \n 122.45 | \n
\n \n 2880 | \n 66.70 | \n 77.41 | \n 83.68 | \n 91.57 | \n 102.29 | \n 113.00 | \n 116.45 | \n 119.26 | \n 127.16 | \n 133.42 | \n 137.87 | \n
\n \n 4320 | \n 71.93 | \n 85.72 | \n 93.78 | \n 103.95 | \n 117.73 | \n 131.52 | \n 135.96 | \n 139.58 | \n 149.75 | \n 157.81 | \n 163.53 | \n
\n \n 5760 | \n 78.95 | \n 95.65 | \n 105.42 | \n 117.72 | \n 134.43 | \n 151.13 | \n 156.50 | \n 160.89 | \n 173.20 | \n 182.97 | \n 189.90 | \n
\n \n 7200 | \n 83.53 | \n 101.38 | \n 111.82 | \n 124.98 | \n 142.83 | \n 160.68 | \n 166.43 | \n 171.12 | \n 184.28 | \n 194.72 | \n 202.13 | \n
\n \n 8640 | \n 85.38 | \n 104.95 | \n 116.40 | \n 130.82 | \n 150.38 | \n 169.95 | \n 176.25 | \n 181.40 | \n 195.82 | \n 207.27 | \n 215.39 | \n
\n \n
\n
"
},
"execution_count": 18,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"idf.result_table(add_names=True).round(2)"
]
},
{
"cell_type": "markdown",
"source": [
"To save the table as a csv:"
],
"metadata": {
"collapsed": false
}
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {},
"outputs": [],
"source": [
"idf.result_table(add_names=True).round(2).to_csv(path.join(output_directory, 'idf_table_UNIX.csv'), sep=',', decimal='.', float_format='%0.2f')"
]
},
{
"cell_type": "code",
"execution_count": 20,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"return period (a) 1 2 3 5 10 20 25 30 50 75 100\n",
"frequency (1/a) 1.000 0.500 0.333 0.200 0.100 0.050 0.040 0.033 0.020 0.013 0.010\n",
"duration (min) \n",
"5 9.39 10.97 11.89 13.04 14.61 16.19 16.69 17.11 18.26 19.18 19.83\n",
"10 15.15 17.62 19.06 20.88 23.35 25.82 26.62 27.27 29.09 30.54 31.56\n",
"15 19.03 22.25 24.13 26.51 29.72 32.94 33.98 34.83 37.20 39.08 40.42\n",
"20 21.83 25.71 27.99 30.85 34.73 38.62 39.87 40.89 43.75 46.02 47.63\n",
"30 25.60 30.66 33.62 37.35 42.41 47.47 49.10 50.43 54.16 57.12 59.22\n",
"45 28.92 35.51 39.37 44.23 50.83 57.42 59.54 61.28 66.14 69.99 72.73\n",
"60 30.93 38.89 43.54 49.40 57.36 65.31 67.88 69.97 75.83 80.49 83.79\n",
"90 33.37 41.74 46.64 52.80 61.17 69.54 72.23 74.43 80.60 85.49 88.96\n",
"120 35.22 43.90 48.97 55.36 64.03 72.70 75.49 77.78 84.17 89.24 92.84\n",
"180 38.01 47.13 52.46 59.18 68.30 77.42 80.36 82.76 89.48 94.81 98.60\n",
"240 40.12 49.57 55.10 62.06 71.51 80.97 84.01 86.49 93.46 98.99 102.91\n",
"360 43.29 53.23 59.04 66.37 76.31 86.25 89.45 92.06 99.39 105.20 109.33\n",
"540 46.71 57.16 63.27 70.98 81.43 91.89 95.25 98.00 105.71 111.82 116.16\n",
"720 49.29 60.13 66.47 74.45 85.29 96.12 99.61 102.46 110.44 116.78 121.28\n",
"1080 54.41 64.97 71.15 78.94 89.50 100.06 103.46 106.24 114.02 120.20 124.58\n",
"1440 58.02 67.72 73.39 80.54 90.24 99.93 103.05 105.61 112.75 118.42 122.45\n",
"2880 66.70 77.41 83.68 91.57 102.29 113.00 116.45 119.26 127.16 133.42 137.87\n",
"4320 71.93 85.72 93.78 103.95 117.73 131.52 135.96 139.58 149.75 157.81 163.53\n",
"5760 78.95 95.65 105.42 117.72 134.43 151.13 156.50 160.89 173.20 182.97 189.90\n",
"7200 83.53 101.38 111.82 124.98 142.83 160.68 166.43 171.12 184.28 194.72 202.13\n",
"8640 85.38 104.95 116.40 130.82 150.38 169.95 176.25 181.40 195.82 207.27 215.39\n"
]
}
],
"source": [
"print(idf.result_table(add_names=True).round(2).to_string())"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"To create a color plot of the IDF curves:"
]
},
{
"cell_type": "code",
"execution_count": 21,
"metadata": {
"scrolled": false
},
"outputs": [
{
"data": {
"text/plain": "