     call get_option("/physical_parameters/gravity/magnitude", gravity_magnitude)     positions => extract_vector_field(state, "Coordinate")      tke => extract_scalar_field(state, "GLSTurbulentKineticEnergy")         psi => extract_scalar_field(state, "GLSGenericSecondQuantity")     