C*MASKTRANA.FOR***************************************************************
C  Program for masking a Fourier transform                                   *
C     Input transform in standard format on stream 1 (IN)                    *
C     Output masked transform in standard format on stream 2 (OUT)           *
C     Control file contains parameters for generation of list of holes in    *
C         mask on stream 5, normal input stream (MASK)                       *
C                                                                            *
C     Mask file in free format                                               *
C                                                                            *
C  -  first record  ISHAPE, IAMPLIMIT, ISPOTFILE                             *
C     ISHAPE=1 Hard edge circular holes                                      *
C          2 Soft edge circular holes (Gaussian weighted to EXP(-2) at edge) *
C          3 Hard edge rectangular holes                                     *
C     IAMPLIMIT  if T then spots allowed through mask are limited in their   *
C                     maximum amplitude (they are reduced to 2x average)     *
C                if F then no amplitude reduction is applied                 *
C     ISPOTFILE  if F then input from stream 5 as in old program MASKTRAN    *
C                if T then input from file SPOTS
C     IHOLEFILE  if F then input holediameter from stream 5                  *
C                if T then input holediameter from file hole1.tmp
C                                                                            *
C  -  second record RAD       (ISHAPE = 1 or 2) radius (circular holes)      *
C                DELX,DELY (ISHAPE = 3)  half edge lengths(rectangular holes)*
C                                                                            *
C  -  third record AX,AY,BX,BY,IHMIN,IHMAX,IKMIN,IKMAX,RMAX,ITYPE            *
C       AX    etc, lattice parameters from refined NNBOX output in grid units*
C       IHMIN etc, set min/max limits on H,K for lattice generator           *
C       RMAX  Cutoff radius in transform for lattice generator in grid units *
C       ITYPE if = 0, then generate all spots                                *
C             if = 1, then include only spots of given index (see below)     *
C                                                                            *
C     Followed by record for each required spot to be allowed through mask   *
C             (only if ITYPE=1; input either from stream 5 or file SPOTS)    *
C         IH,IK                                                              *
C                                                                            *
C     Derived from TRMASK(Version 1.04)22-SEP-82 (RAC)                       *
C                                                                            *
C     Version 1.00    3-JUL-84        RH   See above.                        *
C     Version 1.01   27-JUL-87        RH   Option to limit maximum amplitude *
C     Version 1.02   31-MAR-88        RH   Debug spots for which XPOS=0      *
C          Converted to Alliant 12-NOV-91  RB   no changes                   *
C     Version 2.00   19-MAY-93        RH   read spots from file SPOTS        *
C     Version 2.01   17-APR-95        RH   check overflow of NHOLEMAX        *
C     Version 2.02   18-APR-95        RH   no (0,0) using lattice generator  * 
C     Version 2.03   25-JUL-95        RH   cosmetic removal of statement 99  * 
C     Version 2.04   03-SEP-95        RH   IDEEP=200 plus several checks for * 
C                                          hole geometry within transform    *
C                                                                            *
C*****************************************************************************