def in_domain(x, y):
    """
    Check if a point x, y is within the geometry domain
    """
    l = 0.1
    h = 0.04
    xc = x - h
    yc = y - h

    check = True

    if (x > l) or (x < 0.0) or (y > l) or (y < 0.0):
        check = False
        return check

    if (xc > 0.0) & (yc > 0.0):
        check = False

    return check