C  CTFAPPLY : derived from CTFPLOT giving automatic application of C.T.F. to
C             phases input from MMBOX and output ready for ORIGTILT.
C  PROGRAM PLOTS SPOTS FROM MMBOX IN RECIPROCAL SPACE AND C.T.F. PLOT OF IMAGE.
C
C       V1.01   20.4.87         RH      changed to allow IQMAX = 9 to be passed.
C       V1.02   24.1.89         RH      plots box sizes proportional to 8.1-IQ.
C       V1.03   13.6.90         RH      Minor change to FORMAT statement 105.
C       V1.04   6.10.92         RH      extended o/p format h,k,A,P,IQ,BCK,CTF
C       V1.05   21.5.93         RH      compatible with Alliant - CCPDPN.
C       V1.06   14.4.95         RH      minor change to I/P FORMAT statement 99
C       V2.00   13.8.00         TSH     major change from plot82 to plot2k
C         "     13.6.01         TSH     P2K_FONT needed string terminator
C       V3.00   01.2.97         HS      PHACON,RESMAX, ps-file completed.
C       V3.01   02.2.98         HS      RESHIG,RESLOW
C               remember to change version number and date in title record.
C
C  INPUT PARAMETERS
C      CARD 1    AX AY BX BY ISIZE DSTEP XMAG
C      CARD 2    DFMID1 DFMID2 ANGAST CS KV RESMAX
C      CARD 3    ISER TITLE
C      CARD 4    PHACON
C      CARD 5    RESHIG,RESLOW
C
C         AX,AY   - LATTICE PARAMETERS (FROM NNBOX) OF (1,0) AND (0,1)
C         BX,BY   -                                     IN GRID UNITS.
C         ISIZE   - SIZE OF DENSITOMETERED ARRAY (E.G. 2048)
C         DSTEP   - DENSITOMETER STEPSIZE IN MICRONS
C         XMAG    - PRECISE MAGNIFICATION NORMALLY WORKED OUT FROM LATTICE
C                    PARAMETERS AND KNOWN CELL DIMENSIONS.
C         DFMID1  - DEFOCUS LEVEL (UNDERFOCUS +VE). IF DFMID2=DFMID1, IMAGE
C         DFMID2  -  IS NON-ASTIGMATIC. OTHERWISE AMOUNT OF DEFOCUS IN TWO
C                    ORTHOGONAL DIRECTIONS, DFMID1 BEING DEFOCUS IN DIRECTION
C                    ANGAST (DEGS) RELATIVE TO X AND Y OF THE FOURIER TRANSFORM
C         CS      - SPHERICAL ABERRATION IN MM.
C         KV      - E.M. ACCELERATING VOLTAGE
C         ISER    - SERIAL NUMBER AT HEAD OF OUTPUT FILE.
C         TITLE   - TITLE FOR OUTPUT FILE.
C         RESMAX  - Maximum resolution for output ps-file.
C         PHACON  - PHASECONTRAST, should be approx. 0.93(Cryo),0.65(neg.stain)
C         RESHIG  - HIGHER VALUE FOR RESOLUTION CUTOFF (JUST ON OUTPUT)
C         RESLOW  - LOWER  VALUE FOR RESOLUTION CUTOFF (JUST ON OUTPUT)
C
C   INPUT  DATASTREAM  'IN'
C   OUTPUT DATASTREAM  'OUT'