; ; Analyze the Goods tip/tilt star statistics ; ; mid RA/Dec: 12:34:00 +62:15:00 ;clearwins arcsec = !pi/(180.*3600.) microns = 1.e-6 mm = 1.e-3 degrees = !pi/180. ; dir = 'C:\Documents and Settings\Don Gavel\My Documents\Adaptive Optics\Keck\Keck NGAO\Keck NGAO Extragalactic Science Modeling\' ; file = 'GoodsStarLists_fromJasonMelborne\goods_n_lgs.dat' u = read_ascii(dir+file,comment_symbol=';') u = u.field1 ra = u[1,*]+u[2,*]/60.+u[3,*]/3600. dec = sign(u[4,*])*(u[4,*]+u[5,*]/60.+u[6,*]/3600.) bmag = u[7,*] rmag = u[8,*] degrees = !pi/180. degreesPerHour = 15. nstars = (size(u))[2] raHours = ra ra *= degreesPerHour ; degrees of RA ; window,/free,title='goods N' plot,-ra,dec,psym=1,yrange=[min(dec),max(dec)] dra = 1./(cos(dec[0]*degrees)*600.) ddec = 1./600. n = fix((max(ra)-min(ra))/dra)+1 m = fix((max(dec)-min(dec))/ddec)+1 ra0 = max(ra) dec0 = min(dec) ira = (ra0-ra)/dra < (n-1) idec = (dec-dec0)/ddec < (m-1) a = fltarr(n,m,3) aa = fltarr(n,m) bg = 19 aa -= bg aa[ira,idec] = -(rmag