Check age classes

[luyssaert]

The results of configurations without LCC should be identical irrespective of whether 1 or x age classes are used. A test run with r6980 showed this is no longer the case. When the vegetation is moved to the next age class, for example, GPP between two setups diverges quickly (7N, 81W, PFT2). Suggesting that one or several (new) variables are not properly handled or not handled at all in age_class_distr.

[luyssaert]

A conflict in slowproc that prevented using age classes in combination with impose_veg is fixed in r7333

[jjeong]

The same problem was observed with r7536, 1N 59W, PFT3: Diameter decreased and IND increased when age changed from age1 to age2, which did not happen without age class. This stemmed from that lstress_fac is not moved so recruitment is calculated maximum with age2.

---

Printed variables in stomate_prescribe.f90 without age class

We are calling stomate_prescribe.f90 for sapling recruitment in PFT

0 3
0 Density of stem 0.683279479178339
0 Mean seasonal transmitted light used in stomate_prescribe.f90 (ivm,ipts,nlevel
0 s_tot) 3 1 10
0 10 1.00021016535762
0 9 1.00021016535762
0 8 0.753268309090018
0 7 0.459064701011401
0 6 0.263385988964993
0 5 0.163530114598119
0 4 0.111529136957455
0 3 8.411678289400436E-002
0 2 6.791003804610461E-002
0 1 4.834454242356064E-002
0 lstress_fac from CN-CAN is 4.83445424235606 (%)
0 PFT, 3
0 canopy cover (LAI), 0.000000000000000E+000
0 used canopy_gap (%) is, 4.83445424235606
0 Original ind is, 6832.79479178339
0 new_ind is, 0.966891048471212

---

Printed variables in stomate_prescribe.f90 with age class

0 We are calling stomate_prescribe.f90 for sapling recruitment in PFT
0 7
0 Density of stem 0.683279474336626
0 Mean seasonal transmitted light used in stomate_prescribe.f90 (ivm,ipts,nlevel
0 s_tot) 7 1 10
0 10 1.00000000000000
0 9 1.00000000000000
0 8 1.00000000000000
0 7 1.00000000000000
0 6 1.00000000000000
0 5 1.00000000000000
0 4 1.00000000000000
0 3 1.00000000000000
0 2 1.00000000000000
0 1 1.00000000000000
0 lstress_fac from CN-CAN is 100.000000000000 (%)
0 PFT, 7
0 canopy cover (LAI), 0.000000000000000E+000
0 used canopy_gap (%) is, 100.000000000000
0 Original ind is, 6832.79474336626
0 new_ind is, 20.0000002000000

Note that lstress_fac is 1 for the below block and new_ind is much higher in the end.

[luyssaert]

This issue is solved in r7554 by adding the variable light_tran_to_level_season in the subroutine age_class_distr

[jjeong]

Even after the fix in r7554 there still are differences between simulation with age classes and without age classes.
Also in the simulations with off the recruitment, the differences are less but still, there are.