diff --git a/tangostationcontrol/tangostationcontrol/beam/old_test.py b/tangostationcontrol/tangostationcontrol/beam/old_test.py deleted file mode 100644 index ae4036987f82611fab01958cb763be026276b2f6..0000000000000000000000000000000000000000 --- a/tangostationcontrol/tangostationcontrol/beam/old_test.py +++ /dev/null @@ -1,159 +0,0 @@ -import datetime -from delays import * - -if __name__ == '__main__': - # # create a frame tied to the reference position - reference_itrf = [3826577.066, 461022.948, 5064892.786] # CS002LBA, in ITRF2005 epoch 2012.5 - d = delay_calculator(reference_itrf) - - # # set the timestamp to solve for - timestamp = datetime.datetime(2021,1,1,0,0,5) - d.set_measure_time(timestamp) - - # compute the delays for an antennas w.r.t. the reference position - antenna_itrf = [[3826923.546, 460915.441, 5064643.489]] # CS001LBA, in ITRF2005 epoch 2012.5 - - # # obtain the direction vector for a specific pointing - direction = "J2000","0deg","0deg" - - # calculate the delays based on the set reference position, the set time and now the set direction and antenna positions. - delays = d.convert(direction, antenna_itrf) - - # print the delays - # pprint.pprint(delays) - - - #test changing the time - - print(f"Changing timestamp test\nBase parametres: Direction: {direction}, position: {antenna_itrf}") - old_val = 0 - - for i in range(10): - # # set the timestamp to solve for - timestamp = datetime.datetime(2021,1,1,0,i,5) - d.set_measure_time(timestamp) - - old_val = delays - delays = d.convert(direction, antenna_itrf) - diff = delays[0] - old_val - - # print the delays - print(f"Timestamp: {timestamp}: {delays[0]}, diff: {diff}") - - - # reset time - timestamp = datetime.datetime(2021, 1, 1, 0, 0, 5) - d.set_measure_time(timestamp) - - - #test changing the antenna position - print(f"Changing Antenna position test.\nBase parametres: Time: {timestamp} Direction: {direction}") - old_val = 0 - - for i in range(10): - antenna_itrf = [[3826577.066 + i, 461022.948, 5064892.786]] # CS002LBA, in ITRF2005 epoch 2012.5 - - old_val = delays - delays = d.convert(direction, antenna_itrf) - diff = delays[0] - old_val - - # print the delays - print(f"Antenna position: {antenna_itrf}: {delays}, diff: {diff}") - - # test changing the direction - - antenna_itrf = [[3826923.546, 460915.441, 5064643.489]] # CS001LBA, in ITRF2005 epoch 2012.5 - # antenna_itrf = [[0, 0, 0]] # CS001LBA, in ITRF2005 epoch 2012.5 - - print(f"Changing direction test.\nBase parametres: Time: {timestamp} , position: {antenna_itrf}") - - old_val = 0 - - for i in range(12): - deg = i * 30 - direction = "J2000", f"0deg", f"{deg}deg" - - old_val = delays - delays = d.convert(direction, antenna_itrf) - diff = delays[0] - old_val - - # print the delays - print(f"Direction: {direction}: {delays}, diff: {diff}") - - # # # create a frame tied to the reference position - # reference_itrf = [3826577.066, 461022.948, 5064892.786] # CS002LBA, in ITRF2005 epoch 2012.5 - # d = delay_calculator(reference_itrf) - # - # # set the antenna position identical to the reference position - # speed_of_light = 299792458.0 - # antenna_itrf = [[reference_itrf[0], reference_itrf[1], reference_itrf[2]]] # CS001LBA, in ITRF2005 epoch 2012.5 - # - # for i in range(24): - # # # set the timestamp to solve for - # timestamp = datetime.datetime(2022, 1, 27, i, 28, 0) # sunrise this particular day - # d.set_measure_time(timestamp) - # - # # compute the delays for an antennas w.r.t. the reference position - # - # # # obtain the direction vector for a specific pointing - # direction = "SUN", "0deg", "0deg" - # - # # calculate the delays based on the set reference position, the set time and now the set direction and antenna positions. - # pointing = d.measure.direction(*direction) - # dir = d.get_direction_vector(pointing) - # - # if i == 9: - # print("sunrise") - # print(f"time: {i:2}:28,\tdirection: {dir}" ) - # - # print(d.measure.epoch()) - - """" - ############################################################ - """ - - # # # create a frame tied to the reference position - # reference_itrf = [3826577.066, 461022.948, 5064892.786] # CS002LBA, in ITRF2005 epoch 2012.5 - # d = delay_calculator(reference_itrf) - # - # # # set the timestamp to solve for - # timestamp = datetime.datetime(2021, 1, 1, 0, 0, 5) - # d.set_measure_time(timestamp) - # - # # compute the delays for an antennas w.r.t. the reference position - # antenna_itrf = [[3826923.546, 460915.441, 5064643.489]] # CS001LBA, in ITRF2005 epoch 2012.5 - # - # # # obtain the direction vector for a specific pointing - # direction = "AZEL", "0deg", "0deg" - # - # # calculate the delays based on the set reference position, the set time and now the set direction and antenna positions. - # delays = d.convert(direction, antenna_itrf) - # - # print(delays) - # # create a frame tied to the reference position - reference_itrf = [3826577.066, 461022.948, 5064892.786] # CS002LBA, in ITRF2005 epoch 2012.5 - d = delay_calculator(reference_itrf) - - # # set the timestamp to solve for - timestamp = datetime.datetime(2000, 1, 1, 0, 0, 0) - d.set_measure_time(timestamp) - - # # obtain the direction vector for a specific pointing - direction = "J2000", "0deg", "0deg" - - speed_of_light = 299792458.0 - antenna_itrf = [[reference_itrf[0], reference_itrf[1] - speed_of_light, reference_itrf[2]]] # CS001LBA, in ITRF2005 epoch 2012.5 - - # calculate the delays based on the set reference position, the set time and now the set direction and antenna positions. - delays = d.convert(direction, antenna_itrf) - - print("One second delay: ", delays) - - antenna_itrf = [[reference_itrf[0], reference_itrf[1] - speed_of_light - speed_of_light, reference_itrf[2]]] # CS001LBA, in ITRF2005 epoch 2012.5 - - - # calculate the delays based on the set reference position, the set time and now the set direction and antenna positions. - delays = d.convert(direction, antenna_itrf) - - print("Two second delay: ", delays) -