PURE FUNCTION get_max_layers_3D() INTEGER :: get_max_layers_3D get_max_layers_3D = layers_3D_C_reduction%max_num_layers_3D END FUNCTION get_max_layers_3D