For a Fortran array DIMENSION A(I,J):
Note that this conforms to the existing IDI output routines, but that (i) the current GKS draw pixel array implementation is reversed in y and (ii) DAOPHOT’s coordinate system differs by 0.5 pixels in each direction.