#!/usr/local/bin/python2.6 import pylab as pb import math as m from matplotlib.ticker import MultipleLocator decs = range(-60,90,10) pb.clf() ctr = 0 color = ['k','b','g','c','r','m'] for i in range(len(decs)): dec = decs[i] desc = pb.subplot(5,3,1+i) azRate = [] elRate = [] hourAngle = [] elevation = [] azimuth = [] f = open('verify.out.dec%d'%dec,'r') lines = f.readlines() for line in lines: if (line[0] == '!'): continue tokens = line.split() azimuth.append(float(tokens[6])) elevation.append(float(tokens[7])) hourAngle.append(float(tokens[8])*60) azRate.append(float(tokens[9])) elRate.append(float(tokens[10])) if (len(elRate) > 1): if (elRate[-1]/abs(elRate[-1]) != elRate[-2]/abs(elRate[-2])): zeropoint = 0.5*(hourAngle[-1]+hourAngle[-2]) zeropointAz = 0.5*(azimuth[-1]+azimuth[-2]) zeropointEl = 0.5*(elevation[-1]+elevation[-2]) pb.plot(hourAngle, elRate, 'b-') pb.ylim([-0.36,0.36]) pb.xlim([zeropoint-8,zeropoint+8]) desc.xaxis.grid(True,which='major') desc.yaxis.grid(True,which='major') fontsize = 10 pb.setp(desc.get_xticklabels(), fontsize=fontsize) pb.setp(desc.get_yticklabels(), fontsize=fontsize) majorLocator = MultipleLocator(2) desc.xaxis.set_major_locator(majorLocator) majorLocator = MultipleLocator(0.1) desc.yaxis.set_major_locator(majorLocator) pb.text(0.04,0.38,'az=%+.0f'%zeropointAz, transform=desc.transAxes,size=10) pb.text(0.04,0.22,'el=%+.2f'%zeropointEl, transform=desc.transAxes,size=10) pb.text(0.04,0.06,'dec=%+d'%dec, transform=desc.transAxes,size=10) pb.hold(True) f.close() ctr = ctr+1 if (i==1): pb.title('Sidereal tracking rates from Mauna Kea') if (i==2): yoff = 0.04 pb.text(1.02,1-yoff,'+1e-4', transform=desc.transAxes,size=10) pb.text(1.02,0.5-yoff,'0.0', transform=desc.transAxes,size=10) pb.text(1.02,0-yoff,'-1e-4', transform=desc.transAxes,size=10) pb.text(1.02,1.1,'deg/sec', transform=desc.transAxes,size=10) if (i == 10): pb.xlabel('Hour angle (minutes)') if (i == 6): pb.ylabel('Elevation tracking rate (arcsec/sec)') pb.savefig('elrates.png',dpi=200) pb.draw()