# ---
# jupyter:
#   jupytext:
#     text_representation:
#       extension: .py
#       format_name: light
#       format_version: '1.5'
#       jupytext_version: 1.14.4
#   kernelspec:
#     display_name: Casa
#     language: casa
#     name: casapy
# ---

# ## Flux density values for standard flux calibrators
#
# To convert correlation coefficients to absolute flux densities

# ### J0408-6545    
# RA: 04h08m20.4s     
# Decl: -65d45m09.6s

# Coefficients: a=-0.9790, b=3.3662, c=-1.1216, d=0.0861     
# $log_{10}(S) = a + b*log_{10}(f) +c*log_{10}(f)^2 + d*log_{10}(f)^3$     
# with $S$ in Jy and $f$ in MHz

# +
a=-0.9790
b=3.3662
c=-1.1216
d=0.0861

f = 1284.  # MHz (L-band)

log_S = a + b*np.log10(f) + c*np.log10(f)**2 + d*np.log10(f)**3

print('Calculated Stokes I for J0408-6545 {} [Jy]'.format(10**log_S))
# -

# ```
# setjy(vis=msfile, field='J0408-6545', scalebychan=True, standard='manual', fluxdensity=[17.1,0,0,0])
# ```

# ### J1939-6342    
# RA: 19h39m25.05s     
# Decl: -63d42m43.63s

# Coefficients: a=-30.7667 b=26.4908 c=-7.0977 d=0.605334     
# $log_{10}(S) = a + b*log_{10}(f) +c*log_{10}(f)^2 + d*log_{10}(f)^3$     
# with $S$ in Jy and $f$ in MHz

# +
a=-30.7667
b=26.4908
c=-7.0977
d=0.605334

f = 1284.  # MHz (L-band)

log_S = a + b*np.log10(f) + c*np.log10(f)**2 + d*np.log10(f)**3

print('Calculated Stokes I for J1939-6342 {} [Jy]'.format(10**log_S))
# -

# ```
# setjy(vis=msfile, field='J1939-6342', scalebychan=True, standard='Stevens-Reynolds 2016', fluxdensity=-1)
# ```

# ### J1331+3030    
# RA: 13h31m08.3s     
# Decl: +30d30m32.96s

# ```
# setjy(vis=msfile, field='3C286', scalebychan=True, standard='Perley-Butler 2013', fluxdensity=-1)
# ```



StarcoderdataJupyterScriptsDedupFiltered