-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathAutumn_Weighted_Network_LIMfile.R
135 lines (98 loc) · 2.69 KB
/
Autumn_Weighted_Network_LIMfile.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
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
! Weighted Network
! Network LIM Declaration File
! Composed with autoLIMR
! Author: gemma
! Date: 2023-06-17
! Respiration included: Yes
! U included: Yes
! Living compartments: 3
! Non-living compartments: 1
! External compartments: 7
! Boundary flows: 10
! Internal flows: 6
! Abbreviations
! GPP = Gross Primary Production (Primary Producers only)
! Q = Consumption
! NPP = Net Primary Production (Primary Producers only)
! P = Production
! R = respiration
! U = Passive flows to non-living compartments/Unassimilated material
! AE = Assimilation Efficiency
! IM = Import flow
! EX = Export Flow
! NLNode = Non-living compartment
### COMPARTMENTS
Invert = 900.1
Plant = 600.0
Vert = 33.0
DetNLNode = 8000.0
### END COMPARTMENTS
### EXTERNALS
CO2
DetNLNodeImport
PlantImport
DetNLNodeExport
InvertExport
PlantExport
VertExport
### END EXTERNALS
### VARIABLES
! Consumption (Q) / Gross Primary Production (GPP) Variables
Invert_Q = Flowto(Invert)
Plant_GPP = Flowto(Plant) - Plant_IM
Vert_Q = Flowto(Vert)
! Production (P/NPP) Variables
Invert_P = Flowfrom(Invert) - Invert_R - Invert_U - Invert_EX
Plant_NPP = Flowfrom(Plant) - Plant_R - Plant_U - Plant_EX
Vert_P = Flowfrom(Vert) - Vert_R - Vert_U - Vert_EX
! Unused Energy/Material (U) Variables
Invert_U = Flowto(Invert) - Invert_P - Invert_R - Invert_EX
Plant_U = Flowto(Plant) - Plant_NPP - Plant_R - Plant_EX
Vert_U = Flowto(Vert) - Vert_P - Vert_R - Vert_EX
### END VARIABLES
### FLOWS
! GPP flows
Plant_GPP: CO2 -> Plant
! Respiration flows
Plant_R: Plant -> CO2
Invert_R: Invert -> CO2
Vert_R: Vert -> CO2
! Import flows
DetNLNode_IM: DetNLNodeImport -> DetNLNode
Plant_IM: PlantImport -> Plant
! Export flows
DetNLNode_EX: DetNLNode -> DetNLNodeExport
Plant_EX: Plant -> PlantExport
Invert_EX: Invert -> InvertExport
Vert_EX: Vert -> VertExport
! Adjacency Matrix flows
DetNLNode_Q_Invert: DetNLNode -> Invert
Invert_Q_Vert: Invert -> Vert
Invert_U_DetNLNode: Invert -> DetNLNode
Plant_Q_Invert: Plant -> Invert
Plant_U_DetNLNode: Plant -> DetNLNode
Vert_U_DetNLNode: Vert -> DetNLNode
### END FLOWS
### INEQUALITIES
! Network Data Input Inequalities
Plant_GPP > 850
Plant_GPP < 1700
Plant_NPP > 0.4*Plant_GPP
Invert_P > 0.2*Invert_Q
Vert_P > 2e-04
Plant_NPP < 0.9*Plant_GPP
Invert_P < 0.4*Invert_Q
Plant_R > 0.2*Plant_NPP
Invert_R > 0.8*Invert_P
Vert_R > 0.35
Plant_R < 0.5*Plant_NPP
Invert_R < 3*Invert_P
Vert_R < 0.35
DetNLNode_IM < 0.1*DetNLNode
Plant_IM < 0.1*Plant
DetNLNode_EX < 0.1*DetNLNode
Plant_EX < 0.1*Plant
Invert_EX < 0.1*Invert
Vert_EX < 0.1*Vert
! Adjacency Matrix Inequalities
### END INEQUALITIES