PROGRAM cpc2_chk * - runs interactively * - read columns unformatted of CPC2 catalog, release 2 * - output the minimal and maximal values for each column * MIZ 990303 * f77 cpc2_chk.f -o cpc2_chk IMPLICIT NONE * variables INTEGER dima !maximal number of columns PARAMETER (dima=9) INTEGER val(dima), minv(dima), maxv(dima), . i, j, n_all, . num1, mag, ra, dc, err_ra, . err_dc, n_obs, epoch_ra, epoch_dc CHARACTER fnin*20, answer*20 * initializing n_all = 0 fnin = 'cpc2_2.m' DO j=1,dima maxv(j)=-999000000 minv(j)= 999000000 ENDDO * start program WRITE (*,'(/a/)') 'program to check catalog entries' WRITE (*,'(3a,$)') 'file to check: ',fnin,' ' READ (*,'(a)') answer IF (answer.NE.' ') READ (answer,*) fnin OPEN (11,FILE=fnin,READONLY) * read data DO i=1,999999 READ (11,*,END=10) num1, mag, ra, dc, err_ra, . err_dc, n_obs, epoch_ra, epoch_dc n_all = n_all + 1 val(1) = num1 val(2) = mag val(3) = ra val(4) = dc val(5) = err_ra val(6) = err_dc val(7) = n_obs val(8) = epoch_ra val(9) = epoch_dc DO j=1,dima IF (val(j).GT.maxv(j)) maxv(j) = val(j) IF (val(j).LT.minv(j)) minv(j) = val(j) ENDDO ENDDO 10 CONTINUE * write to screen WRITE (*,'(2a)') 'input file name: ', fnin WRITE (*,'(a,i7)') 'number of lines: ', n_all WRITE (*,'(/a)')'minimal and maximal values for each column:' WRITE (*,'(i7,i4,2(i11),2(i4),i3,2(i6))') . (minv(j), j=1,dima) WRITE (*,'(i7,i4,2(i11),2(i4),i3,2(i6))') . (maxv(j), j=1,dima) END