Clean up rdy for elke (#15)

* added files and directory for Elke, a new BSc student

* Cleaned up the code cells of Thirza, but kept original as well

* everything ready for Elke

Author: MarkMelotto

book/thesis_projects/BSc/2025_Q3_ThirzaVanEsch_CEG/Report/CHAPTER2.ipynb

@@ -258,20 +258,6 @@
     "The lines from fit 1 to fit 5 are plotted as well to give an indication of the influence of certain peak discharges on the extrapolation. This will be useful later in the research, when the significance of the difference in future return periods with different climate change scenarios compared to the current return periods needs to be assessed. \n"
    ]
   },
-  {
-   "cell_type": "code",
-   "execution_count": 8,
-   "id": "9dc2476f-dcba-4ed4-aa72-72c4acb29840",
-   "metadata": {
-    "jupyter": {
-     "source_hidden": true
-    }
-   },
-   "outputs": [],
-   "source": [
-    "# max(discharge)"
-   ]
-  },
   {
    "cell_type": "code",
    "execution_count": 9,
@@ -289,60 +275,6 @@
     "# print(maxdischarge)"
    ]
   },
-  {
-   "cell_type": "code",
-   "execution_count": 10,
-   "id": "e6693540-e767-4a18-af9e-e5b6113d8ae7",
-   "metadata": {
-    "jupyter": {
-     "source_hidden": true
-    }
-   },
-   "outputs": [],
-   "source": [
-    "# import scipy.stats as stats\n",
-    "# import math\n",
-    "# # # discharge data > 20 m³/s\n",
-    "# # data = [value for value in discharge if value > 20]\n",
-    "\n",
-    "# # all discharge data\n",
-    "# cleaned_discharge = [x for x in discharge if not math.isnan(x)]\n",
-    "# data = cleaned_discharge\n",
-    "\n",
-    "# # # maximum annual discharge\n",
-    "# # data = [x for x in maxdischarge if not math.isnan(x)]\n",
-    "# # # print(data)\n",
-    "\n",
-    "# # Generalized Extreme Value (GEV) distribution \n",
-    "# shape, loc, scale = stats.genextreme.fit(data)\n",
-    "\n",
-    "# # Define threshold value\n",
-    "# threshold = 530\n",
-    "\n",
-    "# # calculate exceedance probability in years\n",
-    "# p = 1 - stats.genextreme.cdf(threshold, shape, loc=loc, scale=scale)\n",
-    "# # calculate return period in years\n",
-    "# T = 1 / (p * 365.25)\n",
-    "# # print(T)\n",
-    "\n",
-    "# # calculate change threshold exceedance at least once in 1000 and 1100 years\n",
-    "# P_1000 = 1 - (1 - p) ** 1000\n",
-    "# P_1100 = 1 - (1 - p) ** 1100\n",
-    "\n",
-    "# # Results 1000\n",
-    "# print(f\"Returnperiod for 530 m³/s: {T:.2f} years\")\n",
-    "# print(f\"Chance of at least one exceedance in 1000 years: {P_1000:.5f} ({P_1000 * 100:.2f}%)\")\n",
-    "# print(f\"Chance of at least one exceedance in 1100 years: {P_1100:.5f} ({P_1100 * 100:.2f}%)\")\n",
-    "\n",
-    "# # expected amount of times threshold is exceeded in 1000 years\n",
-    "# expected_exceedances1000 = 1000 / T if T > 0 else 0\n",
-    "# print(f\"The discharge currently exceeds the threshold of 530 m³/s {expected_exceedances1000:.3f} times in 1000 years\")\n",
-    "\n",
-    "# # expected amount threshold is exceeded in 1100 years\n",
-    "# expected_exceedances1100 = 1100 / T if T > 0 else 0\n",
-    "# print(f\"The discharge currently exceeds the threshold of 530 m³/s {expected_exceedances1100:.3f} times in 1100 years\")"
-   ]
-  },
   {
    "cell_type": "code",
    "execution_count": 27,