def brightness_temperature(beam_area, disp):
 beam = beam_area.to(si.sr).value
 nu = disp.to(si.GHz, spectral())
 def convert_Jy_to_K(x_jybm):
  factor = ((((2 * _si.k_B) * si.K) * (nu ** 2)) / (_si.c ** 2)).to(astrophys.Jy).value
  return ((x_jybm / beam) / factor)
 def convert_K_to_Jy(x_K):
  factor = (astrophys.Jy / (((2 * _si.k_B) * (nu ** 2)) / (_si.c ** 2))).to(si.K).value
  return ((x_K * beam) / factor)
 return [(astrophys.Jy, si.K, convert_Jy_to_K, convert_K_to_Jy)]

def brightness_temperature ( beam_area , disp ) : beam = beam_area . to ( si . sr ) . value nu = disp . to ( si . GHz , spectral ( ) ) def convert_Jy_to_K ( x_jybm ) : factor = ( ( ( ( 2 * _si . k_B ) * si . K ) * ( nu ** 2 ) ) / ( _si . c ** 2 ) ) . to ( astrophys . Jy ) . value return ( ( x_jybm / beam ) / factor ) def convert_K_to_Jy ( x_K ) : factor = ( astrophys . Jy / ( ( ( 2 * _si . k_B ) * ( nu ** 2 ) ) / ( _si . c ** 2 ) ) ) . to ( si . K ) . value return ( ( x_K * beam ) / factor ) return [ ( astrophys . Jy , si . K , convert_Jy_to_K , convert_K_to_Jy ) ]
defines the conversion between jy/beam and "brightness temperature" .