def pbsize(freq, d=12.0):
    """Returns the primary beam size in arcseconds for a dish with diameter d (default 12) m at frequency freq GHz."""
    return 1.22 * 2.9978e-1 * 180.0 * 3600.0 / (3.1415927 * freq * d)