class: center, middle, inverse, title-slide .title[ # Découverte de <strong>airGRiwrm</strong> sur la Sèvre-Nantaise ] .subtitle[ ## Gestion intégrée de la ressource en eau avec <strong>airGR</strong> ] .author[ ### David Dorchies ] .institute[ ### UMR G-EAU ] .date[ ### Montpellier / Clisson, le 6 juillet 2023 (MAJ 30/07/2024) ] --- # Programme de la formation - Introduction au modèle semi-distribué (5 mn) - Présentation du package **airGRiwrm** (10 mn) - Installation du package **airGRiwrm** (5 mn) - Présentation du jeu de données de la formation (15 mn) - Créer un réseau semi-distribué avec **airGRiwrm** (15 mn) - Calage d'un réseau semi-distribué de modèles hydrologiques (60 mn) - Modélisation du barrage de la Bultière (30 mn) - Régulation du barrage de la Bultière (30 mn) --- class: inverse, middle, center # Le modèle semi-distribué --- # Un modèle semi-distribué : c'est quoi ? .pull-left[ - Découpage du bassin versant en sous-bassins en fonction de la position des stations hydrométriques exutoires de chaque sous-bassin - Le ruissellement de chaque sous-bassin est modélisé par un modèle global (GR...) - Le modèle semi-distribué relie les stations hydrométrique grâce à un routage hydraulique de type "lag" <img src="img/scheme_bvi.png" width="50%" style="display: block; margin: auto;" /> ] .pull-right[ .vertical-center[  ]] --- # Intérêt du modèle hydrologique semi-distribué .left-column[ Au delà de la simulation de débits ayant une cohérence amont-aval dans le modèle, la semi-distribution et le routage hydraulique permettent l'introduction d'influences humaines&nbsp;: - Prélèvements de débit - Réservoirs - Transferts inter-bassins *Exemple ci-contre avec le bassin versant de l'Aude où sont modélisés 7 réservoirs, des transferts inter-bassins et des canaux de navigation.* ] .right-column[ .center[  ] ] --- class: inverse, middle, center  # *airGRiwrm*: extension de *airGR* pour un réseau de modèles semi-distribués --- # Introduction au package **airGRiwrm** .pull-left[ **airGRiwrm** complète la galaxie **airGR** en automatisant la semi-distribution et en ajoutant des fonctionnalités pour simuler les influences humaines. Principalement utilisé pour les études d'impact et d'adaptation aux changements globaux de bassins versants anthropisés. - Le package airGRiwrm est en développement depuis mai 2020 - Diffusé sur [le dépôt officiel "CRAN" en août 2022 *(v0.6.1)*](https://cloud.r-project.org/package=airGRiwrm) - Logiciel gratuit et open-source (Github : https://github.com/inrae/airGRiwrm) .center[  ] ] .pull-right[ .center[  ] Documentation avec des exemples d'utilisation: - https://airgriwrm.g-eau.fr/ *(v0.6.1)* - https://airgriwrm.g-eau.fr/dev *(v0.6.1.9000 en développement)* ] --- class: hide_logo # *airGRiwrm*: extension de *airGR* pour un réseau de modèles semi-distribués **airGRiwrm** redéfinit les fonctions de **airGR** pour automatiser la gestion d'un réseau semi-distribué: ``` r library(airGRiwrm) ``` Les fonctions originales de **airGR** restent accessibles (Redirection avec [les classes S3](http://www.duclert.org/r-divers/classes-S3-R.php)). --- # Procédure de calage d'un réseau avec *airGRiwrm* Par rapport à *airGR*, on ajoute une étape préalable de description du réseau&nbsp;: <div id="htmlwidget-a660d75d2839daf28b87" style="width:100%;height:741.6px;" class="DiagrammeR html-widget"></div> <script type="application/json" data-for="htmlwidget-a660d75d2839daf28b87">{"x":{"diagram":"\ngraph LR\n subgraph Description du réseau\n CreateGRiwrm --> griwrm((griwrm))\n end\n subgraph Mêmes procédures que airGR\n griwrm --> CreateInputsModel\n CreateInputsModel --> IM((IM))\n IM --> CreateRunOptions\n CreateRunOptions --> RO((RO))\n IM --> CreateInputsCrit\n RO --> CreateInputsCrit\n CreateInputsCrit --> IC((IC))\n CreateCalibOptions --> CO((CO))\n IM --> Calibration\n RO --> Calibration\n IC --> Calibration\n CO --> Calibration\n end\n style griwrm fill:#fcc\n style IM fill:#fcc\n style RO fill:#fcc\n style IC fill:#fcc\n style CO fill:#fcc\n"},"evals":[],"jsHooks":[]}</script> --- # Les fonctionnalités de *airGRiwrm* ## Version officielle diffusée sur le CRAN *(v0.6.1)* - Automatisation des fonctionnalités de *airGR* pour un réseau semi-distribué via la description d'un réseau de "noeuds": - Modèles GR (pas de temps journalier et horaire) et transfert hydraulique (Lag) - Calage et run des simulations - Régularisation des paramètres (Méthode de De Lavenne) - Noeuds "DirectInjection" pour introduire des influences humaines dans le réseau via des chroniques prédéfinies de débits prélevés ou lâchés ## Version en développement *(v0.6.1.9000)* - Ajout de noeuds de type "Ungauged" pour simuler les débits du cours d'eau sur des points non jaugés - Ajout de noeuds de type "Reservoir" pour simuler l'évolution du volume d'un réservoir - Ajout de noeuds de type "Diversion" pour transférer des débits entre les noeuds du réseau --- class: inverse, middle, center # Installation du package *airGRiwrm*  --- # Installation du package *airGRiwrm* ### La version diffusée sur le CRAN ``` r install.packages("airGRiwrm") ``` ### La version en cours de développement Installation de **airGRiwrm** branche "dev"&nbsp;: ``` r install.packages("remotes") remotes::install_github("inrae/airGRiwrm", "dev") ``` --- class: inverse, middle, center # Le jeu de données de la formation: le cas d'étude "Maine" --- # Présentation du cas d'étude "Maine" <div class="leaflet html-widget html-fill-item" id="htmlwidget-2f9b76f96433d6c04656" style="width:100%;height:540px;"></div> <script type="application/json" data-for="htmlwidget-2f9b76f96433d6c04656">{"x":{"options":{"crs":{"crsClass":"L.CRS.EPSG3857","code":null,"proj4def":null,"projectedBounds":null,"options":{}}},"calls":[{"method":"createMapPane","args":["tmap401",401]},{"method":"createMapPane","args":["tmap402",402]},{"method":"createMapPane","args":["tmap403",403]},{"method":"addProviderTiles","args":["OpenStreetMap",null,"OpenStreetMap",{"minZoom":0,"maxZoom":18,"tileSize":256,"subdomains":"abc","errorTileUrl":"","tms":false,"noWrap":false,"zoomOffset":0,"zoomReverse":false,"opacity":1,"zIndex":1,"detectRetina":false,"pane":"tilePane"}]},{"method":"addCircles","args":[[47.05552408903289,46.90370052398146,46.92258793435403,46.91795301499165,46.91795301499165,46.95218178641273,46.93950367953726,46.88800866619889,46.88074155217706,46.84514520224651,46.84733225195411],[-1.370098590461308,-1.157340377029569,-1.274327611501934,-1.214063263499455,-1.214063263499455,-1.282695202306296,-1.243583949451833,-1.086383546421043,-1.07920635863741,-1.237362066078344,-1.226064249122383],[551.380692521076,551.380692521076,551.380692521076,551.380692521076,551.380692521076,551.380692521076,551.380692521076,551.380692521076,551.380692521076,551.380692521076,551.380692521076],["Maine_Remouille","GrandeMaine_StFulgent","PetiteMaine_StGeorges","P_BULTIERE","R_BULTIERE","P_MAINE","R_MAINE","P_GRANDEMAINE","R_GRANDEMAINE","P_PETITEMAINE","R_PETITEMAINE"],"sfReseau",{"interactive":true,"className":"","pane":"tmap401","stroke":true,"color":"#666666","weight":1,"opacity":0.5,"fill":true,"fillColor":["#A0F0A0","#A0F0A0","#A0F0A0","#F0A0A0","#F0A0A0","#F0A0A0","#F0A0A0","#F0A0A0","#F0A0A0","#F0A0A0","#F0A0A0"],"fillOpacity":[1,1,1,1,1,1,1,1,1,1,1]},["<style> div.leaflet-popup-content {width:auto !important;overflow-y:auto; overflow-x:hidden;}<\/style><div style=\"max-height:25em;padding-right:0px;\"><table>\n\t\t\t <thead><tr><th colspan=\"2\"><b>Maine_Remouille<\/b><\/th><\/thead><\/tr><tr><td style=\"color: #888888;\"><nobr>type<\/nobr><\/td><td align=\"right\"><nobr>hydrologie<\/nobr><\/td><\/tr><\/table><\/div>","<style> div.leaflet-popup-content {width:auto !important;overflow-y:auto; overflow-x:hidden;}<\/style><div style=\"max-height:25em;padding-right:0px;\"><table>\n\t\t\t <thead><tr><th colspan=\"2\"><b>GrandeMaine_StFulgent<\/b><\/th><\/thead><\/tr><tr><td style=\"color: #888888;\"><nobr>type<\/nobr><\/td><td align=\"right\"><nobr>hydrologie<\/nobr><\/td><\/tr><\/table><\/div>","<style> div.leaflet-popup-content {width:auto !important;overflow-y:auto; overflow-x:hidden;}<\/style><div style=\"max-height:25em;padding-right:0px;\"><table>\n\t\t\t <thead><tr><th colspan=\"2\"><b>PetiteMaine_StGeorges<\/b><\/th><\/thead><\/tr><tr><td style=\"color: #888888;\"><nobr>type<\/nobr><\/td><td align=\"right\"><nobr>hydrologie<\/nobr><\/td><\/tr><\/table><\/div>","<style> div.leaflet-popup-content {width:auto !important;overflow-y:auto; overflow-x:hidden;}<\/style><div style=\"max-height:25em;padding-right:0px;\"><table>\n\t\t\t <thead><tr><th colspan=\"2\"><b>P_BULTIERE<\/b><\/th><\/thead><\/tr><tr><td style=\"color: #888888;\"><nobr>type<\/nobr><\/td><td align=\"right\"><nobr>influences<\/nobr><\/td><\/tr><\/table><\/div>","<style> div.leaflet-popup-content {width:auto !important;overflow-y:auto; overflow-x:hidden;}<\/style><div style=\"max-height:25em;padding-right:0px;\"><table>\n\t\t\t <thead><tr><th colspan=\"2\"><b>R_BULTIERE<\/b><\/th><\/thead><\/tr><tr><td style=\"color: #888888;\"><nobr>type<\/nobr><\/td><td align=\"right\"><nobr>influences<\/nobr><\/td><\/tr><\/table><\/div>","<style> div.leaflet-popup-content {width:auto !important;overflow-y:auto; overflow-x:hidden;}<\/style><div style=\"max-height:25em;padding-right:0px;\"><table>\n\t\t\t <thead><tr><th colspan=\"2\"><b>P_MAINE<\/b><\/th><\/thead><\/tr><tr><td style=\"color: #888888;\"><nobr>type<\/nobr><\/td><td align=\"right\"><nobr>influences<\/nobr><\/td><\/tr><\/table><\/div>","<style> div.leaflet-popup-content {width:auto !important;overflow-y:auto; overflow-x:hidden;}<\/style><div style=\"max-height:25em;padding-right:0px;\"><table>\n\t\t\t <thead><tr><th colspan=\"2\"><b>R_MAINE<\/b><\/th><\/thead><\/tr><tr><td style=\"color: #888888;\"><nobr>type<\/nobr><\/td><td align=\"right\"><nobr>influences<\/nobr><\/td><\/tr><\/table><\/div>","<style> div.leaflet-popup-content {width:auto !important;overflow-y:auto; overflow-x:hidden;}<\/style><div style=\"max-height:25em;padding-right:0px;\"><table>\n\t\t\t <thead><tr><th colspan=\"2\"><b>P_GRANDEMAINE<\/b><\/th><\/thead><\/tr><tr><td style=\"color: #888888;\"><nobr>type<\/nobr><\/td><td align=\"right\"><nobr>influences<\/nobr><\/td><\/tr><\/table><\/div>","<style> div.leaflet-popup-content {width:auto !important;overflow-y:auto; overflow-x:hidden;}<\/style><div style=\"max-height:25em;padding-right:0px;\"><table>\n\t\t\t <thead><tr><th colspan=\"2\"><b>R_GRANDEMAINE<\/b><\/th><\/thead><\/tr><tr><td style=\"color: #888888;\"><nobr>type<\/nobr><\/td><td align=\"right\"><nobr>influences<\/nobr><\/td><\/tr><\/table><\/div>","<style> div.leaflet-popup-content {width:auto !important;overflow-y:auto; overflow-x:hidden;}<\/style><div style=\"max-height:25em;padding-right:0px;\"><table>\n\t\t\t <thead><tr><th colspan=\"2\"><b>P_PETITEMAINE<\/b><\/th><\/thead><\/tr><tr><td style=\"color: #888888;\"><nobr>type<\/nobr><\/td><td align=\"right\"><nobr>influences<\/nobr><\/td><\/tr><\/table><\/div>","<style> div.leaflet-popup-content {width:auto !important;overflow-y:auto; overflow-x:hidden;}<\/style><div style=\"max-height:25em;padding-right:0px;\"><table>\n\t\t\t <thead><tr><th colspan=\"2\"><b>R_PETITEMAINE<\/b><\/th><\/thead><\/tr><tr><td style=\"color: #888888;\"><nobr>type<\/nobr><\/td><td align=\"right\"><nobr>influences<\/nobr><\/td><\/tr><\/table><\/div>"],null,["Maine_Remouille","GrandeMaine_StFulgent","PetiteMaine_StGeorges","P_BULTIERE","R_BULTIERE","P_MAINE","R_MAINE","P_GRANDEMAINE","R_GRANDEMAINE","P_PETITEMAINE","R_PETITEMAINE"],{"interactive":false,"permanent":false,"direction":"auto","opacity":1,"offset":[0,0],"textsize":"10px","textOnly":false,"className":"","sticky":true},null,null]},{"method":"addLegend","args":[{"colors":["#A0F0A0","#F0A0A0"],"labels":["hydrologie","influences"],"na_color":null,"na_label":"NA","opacity":1,"position":"topright","type":"unknown","title":"type","extra":null,"layerId":"legend401","className":"info legend sfReseau","group":"sfReseau"}]},{"method":"addMarkers","args":[[47.05552408903289,46.90370052398146,46.92258793435403,46.91795301499165,46.91795301499165,46.95218178641273,46.93950367953726,46.88800866619889,46.88074155217706,46.84514520224651,46.84733225195411],[-1.370098590461308,-1.157340377029569,-1.274327611501934,-1.214063263499455,-1.214063263499455,-1.282695202306296,-1.243583949451833,-1.086383546421043,-1.07920635863741,-1.237362066078344,-1.226064249122383],{"iconUrl":{"data":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAAEAAAABCAQAAAC1HAwCAAAAC0lEQVR4nGP6zwAAAgcBApocMXEAAAAASUVORK5CYII=","index":0},"iconWidth":1,"iconHeight":1},["Maine_Remouille.1","GrandeMaine_StFulgent.1","PetiteMaine_StGeorges.1","P_BULTIERE.1","R_BULTIERE.1","P_MAINE.1","R_MAINE.1","P_GRANDEMAINE.1","R_GRANDEMAINE.1","P_PETITEMAINE.1","R_PETITEMAINE.1"],"sfReseau",{"interactive":true,"draggable":false,"keyboard":true,"title":"","alt":"","zIndexOffset":0,"opacity":1,"riseOnHover":false,"riseOffset":250,"pane":"tmap402"},null,null,null,null,["Maine_Remouille","GrandeMaine_StFulgent","PetiteMaine_StGeorges","P_BULTIERE","R_BULTIERE","P_MAINE","R_MAINE","P_GRANDEMAINE","R_GRANDEMAINE","P_PETITEMAINE","R_PETITEMAINE"],{"interactive":false,"permanent":true,"direction":"center","offset":[0,0],"textsize":"12px","textOnly":true,"style":{"color":"#000000"},"className":"","sticky":true},null]},{"method":"addPolygons","args":[[[[{"lng":[-1.370813045430074,-1.369498174999511,-1.369436120191693,-1.370750968675223,-1.370813045430074],"lat":[47.05594683806952,47.05598925544965,47.05509065313351,47.05504823645349,47.05594683806952]}]],[[{"lng":[-1.372065815081659,-1.370750968675223,-1.370626821396708,-1.371941623911479,-1.372065815081659],"lat":[47.0550058047406,47.05504823645349,47.05325103250799,47.0532086021957,47.0550058047406]}]],[[{"lng":[-1.370626821396708,-1.369312016805049,-1.369249968226015,-1.370564750872835,-1.370626821396708],"lat":[47.05325103250799,47.05329344778791,47.05239484475887,47.05235243017899,47.05325103250799]}]],[[{"lng":[-1.402682374241768,-1.400053437632155,-1.399990889614317,-1.402619782358369,-1.402682374241768],"lat":[47.04050414490015,47.04058969366909,47.03969110443916,47.03960555708209,47.04050414490015]}]],[[{"lng":[-1.371879531442871,-1.370564750872835,-1.370502682425745,-1.369187921723034,-1.369249968226015,-1.367935183503142,-1.367563067276722,-1.364933754916242,-1.364747856448377,-1.36211866743172,-1.36218058763715,-1.356922096995959,-1.356798434165974,-1.354169264186543,-1.354107479750042,-1.351478345378908,-1.351416606884017,-1.350102058537493,-1.350040344044703,-1.348725815570499,-1.348664125078348,-1.347349616476673,-1.347164623186073,-1.345850178326226,-1.345726903645164,-1.344412500589828,-1.344289278018872,-1.342974916766271,-1.342913330502668,-1.341598989120912,-1.341537426847311,-1.340223105336609,-1.340161567051548,-1.338847265412105,-1.338785751114125,-1.337471469346147,-1.337348490777984,-1.336034250809281,-1.335972786536805,-1.333344344287009,-1.333282925922747,-1.330654519309312,-1.330593146850654,-1.329278962390601,-1.329217613909634,-1.327903449323117,-1.327842124818382,-1.325213833342729,-1.325152554736812,-1.321210167927638,-1.320965337064723,-1.319651291723727,-1.319590111045729,-1.318276085579708,-1.318153774200741,-1.311583835409881,-1.311522792382933,-1.308894846411497,-1.30883384926301,-1.306205938970538,-1.306023091266595,-1.304709198815465,-1.304648275572692,-1.302020528392382,-1.301959651017635,-1.299331939526648,-1.299271108017316,-1.296643432220847,-1.296521862960092,-1.29520806584139,-1.295147306172531,-1.288578399688821,-1.288639049790348,-1.284697615819956,-1.284637031460048,-1.282009442511961,-1.282130525599454,-1.275561298672167,-1.275621733672583,-1.270366228089386,-1.270426577442034,-1.26779876870576,-1.267617858201362,-1.266304016558728,-1.2662437390051,-1.263616093505372,-1.263555861796696,-1.260928252071538,-1.2608680662052,-1.258240492259797,-1.258180352233182,-1.254239041978343,-1.254118899434443,-1.252805169157845,-1.253105431010216,-1.249163899422002,-1.24898393352358,-1.246356367260053,-1.24641631005786,-1.245102502004799,-1.245042581121807,-1.243728792980212,-1.243788691948289,-1.241161065827838,-1.241101210689892,-1.238473620394931,-1.238354003779112,-1.234412734827269,-1.234352995258646,-1.231725516463464,-1.231785212204931,-1.230471447899617,-1.23053112372248,-1.229217335507484,-1.229157681599179,-1.227843913304349,-1.227784283304368,-1.225156784562597,-1.224978037999358,-1.22235066276608,-1.222469737499783,-1.219842266654504,-1.219901763179631,-1.213332941699686,-1.213154786791971,-1.215782189864577,-1.215663342205702,-1.213036026779867,-1.212739161486106,-1.211425610341705,-1.211366265144087,-1.210052733923072,-1.209934093291098,-1.207307112525859,-1.207247838994428,-1.20593436754476,-1.205875117899141,-1.204561666373653,-1.204502440612388,-1.203189009011279,-1.203070607233286,-1.201757217459445,-1.201698041440431,-1.199071299763949,-1.198953041263351,-1.197639711262538,-1.197462404030647,-1.196149137752729,-1.19609006119049,-1.194776814836251,-1.194717762145677,-1.190778076937053,-1.190660108855011,-1.18934692030189,-1.189170048646322,-1.186543796999608,-1.186484887508416,-1.185171780625447,-1.185054011333679,-1.18111481023763,-1.181055994218241,-1.178429894109833,-1.178488666345358,-1.175862514597291,-1.175803786146024,-1.173177670332645,-1.173118987626964,-1.17180594867025,-1.171747289816424,-1.163869146287508,-1.163810620737476,-1.161184601023659,-1.161243082791646,-1.159930048112075,-1.159871588235259,-1.150680444063272,-1.150622139375637,-1.149309132872004,-1.149250852023022,-1.147937865462552,-1.147879608450767,-1.146566641833653,-1.146391948142467,-1.145079045245378,-1.145137252642283,-1.143824325912464,-1.143882513365527,-1.132065888396184,-1.132123880784251,-1.129497899010451,-1.129555849553932,-1.128242833870885,-1.128126978502146,-1.12681400466553,-1.126756101774396,-1.125443147893083,-1.125385268825937,-1.124072334900124,-1.123898774965207,-1.122585904770488,-1.122528077211216,-1.119902374799764,-1.119844592944535,-1.118531760728819,-1.118474002689469,-1.117161190429587,-1.11710345620466,-1.114477869668892,-1.11430480986744,-1.112992079359135,-1.112934418495048,-1.111621707942717,-1.111275914483763,-1.112588493757672,-1.112530846375066,-1.111218288978347,-1.111103043741856,-1.109790528172713,-1.109675334384176,-1.108362860640812,-1.108305288505884,-1.106992834713389,-1.106935286374447,-1.105622852533012,-1.105565327988602,-1.102940498287741,-1.102883019408353,-1.101570623549942,-1.101513168461543,-1.100200792555593,-1.100258225772583,-1.097633424463934,-1.097576034990051,-1.09363888426854,-1.093581562324524,-1.08833208871537,-1.08838932317378,-1.083139731493812,-1.083082584522537,-1.079145436335386,-1.078974203708568,-1.080286522732697,-1.080229427133451,-1.081541722381546,-1.08148460682055,-1.082796878292801,-1.082739742771513,-1.08536423625598,-1.085307058909116,-1.086619280920893,-1.086504888224031,-1.087817064594812,-1.087702635809513,-1.089014766541861,-1.088957533152087,-1.090269640110559,-1.090212386769851,-1.094148630586237,-1.094091313573996,-1.095403369158467,-1.095460708031591,-1.09939692887715,-1.099454335249128,-1.103390604427057,-1.103333132471685,-1.10595725860778,-1.106014774284933,-1.107326856333466,-1.107269318795567,-1.108581377062211,-1.108523819585532,-1.109835854070479,-1.109720701164874,-1.111032690010041,-1.110917501079298,-1.114853324931714,-1.114795667782944,-1.117419497015993,-1.117304101074847,-1.118615969093306,-1.118442824069924,-1.119754624603925,-1.119812359536754,-1.122435998525261,-1.122378219889145,-1.126313598349125,-1.126197915761335,-1.127509661020089,-1.127393942458133,-1.12870564208823,-1.128647763856807,-1.133894455660471,-1.133952421280848,-1.136575799270326,-1.136517789956053,-1.137829454199739,-1.137771424975702,-1.14039470395937,-1.140336632981864,-1.141648247721413,-1.141590156837871,-1.142901747794159,-1.142843637006035,-1.144155204179257,-1.14403894473693,-1.145350466285418,-1.145292317633298,-1.150538297048278,-1.150421830847643,-1.151733277165025,-1.151675025141423,-1.152986447676277,-1.152928175759078,-1.154239574511599,-1.154122992838224,-1.155434345972602,-1.155376036218505,-1.15799869297845,-1.157823644388742,-1.156512382482953,-1.156279098945296,-1.154967922417993,-1.154851335982907,-1.153540201168667,-1.153423666168771,-1.15473475732769,-1.154676470917291,-1.153365401584944,-1.153190619496775,-1.154501623353754,-1.15426852048996,-1.156890347790537,-1.156948664232951,-1.160881456054493,-1.160939839907948,-1.16487267966466,-1.164931130932726,-1.166242095441355,-1.166300570482303,-1.167611554862819,-1.167670053678101,-1.168981057930308,-1.169039580521373,-1.171661626815495,-1.171720195008375,-1.178275385657337,-1.178334064929382,-1.187511402363986,-1.18762907233883,-1.196806619138743,-1.196865609852158,-1.198176701925049,-1.198294732910015,-1.200916998453915,-1.20097606055524,-1.20228721219898,-1.202346298099176,-1.204968639127476,-1.205027770656647,-1.208961332885002,-1.209079732840062,-1.210390959958584,-1.210450184729738,-1.211761431702378,-1.21182068028091,-1.213131947107469,-1.213191219494831,-1.215813790876704,-1.215754474826265,-1.21706573571705,-1.217243755303064,-1.218555079709857,-1.218614445365878,-1.219925789624976,-1.219985179096256,-1.222607905335946,-1.222667340458558,-1.223978722438129,-1.22403818137953,-1.225349583209985,-1.225290102433434,-1.229224230506502,-1.229283776787903,-1.230595170673642,-1.230654740777788,-1.231966154510913,-1.232025748439255,-1.235960043210534,-1.236019704637684,-1.23733115408332,-1.23739083933882,-1.240013775927623,-1.240073506850609,-1.24138499399275,-1.241444748747617,-1.24275625573367,-1.242816034321867,-1.244127561151636,-1.244247167993644,-1.245558736507113,-1.245618564763201,-1.24693015312061,-1.247109715397909,-1.249733017171595,-1.249792918947027,-1.251104588678472,-1.251164514297598,-1.252476203872811,-1.252536153337082,-1.253847862755862,-1.253907836066731,-1.255219565328875,-1.255279562487795,-1.256591311593105,-1.256651332601528,-1.257963101549802,-1.25802314640918,-1.259334935200219,-1.259515147351682,-1.260826999683059,-1.260947194845562,-1.266194758901899,-1.266134572915658,-1.270070159311788,-1.269949662272944,-1.271261476685443,-1.27120120933522,-1.269889416770441,-1.26976893179877,-1.268457180916309,-1.268336747678974,-1.274895263523799,-1.27477461989027,-1.282644510328888,-1.282523612649699,-1.283835210315535,-1.283774742662211,-1.285086316469645,-1.285025828995193,-1.28633737894443,-1.286276871650303,-1.287588397741539,-1.287406822464676,-1.283472434666931,-1.283351525593559,-1.282040102628829,-1.281919245293546,-1.283230624590256,-1.283170177114629,-1.285792885989838,-1.285732396867335,-1.288355054004199,-1.288294523238024,-1.290917128632999,-1.290674851133613,-1.291986063483611,-1.291864893206534,-1.293176059881501,-1.293115455949651,-1.294426598776926,-1.294365975042423,-1.295677094022205,-1.295616450486499,-1.292994254152172,-1.292933656286341,-1.291622576909371,-1.291562002888029,-1.290250943312768,-1.290190393134469,-1.291501430887667,-1.291319727011435,-1.292630697282221,-1.292388366462905,-1.293699247440897,-1.293456861695986,-1.292146067977897,-1.291964365329641,-1.293275093610409,-1.293153924991417,-1.294464607630217,-1.294404004542394,-1.29571466335042,-1.295654040475314,-1.296964675452769,-1.296904032791827,-1.298214643938912,-1.298153981493583,-1.302085737371838,-1.302025011531592,-1.307267237121252,-1.307206426086435,-1.309827483105836,-1.30976663049122,-1.31107713415323,-1.310955391406751,-1.312265849433435,-1.312022301097199,-1.313332669891943,-1.313149964564418,-1.314460265935129,-1.314277513497109,-1.312967277512944,-1.312845496297303,-1.311535301870386,-1.311474436100587,-1.310164261437052,-1.31004257956885,-1.307422311324221,-1.307057575866281,-1.305747569438054,-1.305626050736546,-1.306936013592259,-1.306814459422866,-1.308124376683003,-1.308063580853664,-1.309373474303601,-1.309312658718036,-1.310622528357974,-1.310500859705109,-1.315740143727529,-1.315618308946205,-1.318237851658205,-1.318298811592114,-1.322228175756411,-1.322350232400679,-1.324969885460118,-1.324725606075232,-1.327345076760382,-1.327283968448636,-1.328593678959419,-1.32822694120976,-1.325607779445055,-1.325546707160111,-1.322927580798187,-1.322805529418236,-1.324115050074261,-1.323992963289488,-1.322683486172785,-1.322622467600244,-1.333098053391625,-1.333036862712405,-1.33434628000093,-1.334285069594385,-1.335594463076977,-1.33565569524974,-1.338274519625319,-1.338458358991712,-1.339767833419107,-1.339890446745808,-1.34119996266841,-1.341383963294684,-1.342693542485739,-1.342754901888191,-1.345374097679089,-1.345435502672065,-1.351983565186528,-1.352045081091728,-1.353354709219962,-1.353416248948949,-1.358654828038343,-1.358593201216153,-1.361212434894088,-1.361274105262263,-1.365203005680921,-1.365326483226461,-1.364016808262735,-1.364201974687951,-1.362892232340307,-1.362953936819534,-1.361644170638093,-1.361705855396153,-1.360396065380706,-1.360704411113852,-1.362014310030133,-1.362076007126385,-1.363385925768277,-1.363694551024472,-1.362384523460575,-1.362446232898719,-1.363756282249186,-1.364003227729717,-1.365313362175221,-1.36537512548137,-1.367995431773913,-1.368057240720667,-1.370677582351111,-1.370739436941055,-1.373359813904234,-1.373421714139951,-1.376042126430702,-1.376104072314775,-1.377414297154465,-1.377476266896811,-1.380096753962032,-1.380158769357385,-1.381469031581881,-1.381655155559513,-1.382965481098995,-1.383027548358531,-1.381717201023074,-1.381903349828256,-1.383213762555979,-1.383337915704345,-1.384648369960277,-1.384710471440004,-1.387331417338326,-1.387393564490698,-1.386083070776711,-1.386393728615636,-1.38770433134413,-1.387766490933856,-1.389077113394169,-1.389139296862842,-1.390449939054766,-1.390512146403837,-1.391822808327164,-1.391947272864727,-1.393257976326806,-1.393382492785016,-1.392071745705658,-1.392133985240146,-1.390823214274909,-1.390885434074995,-1.389574639223671,-1.389636839287908,-1.388326020550292,-1.388388200877231,-1.38576651355543,-1.385828652330479,-1.384517783746707,-1.384766272317187,-1.386077228161184,-1.386139377302005,-1.380895454225004,-1.380957518168299,-1.379646510402815,-1.379708554598049,-1.378397522943786,-1.378521573904663,-1.375899417108583,-1.376147369287536,-1.374836200497387,-1.37508409847105,-1.379017860537715,-1.379079905682378,-1.380391177391655,-1.380453246434314,-1.381764537899833,-1.381826630841937,-1.383137942063491,-1.383324298811445,-1.384635673448344,-1.384759965308344,-1.386071381533574,-1.386320077476606,-1.385008573918249,-1.385070731256727,-1.383759203788666,-1.383821341365827,-1.382509789987856,-1.38257190780225,-1.383883481016975,-1.383945622742214,-1.385257215720225,-1.385319381357759,-1.384007766541646,-1.384069912415376,-1.385381549070315,-1.385692418762086,-1.389627643885184,-1.389752137350708,-1.391063918589985,-1.391188464056796,-1.392500286906269,-1.392811796192098,-1.3967475799616,-1.396997087034366,-1.400933114259522,-1.400995561789208,-1.402307588546629,-1.402370060013501,-1.403682106537097,-1.403744601942608,-1.405056668232173,-1.405119187577774,-1.403807099434869,-1.403932100680064,-1.401307830716069,-1.401370290757859,-1.402682448638753,-1.402994889425072,-1.40430715450854,-1.404369670786263,-1.405681955641919,-1.405807038182794,-1.403182378939605,-1.403244879618866,-1.401932525003724,-1.402369935154897,-1.403682442826262,-1.404057578439129,-1.397494363112424,-1.397556783661281,-1.39624411246475,-1.39611931719735,-1.394806687661753,-1.394744315026133,-1.389493863543077,-1.389431580475504,-1.390744174597578,-1.390681871739854,-1.388056725156543,-1.388118984273002,-1.386806385990799,-1.386868625315445,-1.3842433787715,-1.384305576430426,-1.3829929281691,-1.382930752382729,-1.381618123917062,-1.381680277830774,-1.380367625416173,-1.380429759532668,-1.379117083168928,-1.379179197486751,-1.377866497173665,-1.377928591691361,-1.376615867428717,-1.376740018934728,-1.374114476683129,-1.374176511791717,-1.372863715679161,-1.372925730981803,-1.368987264580166,-1.369111170122613,-1.367798300088344,-1.367984108281357,-1.370609977570474,-1.370733947580842,-1.372046922886139,-1.372108932882548,-1.376047912016995,-1.3761720695383,-1.377485102206069,-1.377547205964549,-1.380173308846211,-1.380297610137827,-1.381610702238101,-1.381735055614335,-1.380421919737034,-1.380670563861554,-1.381983787301643,-1.382108165614382,-1.387361213725927,-1.387423493574789,-1.38611021277358,-1.386234734926537,-1.383608079503383,-1.383670299909445,-1.384983650557086,-1.385232640569502,-1.38391920233077,-1.384105900986614,-1.38279239496658,-1.382916827334135,-1.384230377158667,-1.384541604000102,-1.383227944651572,-1.383290174355643,-1.380662805604029,-1.380787183620645,-1.379473452309218,-1.379535622526769,-1.378221867227294,-1.378470481238521,-1.37978432418365,-1.37984650479506,-1.378532639936693,-1.378656963568208,-1.379970872255236,-1.38003305910421,-1.381346987627269,-1.38140919847169,-1.382723146830567,-1.382785381671895,-1.384099349866384,-1.384161608706076,-1.388103566547142,-1.388165893222197,-1.389479896919632,-1.389542247597219,-1.390856271127972,-1.391043401430395,-1.392357488634171,-1.392419891494694,-1.393733998532618,-1.393796425401722,-1.395110552273586,-1.395235456120441,-1.399177955737064,-1.39949058078304,-1.398176306709523,-1.398301327511984,-1.396987007492676,-1.39711199280056,-1.395797626832851,-1.395860100692037,-1.398488874397458,-1.398551394208176,-1.399865799855225,-1.399990889614317,-1.397361988501355,-1.397486998895942,-1.394858001685648,-1.395045401347777,-1.392416264165815,-1.392603551009673,-1.388659632133301,-1.388721999439938,-1.382148649892886,-1.382210909582801,-1.379581511296202,-1.379643729185645,-1.374384819870964,-1.374446952070881,-1.371817441051777,-1.371879531442871],"lat":[47.052310000567,47.05235243017899,47.0514538276128,47.05149624149266,47.05239484475887,47.05243724430659,47.04704561696083,47.04713036256718,47.04443454151157,47.04451922280189,47.04541783145348,47.04558701646486,47.04378979334206,47.04387429287569,47.04297567956733,47.04306011758705,47.04216150265052,47.04220369841915,47.04130508255185,47.04134726259465,47.04044864579709,47.04049081011455,47.03779495623189,47.03783710343318,47.03603986495389,47.03608199573596,47.03428475493507,47.03432686929892,47.03342824785487,47.03347034649629,47.0325717241254,47.03261380704487,47.03171518374766,47.03175725094564,47.03085862672256,47.03090067819952,47.0291034276716,47.02914546273426,47.02824683643014,47.02833086008972,47.02743223216824,47.02751619433757,47.02661756479999,47.02665952265326,47.02576089219306,47.02580283432836,47.02490420294606,47.02498804075635,47.02408940776073,47.0242150497123,47.02062050714369,47.02066235498111,47.01976371807652,47.01980555019932,47.01800827432475,47.01821720268096,47.01731856095376,47.01740202575293,47.01650338242076,47.01658678574923,47.01389085026064,47.01393252732631,47.01303388102218,47.01311718871059,47.01221854080523,47.01230178702983,47.01140313752449,47.011486322287,47.00968901985338,47.00973058832938,47.00883193608835,47.0090395497162,47.00993820538365,47.01006259535446,47.00916393763411,47.0092467878043,47.01104410575446,47.01125097512673,47.01214963717789,47.01231486499911,47.01321352955163,47.01329605468747,47.01060005626757,47.01064129425902,47.00974262698794,47.00982505654306,47.00892638768673,47.00900875579347,47.00811008535313,47.00819239201323,47.00729371999014,47.00741706528036,47.00561971649071,47.00566080018791,47.01015417387119,47.01027734499871,47.00758131536283,47.00766335027298,47.00856202839653,47.00860302399529,47.00770434519515,47.00774532509535,47.00864400457184,47.00872592065812,47.00782723982969,47.00790909447571,47.00611172944535,47.00623439469917,47.00533570982418,47.00541741020295,47.00631609642637,47.00635692475723,47.00725561143106,47.00729642541369,47.00639773806625,47.00643853635338,47.00553984810922,47.00562139827175,47.00292532816345,47.00300681420434,47.00480419718913,47.00488562583355,47.00578431833566,47.00598763042311,47.00329154218204,47.00321026643726,47.00141287584711,47.00149414890914,46.99700066184716,46.99704127249561,46.99614257374967,46.99618316870851,46.9943857692152,46.99446691139632,46.99356820998019,46.99360875787313,46.99271005556798,46.99275058777318,46.99185188457943,46.99189240109742,46.99009499271463,46.99013549287745,46.9892367876891,46.98931774162506,46.98752032792041,46.98756078102661,46.98486465683057,46.98490509291744,46.98400638374935,46.98404680415259,46.98314809410009,46.98326926321213,46.98147183845667,46.98151219679627,46.97881605604093,46.97889672367898,46.97799800833042,46.97803831896164,46.97624088628717,46.97636172410336,46.97546300544894,46.97554348759281,46.97644220757549,46.97652263098974,46.97562390967979,46.97570427170958,46.97480554885845,46.97484570668906,46.97394698296051,46.9741876106747,46.97328888276049,46.97336897056437,46.97426769980031,46.9743077218419,46.97340899194543,46.97368872122049,46.97278998649337,46.97282988710339,46.97193115150399,46.97197103644177,46.97107229997056,46.9711121692366,46.96841595656993,46.96845580884924,46.96935454727546,46.96939438519942,46.97029312407027,46.97065099569815,46.97154974026222,46.97162910345109,46.97252784911181,46.97256750883984,46.97077001599616,46.97080965940135,46.96991091200577,46.96995053974317,46.96905179148087,46.96909140355095,46.96639515552901,46.96643475062525,46.96553600020727,46.9656151440556,46.96471639211994,46.9647559408727,46.96385718807304,46.96389672116111,46.96299796749794,46.96307698733379,46.96038072116797,46.96042020661308,46.95952145015226,46.95956091993484,46.95416837286155,46.95412890698716,46.9532301490617,46.95326961428478,46.95147209650643,46.95151154541821,46.94971402550624,46.94975345810771,46.9488546971897,46.94889411413235,46.94799535235666,46.948034753641,46.94713599100807,46.94721474725343,46.94631598311403,46.94635533807573,46.94545657308038,46.94549591238539,46.94639467802998,46.94647331291597,46.94557454597366,46.94569238380119,46.94479361470797,46.94495051911808,46.94584929080079,46.94600595767987,46.94510718341165,46.94522452405408,46.94252819420569,46.94248909760142,46.94159032121988,46.94155121025434,46.94065243431331,46.94061330898704,46.93971453348696,46.93963623910798,46.93873746469549,46.93869829564363,46.93690074749868,46.93686156473451,46.93506401706723,46.93502482059174,46.93412604709965,46.93408683626673,46.93318806321849,46.93307034263708,46.9321715713285,46.93213230177491,46.93303107243543,46.9329131718073,46.933811940319,46.93369390270888,46.93279513614513,46.93271637068523,46.93361513594914,46.93357573006345,46.93267696544985,46.93263754521087,46.93173878104461,46.93169934645286,46.92990181881311,46.92986237052012,46.92806484337211,46.9279464123857,46.92704765046278,46.92696862276517,46.92517110092201,46.92513156587529,46.92243528355695,46.92239573546684,46.92329449578815,46.92321535330039,46.9223165942851,46.92219777000898,46.92040025529757,46.92036061851293,46.91856310430769,46.91852345383188,46.91762469708336,46.91746594780682,46.91836470193563,46.91828523599222,46.91738648317474,46.91734672836034,46.91644797599897,46.9163684226865,46.91546967143816,46.91542987294085,46.91453112214973,46.91449130931116,46.91359255897771,46.91355273179839,46.91175523184884,46.91171539098671,46.91081664137128,46.91065713058364,46.90885963602182,46.90881972215015,46.9079209752307,46.9078810470216,46.90698230056315,46.90694235801724,46.90514486582581,46.90510490960287,46.90420616387056,46.90412620775855,46.90142997333851,46.90146995691256,46.89787497229889,46.8979149382408,46.89611744344139,46.89615739307168,46.89435989617385,46.89431994786198,46.89342119978033,46.89346114743304,46.8907649000455,46.89072495437023,46.88712995777977,46.88705002673144,46.88794877484989,46.88782876387094,46.88872750981606,46.88860736194172,46.88950610571102,46.88946602577725,46.89036476869165,46.89032467310562,46.89122341516464,46.89118330392585,46.89208204512899,46.89200177635362,46.89290051603802,46.89269957843702,46.89359831461131,46.89331636768094,46.89511383014022,46.89483113930303,46.89572986557522,46.89568942053572,46.89748687115829,46.89740593343272,46.89830465711449,46.89826416509511,46.89916288791242,46.89908185755921,46.89998057884289,46.8998589188606,46.90165635682236,46.90161577216993,46.90251449018533,46.90247388986968,46.90337260701763,46.90333199103828,46.90423070731829,46.90414942903738,46.90325071409852,46.90321005313815,46.90590619534769,46.90586551738014,46.90676423038213,46.90672353674878,46.90762224888059,46.90754081528664,46.9084395258758,46.90839878591452,46.90929749563225,46.90925674000355,46.9083580309583,46.90823567612048,46.90913438314676,46.90909356753772,46.90999227369107,46.90995144241328,46.91085014769316,46.91072756186593,46.91162626492311,46.91158537231519,46.91248407449743,46.9124022429446,46.91330094357613,46.91326000463057,46.91415870438588,46.91411774976857,46.91501644864717,46.91497547835768,46.91677287415914,46.91673188752066,46.91763058444278,46.91758958213086,46.92028566965558,46.92020361597883,46.92110231006171,46.92106126004937,46.9219599532521,46.92191888756419,46.9228175798863,46.92277649852235,46.92367518996335,46.92363409292291,46.92453278348229,46.92449167076484,46.92539036044214,46.92534923204725,46.92624792084194,46.9262067767691,46.92890283989048,46.92886167878088,46.93065905181283,46.93049425193678,46.9295955682427,46.92947181287457,46.92767444895533,46.92763316852839,46.92673448694272,46.92677576668853,46.92497840154159,46.92501966492638,46.92322229760148,46.92301583750434,46.9212184761784,46.92097023731797,46.91917288337115,46.9191314591038,46.91823278250951,46.91819134392905,46.91729266781588,46.91725121492285,46.91635253929121,46.91631107208612,46.91361504603139,46.91373939650288,46.91194204068796,46.91198345948755,46.91018610150144,46.91014468406866,46.90924600545794,46.90916312697675,46.90826444953309,46.9081815124441,46.90728283616868,46.9071998404741,46.90360513885182,46.90356362125461,46.90176627061761,46.90172473939851,46.90082606446704,46.90078451894155,46.89988584449708,46.89984428466569,46.89894561070859,46.89902871400879,46.89813003848243,46.89817156696059,46.89727289055119,46.89731440335348,46.89641572606151,46.89637421394415,46.89367818293949,46.89363665788716,46.89004194986813,46.89000041256757,46.88640570415604,46.88644723871536,46.88375120330274,46.88370967079923,46.88191231425608,46.88187076813661,46.88097209025958,46.88093052983942,46.8800318524544,46.87999027773409,46.87909160084157,46.87905001182159,46.87815133542198,46.87802648050887,46.87712780597615,46.87696112571608,46.87606245374048,46.87597902507789,46.87508035428606,46.87503861816659,46.87324127738282,46.87319952765678,46.8696048465415,46.86956308458788,46.86686707381021,46.86682529894159,46.86412928851677,46.86417106131782,46.86237371870652,46.86241547514855,46.86151680286977,46.86155854364241,46.85976119714027,46.8598446309909,46.85445257870368,46.85449426903183,46.85269691539254,46.85265522643994,46.85085787342852,46.85081617087383,46.84991749477626,46.84987577793223,46.84897710233662,46.84893537120372,46.84713802083146,46.84697095204242,46.84517360643996,46.84508998494599,46.84598865646041,46.84586310982957,46.84766044815927,46.84757667276129,46.84398200088791,46.84389817111634,46.84299950406761,46.84295756741062,46.83756556539376,46.83764941543252,46.83675074640257,46.83683453516439,46.83503719379021,46.83499530827802,46.83319796755428,46.83323985168454,46.83234118035795,46.83200569362332,46.83110702764893,46.83106502512576,46.83016635966213,46.83012434286042,46.83102300763087,46.83093892772616,46.83363491678268,46.83359285229119,46.83539017602732,46.83534809517518,46.8380440773217,46.83800197941336,46.8389006390668,46.83881639694011,46.83971505501962,46.83950418418621,46.84040283860252,46.8403606188181,46.84125927235323,46.84109024068815,46.84019158994148,46.84010698566035,46.84100563501132,46.84087861419145,46.84267590814714,46.8427182647918,46.84541420642979,46.84545655019563,46.84635519773484,46.84639752722383,46.84729617527496,46.84733849048659,46.85183173142979,46.85178941272762,46.85268805965597,46.85264572527853,46.85713895360519,46.8571812914748,46.85807993727239,46.8580375987043,46.8616321772061,46.86158982086443,46.86248846431631,46.862403705296,46.86330234715918,46.86321752682063,46.86411616709387,46.86403128543538,46.86492992411737,46.86484498113722,46.86574361822673,46.8657011235652,46.86659975976229,46.86651472409496,46.86741335869757,46.86737081769109,46.87006671824142,46.87002416014784,46.87092279257775,46.87096535137344,46.87366124961215,46.87361868871002,46.87541595183144,46.87537337454192,46.8762720051095,46.87618680417587,46.87708543314369,46.87712804180508,46.8816211872394,46.8815785750631,46.88247720286149,46.88243457499728,46.88333320189668,46.88329055834407,46.88418918434402,46.88414652510254,46.88594377510613,46.88590109947122,46.88769834728023,46.88774102432329,46.88863964863639,46.88868231139701,46.8895809362167,46.88962358469438,46.89052221002022,46.89056484421442,46.89146347004595,46.89154869487901,46.89244732191872,46.89248991255678,46.89608442154319,46.89604182809432,46.89694045414132,46.8971107408138,46.89800936947107,46.8980519043684,46.89895053352808,46.89899305413753,46.90079031326099,46.90087531231575,46.90446983377031,46.9045123136158,46.9081068346598,46.90797934173608,46.90887796938971,46.90883544106255,46.90973406781245,46.90969152379107,46.91059014963678,46.91054758992067,46.9132434641353,46.91320088731909,46.91499813437503,46.91495554115989,46.91855003001535,46.91859262604162,46.91949124844752,46.91953383018146,46.9204324530852,46.92047502052627,46.9213736439274,46.92133107578399,46.92222969827763,46.92218711443592,46.92308573602153,46.92312832056587,46.92402694264848,46.92398435740154,46.92847746120928,46.92834960494271,46.93014684079937,46.93010419064133,46.93190142426077,46.93185875769681,46.9363518345831,46.93622373433295,46.93981818362597,46.93968993991592,46.94058854959879,46.94054577098866,46.94144437975547,46.94140158543858,46.94230019328882,46.94225738326467,46.94315599019786,46.94319880092838,46.94499601557455,46.94508159485554,46.94598020327378,46.94593742042832,46.95043045581214,46.95038765443352,46.9512862601669,46.95124344307845,46.95304065249324,46.95312627449213,46.95402488029239,46.95406766949172,46.96035790904406,46.96031511490214,46.96570673618088,46.96592057799928,46.96681918431806,46.96686190836358,46.96506469409967,46.96510740172771,46.9642087935667,46.96437947118635,46.9634808599931,46.96343821380302,46.96253960309784,46.96262487906362,46.96352349117598,46.96356610735157,46.96446471995151,46.9645499086865,46.96544852247631,46.96549109503488,46.9645924805426,46.96463503739102,46.96553365258557,46.96557619512832,46.96647481080885,46.9665173390454,46.96741595521149,46.96745846914132,46.96835708579254,46.96839958541513,46.97019681947003,46.97028177649248,46.97118039459637,46.97122285129414,46.97212146988139,46.97224875201696,46.97404599273952,46.97408839149732,46.97678425304551,46.97669943632062,46.97849667346113,46.97845424118147,46.97935285872313,46.97922546971215,46.98102269993316,46.98098020550339,46.98187881958317,46.98179378428363,46.9835910089766,46.98354846740413,46.98534568981131,46.98538823278782,46.98898267775801,46.98894013197324,46.9907373517244,46.99056701282002,46.99146561955143,46.99150822750283,46.99330544170417,46.99339061537294,46.9942892235449,46.99424664351565,46.99784107115921,46.99788365399923,47.00057947448862,47.00062204441929,47.00241925836084,47.00237668702514,47.00686971442288,47.00691228927141,47.00781089477573,47.0078960008192,47.00969321395729,47.00973574585334,47.01063435278409,47.01067687036161,47.01427129861837,47.01422877823409,47.01512738402723,47.01516990511322,47.01696711741818,47.01692459492872,47.01782320003661,47.01778066182297,47.01867926600001,47.01863671206172,47.01953531530737,47.01949274564399,47.02039134795777,47.02026354672194,47.021162146697,47.02111951553559,47.02201811457729,47.02197546768894,47.02467126132126,47.02462859729733,47.02552719400948,47.02548451425713,47.02638311003395,47.02634041455273,47.02813760400323,47.028009423182,47.03250238217262,47.03254512766501,47.03434231104397,47.03438504292068,47.03618222677707,47.03622494503708,47.03712353731988,47.03703808436242,47.03793667500103,47.0378939252757,47.03969110443916,47.03977659168387,47.04157377273067,47.04165920268253,47.04435497671938,47.04444035078323,47.04713612692628,47.04726408163553,47.04816267532268,47.04837563607795,47.04927423304339,47.04935931353423,47.05025791166723,47.0504278950805,47.05132649578216,47.05141139870108,47.052310000567]},{"lng":[-1.392071745705658,-1.39200950824722,-1.389388051479397,-1.389450245319519,-1.392071745705658],"lat":[46.88774102432329,46.88684239981314,46.88692770753291,46.8878263334505,46.88774102432329]}]],[[{"lng":[-1.07652066134414,-1.072583484572107,-1.072526514482436,-1.071214140290684,-1.071157193974997,-1.069844839752745,-1.069787917209588,-1.067163246804185,-1.067106369902264,-1.064481735639259,-1.064424904375948,-1.061800306260468,-1.061743520633141,-1.06043124059954,-1.060374478737024,-1.057749956720855,-1.05763652615235,-1.056324306038076,-1.056267615461262,-1.054955415320477,-1.054842083577079,-1.053529925275261,-1.053473284106399,-1.052161145777815,-1.052104528364213,-1.050792410009053,-1.050735816349257,-1.049423717967708,-1.049367148060263,-1.045430907173345,-1.045374404741144,-1.044062342533487,-1.044005863849462,-1.040069731499573,-1.040013320284465,-1.037389266354646,-1.037332900744147,-1.03339687131835,-1.033284276901305,-1.031972307052,-1.032028581461162,-1.03071658787339,-1.03066033532385,-1.025412429637348,-1.025356266402446,-1.016172512301024,-1.016116503953672,-1.013492601087951,-1.01338067743595,-1.012068766916788,-1.012012829750571,-1.009389046825162,-1.009333155239843,-1.00277377456895,-1.002717994131115,-1.001406134263755,-1.001350377543982,-1.002662215555477,-1.002271817659661,-1.00096013261887,-1.00079291270654,-1.002104532199359,-1.002048774101724,-1.003360369885593,-1.003025729307752,-1.001714264595295,-1.001602776311633,-1.000291353423502,-1.00034707479306,-0.9924783673836531,-0.9924227770684068,-0.9911113411782831,-0.9909446468874408,-0.9883219005927354,-0.9882663832176232,-0.9856436731920101,-0.9855882013484051,-0.9842768653979083,-0.9842214172442748,-0.9829101012822589,-0.9828546768171449,-0.9815433808437938,-0.9814879800657473,-0.9801767040812479,-0.9801213269888157,-0.9774988131507354,-0.9774434815788404,-0.9761322437262633,-0.9759663255932711,-0.974655151400407,-0.9745998708805214,-0.9732887166767979,-0.9730677003692165,-0.9717566316498062,-0.9715909542638131,-0.9702799491904426,-0.9702247489081974,-0.9689137638213532,-0.968803412435832,-0.9648705770968419,-0.9648154696438414,-0.9635045426792598,-0.9633392968252998,-0.9646501583150048,-0.9645950582151414,-0.9659058960432838,-0.9657405413057567,-0.967051311829499,-0.9669961754403643,-0.9683069223045484,-0.9681414587149575,-0.9694521382829125,-0.9693969656150067,-0.970707621525403,-0.970652428881195,-0.9719630611342196,-0.9719078485151593,-0.9732184571109942,-0.9730527648598716,-0.9782949190466744,-0.9780183659115996,-0.9767079391848102,-0.9766526559000019,-0.9779630608180709,-0.9779077575688222,-0.9792181388347077,-0.9788862272317296,-0.980196475821048,-0.9801411418736062,-0.9814513668151923,-0.9813960129108179,-0.9840164136525499,-0.9839056241777833,-0.9852157781772823,-0.9851603644108013,-0.9864704947633076,-0.9864150610456875,-0.987725167751385,-0.987669714084075,-0.9889797971431491,-0.9888688517614571,-0.98755881229571,-0.9875033641744703,-0.9848833232893394,-0.9846063283443176,-0.985916240742822,-0.985805412109519,-0.9962842960830054,-0.9962287102206578,-0.9975385405976451,-0.9974829347997132,-1.000102546415743,-1.000158195790229,-1.002777843562805,-1.002666459517825,-1.003976237044619,-1.003920526020066,-1.00523027990346,-1.005118819854294,-1.006428528309601,-1.006372779287064,-1.007682464100158,-1.007626695152352,-1.00893635632342,-1.008824780439484,-1.012753622101475,-1.012697771605747,-1.01400735999054,-1.013951489576144,-1.015261054319026,-1.015205163987408,-1.019133781700272,-1.019245694902179,-1.020555273969588,-1.020611255148115,-1.023230451245457,-1.023174426506743,-1.031031839284844,-1.03097568573924,-1.032285226211129,-1.032229052758916,-1.036157597606043,-1.036213836391536,-1.040142429703874,-1.040198735698205,-1.041508284828497,-1.041395631160929,-1.042705134858888,-1.042592445147208,-1.050449166567458,-1.050279754149831,-1.052898515600804,-1.052955028070695,-1.054264427744686,-1.054207893502057,-1.055517269519522,-1.055573825534796,-1.062120786201833,-1.062177452967133,-1.06348686121158,-1.063430172673479,-1.064739557256184,-1.064796267566993,-1.06610567203265,-1.066219141871031,-1.067528587993412,-1.067585347523246,-1.070204277642202,-1.070147474563277,-1.074075790219988,-1.074132658621839,-1.076751569800635,-1.07669465785064,-1.079313517895486,-1.079256564295624,-1.083184774808095,-1.083127757789434,-1.087055885886431,-1.087112968220838,-1.088422362221773,-1.088479468230767,-1.089788882101195,-1.089903143375021,-1.091212598888455,-1.09126975415555,-1.099126577859,-1.099241155409472,-1.100550662899036,-1.100607976313974,-1.101917503669435,-1.102032179783337,-1.100722608874042,-1.100779928019364,-1.099470333421752,-1.099584933885085,-1.098275293819324,-1.098561719491559,-1.099871468466356,-1.09992878111302,-1.101238549961434,-1.101295886303061,-1.102605675024903,-1.102663035062941,-1.103972843658018,-1.104030227393913,-1.105340055862036,-1.105397463297234,-1.106707311638211,-1.106764742774159,-1.108074610987798,-1.108132065825943,-1.112061724335565,-1.112119246453809,-1.113429150577846,-1.113486696402435,-1.116106542472329,-1.116336919279005,-1.117646926595919,-1.117704547389957,-1.120324599844314,-1.120382266144794,-1.121692311281224,-1.121750001297553,-1.123060066303908,-1.123348654214904,-1.127279164614195,-1.127336953361789,-1.131267511811541,-1.131325367879336,-1.13394577413039,-1.134003675723834,-1.135313897752008,-1.135487679515378,-1.136797965010596,-1.136855917929187,-1.14078682821978,-1.140844848475222,-1.142155169838322,-1.142213213830508,-1.143523555060014,-1.143581622790391,-1.144891983886107,-1.144950075356125,-1.146260456317859,-1.146376688678115,-1.145066264107923,-1.145124361389898,-1.142503462827073,-1.142909892739659,-1.14159928789427,-1.141715384532133,-1.144336679533029,-1.144394774381831,-1.143084106037097,-1.143258336767823,-1.144569070555164,-1.144743384170923,-1.147364976828126,-1.147423128878621,-1.14873394411544,-1.14879211992408,-1.150102955038452,-1.150161154606687,-1.15147200959842,-1.151588458199329,-1.152899354888372,-1.152957603923318,-1.15426852048996,-1.154501623353754,-1.153190619496775,-1.153365401584944,-1.154676470917291,-1.15473475732769,-1.153423666168771,-1.153540201168667,-1.154851335982907,-1.154967922417993,-1.156279098945296,-1.156512382482953,-1.157823644388742,-1.15799869297845,-1.155376036218505,-1.155434345972602,-1.154122992838224,-1.154239574511599,-1.152928175759078,-1.152986447676277,-1.151675025141423,-1.151733277165025,-1.150421830847643,-1.150538297048278,-1.145292317633298,-1.145350466285418,-1.14403894473693,-1.144155204179257,-1.142843637006035,-1.142901747794159,-1.141590156837871,-1.141648247721413,-1.140336632981864,-1.14039470395937,-1.137771424975702,-1.137829454199739,-1.136517789956053,-1.136575799270326,-1.133952421280848,-1.133894455660471,-1.128647763856807,-1.12870564208823,-1.127393942458133,-1.127509661020089,-1.126197915761335,-1.126313598349125,-1.122378219889145,-1.122435998525261,-1.119812359536754,-1.119754624603925,-1.118442824069924,-1.118615969093306,-1.117304101074847,-1.117419497015993,-1.114795667782944,-1.114853324931714,-1.110917501079298,-1.111032690010041,-1.109720701164874,-1.109835854070479,-1.108523819585532,-1.108581377062211,-1.107269318795567,-1.107326856333466,-1.106014774284933,-1.10595725860778,-1.103333132471685,-1.103390604427057,-1.099454335249128,-1.09939692887715,-1.095460708031591,-1.095403369158467,-1.094091313573996,-1.094148630586237,-1.090212386769851,-1.090269640110559,-1.088957533152087,-1.089014766541861,-1.087702635809513,-1.087817064594812,-1.086504888224031,-1.086619280920893,-1.085307058909116,-1.08536423625598,-1.082739742771513,-1.082796878292801,-1.08148460682055,-1.081541722381546,-1.080229427133451,-1.080286522732697,-1.078974203708568,-1.079088356886854,-1.076463625638866,-1.07652066134414],"lat":[46.94530267611154,46.94541979164084,46.94452101200764,46.94456001985834,46.94366123937586,46.94370023157554,46.94280145024424,46.94287938833565,46.94198060551302,46.94205848229092,46.94115969797826,46.94123751344442,46.94033872764293,46.94037761222423,46.9394788255767,46.93955654843713,46.93775897196461,46.93779780960346,46.9368990204207,46.93693784241299,46.9351402621561,46.93517906786192,46.93428027678852,46.93431906684901,46.93342027493273,46.93345904934836,46.9325602565897,46.93259901536096,46.93170022176041,46.9318164061295,46.93091761040956,46.93095630788525,46.9300575113249,46.9301735118126,46.92927471313636,46.92935197049986,46.92845317034738,46.92856894195335,46.92677133722436,46.92680989647967,46.9277086995809,46.92774724446914,46.92684844073179,46.9270024677076,46.9261036612278,46.92637262636605,46.92547381524723,46.92555052613405,46.92375290076485,46.92379123243833,46.9228924188218,46.92296903589619,46.922070220816,46.9222614978009,46.92136267936502,46.92140088912429,46.92050206985912,46.92046386073044,46.91417212474301,46.9142103294579,46.91151386632766,46.9114756635042,46.91057684269842,46.91053862550621,46.90514570034624,46.90518391375449,46.90338626921305,46.90342446636283,46.90432328936114,46.9045521610932,46.9036533343179,46.90369142648616,46.90099494311008,46.90107107868474,46.90017224924897,46.90024832358096,46.89934949269647,46.8993875067402,46.89848867503537,46.89852667345561,46.89762784093091,46.89766582372816,46.89676699038408,46.89680495755879,46.89590612339584,46.89598201150555,46.89508317589816,46.89512109683389,46.89242458698419,46.89246249104796,46.89156365342325,46.8916015418671,46.88800618696018,46.88804405790965,46.88534753785846,46.88538539194037,46.88448655091961,46.88452438938431,46.88272670552683,46.88284012722197,46.88194128313841,46.88197905976247,46.87928252451201,46.87924474975811,46.87834590492223,46.87830811580077,46.87561158203955,46.87557377979872,46.87467493546143,46.87463711885452,46.87194058659564,46.87190275687173,46.87100391303738,46.87096606894895,46.87006722555329,46.87002936710078,46.86913052414427,46.86909265132824,46.86639612322212,46.86624448958286,46.86175028489162,46.86178821265308,46.86088937053711,46.86085144340148,46.85995260172793,46.8599146602325,46.85452161010761,46.85448365738439,46.85358481568523,46.8535468486044,46.85264800735005,46.85257202949073,46.85077434894521,46.8507363387947,46.84983749887741,46.84979947437139,46.8489006349007,46.84886259603964,46.84796375701598,46.84792570380042,46.84612802646784,46.8461660784276,46.84526723886332,46.84534329658274,46.84084908979455,46.84081107156235,46.83901338884971,46.83870871376096,46.83780987716448,46.83777172599985,46.83687288985481,46.8367965438477,46.837695378733,46.83761897154956,46.83582130412253,46.83578307932361,46.83488424597415,46.83484600682763,46.83304834085843,46.83301008799602,46.83211125537689,46.83207298816812,46.83117415600358,46.83113587444898,46.82933821085352,46.82922328011945,46.82832444995364,46.82828611038829,46.82738728067908,46.82734892677022,46.82645009751802,46.82633494783521,46.82813260236377,46.82809418791705,46.82899301428353,46.82891613919315,46.82801731409508,46.82778633320344,46.82688751174038,46.82684896314523,46.82595014214225,46.8258344084126,46.82673322750601,46.82661735708686,46.82751617409244,46.82747752003067,46.82567988711901,46.82564121935757,46.82384358701859,46.82361127363935,46.82091483531483,46.82083728156455,46.8217360932349,46.82169729325911,46.82079848222891,46.82075966791967,46.82165847830943,46.82146417895414,46.82236298596305,46.82232408052889,46.82142527416193,46.82138635439604,46.8222851601209,46.82224622473903,46.82404383437756,46.82400488273653,46.82490368664951,46.82482573715853,46.82392693453164,46.82380989991771,46.82470870061361,46.82463060137839,46.82373180197106,46.82365364412814,46.82275484583443,46.82263749870126,46.82173870216792,46.82162122221016,46.82252001680521,46.82248082622602,46.82337961999865,46.82334041379909,46.82513799952232,46.82509877705478,46.82599756900466,46.82576191585428,46.82755949144741,46.8275201622145,46.82841894909646,46.8283796042397,46.83017717617266,46.83021652232775,46.83111530867662,46.8311546405052,46.8329522139663,46.83299153211657,46.83748546588437,46.83744614449037,46.83834493005706,46.83830559303659,46.83920437777434,46.83916502512689,46.84006380903524,46.84002444076037,46.84092322383879,46.84088383993599,46.84178262218402,46.84174322265284,46.84264200406996,46.84260258890997,46.84350136949568,46.84338303219082,46.84428181064271,46.84424233426527,46.8451411118843,46.84506211288868,46.84865721633074,46.84861769175551,46.84951646650769,46.84943737111154,46.85033614437575,46.85029657355413,46.85119534598225,46.85115575952684,46.85564961564594,46.85553075658942,46.8564295252994,46.85631052944165,46.85720929600342,46.85712988921068,46.85802865427719,46.85798892775186,46.86068521987226,46.86064547639716,46.86154423940988,46.86142491711658,46.8623236779743,46.86228387325296,46.86318263326741,46.86314281290546,46.86404157207618,46.86400173607313,46.86490049439965,46.86486064275501,46.86665815753122,46.86669801049094,46.86759676825531,46.86767643053488,46.87396773881539,46.8740075520758,46.87580506834036,46.8757254242045,46.87662418073807,46.87666401095686,46.8793602813881,46.87932044919763,46.88201671594089,46.88193700264939,46.88283575653274,46.88279587674555,46.88369462977924,46.88365473434462,46.88455348652816,46.88451357544562,46.88631107791908,46.8862711505294,46.8871699008186,46.88712995777977,46.89072495437023,46.8907649000455,46.89346114743304,46.89342119978033,46.89431994786198,46.89435989617385,46.89615739307168,46.89611744344139,46.8979149382408,46.89787497229889,46.90146995691256,46.90142997333851,46.90412620775855,46.90420616387056,46.90510490960287,46.90514486582581,46.90694235801724,46.90698230056315,46.9078810470216,46.9079209752307,46.90881972215015,46.90885963602182,46.91065713058364,46.91081664137128,46.91171539098671,46.91175523184884,46.91355273179839,46.91359255897771,46.91449130931116,46.91453112214973,46.91542987294085,46.91546967143816,46.9163684226865,46.91644797599897,46.91734672836034,46.91738648317474,46.91828523599222,46.91836470193563,46.91746594780682,46.91762469708336,46.91852345383188,46.91856310430769,46.92036061851293,46.92040025529757,46.92219777000898,46.9223165942851,46.92321535330039,46.92329449578815,46.92239573546684,46.92243528355695,46.92513156587529,46.92517110092201,46.92696862276517,46.92704765046278,46.9279464123857,46.92806484337211,46.92986237052012,46.92990181881311,46.93169934645286,46.93173878104461,46.93263754521087,46.93267696544985,46.93357573006345,46.93361513594914,46.93271637068523,46.93279513614513,46.93369390270888,46.933811940319,46.9329131718073,46.93303107243543,46.93213230177491,46.9321715713285,46.93307034263708,46.93318806321849,46.93408683626673,46.93412604709965,46.93502482059174,46.93506401706723,46.93686156473451,46.93690074749868,46.93869829564363,46.93873746469549,46.93963623910798,46.93971453348696,46.94061330898704,46.94065243431331,46.94155121025434,46.94159032121988,46.94248909760142,46.94252819420569,46.94432574764348,46.94440389841117,46.94530267611154]}]],[[{"lng":[-1.266194758901899,-1.260947194845562,-1.260826999683059,-1.259515147351682,-1.259334935200219,-1.25802314640918,-1.257963101549802,-1.256651332601528,-1.256591311593105,-1.255279562487795,-1.255219565328875,-1.253907836066731,-1.253847862755862,-1.252536153337082,-1.252476203872811,-1.251164514297598,-1.251104588678472,-1.249792918947027,-1.249733017171595,-1.247109715397909,-1.24693015312061,-1.245618564763201,-1.245558736507113,-1.244247167993644,-1.244127561151636,-1.242816034321867,-1.24275625573367,-1.241444748747617,-1.24138499399275,-1.240073506850609,-1.240013775927623,-1.23739083933882,-1.23733115408332,-1.236019704637684,-1.235960043210534,-1.232025748439255,-1.231966154510913,-1.230654740777788,-1.230595170673642,-1.229283776787903,-1.229224230506502,-1.225290102433434,-1.225349583209985,-1.22403818137953,-1.223978722438129,-1.222667340458558,-1.222607905335946,-1.219985179096256,-1.219925789624976,-1.218614445365878,-1.218555079709857,-1.217243755303064,-1.21706573571705,-1.215754474826265,-1.215813790876704,-1.213191219494831,-1.213131947107469,-1.21182068028091,-1.211761431702378,-1.210450184729738,-1.210390959958584,-1.209079732840062,-1.208961332885002,-1.205027770656647,-1.204968639127476,-1.202346298099176,-1.20228721219898,-1.20097606055524,-1.200916998453915,-1.198294732910015,-1.198176701925049,-1.196865609852158,-1.196806619138743,-1.18762907233883,-1.187511402363986,-1.178334064929382,-1.178275385657337,-1.171720195008375,-1.171661626815495,-1.169039580521373,-1.168981057930308,-1.167670053678101,-1.167611554862819,-1.166300570482303,-1.166242095441355,-1.164931130932726,-1.16487267966466,-1.160939839907948,-1.160881456054493,-1.156948664232951,-1.156890347790537,-1.152957603923318,-1.152899354888372,-1.151588458199329,-1.15147200959842,-1.150161154606687,-1.150102955038452,-1.14879211992408,-1.14873394411544,-1.147423128878621,-1.147364976828126,-1.144743384170923,-1.144569070555164,-1.143258336767823,-1.143084106037097,-1.144394774381831,-1.144336679533029,-1.141715384532133,-1.14159928789427,-1.142909892739659,-1.142503462827073,-1.145124361389898,-1.145066264107923,-1.146376688678115,-1.146260456317859,-1.144950075356125,-1.144891983886107,-1.143581622790391,-1.143523555060014,-1.142213213830508,-1.142155169838322,-1.140844848475222,-1.14078682821978,-1.136855917929187,-1.136797965010596,-1.135487679515378,-1.135313897752008,-1.134003675723834,-1.13394577413039,-1.131325367879336,-1.131267511811541,-1.127336953361789,-1.127279164614195,-1.123348654214904,-1.123060066303908,-1.121750001297553,-1.121692311281224,-1.120382266144794,-1.120324599844314,-1.117704547389957,-1.117646926595919,-1.116336919279005,-1.116106542472329,-1.113486696402435,-1.113429150577846,-1.112119246453809,-1.112061724335565,-1.108132065825943,-1.108074610987798,-1.106764742774159,-1.106707311638211,-1.105397463297234,-1.105340055862036,-1.104030227393913,-1.103972843658018,-1.102663035062941,-1.102605675024903,-1.101295886303061,-1.101238549961434,-1.09992878111302,-1.099871468466356,-1.098561719491559,-1.098275293819324,-1.099584933885085,-1.099470333421752,-1.100779928019364,-1.100722608874042,-1.102032179783337,-1.101917503669435,-1.100607976313974,-1.100550662899036,-1.099241155409472,-1.099012007944622,-1.100321428337702,-1.100264124471497,-1.101573521182804,-1.101516197454562,-1.105444310810586,-1.105329538462941,-1.106638862218009,-1.106581457143384,-1.109200055376834,-1.109027721024573,-1.11033695196881,-1.110279489246779,-1.118134704088319,-1.118077112692485,-1.124622963478876,-1.124219115781243,-1.125528127847367,-1.125470421248393,-1.126779409634117,-1.126606236100736,-1.127915157294002,-1.127799674376109,-1.12910855013677,-1.128877520369446,-1.130186307196846,-1.130128532813685,-1.131437295966359,-1.131495092099945,-1.138038987725667,-1.138096894538563,-1.142023274027624,-1.142081248022829,-1.148625283964528,-1.14868336864837,-1.152609832257449,-1.15266798413158,-1.157903328830431,-1.157961569650508,-1.161888123064585,-1.161946431082492,-1.163255300089713,-1.163313631802978,-1.168549175407115,-1.168607596077947,-1.169916498869423,-1.169974943240372,-1.172592786490376,-1.17265127631753,-1.175269155283478,-1.175152090560886,-1.176460983611534,-1.176343883182283,-1.1776527307777,-1.17759416173696,-1.17890298563025,-1.17884439678972,-1.180153196981076,-1.180035981655335,-1.181344736396043,-1.181286109912897,-1.18259484095288,-1.182477550350075,-1.183786235942976,-1.183727571825383,-1.186344893657751,-1.186227484305109,-1.184918867855236,-1.184801509808548,-1.183492934890883,-1.183434280540739,-1.182125725412006,-1.181773970763086,-1.183082395453607,-1.183023754773797,-1.184332155774626,-1.184273495310971,-1.188198621382664,-1.188139897665812,-1.190756595134574,-1.190639066629903,-1.193255669336392,-1.193196864553092,-1.195813415962725,-1.195460367611627,-1.196768509998076,-1.196650799430399,-1.197958896403478,-1.197900022332526,-1.196591947085367,-1.196356557295005,-1.197664545645685,-1.197605681372728,-1.200221608750814,-1.200162702994255,-1.201470642021235,-1.201352792949511,-1.202660686574867,-1.202601743261764,-1.203909613206278,-1.203791689029614,-1.205099513576067,-1.204981553816009,-1.20759711015434,-1.207361040475325,-1.208668728842907,-1.208550662363684,-1.20985830534341,-1.209681155227955,-1.212296305033221,-1.212237215508065,-1.217467404685743,-1.217408230297863,-1.220023269661382,-1.220082487463126,-1.224005096848693,-1.223945813926893,-1.229175845328466,-1.22923521507563,-1.230542739695341,-1.230483348242081,-1.233098348140143,-1.233157783005297,-1.238387845887216,-1.238447369555539,-1.239754902032531,-1.239873998725668,-1.242489144563006,-1.242548739298161,-1.249086677715043,-1.249206090223459,-1.25051371536671,-1.250573446310884,-1.251881091176054,-1.251940845818226,-1.259786803486157,-1.259966470148826,-1.261274187922353,-1.261334102508554,-1.262641840000549,-1.262701778293802,-1.264009535504067,-1.264129461503177,-1.262821660864812,-1.262881605142768,-1.261573780794174,-1.26163370535128,-1.262941551415822,-1.263001499684077,-1.261693631902867,-1.261873423525499,-1.260565488595088,-1.260924995376059,-1.26223306063515,-1.262472858616423,-1.263781008779826,-1.263900963199519,-1.265209154820347,-1.265329160684107,-1.264020925606382,-1.264080909805252,-1.265389166612506,-1.265449174538716,-1.266757451078308,-1.266877516387023,-1.268185834390244,-1.268305951159428,-1.27092266809618,-1.271042879799821,-1.269734478864055,-1.269914746223368,-1.271223212365975,-1.271343444085037,-1.270034934467599,-1.270155130716515,-1.274080783994978,-1.274140950340957,-1.278066650806241,-1.278126884376632,-1.282052632009636,-1.28199233321937,-1.284609444789535,-1.284669787059219,-1.285978361570121,-1.286099095619366,-1.287407711601125,-1.287528497176018,-1.292763114918264,-1.292702634156838,-1.294011261816797,-1.294132268838277,-1.295440937967586,-1.295501466244028,-1.296810155100648,-1.29687070713676,-1.300796826834993,-1.300857446120756,-1.302166170394205,-1.302287458513098,-1.303596224257069,-1.303656893092464,-1.304965678561915,-1.305026371164811,-1.307643979532033,-1.307765457793415,-1.306456611125565,-1.306578053981367,-1.307886944147463,-1.308069188854052,-1.306760233435176,-1.306942431095536,-1.309560470427858,-1.309682032670237,-1.310991092809101,-1.311051898723362,-1.312360978592278,-1.312421808288656,-1.313730907887424,-1.31379176136736,-1.316409997996109,-1.316470897018582,-1.317780034047535,-1.317840956857074,-1.31915011361445,-1.319211060212497,-1.320520236698092,-1.320764130438785,-1.322073391934606,-1.322317405235312,-1.324936096233912,-1.32518031618762,-1.323870884641292,-1.323931922948369,-1.322622467600244,-1.322683486172785,-1.323992963289488,-1.324115050074261,-1.322805529418236,-1.322927580798187,-1.325546707160111,-1.325607779445055,-1.32822694120976,-1.328593678959419,-1.327283968448636,-1.327345076760382,-1.324725606075232,-1.324969885460118,-1.322350232400679,-1.322228175756411,-1.318298811592114,-1.318237851658205,-1.315618308946205,-1.315740143727529,-1.310500859705109,-1.310622528357974,-1.309312658718036,-1.309373474303601,-1.308063580853664,-1.308124376683003,-1.306814459422866,-1.306936013592259,-1.305626050736546,-1.305747569438054,-1.307057575866281,-1.307422311324221,-1.31004257956885,-1.310164261437052,-1.311474436100587,-1.311535301870386,-1.312845496297303,-1.312967277512944,-1.314277513497109,-1.314460265935129,-1.313149964564418,-1.313332669891943,-1.312022301097199,-1.312265849433435,-1.310955391406751,-1.31107713415323,-1.30976663049122,-1.309827483105836,-1.307206426086435,-1.307267237121252,-1.302025011531592,-1.302085737371838,-1.298153981493583,-1.298214643938912,-1.296904032791827,-1.296964675452769,-1.295654040475314,-1.29571466335042,-1.294404004542394,-1.294464607630217,-1.293153924991417,-1.293275093610409,-1.291964365329641,-1.292146067977897,-1.293456861695986,-1.293699247440897,-1.292388366462905,-1.292630697282221,-1.291319727011435,-1.291501430887667,-1.290190393134469,-1.290250943312768,-1.291562002888029,-1.291622576909371,-1.292933656286341,-1.292994254152172,-1.295616450486499,-1.295677094022205,-1.294365975042423,-1.294426598776926,-1.293115455949651,-1.293176059881501,-1.291864893206534,-1.291986063483611,-1.290674851133613,-1.290917128632999,-1.288294523238024,-1.288355054004199,-1.285732396867335,-1.285792885989838,-1.283170177114629,-1.283230624590256,-1.281919245293546,-1.282040102628829,-1.283351525593559,-1.283472434666931,-1.287406822464676,-1.287588397741539,-1.286276871650303,-1.28633737894443,-1.285025828995193,-1.285086316469645,-1.283774742662211,-1.283835210315535,-1.282523612649699,-1.282644510328888,-1.27477461989027,-1.274895263523799,-1.268336747678974,-1.268457180916309,-1.26976893179877,-1.269889416770441,-1.27120120933522,-1.271261476685443,-1.269949662272944,-1.270070159311788,-1.266134572915658,-1.266194758901899],"lat":[46.93049425193678,46.93065905181283,46.92886167878088,46.92890283989048,46.9262067767691,46.92624792084194,46.92534923204725,46.92539036044214,46.92449167076484,46.92453278348229,46.92363409292291,46.92367518996335,46.92277649852235,46.9228175798863,46.92191888756419,46.9219599532521,46.92106126004937,46.92110231006171,46.92020361597883,46.92028566965558,46.91758958213086,46.91763058444278,46.91673188752066,46.91677287415914,46.91497547835768,46.91501644864717,46.91411774976857,46.91415870438588,46.91326000463057,46.91330094357613,46.9124022429446,46.91248407449743,46.91158537231519,46.91162626492311,46.91072756186593,46.91085014769316,46.90995144241328,46.90999227369107,46.90909356753772,46.90913438314676,46.90823567612048,46.9083580309583,46.90925674000355,46.90929749563225,46.90839878591452,46.9084395258758,46.90754081528664,46.90762224888059,46.90672353674878,46.90676423038213,46.90586551738014,46.90590619534769,46.90321005313815,46.90325071409852,46.90414942903738,46.90423070731829,46.90333199103828,46.90337260701763,46.90247388986968,46.90251449018533,46.90161577216993,46.90165635682236,46.8998589188606,46.89998057884289,46.89908185755921,46.89916288791242,46.89826416509511,46.89830465711449,46.89740593343272,46.89748687115829,46.89568942053572,46.89572986557522,46.89483113930303,46.89511383014022,46.89331636768094,46.89359831461131,46.89269957843702,46.89290051603802,46.89200177635362,46.89208204512899,46.89118330392585,46.89122341516464,46.89032467310562,46.89036476869165,46.88946602577725,46.88950610571102,46.88860736194172,46.88872750981606,46.88782876387094,46.88794877484989,46.88705002673144,46.8871699008186,46.8862711505294,46.88631107791908,46.88451357544562,46.88455348652816,46.88365473434462,46.88369462977924,46.88279587674555,46.88283575653274,46.88193700264939,46.88201671594089,46.87932044919763,46.8793602813881,46.87666401095686,46.87662418073807,46.8757254242045,46.87580506834036,46.8740075520758,46.87396773881539,46.86767643053488,46.86759676825531,46.86669801049094,46.86665815753122,46.86486064275501,46.86490049439965,46.86400173607313,46.86404157207618,46.86314281290546,46.86318263326741,46.86228387325296,46.8623236779743,46.86142491711658,46.86154423940988,46.86064547639716,46.86068521987226,46.85798892775186,46.85802865427719,46.85712988921068,46.85720929600342,46.85631052944165,46.8564295252994,46.85553075658942,46.85564961564594,46.85115575952684,46.85119534598225,46.85029657355413,46.85033614437575,46.84943737111154,46.84951646650769,46.84861769175551,46.84865721633074,46.84506211288868,46.8451411118843,46.84424233426527,46.84428181064271,46.84338303219082,46.84350136949568,46.84260258890997,46.84264200406996,46.84174322265284,46.84178262218402,46.84088383993599,46.84092322383879,46.84002444076037,46.84006380903524,46.83916502512689,46.83920437777434,46.83830559303659,46.83834493005706,46.83744614449037,46.83748546588437,46.83299153211657,46.8329522139663,46.8311546405052,46.83111530867662,46.83021652232775,46.83017717617266,46.8283796042397,46.82841894909646,46.8275201622145,46.82755949144741,46.82396433955395,46.82392501291653,46.82302622515061,46.82298688418853,46.82208809689553,46.82196998611744,46.8201724149012,46.82013301599567,46.81923423077421,46.81915538934622,46.81645903653304,46.81641959531287,46.81552081134288,46.81528385353088,46.81438507329585,46.8141872000128,46.80789575635002,46.80785614135662,46.80695736365104,46.80691773434257,46.80422140215244,46.80418175983745,46.80238420549192,46.80234454951758,46.79874944138689,46.79870977306309,46.79781099625768,46.79777131362226,46.79867008977299,46.79847144882928,46.79937022153221,46.79925085540332,46.80014962596606,46.7999503798064,46.80084914691082,46.80072941764364,46.80162818260128,46.80146833141075,46.8023670935595,46.80224704603683,46.80314580603304,46.80310575959799,46.80400451876125,46.80384418070668,46.80474293705225,46.80470281445894,46.8056015699698,46.80552127855681,46.80642003257002,46.8063396799616,46.80454217441332,46.80450197698395,46.80270447207057,46.8026642610009,46.80176550894885,46.80172528357618,46.80082653201553,46.8007862923404,46.798988790031,46.79894853671799,46.79804978596992,46.79800951835571,46.79621201767425,46.79617173642391,46.79527298649119,46.79519238042624,46.79339488270836,46.79343519189307,46.79163769216408,46.79167798505541,46.79077923427136,46.79081951153481,46.78542699928998,46.785386726013,46.78448797404635,46.78444768647225,46.78354893500214,46.78342798450651,46.78252923486378,46.7824485277258,46.78065103060033,46.78057026628381,46.77967151880242,46.77959069597821,46.7741982155889,46.77415778573871,46.77236029228116,46.77231984880689,46.7714211024972,46.77146154530433,46.76786655574826,46.76782611560962,46.76692736847664,46.76684644466431,46.76594769870223,46.76590721502934,46.76410972394959,46.76406922665673,46.76317048153972,46.76312996995956,46.76133248057279,46.76129195537437,46.75949446667367,46.75941337408801,46.75581840009127,46.75577783404384,46.75398034741426,46.75393976775245,46.75124353860917,46.75116233844597,46.7502635964166,46.75010101935023,46.74920227984238,46.74912090294266,46.75001964110815,46.74989746160978,46.74899872545957,46.74883561283442,46.74973434629405,46.74969353008715,46.74879479730077,46.74871312138094,46.74961185282006,46.7494483188735,46.7503470474548,46.75030612591578,46.75210358124629,46.75202169061344,46.75292041668639,46.75271542507139,46.75451286997039,46.75447182543833,46.75537054696817,46.75532948680676,46.75622820749725,46.75598152846279,46.75867767735242,46.75863650980951,46.75953522510122,46.75949404192596,46.76039275637508,46.76035155756701,46.76214898461532,46.76219018478255,46.76308889874051,46.76313008463313,46.76402879910626,46.76398761253429,46.76488632616365,46.76492751341503,46.76762365535281,46.76766482968688,46.77305711316962,46.77301593476025,46.77661078436964,46.77656958828518,46.77836701072675,46.77832579832474,46.78012321873491,46.7801644324966,46.78106314312905,46.78102192868751,46.7819206384715,46.78187940839076,46.78367682609175,46.78363557969088,46.78543299535393,46.78535045494974,46.78714786721063,46.7871891462547,46.7898852654117,46.78984398432478,46.79164139488019,46.79168267732905,46.79348008856119,46.79335619224126,46.79425489555601,46.79413086252868,46.79502956362535,46.79490539388709,46.79400669483873,46.79392384156662,46.79482253924823,46.79478108948489,46.79657848296338,46.79653701686946,46.7983344082887,46.79816838880599,46.79726969592154,46.79722815432756,46.7990255385527,46.79898398062463,46.79988267179147,46.79984109821405,46.8007397885213,46.80061497594774,46.80151366402209,46.80147202921625,46.80326940346826,46.80322775232418,46.80412643850115,46.80408477170514,46.80498345701982,46.80490007715807,46.8066974445107,46.8067391432998,46.8085365113256,46.80849481116076,46.8111908598145,46.8112325620431,46.81392861116987,46.81384518761757,46.81564255006297,46.81560081445877,46.81649949472625,46.81645774346546,46.81735642286608,46.81731465594821,46.81821333448147,46.8181297543622,46.81902843133808,46.81898661813599,46.81988529424328,46.8198434653825,46.82074214062067,46.82070029610075,46.82429499249693,46.82425313024568,46.8278478209962,46.82776404605578,46.83135872838042,46.83140062610077,46.83229929691864,46.83234118035795,46.83323985168454,46.83319796755428,46.83499530827802,46.83503719379021,46.83683453516439,46.83675074640257,46.83764941543252,46.83756556539376,46.84295756741062,46.84299950406761,46.84389817111634,46.84398200088791,46.84757667276129,46.84766044815927,46.84586310982957,46.84598865646041,46.84508998494599,46.84517360643996,46.84697095204242,46.84713802083146,46.84893537120372,46.84897710233662,46.84987577793223,46.84991749477626,46.85081617087383,46.85085787342852,46.85265522643994,46.85269691539254,46.85449426903183,46.85445257870368,46.8598446309909,46.85976119714027,46.86155854364241,46.86151680286977,46.86241547514855,46.86237371870652,46.86417106131782,46.86412928851677,46.86682529894159,46.86686707381021,46.86956308458788,46.8696048465415,46.87319952765678,46.87324127738282,46.87503861816659,46.87508035428606,46.87597902507789,46.87606245374048,46.87696112571608,46.87712780597615,46.87802648050887,46.87815133542198,46.87905001182159,46.87909160084157,46.87999027773409,46.8800318524544,46.88093052983942,46.88097209025958,46.88187076813661,46.88191231425608,46.88370967079923,46.88375120330274,46.88644723871536,46.88640570415604,46.89000041256757,46.89004194986813,46.89363665788716,46.89367818293949,46.89637421394415,46.89641572606151,46.89731440335348,46.89727289055119,46.89817156696059,46.89813003848243,46.89902871400879,46.89894561070859,46.89984428466569,46.89988584449708,46.90078451894155,46.90082606446704,46.90172473939851,46.90176627061761,46.90356362125461,46.90360513885182,46.9071998404741,46.90728283616868,46.9081815124441,46.90826444953309,46.90916312697675,46.90924600545794,46.91014468406866,46.91018610150144,46.91198345948755,46.91194204068796,46.91373939650288,46.91361504603139,46.91631107208612,46.91635253929121,46.91725121492285,46.91729266781588,46.91819134392905,46.91823278250951,46.9191314591038,46.91917288337115,46.92097023731797,46.9212184761784,46.92301583750434,46.92322229760148,46.92501966492638,46.92497840154159,46.92677576668853,46.92673448694272,46.92763316852839,46.92767444895533,46.92947181287457,46.9295955682427,46.93049425193678]}]]],["X1","X1.1","X1.2","X1.3","X1.4","X3","X2"],"sfBVI",{"interactive":true,"className":"","pane":"tmap403","stroke":true,"color":"#008C8E","weight":2,"opacity":1,"fill":true,"fillColor":["#008C8E","#008C8E","#008C8E","#008C8E","#008C8E","#008C8E","#008C8E"],"fillOpacity":[0.09803921568627451,0.09803921568627451,0.09803921568627451,0.09803921568627451,0.09803921568627451,0.09803921568627451,0.09803921568627451],"dashArray":"none","smoothFactor":1,"noClip":false},["<style> div.leaflet-popup-content {width:auto !important;overflow-y:auto; overflow-x:hidden;}<\/style><div style=\"max-height:25em;padding-right:0px;\"><table>\n\t\t\t <thead><tr><th colspan=\"2\"><b>1<\/b><\/th><\/thead><\/tr><tr><td style=\"color: #888888;\"><nobr>Nmr_Mll<\/nobr><\/td><td align=\"right\"><nobr>1<\/nobr><\/td><\/tr><\/table><\/div>","<style> div.leaflet-popup-content {width:auto !important;overflow-y:auto; overflow-x:hidden;}<\/style><div style=\"max-height:25em;padding-right:0px;\"><table>\n\t\t\t <thead><tr><th colspan=\"2\"><b>1<\/b><\/th><\/thead><\/tr><tr><td style=\"color: #888888;\"><nobr>Nmr_Mll<\/nobr><\/td><td align=\"right\"><nobr>1<\/nobr><\/td><\/tr><\/table><\/div>","<style> div.leaflet-popup-content {width:auto !important;overflow-y:auto; overflow-x:hidden;}<\/style><div style=\"max-height:25em;padding-right:0px;\"><table>\n\t\t\t <thead><tr><th colspan=\"2\"><b>1<\/b><\/th><\/thead><\/tr><tr><td style=\"color: #888888;\"><nobr>Nmr_Mll<\/nobr><\/td><td align=\"right\"><nobr>1<\/nobr><\/td><\/tr><\/table><\/div>","<style> div.leaflet-popup-content {width:auto !important;overflow-y:auto; overflow-x:hidden;}<\/style><div style=\"max-height:25em;padding-right:0px;\"><table>\n\t\t\t <thead><tr><th colspan=\"2\"><b>1<\/b><\/th><\/thead><\/tr><tr><td style=\"color: #888888;\"><nobr>Nmr_Mll<\/nobr><\/td><td align=\"right\"><nobr>1<\/nobr><\/td><\/tr><\/table><\/div>","<style> div.leaflet-popup-content {width:auto !important;overflow-y:auto; overflow-x:hidden;}<\/style><div style=\"max-height:25em;padding-right:0px;\"><table>\n\t\t\t <thead><tr><th colspan=\"2\"><b>1<\/b><\/th><\/thead><\/tr><tr><td style=\"color: #888888;\"><nobr>Nmr_Mll<\/nobr><\/td><td align=\"right\"><nobr>1<\/nobr><\/td><\/tr><\/table><\/div>","<style> div.leaflet-popup-content {width:auto !important;overflow-y:auto; overflow-x:hidden;}<\/style><div style=\"max-height:25em;padding-right:0px;\"><table>\n\t\t\t <thead><tr><th colspan=\"2\"><b>3<\/b><\/th><\/thead><\/tr><tr><td style=\"color: #888888;\"><nobr>Nmr_Mll<\/nobr><\/td><td align=\"right\"><nobr>3<\/nobr><\/td><\/tr><\/table><\/div>","<style> div.leaflet-popup-content {width:auto !important;overflow-y:auto; overflow-x:hidden;}<\/style><div style=\"max-height:25em;padding-right:0px;\"><table>\n\t\t\t <thead><tr><th colspan=\"2\"><b>2<\/b><\/th><\/thead><\/tr><tr><td style=\"color: #888888;\"><nobr>Nmr_Mll<\/nobr><\/td><td align=\"right\"><nobr>2<\/nobr><\/td><\/tr><\/table><\/div>"],null,["1","1","1","1","1","3","2"],{"interactive":false,"permanent":false,"direction":"auto","opacity":1,"offset":[0,0],"textsize":"10px","textOnly":false,"className":"","sticky":true},null]},{"method":"addLayersControl","args":[[],["sfReseau","sfBVI"],{"collapsed":true,"autoZIndex":true,"position":"topleft"}]}],"limits":{"lat":[46.74871312138094,47.05598925544965],"lng":[-1.405807038182794,-0.9633392968252998]},"fitBounds":[46.84514520224651,-1.370098590461308,47.05552408903289,-1.07920635863741,[]]},"evals":[],"jsHooks":[]}</script> --- # Cas d'étude: chargement des données ## Exercice - Créer un R Notebook sous RStudio - Enregistrer le dossier `data` et le R Notebook dans le même dossier - Définir le working directory (Menu Session > Set working directory > To source file location) - Créer un premier chunk pour - effacer l'environnement de travail - charger les données du fichier `data/Maine_SevreNantaise.RData` - afficher la liste des variables présentes dans l'environnement de travail -- ``` r rm(list = ls()) load("data/Maine_SevreNantaise.RData") ls() #> [1] "BasinObs" "connections_reservoir_use" #> [3] "DatesR" "Q_water_use" #> [5] "reseau_maine" ``` --- # Cas d'étude: chargement des données ## Description des données fournies ``` #> [1] "BasinObs" "connections_reservoir_use" #> [3] "DatesR" "Q_water_use" #> [5] "reseau_maine" ``` - `BasinObs` (*list*): Données hydroclimatiques de chaque bassin - `connections_reservoir_use` (*data.frame*): Description du réseau des influences - `DatesR` (*POSIXct*): Dates des pas de temps d'observation - `Q_water_use` (*matrix*): Débits observés à la prise et à la restitution des réservoirs - `reseau_maine` (*data.frame*): Description du réseau de stations hydrométriques --- # Cas d'étude: description du réseau de stations hydrométriques Exercice: *Afficher le tableau décrivant le réseau de stations* -- ``` r reseau_maine ``` |id | area|nom |id_aval | distance_aval| |:---------------------|------:|:-------------------------------------------|:---------------|-------------:| |Maine_Remouille | 594.58|La Maine a Remouille |NA | NA| |GrandeMaine_StFulgent | 132.17|La Grande Maine a Saint-Fulgent |Maine_Remouille | 33.41| |PetiteMaine_StGeorges | 191.82|La Petite Maine a Saint-Georges-de-Montaigu |Maine_Remouille | 23.07| --- # Cas d'étude: description des observations hydroclimatiques Exercice: *Afficher des informations concernant les dates de chroniques (début/fin) et la structure des données d'observation* -- ``` r head(DatesR, 3) #> [1] "2008-01-01 UTC" "2008-01-02 UTC" "2008-01-03 UTC" ``` ``` r tail(DatesR, 3) #> [1] "2020-12-29 UTC" "2020-12-30 UTC" "2020-12-31 UTC" ``` ``` r str(BasinObs) #> List of 3 #> $ Maine_Remouille :'data.frame': 4749 obs. of 3 variables: #> ..$ P: num [1:4749] 0 0.121 4.842 11.442 9.305 ... #> ..$ E: num [1:4749] 0.4 0.3 0.3 0.2 0.3 ... #> ..$ Q: num [1:4749] 0.278 0.17 0.189 0.219 0.722 ... #> $ GrandeMaine_StFulgent:'data.frame': 4749 obs. of 3 variables: #> ..$ P: num [1:4749] 0 0.2 4.7 12.2 9.4 17.7 0.5 1.9 2 4 ... #> ..$ E: num [1:4749] 0.4 0.3 0.3 0.2 0.3 0.3 0.5 0.3 0.3 0.2 ... #> ..$ Q: num [1:4749] 0.244 0.243 0.268 0.307 1.405 ... #> $ PetiteMaine_StGeorges:'data.frame': 4749 obs. of 3 variables: #> ..$ P: num [1:4749] 0 0.2 5 11.6 9.7 16.7 0.4 2 2.1 4.4 ... #> ..$ E: num [1:4749] 0.4 0.3 0.3 0.2 0.3 0.3 0.4 0.3 0.3 0.2 ... #> ..$ Q: num [1:4749] 0.146 0.149 0.159 0.21 1.221 ... ``` --- # Cas d'étude: description des influences ## Exercice *Afficher le tableau de description du réseau des influences* -- ```r connections_reservoir_use ``` |id |id_aval | distance_aval| lambert2.x| lambert2.y| |:-------------|:---------------------|-------------:|----------:|----------:| |P_BULTIERE |Maine_Remouille | 28.035979| 329665.5| 2219226| |R_BULTIERE |Maine_Remouille | 28.035979| 329665.5| 2219226| |P_MAINE |Maine_Remouille | 19.115837| 324618.9| 2223265| |R_MAINE |Maine_Remouille | 21.844591| 327528.5| 2221721| |P_GRANDEMAINE |GrandeMaine_StFulgent | 9.025684| 339235.6| 2215470| |R_GRANDEMAINE |GrandeMaine_StFulgent | 10.418234| 339746.9| 2214639| |P_PETITEMAINE |PetiteMaine_StGeorges | 15.472812| 327524.9| 2211222| |R_PETITEMAINE |PetiteMaine_StGeorges | 15.472703| 328396.7| 2211426| --- # Cas d'étude: description des influences ## Exercice *Afficher les données des influences de débit (`str`, `summary`, `plot`)* (Attention aux unités) -- ``` r str(Q_water_use) #> num [1:4749, 1:8] -0.438 -0.438 -0.829 -2.223 -4.083 ... #> - attr(*, "dimnames")=List of 2 #> ..$ : chr [1:4749] "732" "733" "734" "735" ... #> ..$ : chr [1:8] "P_BULTIERE" "R_BULTIERE" "P_MAINE" "R_MAINE" ... ``` ``` r summary(Q_water_use, echo=Sys.getenv("FORMATION_HIDE_SOLUTIONS")=="") #> P_BULTIERE R_BULTIERE P_MAINE R_MAINE #> Min. :-23.87281 Min. : 0.01599 Min. :-1.098290 Min. :0.02074 #> 1st Qu.: -1.40618 1st Qu.: 0.15800 1st Qu.:-0.139724 1st Qu.:0.06151 #> Median : -0.35891 Median : 0.16000 Median :-0.057601 Median :0.07197 #> Mean : -1.24848 Mean : 0.99701 Mean :-0.092263 Mean :0.07599 #> 3rd Qu.: -0.06629 3rd Qu.: 0.95764 3rd Qu.:-0.030761 3rd Qu.:0.08801 #> Max. : 0.00000 Max. :25.00000 Max. :-0.009128 Max. :0.26807 #> P_GRANDEMAINE R_GRANDEMAINE P_PETITEMAINE R_PETITEMAINE #> Min. :-0.386404 Min. :0.01875 Min. :-0.92594 Min. :0.01175 #> 1st Qu.:-0.037282 1st Qu.:0.03265 1st Qu.:-0.06940 1st Qu.:0.02482 #> Median :-0.020975 Median :0.03892 Median :-0.03683 Median :0.02967 #> Mean :-0.029479 Mean :0.04408 Mean :-0.05502 Mean :0.03367 #> 3rd Qu.:-0.012459 3rd Qu.:0.05043 3rd Qu.:-0.02440 3rd Qu.:0.03916 #> Max. :-0.006554 Max. :0.15985 Max. :-0.01015 Max. :0.14649 ``` --- # Cas d'étude: description des influences ## Exercice *Afficher les données des influences de débit (`str`, `summary`, `plot`)* -- ``` r matplot(DatesR, Q_water_use[, c("P_BULTIERE", "R_BULTIERE")], type = "l") ``` <img src="index_files/figure-html/plot-bultiere-1.png" width="65%" style="display: block; margin: auto;" /> --- class: inverse, middle, center # Créer un réseau de modèles semi-distribués avec *airGRiwrm* --- # Créer un réseau de modèles semi-distribués avec *airGRiwrm* Le réseau est décrit par un objet `GRiwrm` qu'il faut créer&nbsp;: ```r ?CreateGRiwrm ``` ## Exercice *Quelles colonnes manquent à nos tableaux ?* -- Il faut définir quel modèle on souhaite utiliser pour chaque nœud du réseau&nbsp;: - Un modèle hydrologique pour les nœuds correspondants à une station hydrométrique; - ... autre chose pour les nœuds utilisés pour les influences humaines observées. -- ``` r reseau_maine$model <- "RunModel_GR4J" connections_reservoir_use$model <- NA ``` --- # Créer un réseau de modèles semi-distribués avec *airGRiwrm* ## Exercice *Fusionner les deux descriptions de réseau (les colonnes doivent être identiques)* -- ``` r names(reseau_maine) #> [1] "id" "lambert2.x" "lambert2.y" "area" #> [5] "nom" "id_aval" "distance_aval" "model" ``` ``` r names(connections_reservoir_use) #> [1] "id" "id_aval" "distance_aval" "lambert2.x" #> [5] "lambert2.y" "model" ``` ``` r setdiff(names(reseau_maine), names(connections_reservoir_use)) #> [1] "area" "nom" ``` Il faut ajouter les colonnes manquantes avec les bonnes valeurs et/ou supprimer les colonnes superflues. -- ``` r connections_reservoir_use$area <- NA connections_reservoir_use$nom <- connections_reservoir_use$id ``` --- # Créer un réseau de modèles semi-distribués avec *airGRiwrm* ## Exercice *Fusionner les deux descriptions de réseau* (`rbind`) -- ``` r reseau_inf <- rbind(reseau_maine, connections_reservoir_use) ``` |id | lambert2.x| lambert2.y| area|nom |id_aval | distance_aval|model | |:---------------------|----------:|----------:|------:|:-------------------------------------------|:---------------------|-------------:|:-------------| |Maine_Remouille | 318515.0| 2235050| 594.58|La Maine a Remouille |NA | NA|RunModel_GR4J | |GrandeMaine_StFulgent | 333911.0| 2217450| 132.17|La Grande Maine a Saint-Fulgent |Maine_Remouille | 33.410000|RunModel_GR4J | |PetiteMaine_StGeorges | 325104.0| 2219950| 191.82|La Petite Maine a Saint-Georges-de-Montaigu |Maine_Remouille | 23.070000|RunModel_GR4J | |P_BULTIERE | 329665.5| 2219226| NA|P_BULTIERE |Maine_Remouille | 28.035979|NA | |R_BULTIERE | 329665.5| 2219226| NA|R_BULTIERE |Maine_Remouille | 28.035979|NA | |P_MAINE | 324618.9| 2223265| NA|P_MAINE |Maine_Remouille | 19.115837|NA | |R_MAINE | 327528.5| 2221721| NA|R_MAINE |Maine_Remouille | 21.844591|NA | |P_GRANDEMAINE | 339235.6| 2215470| NA|P_GRANDEMAINE |GrandeMaine_StFulgent | 9.025684|NA | |R_GRANDEMAINE | 339746.9| 2214639| NA|R_GRANDEMAINE |GrandeMaine_StFulgent | 10.418234|NA | |P_PETITEMAINE | 327524.9| 2211222| NA|P_PETITEMAINE |PetiteMaine_StGeorges | 15.472812|NA | |R_PETITEMAINE | 328396.7| 2211426| NA|R_PETITEMAINE |PetiteMaine_StGeorges | 15.472703|NA | --- # Création de l'objet `GRiwrm` On peut faire correspondre les noms de colonne avec le paramètre `cols`&nbsp;: ## Exercice *Créer l'objet GRiwrm à partir du réseau fusionné* -- ``` r griwrm <- CreateGRiwrm( reseau_inf, cols = list(down = "id_aval", length = "distance_aval")) str(griwrm) #> GRiwrm [11 × 6] (S3: GRiwrm/tbl_df/tbl/data.frame) #> $ id : chr [1:11] "GrandeMaine_StFulgent" "PetiteMaine_StGeorges" "Maine_Remouille" "P_BULTIERE" ... #> $ down : chr [1:11] "Maine_Remouille" "Maine_Remouille" NA "Maine_Remouille" ... #> $ length: num [1:11] 33.4 23.1 NA 28 28 ... #> $ area : num [1:11] 132 192 595 NA NA ... #> $ model : chr [1:11] "RunModel_GR4J" "RunModel_GR4J" "RunModel_GR4J" NA ... #> $ donor : chr [1:11] "GrandeMaine_StFulgent" "PetiteMaine_StGeorges" "Maine_Remouille" "P_BULTIERE" ... ``` --- # Créer un réseau de modèles semi-distribués avec **airGRiwrm** La fonction `plot.GRiwrm` permet de tracer une représentation schématique du réseau. ```r ?plot.GRiwrm ``` Exercice: *afficher le schéma du réseau en mode "top-bottom"* -- ``` r plot(griwrm, orientation = "BT") ``` <img src="index_files/figure-html/plot.GRiwrm-1.png" width="100%" style="display: block; margin: auto;" /> --- class: inverse, middle, center # Calage d'un réseau de modèles hydrologiques semi-distribués avec *airGRiwrm* --- # `CreateInputsModel` dans *airGRiwrm* ```r ?CreateInputsModel.GRiwrm ``` -- ## Exercice *Créer des matrices pour les "inputs" avec le code de la station en nom de colonne* -- ``` r P <- cbind(BasinObs[[1]]$P, BasinObs[[2]]$P, BasinObs[[3]]$P) str(P) #> num [1:4749, 1:3] 0 0.121 4.842 11.442 9.305 ... ``` ``` r colnames(P) <- names(BasinObs) str(P) #> num [1:4749, 1:3] 0 0.121 4.842 11.442 9.305 ... #> - attr(*, "dimnames")=List of 2 #> ..$ : NULL #> ..$ : chr [1:3] "Maine_Remouille" "GrandeMaine_StFulgent" "PetiteMaine_StGeorges" ``` Heureusement qu'il n'y a que 3 colonnes ! Mais ne peut-on pas mieux faire ? --- # `CreateInputsModel` dans *airGRiwrm* ## Exercice *Créer des matrices pour les "inputs" avec le code de la station en nom de colonne (master class)* -- ``` r l <- lapply(BasinObs, function(x) {x$P}) str(l) #> List of 3 #> $ Maine_Remouille : num [1:4749] 0 0.121 4.842 11.442 9.305 ... #> $ GrandeMaine_StFulgent: num [1:4749] 0 0.2 4.7 12.2 9.4 17.7 0.5 1.9 2 4 ... #> $ PetiteMaine_StGeorges: num [1:4749] 0 0.2 5 11.6 9.7 16.7 0.4 2 2.1 4.4 ... ``` ``` r P <- do.call(cbind, l) str(P) #> num [1:4749, 1:3] 0 0.121 4.842 11.442 9.305 ... #> - attr(*, "dimnames")=List of 2 #> ..$ : NULL #> ..$ : chr [1:3] "Maine_Remouille" "GrandeMaine_StFulgent" "PetiteMaine_StGeorges" ``` -- En une seule ligne, avec la fonction `"[["` et un argument&nbsp;: ``` r E <- do.call(cbind, lapply(BasinObs, "[[", "E")) ``` --- # `CreateInputsModel` dans *airGRiwrm* ## Exercice *Créer l'objet `InputsModel` avec le réseau, les dates, les entrées climatiques et les influences observées (Attention aux unités des prélèvements/rejets)* -- ``` r IM <- CreateInputsModel( griwrm, DatesR = DatesR, Precip = P, PotEvap = E, Qobs = Q_water_use * 86400 # m3/s -> m3/j ) #> Warning in CreateInputsModel.GRiwrm(griwrm, DatesR = DatesR, Precip = P, : The #> usage of 'Qobs' is deprecated, use 'Qinf' instead #> CreateInputsModel.GRiwrm: Processing sub-basin GrandeMaine_StFulgent... #> CreateInputsModel.GRiwrm: Processing sub-basin PetiteMaine_StGeorges... #> CreateInputsModel.GRiwrm: Processing sub-basin Maine_Remouille... ``` --- # L'objet `GRiwrmInputsModel` L'objet `GRiwrmInputsModel` est une liste d'`InputsModel` : 1 élément de liste = 1 modèle *GR ( + Lag )* ``` r str(IM) #> List of 3 #> $ GrandeMaine_StFulgent:List of 19 #> ..$ DatesR : POSIXlt[1:4749], format: "2008-01-01" "2008-01-02" ... #> ..$ Precip : num [1:4749] 0 0.2 4.7 12.2 9.4 17.7 0.5 1.9 2 4 ... #> ..$ PotEvap : num [1:4749] 0.4 0.3 0.3 0.2 0.3 0.3 0.5 0.3 0.3 0.2 ... #> ..$ Qupstream : num [1:4749, 1:2] -943 -1009 -1009 -989 -947 ... #> .. ..- attr(*, "dimnames")=List of 2 #> .. .. ..$ : NULL #> .. .. ..$ : chr [1:2] "P_GRANDEMAINE" "R_GRANDEMAINE" #> ..$ LengthHydro : Named num [1:2] 9.03 10.42 #> .. ..- attr(*, "names")= chr [1:2] "P_GRANDEMAINE" "R_GRANDEMAINE" #> ..$ BasinAreas : Named num [1:3] NA NA 132 #> .. ..- attr(*, "names")= chr [1:3] "P_GRANDEMAINE" "R_GRANDEMAINE" "GrandeMaine_StFulgent" #> ..$ id : chr "GrandeMaine_StFulgent" #> ..$ down : chr "Maine_Remouille" #> ..$ UpstreamNodes : chr [1:2] "P_GRANDEMAINE" "R_GRANDEMAINE" #> ..$ UpstreamIsModeled: Named logi [1:2] FALSE FALSE #> .. ..- attr(*, "names")= chr [1:2] "P_GRANDEMAINE" "R_GRANDEMAINE" #> ..$ UpstreamVarQ : Named chr [1:2] "Qsim_m3" "Qsim_m3" #> .. ..- attr(*, "names")= chr [1:2] "P_GRANDEMAINE" "R_GRANDEMAINE" #> ..$ FUN_MOD : chr "RunModel_GR4J" #> ..$ inUngaugedCluster: logi FALSE #> ..$ isReceiver : logi FALSE #> ..$ gaugedId : chr "GrandeMaine_StFulgent" #> ..$ hasUngaugedNodes : logi FALSE #> ..$ model :List of 4 #> .. ..$ indexParamUngauged: num [1:5] 1 2 3 4 5 #> .. ..$ hasX4 : logi TRUE #> .. ..$ iX4 : num 5 #> .. ..$ IsHyst : logi FALSE #> ..$ hasDiversion : logi FALSE #> ..$ isReservoir : logi FALSE #> ..- attr(*, "class")= chr [1:5] "RunModel_GR4J" "InputsModel" "daily" "GR" ... #> $ PetiteMaine_StGeorges:List of 19 #> ..$ DatesR : POSIXlt[1:4749], format: "2008-01-01" "2008-01-02" ... #> ..$ Precip : num [1:4749] 0 0.2 5 11.6 9.7 16.7 0.4 2 2.1 4.4 ... #> ..$ PotEvap : num [1:4749] 0.4 0.3 0.3 0.2 0.3 0.3 0.4 0.3 0.3 0.2 ... #> ..$ Qupstream : num [1:4749, 1:2] -1771 -1891 -1891 -1859 -1424 ... #> .. ..- attr(*, "dimnames")=List of 2 #> .. .. ..$ : NULL #> .. .. ..$ : chr [1:2] "P_PETITEMAINE" "R_PETITEMAINE" #> ..$ LengthHydro : Named num [1:2] 15.5 15.5 #> .. ..- attr(*, "names")= chr [1:2] "P_PETITEMAINE" "R_PETITEMAINE" #> ..$ BasinAreas : Named num [1:3] NA NA 192 #> .. ..- attr(*, "names")= chr [1:3] "P_PETITEMAINE" "R_PETITEMAINE" "PetiteMaine_StGeorges" #> ..$ id : chr "PetiteMaine_StGeorges" #> ..$ down : chr "Maine_Remouille" #> ..$ UpstreamNodes : chr [1:2] "P_PETITEMAINE" "R_PETITEMAINE" #> ..$ UpstreamIsModeled: Named logi [1:2] FALSE FALSE #> .. ..- attr(*, "names")= chr [1:2] "P_PETITEMAINE" "R_PETITEMAINE" #> ..$ UpstreamVarQ : Named chr [1:2] "Qsim_m3" "Qsim_m3" #> .. ..- attr(*, "names")= chr [1:2] "P_PETITEMAINE" "R_PETITEMAINE" #> ..$ FUN_MOD : chr "RunModel_GR4J" #> ..$ inUngaugedCluster: logi FALSE #> ..$ isReceiver : logi FALSE #> ..$ gaugedId : chr "PetiteMaine_StGeorges" #> ..$ hasUngaugedNodes : logi FALSE #> ..$ model :List of 4 #> .. ..$ indexParamUngauged: num [1:5] 1 2 3 4 5 #> .. ..$ hasX4 : logi TRUE #> .. ..$ iX4 : num 5 #> .. ..$ IsHyst : logi FALSE #> ..$ hasDiversion : logi FALSE #> ..$ isReservoir : logi FALSE #> ..- attr(*, "class")= chr [1:5] "RunModel_GR4J" "InputsModel" "daily" "GR" ... #> $ Maine_Remouille :List of 19 #> ..$ DatesR : POSIXlt[1:4749], format: "2008-01-01" "2008-01-02" ... #> ..$ Precip : num [1:4749] 0 0.121 4.842 11.442 9.305 ... #> ..$ PotEvap : num [1:4749] 0.4 0.3 0.3 0.2 0.3 ... #> ..$ Qupstream : num [1:4749, 1:6] 0 0 0 0 0 0 0 0 0 0 ... #> .. ..- attr(*, "dimnames")=List of 2 #> .. .. ..$ : NULL #> .. .. ..$ : chr [1:6] "GrandeMaine_StFulgent" "PetiteMaine_StGeorges" "P_BULTIERE" "R_BULTIERE" ... #> ..$ LengthHydro : Named num [1:6] 33.4 23.1 28 28 19.1 ... #> .. ..- attr(*, "names")= chr [1:6] "GrandeMaine_StFulgent" "PetiteMaine_StGeorges" "P_BULTIERE" "R_BULTIERE" ... #> ..$ BasinAreas : Named num [1:7] 132 192 NA NA NA ... #> .. ..- attr(*, "names")= chr [1:7] "GrandeMaine_StFulgent" "PetiteMaine_StGeorges" "P_BULTIERE" "R_BULTIERE" ... #> ..$ id : chr "Maine_Remouille" #> ..$ down : chr NA #> ..$ UpstreamNodes : chr [1:6] "GrandeMaine_StFulgent" "PetiteMaine_StGeorges" "P_BULTIERE" "R_BULTIERE" ... #> ..$ UpstreamIsModeled: Named logi [1:6] TRUE TRUE FALSE FALSE FALSE FALSE #> .. ..- attr(*, "names")= chr [1:6] "GrandeMaine_StFulgent" "PetiteMaine_StGeorges" "P_BULTIERE" "R_BULTIERE" ... #> ..$ UpstreamVarQ : Named chr [1:6] "Qsim_m3" "Qsim_m3" "Qsim_m3" "Qsim_m3" ... #> .. ..- attr(*, "names")= chr [1:6] "GrandeMaine_StFulgent" "PetiteMaine_StGeorges" "P_BULTIERE" "R_BULTIERE" ... #> ..$ FUN_MOD : chr "RunModel_GR4J" #> ..$ inUngaugedCluster: logi FALSE #> ..$ isReceiver : logi FALSE #> ..$ gaugedId : chr "Maine_Remouille" #> ..$ hasUngaugedNodes : logi FALSE #> ..$ model :List of 4 #> .. ..$ indexParamUngauged: num [1:5] 1 2 3 4 5 #> .. ..$ hasX4 : logi TRUE #> .. ..$ iX4 : num 5 #> .. ..$ IsHyst : logi FALSE #> ..$ hasDiversion : logi FALSE #> ..$ isReservoir : logi FALSE #> ..- attr(*, "class")= chr [1:5] "RunModel_GR4J" "InputsModel" "daily" "GR" ... #> - attr(*, "class")= chr [1:2] "GRiwrmInputsModel" "list" #> - attr(*, "GRiwrm")= GRiwrm [11 × 6] (S3: GRiwrm/tbl_df/tbl/data.frame) #> ..$ id : chr [1:11] "GrandeMaine_StFulgent" "PetiteMaine_StGeorges" "Maine_Remouille" "P_BULTIERE" ... #> ..$ down : chr [1:11] "Maine_Remouille" "Maine_Remouille" NA "Maine_Remouille" ... #> ..$ length: num [1:11] 33.4 23.1 NA 28 28 ... #> ..$ area : num [1:11] 132 192 595 NA NA ... #> ..$ model : chr [1:11] "RunModel_GR4J" "RunModel_GR4J" "RunModel_GR4J" NA ... #> ..$ donor : chr [1:11] "GrandeMaine_StFulgent" "PetiteMaine_StGeorges" "Maine_Remouille" "P_BULTIERE" ... #> - attr(*, "TimeStep")= num 86400 ``` --- # Calage avec *airGRiwrm*: CreateRunOptions Définition des options de simulation: - Pas de temps de simulation (*warm-up* + run) - *états initiaux du modèle* - *Liste des variables en sortie du modèles* ```r ?CreateRunOptions ``` ## Exercice Définir les indices des pas de temps des périodes de : - initialisation du modèle (Warm-up) sur un an en début de chronique - simulation (Run) sur le reste de la chronique -- ``` r I_WarmUp <- 1:365 I_Run <- seq(366, length(DatesR)) ``` --- # Calage avec *airGRiwrm*: CreateRunOptions ## Exercice *Créer l'objet `RunOptions` à partir de l'objet `InputsModel` et des indices des périodes de chauffe et simulation* -- ``` r RO <- CreateRunOptions( IM, IndPeriod_WarmUp = I_WarmUp, IndPeriod_Run = I_Run ) ``` -- ``` r str(RO) #> List of 3 #> $ GrandeMaine_StFulgent:List of 9 #> ..$ IndPeriod_WarmUp: int [1:365] 1 2 3 4 5 6 7 8 9 10 ... #> ..$ IndPeriod_Run : int [1:4384] 366 367 368 369 370 371 372 373 374 375 ... #> ..$ IniStates : num [1:67] 0 0 0 0 0 0 0 0 0 0 ... #> ..$ IniResLevels : num [1:4] 0.3 0.5 NA NA #> ..$ Outputs_Cal : chr [1:2] "Qsim" "Param" #> ..$ Outputs_Sim : Named chr [1:24] "DatesR" "PotEvap" "Precip" "Prod" ... #> .. ..- attr(*, "names")= chr [1:24] "" "GR1" "GR2" "GR3" ... #> ..$ FortranOutputs :List of 2 #> .. ..$ GR: chr [1:18] "PotEvap" "Precip" "Prod" "Pn" ... #> .. ..$ CN: NULL #> ..$ FeatFUN_MOD :List of 12 #> .. ..$ CodeMod : chr "GR4J" #> .. ..$ NameMod : chr "GR4J" #> .. ..$ NbParam : num 5 #> .. ..$ TimeUnit : chr "daily" #> .. ..$ Id : logi NA #> .. ..$ Class : chr [1:2] "daily" "GR" #> .. ..$ Pkg : chr "airGR" #> .. ..$ NameFunMod : chr "RunModel_GR4J" #> .. ..$ TimeStep : num 86400 #> .. ..$ TimeStepMean: int 86400 #> .. ..$ CodeModHydro: chr "GR4J" #> .. ..$ IsSD : logi TRUE #> ..$ id : chr "GrandeMaine_StFulgent" #> ..- attr(*, "class")= chr [1:3] "RunOptions" "daily" "GR" #> $ PetiteMaine_StGeorges:List of 9 #> ..$ IndPeriod_WarmUp: int [1:365] 1 2 3 4 5 6 7 8 9 10 ... #> ..$ IndPeriod_Run : int [1:4384] 366 367 368 369 370 371 372 373 374 375 ... #> ..$ IniStates : num [1:67] 0 0 0 0 0 0 0 0 0 0 ... #> ..$ IniResLevels : num [1:4] 0.3 0.5 NA NA #> ..$ Outputs_Cal : chr [1:2] "Qsim" "Param" #> ..$ Outputs_Sim : Named chr [1:24] "DatesR" "PotEvap" "Precip" "Prod" ... #> .. ..- attr(*, "names")= chr [1:24] "" "GR1" "GR2" "GR3" ... #> ..$ FortranOutputs :List of 2 #> .. ..$ GR: chr [1:18] "PotEvap" "Precip" "Prod" "Pn" ... #> .. ..$ CN: NULL #> ..$ FeatFUN_MOD :List of 12 #> .. ..$ CodeMod : chr "GR4J" #> .. ..$ NameMod : chr "GR4J" #> .. ..$ NbParam : num 5 #> .. ..$ TimeUnit : chr "daily" #> .. ..$ Id : logi NA #> .. ..$ Class : chr [1:2] "daily" "GR" #> .. ..$ Pkg : chr "airGR" #> .. ..$ NameFunMod : chr "RunModel_GR4J" #> .. ..$ TimeStep : num 86400 #> .. ..$ TimeStepMean: int 86400 #> .. ..$ CodeModHydro: chr "GR4J" #> .. ..$ IsSD : logi TRUE #> ..$ id : chr "PetiteMaine_StGeorges" #> ..- attr(*, "class")= chr [1:3] "RunOptions" "daily" "GR" #> $ Maine_Remouille :List of 9 #> ..$ IndPeriod_WarmUp: int [1:365] 1 2 3 4 5 6 7 8 9 10 ... #> ..$ IndPeriod_Run : int [1:4384] 366 367 368 369 370 371 372 373 374 375 ... #> ..$ IniStates : num [1:67] 0 0 0 0 0 0 0 0 0 0 ... #> ..$ IniResLevels : num [1:4] 0.3 0.5 NA NA #> ..$ Outputs_Cal : chr [1:2] "Qsim" "Param" #> ..$ Outputs_Sim : Named chr [1:24] "DatesR" "PotEvap" "Precip" "Prod" ... #> .. ..- attr(*, "names")= chr [1:24] "" "GR1" "GR2" "GR3" ... #> ..$ FortranOutputs :List of 2 #> .. ..$ GR: chr [1:18] "PotEvap" "Precip" "Prod" "Pn" ... #> .. ..$ CN: NULL #> ..$ FeatFUN_MOD :List of 12 #> .. ..$ CodeMod : chr "GR4J" #> .. ..$ NameMod : chr "GR4J" #> .. ..$ NbParam : num 5 #> .. ..$ TimeUnit : chr "daily" #> .. ..$ Id : logi NA #> .. ..$ Class : chr [1:2] "daily" "GR" #> .. ..$ Pkg : chr "airGR" #> .. ..$ NameFunMod : chr "RunModel_GR4J" #> .. ..$ TimeStep : num 86400 #> .. ..$ TimeStepMean: int 86400 #> .. ..$ CodeModHydro: chr "GR4J" #> .. ..$ IsSD : logi TRUE #> ..$ id : chr "Maine_Remouille" #> ..- attr(*, "class")= chr [1:3] "RunOptions" "daily" "GR" #> - attr(*, "class")= chr [1:2] "list" "GRiwrmRunOptions" ``` --- # Calage avec *airGRiwrm*: `CreateInputsCrit` Définition des critères de calage&nbsp;: ```r ?CreateInputsCrit.GRiwrmInputsModel ``` ## Exercice *Créer un critère de calage utilisant les débit observées au stations hydrométriques, un critère **KGE2** et une transformation "Racine carrée" des débits* -- ``` r Q <- do.call(cbind, lapply(BasinObs, "[[", "Q")) IC <- CreateInputsCrit( IM, FUN_CRIT = airGR::ErrorCrit_KGE2, RunOptions = RO, Obs = Q[I_Run,], transfo = "sqrt" ) ``` -- ``` r str(IC) #> List of 3 #> $ GrandeMaine_StFulgent:List of 8 #> ..$ FUN_CRIT:function (InputsCrit, OutputsModel, warnings = TRUE, verbose = TRUE) #> .. ..- attr(*, "class")= chr [1:2] "FUN_CRIT" "function" #> ..$ Obs : num [1:4384] 0.409 0.386 0.377 0.352 0.337 0.335 0.333 0.279 0.291 0.293 ... #> ..$ VarObs : chr "Q" #> ..$ BoolCrit: logi [1:4384] TRUE TRUE TRUE TRUE TRUE TRUE ... #> ..$ idLayer : logi NA #> ..$ transfo : chr "sqrt" #> ..$ epsilon : NULL #> ..$ Weights : NULL #> ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit" #> $ PetiteMaine_StGeorges:List of 8 #> ..$ FUN_CRIT:function (InputsCrit, OutputsModel, warnings = TRUE, verbose = TRUE) #> .. ..- attr(*, "class")= chr [1:2] "FUN_CRIT" "function" #> ..$ Obs : num [1:4384] 0.278 0.273 0.246 0.245 0.227 0.211 0.222 0.203 0.191 0.118 ... #> ..$ VarObs : chr "Q" #> ..$ BoolCrit: logi [1:4384] TRUE TRUE TRUE TRUE TRUE TRUE ... #> ..$ idLayer : logi NA #> ..$ transfo : chr "sqrt" #> ..$ epsilon : NULL #> ..$ Weights : NULL #> ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit" #> $ Maine_Remouille :List of 8 #> ..$ FUN_CRIT:function (InputsCrit, OutputsModel, warnings = TRUE, verbose = TRUE) #> .. ..- attr(*, "class")= chr [1:2] "FUN_CRIT" "function" #> ..$ Obs : num [1:4384] 0.43 0.43 0.372 0.224 0.246 0.34 0.212 0.201 0.179 0.176 ... #> ..$ VarObs : chr "Q" #> ..$ BoolCrit: logi [1:4384] TRUE TRUE TRUE TRUE TRUE TRUE ... #> ..$ idLayer : logi NA #> ..$ transfo : chr "sqrt" #> ..$ epsilon : NULL #> ..$ Weights : NULL #> ..- attr(*, "class")= chr [1:2] "Single" "InputsCrit" #> - attr(*, "class")= chr [1:2] "GRiwrmInputsCrit" "list" ``` --- # Calage avec *airGRiwrm*: `CreateCalibOptions` ```r ?CreateInputsCrit.GRiwrmInputsModel ``` ## Exercice *Créer l'objet `CalibOptions` avec les options par défaut* ``` r CO <- CreateCalibOptions(IM) ``` ``` r str(CO) #> List of 3 #> $ GrandeMaine_StFulgent:List of 4 #> ..$ FixedParam : logi [1:5] NA NA NA NA NA #> ..$ SearchRanges : num [1:2, 1:5] 1.00e-02 2.00e+01 4.59e-05 2.18e+04 -1.09e+04 ... #> ..$ FUN_TRANSFO :function (ParamIn, Direction) #> ..$ StartParamDistrib: num [1:3, 1:5] 1.25 2.5 5 169.02 247.15 ... #> ..- attr(*, "class")= chr [1:5] "CalibOptions" "daily" "GR" "SD" ... #> $ PetiteMaine_StGeorges:List of 4 #> ..$ FixedParam : logi [1:5] NA NA NA NA NA #> ..$ SearchRanges : num [1:2, 1:5] 1.00e-02 2.00e+01 4.59e-05 2.18e+04 -1.09e+04 ... #> ..$ FUN_TRANSFO :function (ParamIn, Direction) #> ..$ StartParamDistrib: num [1:3, 1:5] 1.25 2.5 5 169.02 247.15 ... #> ..- attr(*, "class")= chr [1:5] "CalibOptions" "daily" "GR" "SD" ... #> $ Maine_Remouille :List of 4 #> ..$ FixedParam : logi [1:5] NA NA NA NA NA #> ..$ SearchRanges : num [1:2, 1:5] 1.00e-02 2.00e+01 4.59e-05 2.18e+04 -1.09e+04 ... #> ..$ FUN_TRANSFO :function (ParamIn, Direction) #> ..$ StartParamDistrib: num [1:3, 1:5] 1.25 2.5 5 169.02 247.15 ... #> ..- attr(*, "class")= chr [1:5] "CalibOptions" "daily" "GR" "SD" ... #> - attr(*, "class")= chr [1:2] "GRiwrmCalibOptions" "list" ``` --- # `Calibration` dans *airGRiwrm*: `Calibration` **airGRiwrm** opère un calage séquentiel amont-aval des modèles hydrologiques des sous-bassins. ```r ?Calibration.GRiwrmInputsModel ``` ## Exercice *Lancer la calibration du modèle avec les inputs (`InputsModel`), les options de simulations (`RunOptions`), les critères de calage (`InputsCrit`) et les options de calage (`CalibOptions`)* -- ``` r OC <- suppressWarnings( Calibration(IM, RunOptions = RO, InputsCrit = IC, CalibOptions = CO) ) #> Calibration.GRiwrmInputsModel: Processing sub-basin GrandeMaine_StFulgent... #> Grid-Screening in progress (0% 20% 40% 60% 80% 100%) #> Screening completed (243 runs) #> Param = 1.250, 169.017, -2.376, 42.098, 1.944 #> Crit. KGE2[sqrt(Q)] = 0.9173 #> Steepest-descent local search in progress #> Calibration completed (70 iterations, 964 runs) #> Param = 0.676, 231.109, -0.745, 15.290, 2.094 #> Crit. KGE2[sqrt(Q)] = 0.9421 #> Calibration.GRiwrmInputsModel: Processing sub-basin PetiteMaine_StGeorges... #> Grid-Screening in progress (0% 20% 40% 60% 80% 100%) #> Screening completed (243 runs) #> Param = 1.250, 169.017, -2.376, 20.697, 2.384 #> Crit. KGE2[sqrt(Q)] = 0.8888 #> Steepest-descent local search in progress #> Calibration completed (67 iterations, 931 runs) #> Param = 0.930, 198.835, -0.892, 10.753, 2.193 #> Crit. KGE2[sqrt(Q)] = 0.9417 #> Calibration.GRiwrmInputsModel: Processing sub-basin Maine_Remouille... #> Grid-Screening in progress (0% 20% 40% 60% 80% 100%) #> Screening completed (243 runs) #> Param = 5.000, 169.017, -0.649, 20.697, 1.944 #> Crit. KGE2[sqrt(Q)] = 0.9522 #> Steepest-descent local search in progress #> Calibration completed (40 iterations, 630 runs) #> Param = 5.454, 187.150, -0.528, 14.447, 2.329 #> Crit. KGE2[sqrt(Q)] = 0.9550 ``` --- # `Calibration` dans *airGRiwrm*: `Calibration` ``` r str(OC) #> List of 3 #> $ GrandeMaine_StFulgent:List of 9 #> ..$ ParamFinalR : num [1:5] 0.676 231.109 -0.745 15.29 2.094 #> ..$ CritFinal : num 0.942 #> ..$ NIter : num 70 #> ..$ NRuns : num 964 #> ..$ HistParamR : num [1:70, 1:5] 1.25 1.25 0.93 0.93 0.93 0.93 0.77 0.77 0.77 0.77 ... #> .. ..- attr(*, "dimnames")=List of 2 #> .. .. ..$ : NULL #> .. .. ..$ : chr [1:5] "Param1" "Param2" "Param3" "Param4" ... #> ..$ HistCrit : num [1:70, 1] 0.917 0.917 0.917 0.917 0.92 ... #> ..$ MatBoolCrit : logi [1:4384, 1:2] TRUE TRUE TRUE TRUE TRUE TRUE ... #> .. ..- attr(*, "dimnames")=List of 2 #> .. .. ..$ : NULL #> .. .. ..$ : chr [1:2] "BoolCrit_Requested" "BoolCrit_Actual" #> ..$ CritName : chr "KGE2[sqrt(Q)]" #> ..$ CritBestValue: num 1 #> ..- attr(*, "class")= chr [1:2] "OutputsCalib" "HBAN" #> $ PetiteMaine_StGeorges:List of 9 #> ..$ ParamFinalR : num [1:5] 0.93 198.835 -0.892 10.753 2.193 #> ..$ CritFinal : num 0.942 #> ..$ NIter : num 67 #> ..$ NRuns : num 931 #> ..$ HistParamR : num [1:67, 1:5] 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 ... #> .. ..- attr(*, "dimnames")=List of 2 #> .. .. ..$ : NULL #> .. .. ..$ : chr [1:5] "Param1" "Param2" "Param3" "Param4" ... #> ..$ HistCrit : num [1:67, 1] 0.889 0.889 0.92 0.931 0.931 ... #> ..$ MatBoolCrit : logi [1:4384, 1:2] TRUE TRUE TRUE TRUE TRUE TRUE ... #> .. ..- attr(*, "dimnames")=List of 2 #> .. .. ..$ : NULL #> .. .. ..$ : chr [1:2] "BoolCrit_Requested" "BoolCrit_Actual" #> ..$ CritName : chr "KGE2[sqrt(Q)]" #> ..$ CritBestValue: num 1 #> ..- attr(*, "class")= chr [1:2] "OutputsCalib" "HBAN" #> $ Maine_Remouille :List of 9 #> ..$ ParamFinalR : num [1:5] 5.454 187.15 -0.528 14.447 2.329 #> ..$ CritFinal : num 0.955 #> ..$ NIter : num 40 #> ..$ NRuns : num 630 #> ..$ HistParamR : num [1:40, 1:5] 5 5.64 6.28 6.28 6.28 6.28 6.28 6.28 6.12 5.96 ... #> .. ..- attr(*, "dimnames")=List of 2 #> .. .. ..$ : NULL #> .. .. ..$ : chr [1:5] "Param1" "Param2" "Param3" "Param4" ... #> ..$ HistCrit : num [1:40, 1] 0.952 0.952 0.952 0.952 0.952 ... #> ..$ MatBoolCrit : logi [1:4384, 1:2] TRUE TRUE TRUE TRUE TRUE TRUE ... #> .. ..- attr(*, "dimnames")=List of 2 #> .. .. ..$ : NULL #> .. .. ..$ : chr [1:2] "BoolCrit_Requested" "BoolCrit_Actual" #> ..$ CritName : chr "KGE2[sqrt(Q)]" #> ..$ CritBestValue: num 1 #> ..- attr(*, "class")= chr [1:2] "OutputsCalib" "HBAN" #> - attr(*, "class")= chr [1:2] "GRiwrmOutputsCalib" "list" ``` --- # `RunModel` dans *airGRiwrm* ```r ?RunModel.GRiwrmInputsModel ``` ## Exercice Récupérer les paramètres de calage et lancer le modèle avec les inputs et les options de simulation -- ``` r Params <- lapply(OC, "[[", "ParamFinalR") str(Params) #> List of 3 #> $ GrandeMaine_StFulgent: num [1:5] 0.676 231.109 -0.745 15.29 2.094 #> $ PetiteMaine_StGeorges: num [1:5] 0.93 198.835 -0.892 10.753 2.193 #> $ Maine_Remouille : num [1:5] 5.454 187.15 -0.528 14.447 2.329 ``` -- ``` r OM <- RunModel(IM, RunOptions = RO, Param = Params) #> RunModel.GRiwrmInputsModel: Processing sub-basin GrandeMaine_StFulgent... #> Warning in RunModel_Lag(InputsModel, RunOptions, Param[1], OutputsModel): 174 #> time steps with negative flow, set to zero. #> RunModel.GRiwrmInputsModel: Processing sub-basin PetiteMaine_StGeorges... #> Warning in RunModel_Lag(InputsModel, RunOptions, Param[1], OutputsModel): 459 #> time steps with negative flow, set to zero. #> RunModel.GRiwrmInputsModel: Processing sub-basin Maine_Remouille... #> Warning in RunModel_Lag(InputsModel, RunOptions, Param[1], OutputsModel): 83 #> time steps with negative flow, set to zero. ``` --- # Visualisation des données de simulation ```r ?plot.GRiwrmOutputsModel ``` ## Exercice Afficher les plots d'analyse de **airGR** avec les débits simulés/observés -- ``` r plot(OM, Q[I_Run,]) ``` <img src="index_files/figure-html/plot.GRiwrmOutputsModel-1.png" width="65%" style="display: block; margin: auto;" /><img src="index_files/figure-html/plot.GRiwrmOutputsModel-2.png" width="65%" style="display: block; margin: auto;" /><img src="index_files/figure-html/plot.GRiwrmOutputsModel-3.png" width="65%" style="display: block; margin: auto;" /> --- # Visualisation des données de simulation Le plot des débits simulés sur tous les noeuds du modèle en m3/s ```r ?plot.Qm3s ``` ## Exercice Afficher tous les débits simulés pour la première année de simulation -- ``` r Qm3s <- attr(OM, "Qm3s") plot(Qm3s[1:365,]) ``` <img src="index_files/figure-html/plot.Qm3s-1.png" width="65%" style="display: block; margin: auto;" /> --- class: inverse, middle, center # Modélisation du barrage de la Bultière <img src="img/Barrage_de_la_Bultiere.jpg" width="65%" style="display: block; margin: auto;" /> --- # Modélisation du barrage de la Bultière Nous implémenterons le barrage de la Bultière à l'aide d'un modèle de type "Reservoir" (`RunModel_Reservoir`) qui simule le volume stocké en fonction des entrées provenant des nœuds amont et d'un débit lâcher fourni par l'utilisateur. Le barrage peut stocker jusqu'à 5 millions de mètres cubes et est utilisé comme réserve pour un prélèvement AEP qui sera lui aussi modélisé. Il est réglé pour lâcher en permanence un débit réservé de 160 L/s. ``` #> Ungauged node 'Barrage_Bultiere' automatically gets the node 'Maine_Remouille' as parameter donor ``` <img src="index_files/figure-html/plot.GRiwrm.bultiere-1.png" width="65%" style="display: block; margin: auto;" /> --- # Installation de la version de développement de airGRiwrm Le noeud de type "Reservoir" n'est présent que dans la version de développement de **airGRiwrm**. ``` r install.packages("remotes") # Si nécessaire remotes::install_github("inrae/airGRiwrm", "dev") ``` --- # Chargement des données du barrage ## Exercice *Charger et explorer les variables `reseau_bultiere` et `Q_AEP_Bultiere` situées dans le fichier "`data/Barrage_Bultiere.RData`".* -- ``` r load("data/Barrage_Bultiere.RData") summary(Q_AEP_Bultiere) #> Min. 1st Qu. Median Mean 3rd Qu. Max. #> 2512 12328 14947 14569 16830 26128 ``` Le débit de prélèvement AEP est ici en m³/jour. --- # Chargement des données du barrage ## Exercice *Charger et explorer les variables `reseau_bultiere` et `Q_AEP_Bultiere` situées dans le fichier "`data/Barrage_Bultiere.RData`".* -- ``` r reseau_bultiere ``` |id |down | length|model | area| |:---------------------|:---------------------|---------:|:------------------|------:| |Maine_Remouille |NA | NA|RunModel_GR4J | 594.58| |GrandeMaine_StFulgent |Barrage_Bultiere | 5.370000|RunModel_GR4J | 132.17| |PetiteMaine_StGeorges |Maine_Remouille | 23.070000|RunModel_GR4J | 191.82| |P_MAINE |Maine_Remouille | 19.115837|NA | NA| |R_MAINE |Maine_Remouille | 21.844591|NA | NA| |P_GRANDEMAINE |GrandeMaine_StFulgent | 9.025684|NA | NA| |R_GRANDEMAINE |GrandeMaine_StFulgent | 10.418234|NA | NA| |P_PETITEMAINE |PetiteMaine_StGeorges | 15.472812|NA | NA| |R_PETITEMAINE |PetiteMaine_StGeorges | 15.472703|NA | NA| |Barrage_Bultiere |Maine_Remouille | 28.040000|RunModel_Reservoir | NA| |AEP_Bultiere |Barrage_Bultiere | 0.000000|NA | NA| --- # Création du réseau : `CreateGRiwrm` ## Exercice *Créer le réseau GRiwrm et afficher le schéma du réseau* -- ``` r g_bultiere <- CreateGRiwrm(reseau_bultiere) #> Ungauged node 'Barrage_Bultiere' automatically gets the node 'Maine_Remouille' as parameter donor ``` ``` r plot(g_bultiere) ``` <img src="index_files/figure-html/g_bultiere-1.png" width="65%" style="display: block; margin: auto;" /> --- # Initialiser les données d'entrée du modèle Les entrées du modèle comprennent les données météo (`Precip`, `PotEvap`), les prélèvements et rejets et le débit de lâcher du barrage (`Qobs`). Le barrage de la Bultière relâche 160 L/s toute l'année (plus en cas de surverse). Ce sera le débit imposé sur le modèle `RunModel_Reservoir` qui relâchera ce débit si le volume contenu dans le barrage le permet. ## Exercice *Créer une chronique constante de débit de lâchers programmés sur le barrage de la Bultière de 160 L/s pour une injection dans les données `Qobs` (en m³/jour)* -- ``` r Q_Barrage_Bultiere <- rep(160 / 1000 * 86400, length(DatesR)) str(Q_Barrage_Bultiere) #> num [1:4749] 13824 13824 13824 13824 13824 ... ``` --- # Initialiser les données d'entrée du modèle Les entrées du modèle comprennent les données météo (`Precip`, `PotEvap`), les prélèvements et rejets et le débit de lâcher du barrage (`Qobs`). ## Exercice *Mettre à jour la matrice des prélèvements/rejets qui sera injectée dans `CreateInputsModel` via le paramètre `Qobs`* -- ``` r str(Q_water_use) #> num [1:4749, 1:8] -0.438 -0.438 -0.829 -2.223 -4.083 ... #> - attr(*, "dimnames")=List of 2 #> ..$ : chr [1:4749] "732" "733" "734" "735" ... #> ..$ : chr [1:8] "P_BULTIERE" "R_BULTIERE" "P_MAINE" "R_MAINE" ... ``` ``` r Q_water_use2 <- Q_water_use[, !grepl("BULTIERE", colnames(Q_water_use))] Q_water_use2 <- cbind(Q_water_use2 * 86400, AEP_Bultiere = -Q_AEP_Bultiere, Barrage_Bultiere = Q_Barrage_Bultiere) str(Q_water_use2) #> num [1:4749, 1:8] -1966 -2116 -2115 -2090 -1417 ... #> - attr(*, "dimnames")=List of 2 #> ..$ : chr [1:4749] "732" "733" "734" "735" ... #> ..$ : chr [1:8] "P_MAINE" "R_MAINE" "P_GRANDEMAINE" "R_GRANDEMAINE" ... ``` --- # Initialiser les données d'entrée du modèle ## Exercice *Créer l'objet `InputsModel`* -- ``` r IM2 <- CreateInputsModel( g_bultiere, DatesR = DatesR, Precip = P, PotEvap = E, Qobs = Q_water_use2 ) #> Warning in CreateInputsModel.GRiwrm(g_bultiere, DatesR = DatesR, Precip = P, : #> The usage of 'Qobs' is deprecated, use 'Qinf' instead #> Warning in updateQinfQrelease(g = x, Qinf = Qinf, Qrelease = Qrelease): Use of the `Qinf` parameter for reservoir releases is deprecated, please use `Qrelease` instead. #> Processing `Qrelease <- cbind(Qrelease, Qinf[, c("Barrage_Bultiere"])`... #> CreateInputsModel.GRiwrm: Processing sub-basin GrandeMaine_StFulgent... #> CreateInputsModel.GRiwrm: Processing sub-basin PetiteMaine_StGeorges... #> CreateInputsModel.GRiwrm: Processing sub-basin Barrage_Bultiere... #> CreateInputsModel.GRiwrm: Processing sub-basin Maine_Remouille... ``` ``` r str(IM2) #> List of 4 #> $ GrandeMaine_StFulgent:List of 19 #> ..$ DatesR : POSIXlt[1:4749], format: "2008-01-01" "2008-01-02" ... #> ..$ Precip : num [1:4749] 0 0.2 4.7 12.2 9.4 17.7 0.5 1.9 2 4 ... #> ..$ PotEvap : num [1:4749] 0.4 0.3 0.3 0.2 0.3 0.3 0.5 0.3 0.3 0.2 ... #> ..$ Qupstream : num [1:4749, 1:2] -943 -1009 -1009 -989 -947 ... #> .. ..- attr(*, "dimnames")=List of 2 #> .. .. ..$ : NULL #> .. .. ..$ : chr [1:2] "P_GRANDEMAINE" "R_GRANDEMAINE" #> ..$ LengthHydro : Named num [1:2] 9.03 10.42 #> .. ..- attr(*, "names")= chr [1:2] "P_GRANDEMAINE" "R_GRANDEMAINE" #> ..$ BasinAreas : Named num [1:3] NA NA 132 #> .. ..- attr(*, "names")= chr [1:3] "P_GRANDEMAINE" "R_GRANDEMAINE" "GrandeMaine_StFulgent" #> ..$ id : chr "GrandeMaine_StFulgent" #> ..$ down : chr "Barrage_Bultiere" #> ..$ UpstreamNodes : chr [1:2] "P_GRANDEMAINE" "R_GRANDEMAINE" #> ..$ UpstreamIsModeled: Named logi [1:2] FALSE FALSE #> .. ..- attr(*, "names")= chr [1:2] "P_GRANDEMAINE" "R_GRANDEMAINE" #> ..$ UpstreamVarQ : Named chr [1:2] "Qsim_m3" "Qsim_m3" #> .. ..- attr(*, "names")= chr [1:2] "P_GRANDEMAINE" "R_GRANDEMAINE" #> ..$ FUN_MOD : chr "RunModel_GR4J" #> ..$ inUngaugedCluster: logi FALSE #> ..$ isReceiver : logi FALSE #> ..$ gaugedId : chr "GrandeMaine_StFulgent" #> ..$ hasUngaugedNodes : logi FALSE #> ..$ model :List of 4 #> .. ..$ indexParamUngauged: num [1:5] 1 2 3 4 5 #> .. ..$ hasX4 : logi TRUE #> .. ..$ iX4 : num 5 #> .. ..$ IsHyst : logi FALSE #> ..$ hasDiversion : logi FALSE #> ..$ isReservoir : logi FALSE #> ..- attr(*, "class")= chr [1:5] "RunModel_GR4J" "InputsModel" "daily" "GR" ... #> $ PetiteMaine_StGeorges:List of 19 #> ..$ DatesR : POSIXlt[1:4749], format: "2008-01-01" "2008-01-02" ... #> ..$ Precip : num [1:4749] 0 0.2 5 11.6 9.7 16.7 0.4 2 2.1 4.4 ... #> ..$ PotEvap : num [1:4749] 0.4 0.3 0.3 0.2 0.3 0.3 0.4 0.3 0.3 0.2 ... #> ..$ Qupstream : num [1:4749, 1:2] -1771 -1891 -1891 -1859 -1424 ... #> .. ..- attr(*, "dimnames")=List of 2 #> .. .. ..$ : NULL #> .. .. ..$ : chr [1:2] "P_PETITEMAINE" "R_PETITEMAINE" #> ..$ LengthHydro : Named num [1:2] 15.5 15.5 #> .. ..- attr(*, "names")= chr [1:2] "P_PETITEMAINE" "R_PETITEMAINE" #> ..$ BasinAreas : Named num [1:3] NA NA 192 #> .. ..- attr(*, "names")= chr [1:3] "P_PETITEMAINE" "R_PETITEMAINE" "PetiteMaine_StGeorges" #> ..$ id : chr "PetiteMaine_StGeorges" #> ..$ down : chr "Maine_Remouille" #> ..$ UpstreamNodes : chr [1:2] "P_PETITEMAINE" "R_PETITEMAINE" #> ..$ UpstreamIsModeled: Named logi [1:2] FALSE FALSE #> .. ..- attr(*, "names")= chr [1:2] "P_PETITEMAINE" "R_PETITEMAINE" #> ..$ UpstreamVarQ : Named chr [1:2] "Qsim_m3" "Qsim_m3" #> .. ..- attr(*, "names")= chr [1:2] "P_PETITEMAINE" "R_PETITEMAINE" #> ..$ FUN_MOD : chr "RunModel_GR4J" #> ..$ inUngaugedCluster: logi FALSE #> ..$ isReceiver : logi FALSE #> ..$ gaugedId : chr "PetiteMaine_StGeorges" #> ..$ hasUngaugedNodes : logi FALSE #> ..$ model :List of 4 #> .. ..$ indexParamUngauged: num [1:5] 1 2 3 4 5 #> .. ..$ hasX4 : logi TRUE #> .. ..$ iX4 : num 5 #> .. ..$ IsHyst : logi FALSE #> ..$ hasDiversion : logi FALSE #> ..$ isReservoir : logi FALSE #> ..- attr(*, "class")= chr [1:5] "RunModel_GR4J" "InputsModel" "daily" "GR" ... #> $ Barrage_Bultiere :List of 18 #> ..$ DatesR : POSIXlt[1:4749], format: "2008-01-01" "2008-01-02" ... #> ..$ Qupstream : num [1:4749, 1:2] 0 0 0 0 0 0 0 0 0 0 ... #> .. ..- attr(*, "dimnames")=List of 2 #> .. .. ..$ : NULL #> .. .. ..$ : chr [1:2] "GrandeMaine_StFulgent" "AEP_Bultiere" #> ..$ LengthHydro : Named num [1:2] 5.37 0 #> .. ..- attr(*, "names")= chr [1:2] "GrandeMaine_StFulgent" "AEP_Bultiere" #> ..$ BasinAreas : Named num [1:3] 132 NA NA #> .. ..- attr(*, "names")= chr [1:3] "GrandeMaine_StFulgent" "AEP_Bultiere" "Barrage_Bultiere" #> ..$ id : chr "Barrage_Bultiere" #> ..$ down : chr "Maine_Remouille" #> ..$ UpstreamNodes : chr [1:2] "GrandeMaine_StFulgent" "AEP_Bultiere" #> ..$ UpstreamIsModeled: Named logi [1:2] TRUE FALSE #> .. ..- attr(*, "names")= chr [1:2] "GrandeMaine_StFulgent" "AEP_Bultiere" #> ..$ UpstreamVarQ : Named chr [1:2] "Qsim_m3" "Qsim_m3" #> .. ..- attr(*, "names")= chr [1:2] "GrandeMaine_StFulgent" "AEP_Bultiere" #> ..$ FUN_MOD : chr "RunModel_Lag" #> ..$ inUngaugedCluster: logi FALSE #> ..$ isReceiver : logi FALSE #> ..$ gaugedId : chr "Barrage_Bultiere" #> ..$ hasUngaugedNodes : logi FALSE #> ..$ model :List of 4 #> .. ..$ indexParamUngauged: num [1:2] 1 2 #> .. ..$ hasX4 : logi FALSE #> .. ..$ iX4 : num 5 #> .. ..$ IsHyst : logi FALSE #> ..$ hasDiversion : logi FALSE #> ..$ isReservoir : logi TRUE #> ..$ Qrelease : Named num [1:4749] 13824 13824 13824 13824 13824 ... #> .. ..- attr(*, "names")= chr [1:4749] "732" "733" "734" "735" ... #> ..- attr(*, "class")= chr [1:5] "RunModel_Reservoir" "InputsModel" "daily" "SD" ... #> $ Maine_Remouille :List of 19 #> ..$ DatesR : POSIXlt[1:4749], format: "2008-01-01" "2008-01-02" ... #> ..$ Precip : num [1:4749] 0 0.121 4.842 11.442 9.305 ... #> ..$ PotEvap : num [1:4749] 0.4 0.3 0.3 0.2 0.3 ... #> ..$ Qupstream : num [1:4749, 1:4] 0 0 0 0 0 0 0 0 0 0 ... #> .. ..- attr(*, "dimnames")=List of 2 #> .. .. ..$ : NULL #> .. .. ..$ : chr [1:4] "PetiteMaine_StGeorges" "Barrage_Bultiere" "P_MAINE" "R_MAINE" #> ..$ LengthHydro : Named num [1:4] 23.1 28 19.1 21.8 #> .. ..- attr(*, "names")= chr [1:4] "PetiteMaine_StGeorges" "Barrage_Bultiere" "P_MAINE" "R_MAINE" #> ..$ BasinAreas : Named num [1:5] 192 132 NA NA 271 #> .. ..- attr(*, "names")= chr [1:5] "PetiteMaine_StGeorges" "Barrage_Bultiere" "P_MAINE" "R_MAINE" ... #> ..$ id : chr "Maine_Remouille" #> ..$ down : chr NA #> ..$ UpstreamNodes : chr [1:4] "PetiteMaine_StGeorges" "Barrage_Bultiere" "P_MAINE" "R_MAINE" #> ..$ UpstreamIsModeled: Named logi [1:4] TRUE TRUE FALSE FALSE #> .. ..- attr(*, "names")= chr [1:4] "PetiteMaine_StGeorges" "Barrage_Bultiere" "P_MAINE" "R_MAINE" #> ..$ UpstreamVarQ : Named chr [1:4] "Qsim_m3" "Qsim_m3" "Qsim_m3" "Qsim_m3" #> .. ..- attr(*, "names")= chr [1:4] "PetiteMaine_StGeorges" "Barrage_Bultiere" "P_MAINE" "R_MAINE" #> ..$ FUN_MOD : chr "RunModel_GR4J" #> ..$ inUngaugedCluster: logi FALSE #> ..$ isReceiver : logi FALSE #> ..$ gaugedId : chr "Maine_Remouille" #> ..$ hasUngaugedNodes : logi FALSE #> ..$ model :List of 4 #> .. ..$ indexParamUngauged: num [1:5] 1 2 3 4 5 #> .. ..$ hasX4 : logi TRUE #> .. ..$ iX4 : num 5 #> .. ..$ IsHyst : logi FALSE #> ..$ hasDiversion : logi FALSE #> ..$ isReservoir : logi FALSE #> ..- attr(*, "class")= chr [1:5] "RunModel_GR4J" "InputsModel" "daily" "GR" ... #> - attr(*, "class")= chr [1:2] "GRiwrmInputsModel" "list" #> - attr(*, "GRiwrm")= GRiwrm [11 × 6] (S3: GRiwrm/tbl_df/tbl/data.frame) #> ..$ id : chr [1:11] "GrandeMaine_StFulgent" "PetiteMaine_StGeorges" "Barrage_Bultiere" "Maine_Remouille" ... #> ..$ down : chr [1:11] "Barrage_Bultiere" "Maine_Remouille" "Maine_Remouille" NA ... #> ..$ length: num [1:11] 5.37 23.07 28.04 NA 19.12 ... #> ..$ area : num [1:11] 132 192 NA 595 NA ... #> ..$ model : chr [1:11] "RunModel_GR4J" "RunModel_GR4J" "RunModel_Reservoir" "RunModel_GR4J" ... #> ..$ donor : chr [1:11] "GrandeMaine_StFulgent" "PetiteMaine_StGeorges" "Barrage_Bultiere" "Maine_Remouille" ... #> - attr(*, "TimeStep")= num 86400 ``` --- # Paramètres de simulation Nous utilisons les paramètres calés précédemment mais il faut ajouter les paramètres du réservoir: - `Vmax`: la capacité de stockage du réservoir (m³) - `celerity`: la vitesse de propagation des débits amont vers le réservoir ## Exercice *Compléter la liste des paramètres avec ceux du réservoir. Pour la célérité, utiliser celle calée sur le bassin versant intermédiaire aval.* -- ``` r str(Params) #> List of 3 #> $ GrandeMaine_StFulgent: num [1:5] 0.676 231.109 -0.745 15.29 2.094 #> $ PetiteMaine_StGeorges: num [1:5] 0.93 198.835 -0.892 10.753 2.193 #> $ Maine_Remouille : num [1:5] 5.454 187.15 -0.528 14.447 2.329 ``` ``` r Params$Barrage_Bultiere <- c(Vmax = 5E6, celerity = Params$Maine_Remouille[1]) ``` --- # Options de simulation du réservoir Le lancement de `CreateRunOptions` est le même que précédemment mais l'objet `InputsModel` a changé... ## Exercice *Définissons les options de simulation* -- ``` r RO2 <- CreateRunOptions( IM2, IndPeriod_WarmUp = I_WarmUp, IndPeriod_Run = I_Run ) ``` Les options de simulations contiennent les états initiaux du model. Le réservoir a pour état initial le volume contenu au premier pas de temps de simulation. Il est fixé par défaut à la moitié de la capacité du réservoir mais peut être redéfini (Voir `?RunModel_Reservoir`). --- # Simulation du réservoir ## Exercice *Lancer la simulation* -- ``` r OM2 <- RunModel(IM2, RunOptions = RO2, Param = Params) #> RunModel.GRiwrmInputsModel: Processing sub-basin GrandeMaine_StFulgent... #> Warning in RunModel_Lag(InputsModel, RunOptions, Param[1], OutputsModel): 174 #> time steps with negative flow, set to zero. #> RunModel.GRiwrmInputsModel: Processing sub-basin PetiteMaine_StGeorges... #> Warning in RunModel_Lag(InputsModel, RunOptions, Param[1], OutputsModel): 459 #> time steps with negative flow, set to zero. #> RunModel.GRiwrmInputsModel: Processing sub-basin Barrage_Bultiere... #> RunModel.GRiwrmInputsModel: Processing sub-basin Maine_Remouille... #> Warning in RunModel_Lag(InputsModel, RunOptions, Param[1], OutputsModel): 53 #> time steps with negative flow, set to zero. ``` --- # Simulation du réservoir ## Exercice *Afficher le graphique des débits et volumes simulés sur le barrage* ```r ?plot.OutputsModelReservoir ``` -- ``` r plot(OM2$Barrage_Bultiere) ``` <img src="index_files/figure-html/plot.OutputModelReservoir-1.png" width="65%" style="display: block; margin: auto;" /> --- class: inverse, middle, center # Régulation du barrage de la Bultière Et si on régulait le barrage pour avoir un débit minimal de 400 L/s à Rémouillé ? --- # Création du superviseur Le superviseur est un ordonnanceur de la simulation qui va : - stopper la simulation tous les *n* pas de temps - lancer un ou plusieurs contrôleurs qui vont: - récupérer des données simulées (mesures) aux pas de temps précédents - calculer des commandes (ex: prélèvement, lâcher de réservoir...) à partir d'une fonction dite de logique de contrôle - appliquer les commandes calculées pour les prochains pas de temps de simulation ```r ?CreateSupervisor ``` ## Exercice *Créer un superviseur avec un pas de temps de décision de 1 jour* -- ``` r sv <- CreateSupervisor(IM2) ``` --- # La fonction de logique de contrôle Elle prend en entrée, les mesures (`Y`) qui seront effectuées durant la simulation. Dans notre cas, la mesure principale correspondra au débit simulé à Rémouillé. Dans le cas où un soutien d'étiage est nécessaire, il faut savoir quelle part du débit de Rémouillé est déjà soutenu par le barrage. Le débit simulé au niveau du barrage sera nécessaire pour naturaliser le débit de Rémouillé et ainsi savoir combien ajouter pour le soutien. Le débit réservé est de 160 L/s et il faut soutenir l'étiage à Rémouillé de façon à avoir un débit minimum de 400 L/s. Le débit à lâcher peut être calculé comme suit. ``` r #' Fonction de logique pour un DOE de 400 L/s à Rémouillé #' @param Y [numeric] avec: #' - Y[1]: le débit simulé à Rémouillé (m3/j) #' - Y[2]: le débit simulé au barrage (m3/j) #' @return Le débit à lâcher au barrage (m3/j) logicFun <- function(Y) { Qnat <- Y[1] - Y[2] # Débit naturel à Rémouillé sans soutien U <- 0 # Commençons avec un soutien nul # Le soutien pour 400 L/s if (Qnat < 0.4 * 86400) U <- 0.4 * 86400 - Qnat # Respectons le débit réservé if (U < 0.16 * 86400) U <- 0.16 * 86400 return(U) } ``` --- # Création du contrôleur ```r ?CreateController ``` Le contrôleur est ajouté au superviseur et doit définir: - Son nom - La localisation des mesures `Y` - La localisation des commandes `U` - La fonction de logique de contrôle utilisée `FUN` ## Exercice *Créer le contrôleur de soutien d'étiage* -- ``` r CreateController(sv, "DOE", Y = c("Maine_Remouille", "Barrage_Bultiere"), U = "Barrage_Bultiere", FUN = logicFun) #> The controller 'DOE' has been added to the supervisor ``` --- ## Simulation avec superviseur ```r ?RunModel.Supervisor ``` ## Exercice *Lancer la simulation supervisée :)* -- ``` r OM3 <- RunModel(sv, RunOptions = RO2, Param = Params) #> Processing: 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% ``` --- # Visualisation du résultat au barrage ``` r plot(OM3$Barrage_Bultiere) ``` <img src="index_files/figure-html/plotBarrageSupervised-1.png" width="100%" style="display: block; margin: auto;" /> --- # Visualisation du résultat à Rémouillé *Traçons la chronique des débits simulés et comparons-la aux débits observés* ``` r plot(OM3$Maine_Remouille, Qobs = Q[I_Run, "Maine_Remouille"], #OM$Maine_Remouille$Qsim, BasinArea = g_bultiere$area[g_bultiere$id == "Maine_Remouille"], which = "Flows", log_scale = TRUE) ``` <img src="index_files/figure-html/plotFlowsRemouille-1.png" width="100%" style="display: block; margin: auto;" /> --- class: inverse, middle, center # Merci pour votre attention **airGRiwrm: 'airGR' Integrated Water Resource Management** **Documentation**: `https://airgriwrm.g-eau.fr` **Contact**: David Dorchies (david.dorchies@inrae.fr)