if (zeta > -2. * pi) and (zeta < 2. * pi): if self.site_latitude.R <= 0.: zeta += pi else: zeta = telaz if zeta < 0: zeta = zeta + 2 * pi elif zeta > 2 * pi: zeta - 2 * pi return zeta * 180 / pi if __name__ == '__main__': import numpy as np dome_radius, mount_dec_height, mount_dec_length, mount_dec_offset = 147, 0, 49.2, 0 site = Site() Model = AzimuthModel(site['latitude'], dome_radius, mount_dec_height, mount_dec_length, mount_dec_offset) # for dra in np.arange(10, 200, 36): # for ddec in np.arange(1, 360, 10): for az, alt in [(ii, jj) for ii in np.arange(5, 360, 10) for jj in np.arange(25, 90, 20)]: tel_pos = Position.altAzToRaDec( Position.fromAltAz(Coord.fromD(alt), Coord.fromD(az)), site['latitude'], site.LST()) model = Model.solve_dome_azimuth(tel_pos, site.LST_inRads()) print 'here', alt, az, model, model - az