Skip to content
GitLab
Select Git revision
  • master default
Find file BlameHistoryPermalink
advanced_settings.R 5.48 KiB
1 
# \\\
2 
# Copyright 2021-2022 Louis Héraut*1,
3 
#                     Éric Sauquet*2,
4 
#                     Valentin Mansanarez
5 
#
6 
# *1   INRAE, France
7 
#      louis.heraut@inrae.fr
8 
# *2   INRAE, France
9 
#      eric.sauquet@inrae.fr
10 
#
11 
# This file is part of ash R toolbox.
12 
#
13 
# Ash R toolbox is free software: you can redistribute it and/or
14 
# modify it under the terms of the GNU General Public License as
15 
# published by the Free Software Foundation, either version 3 of the
16 
# License, or (at your option) any later version.
17 
#
18 
# Ash R toolbox is distributed in the hope that it will be useful, but
19 
# WITHOUT ANY WARRANTY; without even the implied warranty of
20 
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
21 
# General Public License for more details.
22 
#
23 
# You should have received a copy of the GNU General Public License
24 
# along with ash R toolbox.
25 
# If not, see <https://www.gnu.org/licenses/>.
26 
# ///
27 
#
28 
#
29 
# Rcode/advanced_settings.R
30 
#
31 
# Configuration file that regroups different advanced settings about
32 
# file structure, analyses and plotting.
33
34
35 
## 1. FILE STRUCTURE _________________________________________________
36 
### 1.1. Input directories ___________________________________________
37 
# Path to the data
38 
computer_data_path = file.path(computer_work_path, 'data')
39
40 
# Resources directory
41 
resources_path = file.path(computer_work_path, 'resources')
42 
if (!(file.exists(resources_path))) {
43 
  dir.create(resources_path)
44 
}
45 
print(paste('resources_path :', resources_path))
46
47 
# Logo filename
48 
logo_dir = 'logo'
49 
PRlogo_file = 'logo_Prefet_bassin.png'
50 
AEAGlogo_file = 'agence-de-leau-adour-garonne_logo.png'
51 
INRAElogo_file = 'Logo-INRAE_Transparent.png'
52 
FRlogo_file = 'Republique_Francaise_RVB.png'
53
54 
shp_dir = 'map'
55 
# Path to the shapefile for france contour from 'computer_data_path'
56 
fr_shpdir = file.path(shp_dir, 'france')
57 
fr_shpname = 'gadm36_FRA_0.shp'
58
59 
# Path to the shapefile for basin shape from 'computer_data_path'
60 
bs_shpdir = file.path(shp_dir, 'bassin')
61 
bs_shpname = 'BassinHydrographique.shp'
62
63 
# Path to the shapefile for sub-basin shape from 'computer_data_path'
64 
sbs_shpdir = file.path(shp_dir, 'sous_bassin')
65 
sbs_shpname = 'SousBassinHydrographique.shp'
66
67 
# Path to the shapefile for station basins shape from 'computer_data_path'
68 
cbs_shpdir = file.path(shp_dir, 'bassin_station')
69 
cbs_shpname = c('BV_4207_stations.shp', '3BVs_FRANCE_L2E_2018.shp')
70 
cbs_coord = c("L93", "L2")
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140
# Path to the shapefile for river shape from 'computer_data_path' rv_shpdir = file.path('map', 'river') rv_shpname = 'CoursEau_FXX.shp' ### 1.2. Output directories __________________________________________ # Result directory resdir = file.path(computer_work_path, 'results') if (!(file.exists(resdir))) { dir.create(resdir) } print(paste('resdir :', resdir)) # Figure directory figdir = file.path(computer_work_path, 'figures') if (!(file.exists(figdir))) { dir.create(figdir) } print(paste('figdir :', figdir)) ## 2. STATION SELECTION BY LIST ______________________________________ ### 2.1. Selection with '.docx' file _________________________________ # Path to a '.docx' list file of station that will be analysed DOCXlistdir = "" DOCXlistname = "" # "Liste-station_RRSE.docx" ### 2.2. Selection with '.txt' file _________________________________ # Path to the '.txt' list file of station that will be analysed # It can be generated with : # create_selection(computer_data_path, 'dirname', 'selection.txt') TXTlistdir = "" TXTlistname = "" # "selection.txt" ## 3. ANALYSIS PARAMETERS ____________________________________________ ### 3.1. Station analysis ____________________________________________ to_analyse = c( 'QA', 'QMNA', 'VCN10', 'tDEB', 'tCEN' ) ### 3.2. Climate analysis ____________________________________________ to_analyse_climate = c( 'PA', 'TA', 'ETPA' ) ### 3.3. Data modification ___________________________________________ # Local corrections of the data df_flag = data.frame( code=c('O3141010', 'O7635010', 'O7635010', 'O7635010', 'O7635010' ), Date=c('1974-07-04',
141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206
'1948-09-06', '1949-02-08', '1950-07-20', '1953-07-22' ), newValue=c(9.5, 4, 3, 1, 3) # /!\ Unit ) # Number of missing continuous days per year before removing the year dayLac_lim = 3 # Maximum continuously missing years before removing everything # before it yearNA_lim = 10 # Sampling span of the data sampleSpan = c('05-01', '11-30') ### 3.4. Saving option _______________________________________________ saving = c( # 'data', # 'meta', # 'analyse' ) fast_format = FALSE ### 3.5. Statistical option __________________________________________ # The risk of the Mann-Kendall trend detection test alpha = 0.1 allMonth = FALSE ## 4. PLOTTING OPTIONS _______________________________________________ # If the hydrological network needs to be plot show_river = FALSE # Graphical selection of period axis_xlim = NULL # c("1982-01-01", "1983-01-01") # What you want to be plotted for station analyses. For example if 'datasheet' is wrote, datasheet about each stations will be drawn. # All the option are : # 'datasheet' : datasheet of trend analyses for each stations # 'table' : summarizing table about trend analyses # 'map' : map about trend analyses to_plot_station = c( 'datasheet' # 'table', # 'map' ) map_to_plot = c( 'regime', 'trend', 'mean' )