program apercenters c program to derive centers of aperture plates of Sung fields c character*6 name integer rah, ram, decd, decm, ram_new c 1 format(a12,69x,2f8.3) 2 format(2f6.1,a6) 3 format(1x,I3,I3,f5.1, 2x, I3, f5.1,2x, a6) 10 format(a4) 100 open(unit=19,file='inp.field',status='old') open(unit=20,file='centers.out',status='new') c read(19,*) rah,ram, ras, decd, decm, decs do I=1,10 read(19,2, end=1000) delta_ra, delta_dec, name dec=decd + decm/60.0 + decs/3600. dec=dec*3.1416/180.0 ranew=60.0*ram + ras +delta_ra*60.0/(15.0*cos(dec)) ram_new = ranew/60.0 ras_new = ranew-ram_new*60.0 decnew=decm*60.0+decs +delta_dec*60.0 decnew=decnew/60.0 c write (20,3) rah,ram_new,ras_new,decd,decnew,name type 3, rah,ram_new,ras_new, decd, decnew,name enddo 1000 end