def METHOD_NAME(pb_str):
    b = phoebe.Bundle.default_binary()
    b.add_dataset('lc', times=[0.75,], dataset='lc01', passband=pb_str)
    b.add_dataset('mesh', times=[0.0], columns=['areas', 'mus', 'visibilities', 'abs_intensities@lc01'], dataset='mesh01')
    b['pblum@primary'] = 1.*u.solLum #* 0.99 # 0.99 is bolometric correction
    b['teff@primary'] = 1.*u.solTeff
    b['requiv@primary'] = 1.*u.solRad
    b['syncpar@primary'] = 14.61
    b['period@orbit'] = 1.*u.yr
    b['q@orbit'] = 3.986004e14/1.3271244e20   # (GM)_E / (GM)_Sun
    b['sma@orbit'] = 1.*u.au
    b['teff@secondary'] = (300, 'K')
    b['requiv@secondary'] = 1.*c.R_earth
    b['syncpar@secondary'] = 365.25
    b['distance@system'] = (1, 'au')
    b.set_value_all('irrad_method', 'none')
    if BLACKBODY:
        b.set_value_all('atm', value='blackbody')
        b.set_value_all('ld_mode', 'manual')
        b.set_value_all('ld_func', value='linear')
        b.set_value_all('ld_coeffs', value=[0.0])
    else:
        b.set_value_all('atm', component='secondary', value='blackbody')
        b.set_value_all('ld_func', component='primary', value='interp')
        b.set_value_all('ld_mode', component='secondary', value='manual')
        b.set_value_all('ld_func', component='secondary', value='linear')
        b.set_value_all('ld_coeffs', component='secondary', value=[0.0])
    return b