Secure your code as it's written. Use Snyk Code to scan source code in minutes - no build needed - and fix issues immediately.
5.099, 3.162, 1.414, 3.162, 5.099, 7.071, 7.280,
5.385, 3.606, 2.828, 2.236, 2.236, 2.828, 3.606, 6.325])
alts = 90. - zas
Nsrcs = alts.size
fluxes = np.ones_like(alts)
catalog = []
if isinstance(array_location, MoonLocation):
localframe = 'lunartopo'
mock_keywords['world'] = 'moon'
else:
localframe = 'altaz'
mock_keywords['world'] = 'earth'
source_coord = SkyCoord(alt=Angle(alts, unit=units.deg), az=Angle(azs, unit=units.deg),
obstime=time, frame=localframe, location=array_location)
icrs_coord = source_coord.transform_to('icrs')
ra = icrs_coord.ra
dec = icrs_coord.dec
names = np.array(['src' + str(si) for si in range(Nsrcs)])
stokes = np.zeros((4, 1, Nsrcs))
stokes[0, :] = fluxes
catalog = pyradiosky.SkyModel(names, ra, dec, stokes, 'flat')
if return_data:
catalog = SkyModelData(catalog)
if get_rank() == 0 and save:
np.savez('mock_catalog_' + arrangement, ra=ra.rad, dec=dec.rad, alts=alts, azs=azs,
fluxes=fluxes)
return catalog, mock_keywords