program stuff2hamana implicit none real:: npix,pixsec,x,y,tx,ty,z,r,var(13) character(len=512)::filename,cnpix,cpixsec call getarg(1,filename) call getarg(2,cnpix) call getarg(3,cpixsec) read(cnpix,*)npix read(cpixsec,*)pixsec open(10,file=trim(filename)) do read(10,*,end=100)var(1:13) x = var(2) y = var(3) if((x>=0.).and.(y>=0.).and.(x<=npix).and.(y<=npix))then x = var(2)-npix/2.0 y = var(3)-npix/2.0 z = var(12) r = sqrt(x*x+y*y) !tx = atan(x/r)*r2d !ty = atan(y/r)*r2d tx = x*pixsec/3600. ty = y*pixsec/3600. write(*,*)tx,ty,z end if end do 100 close(10) end program stuff2hamana