-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathresponse_time.R
106 lines (103 loc) · 2.66 KB
/
response_time.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
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
get_sensitivities <- function() {
Permafrost <- list(
states = list(
list(ett = c(1e1, 1e2), cy = 74), # ett = response time interval, here 10-100 years
list(ett = c(1e1, 1e2), cy = 57), # cy = southern limit, here 57°N
list(ett = c(1e1, 1e2), cy = 54),
list(ett = c(1e1, 1e3), cy = 51)
)
)
# ==========================================================================
peatlands <- list(
states = list(
list(ett = c(1e1, 1e2), cy = 57),
list(ett = c(1e2, 1e3), cy = 56),
list(ett = c(1e1, 1e2), cy = 54)
)
)
Lakes <- list(
states = list(
list(ett = c(1e0, 1e1), cy = 83),
list(ett = c(1e0, 1e1), cy = 74),
list(ett = c(1e1, 1e3), cy = 58),
list(ett = c(1e1, 1e2), cy = 52)
)
)
# ==========================================================================
Snowpack <- list(
states = list(
list(ett = c(1e1, 1e2), cy = 74),
list(ett = c(1e1, 1e2), cy = 61),
list(ett = c(1e2, 1e3), cy = 58)
)
)
# ==========================================================================
Vegetation <- list(
states = list(
list(ett = c(1e1, 1e2), cy = 74),
list(ett = c(1e1, 1e2), cy = 61),
list(ett = c(1e2, 1e3), cy = 58),
list(ett = c(1e2, 1e3), cy = 51),
list(ett = c(1e2, 1e3), cy = 49)
)
)
# ==========================================================================
vertebrates <- list(
states = list(
list(ett = c(1e0, 1e2), cy = 83),
list(ett = c(1e1, 1e2), cy = 61),
list(ett = c(1e2, 1e3), cy = 58)
)
)
coords <- list(
Permafrost = Permafrost,
peatlands = peatlands,
Lakes = Lakes,
Snowpack = Snowpack,
Vegetation = Vegetation,
vertebrates = vertebrates
)
ls_nm <- list(
per = c(
"Polar desert continuous",
"Polar tundra continuous",
"Extensive discontinuous",
"Sporadic discontinuous",
"Relict"
),
pea = c(
"Polygonal peatlands",
"Complex of tundra peatlands",
"Mosaic palsa and plateau",
"Non-permafrost peatlands"
),
lak = c(
"Perrenial ice cover",
"Intermittent ice-free",
"Season. ice cover/clear",
"Season. ice cover/low colour",
"Season. ice cover/high colour"
),
sno = c(
"Polar desert",
"Polar tundra",
"Shrub tundra",
"Boreal forest"
),
veg = c(
"Polar desert",
"Polar tundra",
"Shrub tundra",
"Open-crown coniferous",
"Closed-crown coniferous",
"Mixed wood boreal"
),
ver = c(
"High-Arctic",
"Arctic",
"Low arctic",
"Boreal"
)
)
list(coords, ls_nm)
}