Time-dependent land cover: boundary condition issue

[BSchilperoort]

In the following file, a boundary condition is set during initialization:

STEMMUS_SCOPE/src/+init/setBoundaryCondition.m

Lines 62 to 69 in a6cda26

	if strcmp(IGBP_veg_long(1:9)', 'Croplands') % ['Croplands']
	% NBChB: Moisture Bottom B.C.:
	% 1 --Specified matric head (BChB);
	% 2 --Specified flux(BChB);
	% 3 --Zero matric head gradient (Gravitiy drainage);
	NBChB = 1;
	else
	NBChB = 3;

This is only different for Croplands. However, as this is done during initialization, it will not be updated when land cover changes.

This is not fixed yet in #182

[Crystal-szj]

The boundary condition might also depend on specific site conditions. Would it be possible to consider NBChB as a site-specific property?