Module vegmod


VARIABLES

optmodeinteger public 
Indicator of optimisation mode
125
timereal*8 public 
Seconds of hour
126
errorreal*8 public 
Cumulative error in water balance
127
finishinteger public 
flag to finish all loops
128
output_matreal*8(:,:)public 
?  
130
p_areal*8 publicparameter
Ratio of diffusivities of water vapour to CO2 in air
132
p_pireal*8 publicparameter
Pi-constant
133
p_mpbarreal*8 publicparameter
Conversion factor from MPa to bar
134
p_Ereal*8 publicparameter
Eurler's number
135
p_R_real*8 publicparameter
Molar gas konstant
136
l_E_real*8 publicparameter
Latent heat of vaporization (MJ/kg)
137
srad2par_hreal*8 publicparameter
Conversion from srad to par hourly (mol/MJ)
138
srad2par_dreal*8 publicparameter
Conversion from srad to par daily (mol/MJ)
139
rho_watreal*8 publicparameter
Density of water (kg/m3)
140
nyearinteger public 
Year
142
ndayinteger public 
Day since start of run
143
nhourinteger public 
Hour of day
144
th_integer public 
Hour since start of run
145
c_testdayinteger public 
Number of days for initial check if netass>0
146
fyearinteger(:)public 
Year for each day
148
fmonthinteger(:)public 
Month for each day
149
fdayinteger(:)public 
Day of month
150
dayyearinteger(:)public 
Day of year
151
tair_hreal*8(:)public 
Hourly air temperature (K)
155
tairmin_dreal*8(:)public 
Daily minimum temperature (K)
156
tairmax_dreal*8(:)public 
Daily maximum temperature (K)
157
topt_real*8 public 
Optimal temperature in temperature response curve
159
press_dreal*8(:)public 
Daily air pressure (Pa)
161
par_hreal*8(:)public 
Hourly photosynthetically active radiation (mol/m2/s)
163
par_dreal*8(:)public 
Daily photosynthetically active radiation
164
par_yreal*8 public 
Annual photosynthetically active radiation
165
srad_dreal*8(:)public 
Daily shortwave radiation
167
srad_yreal*8 public 
Annual shortwave radiation
168
ca_hreal*8(:)public 
Hourly atmospheric CO2 mole fraction
170
ca_dreal*8(:)public 
Daily atmospheric CO2 mole fraction
171
vp_dreal*8(:)public 
Daily absolute vapour pressure (Pa)
173
vd_hreal*8(:)public 
Hourly atmospheric vapour deficit (VPD/air pressure)
175
vd_dreal*8 public 
Mean daily atmospheric vapour deficit
176
vd_yreal*8 public 
Mean annual atmospheric vapour deficit
177
rain_hreal*8(:)public 
Hourly rainfall rate (m/s)
179
rain_dreal*8(:)public 
Daily rainfall
180
rain_yreal*8 public 
Annual rainfall
181
c_hhydrstreal*8(:)public 
Hydrostatic head in each layer relative to soil surface
185
gammastarreal*8 public 
CO2 compensation point
187
wsnewreal*8 public 
Total soil water store at next time step
189
wsoldreal*8 public 
Previous total soil water storage
190
o_caireal*8 public 
Projected cover perennial vegetation (0-1)
192
pcg_dreal*8(3)public 
Projected cover seasonal vegetation (pcg_d(2) is actual value)
193
c_pcgminreal*8 public 
Minimum grass pc; initial point for growth
194
o_wstexpreal*8 public 
Exponent for calculating lambdat_d
198
o_wsgexpreal*8 public 
Exponent for calculating lambdag
199
o_lambdatfreal*8 public 
Factor for calculating lambdat_d
200
o_lambdagfreal*8 public 
Factor for calculating lambdag_d
201
lambdat_dreal*8 public 
Target dE/dA for calculating gstomt
202
lambdag_dreal*8 public 
Target dE/dA for calculating gstomg
203
gstomtreal*8 public 
Tree stomatal conductance
204
gstomgreal*8(3,3,3)public 
Grass stomatal conductance
205
rlt_hreal*8(3,3)public 
Tree leaf respiration for different values of Jmax (rlt_h(2) is actual value)
207
rlt_dreal*8 public 
Daily tree leaf respiration
208
rlt_yreal*8 public 
Annual tree leaf respiration
209
rlg_hreal*8(3,3,3)public 
Grass leaf respiration
210
rlg_dreal*8 public 
Daily grass leaf respiration
211
rlg_yreal*8 public 
Annual grass leaf respiration
212
transptreal*8 public 
Tree transpiration rate
214
transpgreal*8(3,3,3)public 
Grass transpiration rate (mol/m2/s)
215
q_tct_dreal*8(3)public 
Tree foliage turnover costs
217
tct_yreal*8 public 
Annual tree foliage turnover costs
218
tcg_dreal*8(3, 3)public 
Grass foliage turnover costs
219
tcg_yreal*8 public 
Annual grass foliage turnover costs
220
jacttreal*8(3,3)public 
Electron transport rates for different values of Jmax (jactt(2) is actual value)
222
jactgreal*8(3,3,3)public 
Grass electron transport rate
223
jmaxt_hreal*8(3)public 
Tree photosynthetic electron transport capacity
225
jmaxg_hreal*8(3)public 
Grass electron transport capacity
226
jmax25t_dreal*8(3)public 
Tree photosynthetic electron transport capacity at 25oC
228
jmax25g_dreal*8(3)public 
Grass photosynthetic electron transport capacity at 25oC
229
lai_ltreal*8(3)public 
Local leaf area index trees
231
lai_lgreal*8(3)public 
Local leaf area index grasses
232
asst_hreal*8(3,3)public 
Tree hourly assimilation rate for different values of Jmax (asst_h(2) is actual value)
236
asst_dreal*8(3,3)public 
Daily tree assimilation
237
asst_yreal*8 public 
Annual tree assimilation
238
assg_hreal*8(3,3,3)public 
Hourly grass assimilation
239
assg_dreal*8(3,3,3)public 
Daily grass assimilation
240
assg_yreal*8 public 
Annual grass assimilation
241
q_cpcct_dreal*8 public 
Tree water transport costs as a function of projected cover and rooting depth (mol/m2/s)
243
cpcct_yreal*8 public 
Annual tree water transport costs
244
cpccg_dreal*8(3)public 
Grass water transport costs
245
cpccg_yreal*8 public 
Annual grass water transport costs
246
etmt__real*8 public 
Transpiration rate (m/s)
248
etmt_hreal*8 public 
Hourly transpiration
249
etmt_dreal*8 public 
Daily transpiration rate
250
etmt_yreal*8 public 
Annual tree transpiration
251
etmg__real*8(3,3,3)public 
Grass transpiration rate (m/s)
252
etmg_hreal*8 public 
Hourly grass transpiration
253
etmg_dreal*8 public 
Daily grass transpiration
254
etmg_yreal*8 public 
Annual grass transpiration
255
etm_yreal*8 public 
Annual total transpiration
256
mqt_real*8 public 
Tree water content
258
mqtnewreal*8 public 
Tree water content in next time step
259
mqtoldreal*8 public 
Previous tree water content
260
dmqtreal*8 public 
Rate of change in tree water content
261
q_mqxreal*8 public 
Tree maximum water content per ground area
262
mqsst_real*8 public 
Tree water content at steady state
263
mqsstminreal*8 public 
Tree water content at turgor loss point
264
q_mdreal*8 public 
Tree dry mass per unit ground area
266
o_mdstorereal*8 public 
Wood water storage parameter of trees
267
o_rtdepthreal*8 public 
Tree rooting depth (m)
271
o_rgdepthreal*8 public 
Grass rooting depth
272
pos_sltinteger public 
Lowest soil layer containing tree roots
274
pos_slginteger public 
Lowest soil layer containing grass roots
275
pos_ultinteger public 
Lowest soil layer containing tree roots within unsaturated zone
276
pos_ulginteger public 
Lowest soil layer containing grass roots within unsaturated zone
277
changefreal*8 public 
Change factor for adjusting root surface area
279
rsurft_real*8(:)public 
Root surface area of trees in each layer
281
rsurftnewreal*8(:)public 
Adjusted root surface area of trees in each layer for next day
282
rsurfg_real*8(:)public 
Root surface area of grasses in each layer
283
rsurfgnewreal*8(:)public 
Adjusted root surface area of grasses in each layer for next day
284
rootlimreal*8(3,3,3)public 
Indicator whether root surface are was limiting root water uptake
286
rsoilreal*8(:)public 
Resistance to water flow towards roots in each soil layer
288
refftreal*8(:)public 
Relative root water uptake efficiency for trees in each layer
290
reffgreal*8(:)public 
Relative root water uptake efficiency for grasses in each layer
291
posmnainteger(3)public 
Pointer to variable values that achieved maximum net assimilation
292
rrt_dreal*8 public 
Tree root respiration rate (mol/m2/s)
294
rrt_yreal*8 public 
Annual tree root respiration
295
rrg_dreal*8 public 
Grass root respiration
296
rrg_yreal*8 public 
Annual grass root respiration
297
prootmreal*8(:)public 
Root hydraulic head in each layer
299
sumruptkt_hreal*8 public 
Hourly total tree root water uptake
301
ruptkt__real*8(:)public 
Root water uptake rate perennial veg (m/s)
302
ruptkt_hreal*8(:)public 
Hourly root water uptake by trees in each layer
303
ruptkt_dreal*8(:)public 
Daily root water uptake by trees in each layer
304
ruptkg__real*8(:)public 
Root water uptake rate seasonal veg (m/s)
305
ruptkg_hreal*8(:)public 
Hourly root water uptake by grasses in each layer
306
ruptkg_dreal*8(:)public 
Daily root water uptake by grasses in each layer
307
perc_cov_vegreal*8(:)public 
Daily coverage of vegetation
308
i_alphareal*8 public 
Initial slope of electron transport curve
315
i_cpccfreal*8 public 
Water transport costs per m root depth and m^2 cover
316
i_tcfreal*8 public 
Turnover cost factor for foliage (tc=i_tcf*LAI)
317
i_maxyearinteger public 
Number of years to process
318
i_testyearinteger public 
Number of years after which to perform initial test of netass
319
i_hareal*8 public 
Temperature response parameter
320
i_hdreal*8 public 
Temperature response parameter
321
i_toptfreal*8 public 
Parameter to calculate adaptation of topt (range 0-1 for no to full adaptation)
322
i_toptstartreal*8 public 
Start parameter for topt to calculate jmax(temp in K)
323
i_rlratioreal*8 public 
Ratio of leaf respiration to photosynthetic capacity
324
i_latreal*8 public 
geogr. latitude
328
i_mdtfreal*8 public 
Total dry mass of living tissues of trees per unit pc (g/m^2)
337
i_mqxtfreal*8 public 
Total water storage capacity in living tissues of trees per unit md
338
i_rrootmreal*8 public 
Root water uptake resistivity in soil
339
i_rsurfminreal*8 public 
Minimum root area per m^3 to be maintained
340
i_rsurf_real*8 public 
Initial root surface area per m^3
341
i_rootradreal*8 public 
Average fine root radius
342
i_prootmgreal*8 public 
Constant root balance pressure of 1.5 MPa in grasses
343
i_growthmaxreal*8 public 
Parameter determining maximum daily growth increment of root surface area
344
i_incrcovgreal*8 public 
parameter determining maximum increment percentage of grass cover
345
i_incrjmaxreal*8 public 
parameter determining maximum increment percentage of jmax25
346
i_incrlaitreal*8 public 
parameter determining maximum increment percentage of lai trees
347
i_incrlaigreal*8 public 
parameter determining maximum increment percentage of lai grasses
348
i_extcoeffgreal*8 public 
extinction coefficient beer's law grasses
349
i_extcoefftreal*8 public 
extinction coefficient beer's law trees
350
i_trans_vegcovreal*8 public 
fraction of radiative energy reaching soil under full cover (0-1)
351
i_firstyearinteger public 
First year for the generation of hourly output in computation mode
353
i_lastyearinteger public 
Last year for the generation of hourly output in computation mode
354
i_write_hinteger public 
Flag to write out hourly input values after conversation from daily values
356
i_read_pcinteger public 
Flag to write out hourly input values after conversation from daily values
357
i_lai_functioninteger public 
Switch to use 1) linear or 2) exponential LAI estimate, as function of cover
358
i_no_veginteger public 
Flag to switch vegetation off (1=no vegetation)
359
c_epslnreal*8 public 
Soil porosity
364
c_maxhourinteger public 
Number of hours to process
366
c_maxdayinteger public 
Number of days to process
367