-
Notifications
You must be signed in to change notification settings - Fork 0
/
simulation.R
89 lines (78 loc) · 3.74 KB
/
simulation.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
### hisafe-calibration
### SIMULATION
### Author: Kevin J. Wolz
A2.WEATHER <- "./raw_data/restinclieres_A2-1994-2018.wth"
A3.WEATHER <- "./raw_data/restinclieres_A3-1994-2018.wth"
A4.WEATHER <- "./raw_data/restinclieres_A4-1994-2018.wth"
CASTRIES.WEATHER <- "./raw_data/castries.wth"
tree.params <- list()
winner.tree.params <- readr::read_csv("/Users/kevinwolz/Desktop/RESEARCH/ACTIVE_PROJECTS/HI-SAFE/hisafe-GA/analysis/WINNER_GA_PARAMS.csv",
col_types = readr::cols()) %>% as.list()
tree.params <- c(tree.params, winner.tree.params)
#crop.params <- list(tminremp = 0, tmaxremp = 45)
crop.params <- list()
## DEFINE
A2.hip <- define_hisafe(path = PATH,
profiles = PROFILES,
template = "restinclieres_agroforestry_A2",
SimulationName = "Restinclieres-A2",
bulk.pass = c(tree.params, crop.params),
weatherFile = A2.WEATHER)
A3.hip <- define_hisafe(path = PATH,
profiles = PROFILES,
template = "restinclieres_agroforestry_A3",
SimulationName = "Restinclieres-A3",
bulk.pass = c(tree.params, crop.params),
weatherFile = A3.WEATHER)
A4.hip <- define_hisafe(path = PATH,
profiles = PROFILES,
template = "restinclieres_forestry_A4",
SimulationName = "Restinclieres-A4",
bulk.pass = tree.params,
weatherFile = A4.WEATHER)
A2.CC.hip <- define_hisafe(path = PATH,
profiles = PROFILES[!str_detect(PROFILES, "trees")],
template = "restinclieres_monocrop_A2",
SimulationName = "Monocrop-A2",
bulk.pass = crop.params,
weatherFile = A2.WEATHER)
A3.CC.hip <- define_hisafe(path = PATH,
profiles = PROFILES[!str_detect(PROFILES, "trees")],
template = "restinclieres_monocrop_A3",
SimulationName = "Monocrop-A3",
bulk.pass = crop.params,
weatherFile = A3.WEATHER)
castries.hip <- define_hisafe(path = PATH,
profiles = PROFILES,
template = "castries_agroforestry",
SimulationName = "Castries",
mainCropSpecies = "fescue.plt",
interCropSpecies = "fescue.plt",
mainCropItk = "fescue.tec",
interCropItk = "fescue.tec",
bulk.pass = tree.params,
weatherFile = CASTRIES.WEATHER)
## BUILD
if(BUILD.SIMU) {
build_hisafe(A2.hip)
build_hisafe(A3.hip)
build_hisafe(A4.hip)
build_hisafe(A2.CC.hip)
build_hisafe(A3.CC.hip)
build_hisafe(castries.hip)
}
if(CLUSTER) {
build_cluster_script(simu.names = ALL.SIMULATIONS,
hip = NULL,
script.path = PATH,
cluster.path = paste0("/lustre/lecomtei/calibration/", NAME, "/"),
model.version = MODEL.VERSION,
email = "[email protected]")
}
## RUN
if(RUN.SIMU) {
run_hisafe(path = PATH,
parallel = TRUE,
num.cores = 4,
capsis.path = "/Applications/Capsis/")
}