Check free surface canopy layers
In palm_dynamic.py
the relative plant canopy height calculation needs to be checked.
# plant canopy height
if 'lad' in ncs.variables.keys():
lad3 = ncs.variables['lad'][:]
# replace non-zero values with 1
lad3[lad3 != 0] = 1
# minimum index of nonzeo value along inverted z
lad = np.argmax(lad3[::-1], axis=0)
# inversion back and masking grids with no buildings
lad = lad3.shape[0] - lad
lad[np.max(lad3, axis=0) == 0] = 0
else:
lad = np.zeros([ny,nx])