Add max statement to ensure at least one axial layer per subregion

Minor bug fix

Merge request reports

Loading