         nbcs=option_count(trim("/material_phase[0]/subgridscale_parameterisations/GLS/scalar_field::GLSTurbulentKineticEnergy/prognostic/boundary_conditions"))         if (nbcs > 0) then             FLExit("You have automatic boundary conditions on, but some boundary conditions on the GLS 