def get_gbsa_force(system):
    """Find and return the first force that contains 'GBSA' in its name"""
    for force in system.getForces():
        if "GBSA" in force.__class__.__name__:
            return force
    raise (Exception('No GBSA force found'))