-
Notifications
You must be signed in to change notification settings - Fork 0
/
Spring_Weighted_Network_LIMfile.R
129 lines (92 loc) · 2.54 KB
/
Spring_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
! 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 = 1000.0
Plant = 700.5
Vert = 44.0
DetNLNode = 9000.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 > 950
Plant_GPP < 1850
Plant_NPP > 0.5*Plant_GPP
Invert_P > 0.3*Invert_Q
Vert_P > 3e-04
Plant_NPP < 1*Plant_GPP
Invert_P < 0.5*Invert_Q
Plant_R > 0.3*Plant_NPP
Invert_R > 0.9*Invert_P
Vert_R > 0.65
Plant_R < 0.6*Plant_NPP
Invert_R < 4*Invert_P
Vert_R < 0.65
! Adjacency Matrix Inequalities
### END INEQUALITIES