Skip to content

GitLab

Commit 740e58fc authored by Oscar Corvi's avatar Oscar Corvi
Browse files

sensitivity analysis regarding the stress factor shape parameter

parent 96db39f4
Pipeline #230378 passed with stage
in 23 seconds
Hide whitespace changes
Inline Side-by-side
.gitattributes
View file @ 740e58fc
......@@ -72,3 +72,4 @@ notebooks/Finished_project/meaning_stress_factor/complete_VS_incomplete.png filt
notebooks/Finished_project/meaning_stress_factor/stress_factor_reconstruct.png filter=lfs diff=lfs merge=lfs -text
notebooks/Finished_project/meaning_stress_factor/constant_VS_fPAR.png filter=lfs diff=lfs merge=lfs -text
notebooks/Finished_project/simpler_model/different_sites.png filter=lfs diff=lfs merge=lfs -text
notebooks/Finished_project/stress_factor_shape/sensitivity_parameters.png filter=lfs diff=lfs merge=lfs -text
notebooks/Finished_project/simpler_model/simpler_model.ipynb
View file @ 740e58fc
......@@ -2,7 +2,7 @@
"cells": [
{
"cell_type": "markdown",
"id": "covered-filling",
"id": "apparent-discretion",
"metadata": {},
"source": [
"# **Location of the stress factor in potential evapo-transpiration models**"
......@@ -10,7 +10,7 @@
},
{
"cell_type": "markdown",
"id": "differential-handling",
"id": "understanding-hacker",
"metadata": {
"slideshow": {
"slide_type": "slide"
......@@ -62,7 +62,7 @@
},
{
"cell_type": "markdown",
"id": "planned-protest",
"id": "extended-romania",
"metadata": {},
"source": [
"# Part II - Functions set up"
......@@ -70,7 +70,7 @@
},
{
"cell_type": "markdown",
"id": "breeding-confusion",
"id": "satisfactory-editor",
"metadata": {},
"source": [
"## Importing relevant packages"
......@@ -79,7 +79,7 @@
{
"cell_type": "code",
"execution_count": 1,
"id": "above-infrared",
"id": "thrown-gamma",
"metadata": {},
"outputs": [
{
......@@ -144,7 +144,7 @@
},
{
"cell_type": "markdown",
"id": "secret-belfast",
"id": "skilled-montana",
"metadata": {},
"source": [
"## Path of the different files (pre-defined python functions, sympy equations, sympy variables)"
......@@ -153,7 +153,7 @@
{
"cell_type": "code",
"execution_count": 2,
"id": "underlying-couple",
"id": "premier-moscow",
"metadata": {
"tags": [
"parameters"
......@@ -180,7 +180,7 @@
},
{
"cell_type": "markdown",
"id": "jewish-malta",
"id": "pediatric-basket",
"metadata": {},
"source": [
"## Importing the sympy variables and equations defined in the theory.ipynb notebook"
......@@ -189,7 +189,7 @@
{
"cell_type": "code",
"execution_count": 3,
"id": "biological-parameter",
"id": "educated-metropolitan",
"metadata": {},
"outputs": [
{
......@@ -286,7 +286,7 @@
},
{
"cell_type": "markdown",
"id": "interested-sperm",
"id": "interested-technical",
"metadata": {},
"source": [
"## Importing the performance assessment functions defined in the analysis_function.py file"
......@@ -295,7 +295,7 @@
{
"cell_type": "code",
"execution_count": 4,
"id": "rubber-consortium",
"id": "colonial-sweden",
"metadata": {},
"outputs": [
{
......@@ -345,7 +345,7 @@
},
{
"cell_type": "markdown",
"id": "photographic-upper",
"id": "material-assumption",
"metadata": {},
"source": [
"## Data import, preprocess and shape for the computations"
......@@ -353,7 +353,7 @@
},
{
"cell_type": "markdown",
"id": "operational-childhood",
"id": "jewish-channels",
"metadata": {},
"source": [
"### Get the different files where data are stored\n",
......@@ -364,7 +364,7 @@
{
"cell_type": "code",
"execution_count": 5,
"id": "impaired-creek",
"id": "funded-telling",
"metadata": {},
"outputs": [
{
......@@ -394,7 +394,7 @@
},
{
"cell_type": "markdown",
"id": "meaning-lighting",
"id": "understanding-productivity",
"metadata": {},
"source": [
"### Define and test a function that process the fPAR data\n",
......@@ -404,7 +404,7 @@
{
"cell_type": "code",
"execution_count": 6,
"id": "hindu-medicare",
"id": "identical-bahrain",
"metadata": {},
"outputs": [],
"source": [
......@@ -443,7 +443,7 @@
{
"cell_type": "code",
"execution_count": 7,
"id": "vital-dinner",
"id": "removable-intranet",
"metadata": {},
"outputs": [
{
......@@ -569,7 +569,7 @@
},
{
"cell_type": "markdown",
"id": "needed-heaven",
"id": "improving-portsmouth",
"metadata": {},
"source": [
"### fPARSet function\n",
......@@ -579,7 +579,7 @@
{
"cell_type": "code",
"execution_count": 8,
"id": "legitimate-academy",
"id": "through-transportation",
"metadata": {},
"outputs": [],
"source": [
......@@ -610,7 +610,7 @@
},
{
"cell_type": "markdown",
"id": "geological-columbia",
"id": "southwest-antibody",
"metadata": {},
"source": [
"### DataChose function\n",
......@@ -631,7 +631,7 @@
{
"cell_type": "code",
"execution_count": 9,
"id": "champion-auckland",
"id": "focused-tracker",
"metadata": {},
"outputs": [],
"source": [
......@@ -724,7 +724,7 @@
},
{
"cell_type": "markdown",
"id": "successful-module",
"id": "rapid-terrorism",
"metadata": {},
"source": [
"## Compile the different functions defined in the symbolic domain\n",
......@@ -733,7 +733,7 @@
},
{
"cell_type": "markdown",
"id": "dated-census",
"id": "annoying-shape",
"metadata": {},
"source": [
"### Water stress functions"
......@@ -742,7 +742,7 @@
{
"cell_type": "code",
"execution_count": 10,
"id": "executive-blackjack",
"id": "excited-globe",
"metadata": {},
"outputs": [],
"source": [
......@@ -763,7 +763,7 @@
},
{
"cell_type": "markdown",
"id": "terminal-front",
"id": "proof-kuwait",
"metadata": {},
"source": [
"### Soil water potential"
......@@ -772,7 +772,7 @@
{
"cell_type": "code",
"execution_count": 11,
"id": "elect-aquarium",
"id": "thermal-accessory",
"metadata": {},
"outputs": [],
"source": [
......@@ -788,7 +788,7 @@
{
"cell_type": "code",
"execution_count": 12,
"id": "laughing-tactics",
"id": "regulated-filing",
"metadata": {},
"outputs": [],
"source": [
......@@ -805,7 +805,7 @@
},
{
"cell_type": "markdown",
"id": "willing-thumb",
"id": "municipal-instruction",
"metadata": {},
"source": [
"### Penman-Monteith"
......@@ -814,7 +814,7 @@
{
"cell_type": "code",
"execution_count": 13,
"id": "generous-fireplace",
"id": "impressive-latex",
"metadata": {},
"outputs": [],
"source": [
......@@ -833,7 +833,7 @@
{
"cell_type": "code",
"execution_count": 14,
"id": "engaging-purse",
"id": "indirect-confidentiality",
"metadata": {},
"outputs": [],
"source": [
......@@ -861,7 +861,7 @@
{
"cell_type": "code",
"execution_count": 15,
"id": "lyric-signature",
"id": "rolled-poland",
"metadata": {},
"outputs": [],
"source": [
......@@ -883,7 +883,7 @@
{
"cell_type": "code",
"execution_count": 16,
"id": "brief-appreciation",
"id": "increasing-practice",
"metadata": {},
"outputs": [],
"source": [
......@@ -905,7 +905,7 @@
{
"cell_type": "code",
"execution_count": 17,
"id": "protecting-exposure",
"id": "sharp-studio",
"metadata": {},
"outputs": [],
"source": [
......@@ -927,7 +927,7 @@
{
"cell_type": "code",
"execution_count": 18,
"id": "convertible-saint",
"id": "underlying-bikini",
"metadata": {},
"outputs": [],
"source": [
......@@ -955,7 +955,7 @@
{
"cell_type": "code",
"execution_count": 19,
"id": "ranking-minneapolis",
"id": "alpine-click",
"metadata": {},
"outputs": [],
"source": [
......@@ -986,7 +986,7 @@
{
"cell_type": "code",
"execution_count": 20,
"id": "israeli-honduras",
"id": "black-postcard",
"metadata": {},
"outputs": [],
"source": [
......@@ -1015,7 +1015,7 @@
{
"cell_type": "code",
"execution_count": 21,
"id": "broken-instruction",
"id": "standard-clearing",
"metadata": {},
"outputs": [],
"source": [
......@@ -1049,7 +1049,7 @@
{
"cell_type": "code",
"execution_count": 22,
"id": "advanced-finnish",
"id": "spanish-marshall",
"metadata": {},
"outputs": [],
"source": [
......@@ -1072,7 +1072,7 @@
{
"cell_type": "code",
"execution_count": 23,
"id": "lined-corps",
"id": "gothic-residence",
"metadata": {},
"outputs": [],
"source": [
......@@ -1090,7 +1090,7 @@
},
{
"cell_type": "markdown",
"id": "pleased-performer",
"id": "incorporated-consolidation",
"metadata": {},
"source": [
"### Assign the different compiled functions to variables functions (create the functions in python)"
......@@ -1099,7 +1099,7 @@
{
"cell_type": "code",
"execution_count": 24,
"id": "documentary-peoples",
"id": "thrown-whale",
"metadata": {},
"outputs": [],
"source": [
......@@ -1127,7 +1127,7 @@
{
"cell_type": "code",
"execution_count": 25,
"id": "visible-vacation",
"id": "bronze-resort",
"metadata": {},
"outputs": [],
"source": [
......@@ -1156,7 +1156,7 @@
},
{
"cell_type": "markdown",
"id": "legitimate-impact",
"id": "settled-lending",
"metadata": {},
"source": [
"## Functions to run the different models\n",
......@@ -1165,7 +1165,7 @@
},
{
"cell_type": "markdown",
"id": "sensitive-buddy",
"id": "passive-malpractice",
"metadata": {},
"source": [
"### Varying surface resistance model"
......@@ -1174,7 +1174,7 @@
{
"cell_type": "code",
"execution_count": 26,
"id": "virgin-tunisia",
"id": "funky-background",
"metadata": {},
"outputs": [],
"source": [
......@@ -1237,7 +1237,7 @@
},
{
"cell_type": "markdown",
"id": "specialized-intelligence",
"id": "threaded-occasions",
"metadata": {},
"source": [
"### Constant surface conductance model"
......@@ -1246,7 +1246,7 @@
{
"cell_type": "code",
"execution_count": 27,
"id": "protected-hardwood",
"id": "temporal-magnet",
"metadata": {},
"outputs": [],
"source": [
......@@ -1310,7 +1310,7 @@
},
{
"cell_type": "markdown",
"id": "leading-gambling",
"id": "caroline-continuity",
"metadata": {},
"source": [
"### Benchmark Penman-Monteith model"
......@@ -1319,7 +1319,7 @@
{
"cell_type": "code",
"execution_count": 28,
"id": "dirty-portrait",
"id": "atmospheric-equilibrium",
"metadata": {},
"outputs": [],
"source": [
......@@ -1378,7 +1378,7 @@
},
{
"cell_type": "markdown",
"id": "emerging-folks",
"id": "individual-tragedy",
"metadata": {},
"source": [
"### Modified version of the PM equation\n",
......@@ -1394,7 +1394,7 @@
{
"cell_type": "code",
"execution_count": 29,
"id": "private-burton",
"id": "elect-confidentiality",
"metadata": {},
"outputs": [],
"source": [
......@@ -1458,7 +1458,7 @@
{
"cell_type": "code",
"execution_count": 30,
"id": "according-eight",
"id": "distinct-provincial",
"metadata": {},
"outputs": [],
"source": [
......@@ -1522,7 +1522,7 @@
},
{
"cell_type": "markdown",
"id": "welcome-caution",
"id": "sharp-exercise",
"metadata": {},
"source": [
"### Priestley and Taylor model"
......@@ -1531,7 +1531,7 @@
{
"cell_type": "code",
"execution_count": 31,
"id": "stuffed-george",
"id": "brutal-plate",
"metadata": {},
"outputs": [],
"source": [
......@@ -1586,7 +1586,7 @@
},
{
"cell_type": "markdown",
"id": "pressing-thought",
"id": "approximate-hawaii",
"metadata": {},
"source": [
"### Inverse modelling\n",
......@@ -1596,7 +1596,7 @@
{
"cell_type": "code",
"execution_count": 32,
"id": "trained-november",
"id": "engaged-devon",
"metadata": {},
"outputs": [],
"source": [
......@@ -1650,7 +1650,7 @@
},
{
"cell_type": "markdown",
"id": "marked-wellington",
"id": "impressive-spain",
"metadata": {},
"source": [
"## Calibration algorithm"
......@@ -1659,7 +1659,7 @@
{
"cell_type": "code",
"execution_count": 33,
"id": "associate-quilt",
"id": "protected-adolescent",
"metadata": {},
"outputs": [],
"source": [
......@@ -1681,7 +1681,7 @@
{
"cell_type": "code",
"execution_count": 34,
"id": "fluid-wireless",
"id": "finite-survey",
"metadata": {
"jupyter": {
"source_hidden": true
......@@ -1816,7 +1816,7 @@
},
{
"cell_type": "markdown",
"id": "stupid-knife",
"id": "incorporated-reflection",
"metadata": {},
"source": [
"Use the global optimizer from the `scipy.optimize` package. Minimize the squared residual :\n",
......@@ -1828,7 +1828,7 @@
{
"cell_type": "code",
"execution_count": 35,
"id": "varied-burning",
"id": "expressed-tuning",
"metadata": {},
"outputs": [],
"source": [
......@@ -1847,7 +1847,7 @@
},
{
"cell_type": "markdown",
"id": "antique-farmer",
"id": "narrow-coordination",
"metadata": {},
"source": [
"# Part III - Experiments"
......@@ -1855,7 +1855,7 @@
},
{
"cell_type": "markdown",
"id": "serial-survey",
"id": "crucial-premium",
"metadata": {},
"source": [
"## One site, one year\n",
......@@ -1865,7 +1865,7 @@
{
"cell_type": "code",
"execution_count": 36,
"id": "surface-chemistry",
"id": "demographic-bennett",
"metadata": {},
"outputs": [
{
......@@ -10552,7 +10552,7 @@
{
"cell_type": "code",
"execution_count": 37,
"id": "breeding-indianapolis",
"id": "developmental-spanking",
"metadata": {},
"outputs": [],
"source": [
......@@ -10562,7 +10562,7 @@
},
{
"cell_type": "markdown",
"id": "varied-medline",
"id": "interpreted-queensland",
"metadata": {},
"source": [
"### Calibrate all the models using the optimized function defined above"
......@@ -10571,7 +10571,7 @@
{
"cell_type": "code",
"execution_count": 38,
"id": "brilliant-model",
"id": "fabulous-utility",
"metadata": {},
"outputs": [
{
......@@ -10593,7 +10593,7 @@
},
{
"cell_type": "markdown",
"id": "terminal-running",
"id": "intelligent-narrow",
"metadata": {},
"source": [
"### Compute the model results with the optimimum parameters"
......@@ -10602,7 +10602,7 @@
{
"cell_type": "code",
"execution_count": 39,
"id": "beneficial-saturday",
"id": "nearby-rhythm",
"metadata": {},
"outputs": [],
"source": [
......@@ -10616,7 +10616,7 @@
},
{
"cell_type": "markdown",
"id": "entertaining-alexandria",
"id": "military-theorem",
"metadata": {},
"source": [
"### Aggregate the results in a single matrix"
......@@ -10625,7 +10625,7 @@
{
"cell_type": "code",
"execution_count": 40,
"id": "suspected-chemical",
"id": "russian-speech",
"metadata": {},
"outputs": [],
"source": [
......@@ -10660,7 +10660,7 @@
},
{
"cell_type": "markdown",
"id": "unique-cosmetic",
"id": "nutritional-introduction",
"metadata": {},
"source": [
"### Plot the time serie"
......@@ -10669,7 +10669,7 @@
{
"cell_type": "code",
"execution_count": 41,
"id": "opposed-typing",
"id": "responsible-grave",
"metadata": {
"slideshow": {
"slide_type": "slide"
......@@ -10746,7 +10746,7 @@
},
{
"cell_type": "markdown",
"id": "floating-conjunction",
"id": "caroline-pursuit",
"metadata": {},
"source": [
"### Compute the model performance \n",
......@@ -10756,8 +10756,12 @@
{
"cell_type": "code",
"execution_count": 42,
"id": "congressional-seven",
"metadata": {},
"id": "liked-enough",
"metadata": {
"jupyter": {
"source_hidden": true
}
},
"outputs": [
{
"name": "stdout",
......@@ -10799,7 +10803,7 @@
},