Skip to content
GitLab
Select Git revision
  • master default
Find file BlameHistoryPermalink
script.R 3.66 KiB
1
2 
###### A MODIFIER ######
3
4
5 
# Path to the data
6 
computer_data_path =
7 
    "/home/louis/Documents/bouleau/INRAE/CDD_stationnarite/data"
8 
    # "C:\\Users\\louis.heraut\\Documents\\CDD_stationnarite\\data"
9
10 
# Work path
11 
computer_work_path =
12 
    "/home/louis/Documents/bouleau/INRAE/CDD_stationnarite/ASH"
13 
    # "C:\\Users\\louis.heraut\\Documents\\CDD_stationnarite\\ASH"
14
15
16 
### BANQUE HYDRO ###
17 
# Path to the directory where BH data is stored
18 
BHfiledir =
19 
    # "test"
20 
    "BanqueHydro_Export2021"
21
22 
## Manual selection ##
23 
# Name of the file that will be analysed from the BH directory
24 
BHfilename =
25 
    FALSE
26 
    # c("H5920011_HYDRO_QJM.txt", "K4470010_HYDRO_QJM.txt")
27 
    # "all"
28
29 
## Or list selection ##
30 
# Path to the list file of BH data that will be analysed
31 
BHlistdir =
32 
    # FALSE
33 
    ""
34
35 
BHlistname =
36 
    # FALSE
37 
    "Liste-station_RRSE.docx"
38
39
40 
### NIVALE ###
41 
# Path to the directory where NV data is stored
42 
NVfiledir =
43 
    "France207"
44
45 
# Name of the file that will be analysed from the NV directory
46 
NVfilename =
47 
    "all"
48
49 
# Path to the list file of information about station that will be analysed
50 
NVlistdir =
51 
    ""
52
53 
NVlistname =
54 
    "liste_bv_principaux_global.txt"
55
56
57 
########################
58
59
60 
# Set working directory
61 
setwd(computer_work_path)
62
63 
# Sourcing R file
64 
source('processing/extractBH.R')
65 
source('processing/extractNV.R')
66 
source('processing/analyse.R')
67 
source('plotting/panel.R')
68
69 
# Usefull library
70
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138
# 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)) # BANQUE HYDRO # # Get only the selected station from a list station file if (is.character(BHlistdir) & is.character(BHlistname) & is.character(BHfiledir)){ df_selec = get_selection(computer_data_path, BHlistdir, BHlistname, cnames=c('code', 'station', 'BV_km2', 'axe_principal_concerne', 'longueur_serie', 'commentaires', 'choix'), c_num=c('BV_km2', 'longueur_serie')) BHfilename = df_selec[df_selec$ok,]$filename } # Extract information about selected stations df_info_BH = extractBH_info(computer_data_path, BHfiledir, BHfilename) # Extract data about selected stations df_data_BH = extractBH_data(computer_data_path, BHfiledir, BHfilename) # NIVALE # # Extract information about selected stations df_info_NV = extractNVlist_info(computer_data_path, NVfiledir, NVlistdir, NVlistname) ### /!\ missing station info # Extract data about selected stations df_data_NV = extractNV_data(computer_data_path, NVfiledir, NVfilename) common = levels(factor(df_info_NV[df_info_NV$code %in% df_info_BH$code,]$code)) NVadd = levels(factor(df_info_NV[!(df_info_NV$code %in% df_info_BH$code),]$code)) df_info_NVadd = df_info_NV[df_info_NV$code %in% NVadd,] df_info = full_join(df_info_BH, df_info_NVadd, by=c("code", "nom", "L93X", "L93Y", "surface_km2", "file_path")) df_data_NVadd = df_data_NV[df_data_NV$code %in% NVadd,] df_data = full_join(df_data_BH, df_data_NVadd, by=c("Date", "Qm3s", "code")) # Plot time panel of debit by stations # panel(df_data, df_info, figdir, filedir) # panel(df_data, df_info, figdir, filedir, is_sqrt=TRUE) # Compute gap parameters for stations # df_lac = get_lacune(df_data, df_info)