optmode | integer |   | public |   | Indicator of optimisation mode
| 125 |
time | real*8 |   | public |   | Seconds of hour
| 126 |
error | real*8 |   | public |   | Cumulative error in water balance
| 127 |
finish | integer |   | public |   | flag to finish all loops
| 128 |
output_mat | real*8 | (:,:) | public |   | ?  
| 130 |
p_a | real*8 |   | public | parameter | Ratio of diffusivities of water vapour to CO2 in air
| 132 |
p_pi | real*8 |   | public | parameter | Pi-constant
| 133 |
p_mpbar | real*8 |   | public | parameter | Conversion factor from MPa to bar
| 134 |
p_E | real*8 |   | public | parameter | Eurler's number
| 135 |
p_R_ | real*8 |   | public | parameter | Molar gas konstant
| 136 |
l_E_ | real*8 |   | public | parameter | Latent heat of vaporization (MJ/kg)
| 137 |
srad2par_h | real*8 |   | public | parameter | Conversion from srad to par hourly (mol/MJ)
| 138 |
srad2par_d | real*8 |   | public | parameter | Conversion from srad to par daily (mol/MJ)
| 139 |
rho_wat | real*8 |   | public | parameter | Density of water (kg/m3)
| 140 |
nyear | integer |   | public |   | Year
| 142 |
nday | integer |   | public |   | Day since start of run
| 143 |
nhour | integer |   | public |   | Hour of day
| 144 |
th_ | integer |   | public |   | Hour since start of run
| 145 |
c_testday | integer |   | public |   | Number of days for initial check if netass>0
| 146 |
fyear | integer | (:) | public |   | Year for each day
| 148 |
fmonth | integer | (:) | public |   | Month for each day
| 149 |
fday | integer | (:) | public |   | Day of month
| 150 |
dayyear | integer | (:) | public |   | Day of year
| 151 |
tair_h | real*8 | (:) | public |   | Hourly air temperature (K)
| 155 |
tairmin_d | real*8 | (:) | public |   | Daily minimum temperature (K)
| 156 |
tairmax_d | real*8 | (:) | public |   | Daily maximum temperature (K)
| 157 |
topt_ | real*8 |   | public |   | Optimal temperature in temperature response curve
| 159 |
press_d | real*8 | (:) | public |   | Daily air pressure (Pa)
| 161 |
par_h | real*8 | (:) | public |   | Hourly photosynthetically active radiation (mol/m2/s)
| 163 |
par_d | real*8 | (:) | public |   | Daily photosynthetically active radiation
| 164 |
par_y | real*8 |   | public |   | Annual photosynthetically active radiation
| 165 |
srad_d | real*8 | (:) | public |   | Daily shortwave radiation
| 167 |
srad_y | real*8 |   | public |   | Annual shortwave radiation
| 168 |
ca_h | real*8 | (:) | public |   | Hourly atmospheric CO2 mole fraction
| 170 |
ca_d | real*8 | (:) | public |   | Daily atmospheric CO2 mole fraction
| 171 |
vp_d | real*8 | (:) | public |   | Daily absolute vapour pressure (Pa)
| 173 |
vd_h | real*8 | (:) | public |   | Hourly atmospheric vapour deficit (VPD/air pressure)
| 175 |
vd_d | real*8 |   | public |   | Mean daily atmospheric vapour deficit
| 176 |
vd_y | real*8 |   | public |   | Mean annual atmospheric vapour deficit
| 177 |
rain_h | real*8 | (:) | public |   | Hourly rainfall rate (m/s)
| 179 |
rain_d | real*8 | (:) | public |   | Daily rainfall
| 180 |
rain_y | real*8 |   | public |   | Annual rainfall
| 181 |
c_hhydrst | real*8 | (:) | public |   | Hydrostatic head in each layer relative to soil surface
| 185 |
gammastar | real*8 |   | public |   | CO2 compensation point
| 187 |
wsnew | real*8 |   | public |   | Total soil water store at next time step
| 189 |
wsold | real*8 |   | public |   | Previous total soil water storage
| 190 |
o_cai | real*8 |   | public |   | Projected cover perennial vegetation (0-1)
| 192 |
pcg_d | real*8 | (3) | public |   | Projected cover seasonal vegetation (pcg_d(2) is actual value)
| 193 |
c_pcgmin | real*8 |   | public |   | Minimum grass pc; initial point for growth
| 194 |
o_wstexp | real*8 |   | public |   | Exponent for calculating lambdat_d
| 198 |
o_wsgexp | real*8 |   | public |   | Exponent for calculating lambdag
| 199 |
o_lambdatf | real*8 |   | public |   | Factor for calculating lambdat_d
| 200 |
o_lambdagf | real*8 |   | public |   | Factor for calculating lambdag_d
| 201 |
lambdat_d | real*8 |   | public |   | Target dE/dA for calculating gstomt
| 202 |
lambdag_d | real*8 |   | public |   | Target dE/dA for calculating gstomg
| 203 |
gstomt | real*8 |   | public |   | Tree stomatal conductance
| 204 |
gstomg | real*8 | (3,3,3) | public |   | Grass stomatal conductance
| 205 |
rlt_h | real*8 | (3,3) | public |   | Tree leaf respiration for different values of Jmax (rlt_h(2) is actual value)
| 207 |
rlt_d | real*8 |   | public |   | Daily tree leaf respiration
| 208 |
rlt_y | real*8 |   | public |   | Annual tree leaf respiration
| 209 |
rlg_h | real*8 | (3,3,3) | public |   | Grass leaf respiration
| 210 |
rlg_d | real*8 |   | public |   | Daily grass leaf respiration
| 211 |
rlg_y | real*8 |   | public |   | Annual grass leaf respiration
| 212 |
transpt | real*8 |   | public |   | Tree transpiration rate
| 214 |
transpg | real*8 | (3,3,3) | public |   | Grass transpiration rate (mol/m2/s)
| 215 |
q_tct_d | real*8 | (3) | public |   | Tree foliage turnover costs
| 217 |
tct_y | real*8 |   | public |   | Annual tree foliage turnover costs
| 218 |
tcg_d | real*8 | (3, 3) | public |   | Grass foliage turnover costs
| 219 |
tcg_y | real*8 |   | public |   | Annual grass foliage turnover costs
| 220 |
jactt | real*8 | (3,3) | public |   | Electron transport rates for different values of Jmax (jactt(2) is actual value)
| 222 |
jactg | real*8 | (3,3,3) | public |   | Grass electron transport rate
| 223 |
jmaxt_h | real*8 | (3) | public |   | Tree photosynthetic electron transport capacity
| 225 |
jmaxg_h | real*8 | (3) | public |   | Grass electron transport capacity
| 226 |
jmax25t_d | real*8 | (3) | public |   | Tree photosynthetic electron transport capacity at 25oC
| 228 |
jmax25g_d | real*8 | (3) | public |   | Grass photosynthetic electron transport capacity at 25oC
| 229 |
lai_lt | real*8 | (3) | public |   | Local leaf area index trees
| 231 |
lai_lg | real*8 | (3) | public |   | Local leaf area index grasses
| 232 |
asst_h | real*8 | (3,3) | public |   | Tree hourly assimilation rate for different values of Jmax (asst_h(2) is actual value)
| 236 |
asst_d | real*8 | (3,3) | public |   | Daily tree assimilation
| 237 |
asst_y | real*8 |   | public |   | Annual tree assimilation
| 238 |
assg_h | real*8 | (3,3,3) | public |   | Hourly grass assimilation
| 239 |
assg_d | real*8 | (3,3,3) | public |   | Daily grass assimilation
| 240 |
assg_y | real*8 |   | public |   | Annual grass assimilation
| 241 |
q_cpcct_d | real*8 |   | public |   | Tree water transport costs as a function of projected cover and rooting depth (mol/m2/s)
| 243 |
cpcct_y | real*8 |   | public |   | Annual tree water transport costs
| 244 |
cpccg_d | real*8 | (3) | public |   | Grass water transport costs
| 245 |
cpccg_y | real*8 |   | public |   | Annual grass water transport costs
| 246 |
etmt__ | real*8 |   | public |   | Transpiration rate (m/s)
| 248 |
etmt_h | real*8 |   | public |   | Hourly transpiration
| 249 |
etmt_d | real*8 |   | public |   | Daily transpiration rate
| 250 |
etmt_y | real*8 |   | public |   | Annual tree transpiration
| 251 |
etmg__ | real*8 | (3,3,3) | public |   | Grass transpiration rate (m/s)
| 252 |
etmg_h | real*8 |   | public |   | Hourly grass transpiration
| 253 |
etmg_d | real*8 |   | public |   | Daily grass transpiration
| 254 |
etmg_y | real*8 |   | public |   | Annual grass transpiration
| 255 |
etm_y | real*8 |   | public |   | Annual total transpiration
| 256 |
mqt_ | real*8 |   | public |   | Tree water content
| 258 |
mqtnew | real*8 |   | public |   | Tree water content in next time step
| 259 |
mqtold | real*8 |   | public |   | Previous tree water content
| 260 |
dmqt | real*8 |   | public |   | Rate of change in tree water content
| 261 |
q_mqx | real*8 |   | public |   | Tree maximum water content per ground area
| 262 |
mqsst_ | real*8 |   | public |   | Tree water content at steady state
| 263 |
mqsstmin | real*8 |   | public |   | Tree water content at turgor loss point
| 264 |
q_md | real*8 |   | public |   | Tree dry mass per unit ground area
| 266 |
o_mdstore | real*8 |   | public |   | Wood water storage parameter of trees
| 267 |
o_rtdepth | real*8 |   | public |   | Tree rooting depth (m)
| 271 |
o_rgdepth | real*8 |   | public |   | Grass rooting depth
| 272 |
pos_slt | integer |   | public |   | Lowest soil layer containing tree roots
| 274 |
pos_slg | integer |   | public |   | Lowest soil layer containing grass roots
| 275 |
pos_ult | integer |   | public |   | Lowest soil layer containing tree roots within unsaturated zone
| 276 |
pos_ulg | integer |   | public |   | Lowest soil layer containing grass roots within unsaturated zone
| 277 |
changef | real*8 |   | public |   | Change factor for adjusting root surface area
| 279 |
rsurft_ | real*8 | (:) | public |   | Root surface area of trees in each layer
| 281 |
rsurftnew | real*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 |
rsurfgnew | real*8 | (:) | public |   | Adjusted root surface area of grasses in each layer for next day
| 284 |
rootlim | real*8 | (3,3,3) | public |   | Indicator whether root surface are was limiting root water uptake
| 286 |
rsoil | real*8 | (:) | public |   | Resistance to water flow towards roots in each soil layer
| 288 |
refft | real*8 | (:) | public |   | Relative root water uptake efficiency for trees in each layer
| 290 |
reffg | real*8 | (:) | public |   | Relative root water uptake efficiency for grasses in each layer
| 291 |
posmna | integer | (3) | public |   | Pointer to variable values that achieved maximum net assimilation
| 292 |
rrt_d | real*8 |   | public |   | Tree root respiration rate (mol/m2/s)
| 294 |
rrt_y | real*8 |   | public |   | Annual tree root respiration
| 295 |
rrg_d | real*8 |   | public |   | Grass root respiration
| 296 |
rrg_y | real*8 |   | public |   | Annual grass root respiration
| 297 |
prootm | real*8 | (:) | public |   | Root hydraulic head in each layer
| 299 |
sumruptkt_h | real*8 |   | public |   | Hourly total tree root water uptake
| 301 |
ruptkt__ | real*8 | (:) | public |   | Root water uptake rate perennial veg (m/s)
| 302 |
ruptkt_h | real*8 | (:) | public |   | Hourly root water uptake by trees in each layer
| 303 |
ruptkt_d | real*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_h | real*8 | (:) | public |   | Hourly root water uptake by grasses in each layer
| 306 |
ruptkg_d | real*8 | (:) | public |   | Daily root water uptake by grasses in each layer
| 307 |
perc_cov_veg | real*8 | (:) | public |   | Daily coverage of vegetation
| 308 |
i_alpha | real*8 |   | public |   | Initial slope of electron transport curve
| 315 |
i_cpccf | real*8 |   | public |   | Water transport costs per m root depth and m^2 cover
| 316 |
i_tcf | real*8 |   | public |   | Turnover cost factor for foliage (tc=i_tcf*LAI)
| 317 |
i_maxyear | integer |   | public |   | Number of years to process
| 318 |
i_testyear | integer |   | public |   | Number of years after which to perform initial test of netass
| 319 |
i_ha | real*8 |   | public |   | Temperature response parameter
| 320 |
i_hd | real*8 |   | public |   | Temperature response parameter
| 321 |
i_toptf | real*8 |   | public |   | Parameter to calculate adaptation of topt (range 0-1 for no to full adaptation)
| 322 |
i_toptstart | real*8 |   | public |   | Start parameter for topt to calculate jmax(temp in K)
| 323 |
i_rlratio | real*8 |   | public |   | Ratio of leaf respiration to photosynthetic capacity
| 324 |
i_lat | real*8 |   | public |   | geogr. latitude
| 328 |
i_mdtf | real*8 |   | public |   | Total dry mass of living tissues of trees per unit pc (g/m^2)
| 337 |
i_mqxtf | real*8 |   | public |   | Total water storage capacity in living tissues of trees per unit md
| 338 |
i_rrootm | real*8 |   | public |   | Root water uptake resistivity in soil
| 339 |
i_rsurfmin | real*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_rootrad | real*8 |   | public |   | Average fine root radius
| 342 |
i_prootmg | real*8 |   | public |   | Constant root balance pressure of 1.5 MPa in grasses
| 343 |
i_growthmax | real*8 |   | public |   | Parameter determining maximum daily growth increment of root surface area
| 344 |
i_incrcovg | real*8 |   | public |   | parameter determining maximum increment percentage of grass cover
| 345 |
i_incrjmax | real*8 |   | public |   | parameter determining maximum increment percentage of jmax25
| 346 |
i_incrlait | real*8 |   | public |   | parameter determining maximum increment percentage of lai trees
| 347 |
i_incrlaig | real*8 |   | public |   | parameter determining maximum increment percentage of lai grasses
| 348 |
i_extcoeffg | real*8 |   | public |   | extinction coefficient beer's law grasses
| 349 |
i_extcoefft | real*8 |   | public |   | extinction coefficient beer's law trees
| 350 |
i_trans_vegcov | real*8 |   | public |   | fraction of radiative energy reaching soil under full cover (0-1)
| 351 |
i_firstyear | integer |   | public |   | First year for the generation of hourly output in computation mode
| 353 |
i_lastyear | integer |   | public |   | Last year for the generation of hourly output in computation mode
| 354 |
i_write_h | integer |   | public |   | Flag to write out hourly input values after conversation from daily values
| 356 |
i_read_pc | integer |   | public |   | Flag to write out hourly input values after conversation from daily values
| 357 |
i_lai_function | integer |   | public |   | Switch to use 1) linear or 2) exponential LAI estimate, as function of cover
| 358 |
i_no_veg | integer |   | public |   | Flag to switch vegetation off (1=no vegetation)
| 359 |
c_epsln | real*8 |   | public |   | Soil porosity
| 364 |
c_maxhour | integer |   | public |   | Number of hours to process
| 366 |
c_maxday | integer |   | public |   | Number of days to process
| 367 |