procedure qdt1( outv, szoutv, nstr, out1fd, out2fd)

int     szoutv
double  outv[szoutv]
int     nstr
int     out1fd
int     out2fd

int	j
int     index
bool    ok
int     ofd
#
begin

do j = 1, nstr {

   index = (j-1)*22

   call fprintf(out1fd, "%9.2f%9.2f%6d%3d%10.3f%10.3f%13.7e%13.7e%13.7e%13.7e%4d\n")
     call pargd(outv[index+1])
     call pargd(outv[index+2])
     call pargi(nint(outv[index+3]))
     call pargi(nint(outv[index+4]))
     call pargd(outv[index+5])
     call pargd(outv[index+6])
     call pargd(outv[index+7])
     call pargd(outv[index+8])
     call pargd(outv[index+9])
     call pargd(outv[index+10])
     call pargi(nint(outv[index+11]))

   call fprintf(out2fd, "%9.2f%9.2f%6d%3d%10.3f%10.3f%13.7e%13.7e%13.7e%13.7e%4d\n")
     call pargd(outv[index+12])
     call pargd(outv[index+13])
     call pargi(nint(outv[index+14]))
     call pargi(nint(outv[index+15]))
     call pargd(outv[index+16])
     call pargd(outv[index+17])
     call pargd(outv[index+18])
     call pargd(outv[index+19])
     call pargd(outv[index+20])
     call pargd(outv[index+21])
     call pargi(nint(outv[index+22]))

}

do j = 1, 2 {
   ok = true
   if (j == 1) { ofd = out1fd }
   if (j == 2) {
      if (out1fd == out2fd) {
         ok = false
      } else {
         ofd = out2fd
      }
   }
   if (ok) {
      call fprintf (ofd, "#     (1)")
      call fprintf (ofd, "      (2)")
      call fprintf (ofd, "   (3)")
      call fprintf (ofd, " (4)")
      call fprintf (ofd, "      (5)")
      call fprintf (ofd, "       (6)")
      call fprintf (ofd, "          (7)")
      call fprintf (ofd, "          (8)")
      call fprintf (ofd, "          (9)")
      call fprintf (ofd, "         (10)")
      call fprintf (ofd, " (11)")
      call fprintf (ofd, "\n")
      call fprintf (ofd, "#\n")
      call fprintf (ofd, "# column description\n")
      call fprintf (ofd, "#    (1) X coordinate position [px]\n") 
      call fprintf (ofd, "#    (2) Y coordinate position [px]\n") 
      call fprintf (ofd, "#    (3) count\n") 
      call fprintf (ofd, "#    (4) image (1|2)\n") 
      call fprintf (ofd, "#    (5) stellar magnitude [0 mag -> 1 electrons]\n") 
      call fprintf (ofd, "#    (6) error of stellar magnitude [mag]\n") 
      call fprintf (ofd, "#    (7) signal within stellar aperture [electrons]\n") 
      call fprintf (ofd, "#    (8) area of stellar aperture [px]\n") 
      call fprintf (ofd, "#    (9) sky value of background aperture [electrons]\n") 
      call fprintf (ofd, "#   (10) error of sky value [electrons]\n") 
      call fprintf (ofd, "#   (11) number of pixels used to estimate sky value [px]\n") 
   }
}

end
#